POJ 3308 Paratroopers(最小割/KM）

题意：一个矩阵，有的方格里有一些敌人，每行或每列都可以安放一架枪，每架枪都有一个花费，而且能消灭他所在的行或列的所有敌人，最后的花费为所有的枪话费的乘积

分析：如果是至少安放几个枪的话，直接就是最小点覆盖；如果求和最小的话，KM就可以解决，但是球的是product 乘积。。。

       log(x)就可以转成求和了。。。用最小割或KM均可

dinic代码：

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#define inf 99999
#define nMAX 150
using namespace std;
int head[nMAX],level[nMAX],qu[nMAX],s_edge,s,t;
struct Edge
{
    int u,v,nxt;
    double w;
}edge[10000];
void addedge(int u,int v,double w)
{
    s_edge++;
    edge[s_edge].u=u;
    edge[s_edge].v=v;
    edge[s_edge].w=w;
    edge[s_edge].nxt=head[u];
    head[u]=s_edge;

    s_edge++;
    edge[s_edge].u=v;
    edge[s_edge].v=u;
    edge[s_edge].w=0;
    edge[s_edge].nxt=head[v];
    head[v]=s_edge;
}

double dinic()
{
    int beg,tail,e;
    double ans=0.0;
    while(1)
    {
        beg=tail=0;
        qu[tail++]=s;
        memset(level,-1,sizeof(level));
        level[s]=0;
        while(beg<tail)
        {
            int u=qu[beg++];
           for(e=head[u];e;e=edge[e].nxt)
            {
                int v=edge[e].v;
                if(level[v]==-1&&edge[e].w>0)
                {
                    level[v]=level[u]+1;
                    qu[tail++]=v;
                    if(v==t)
                    {
                        beg=tail;break;
                    }
                }
            }
        }
        if(level[t]==-1)break;
        beg=s,tail=0;
        while(1)
        {
            if(beg==t)
            {
                double flow=(double)inf;
                int mark;
                for(int i=0;i<tail;i++)
                    if(edge[qu[i]].w<flow)
                        mark=i,flow=edge[qu[i]].w;
                for(int i=0;i<tail;i++)
                {
                    edge[qu[i]].w-=flow;
                    edge[qu[i]^1].w+=flow;
                }
                ans+=flow;
                tail=mark;
                beg=edge[qu[mark]].u;
            }
            for(e=head[beg];e;e=edge[e].nxt)
            {
                int v=edge[e].v;
                if(level[v]==level[beg]+1&&edge[e].w>0)
                    break;
            }
            if(e)
            {
                qu[tail++]=e;
                beg=edge[e].v;
            }
            else
            {
                if(tail==0)break;
                level[edge[qu[--tail]].v]=-1;
                beg=edge[qu[tail]].u;
            }
        }
    }
    return ans;
}
int main()
{
    int CASE,i,j,k,r,c,num;
    double w;
    scanf("%d",&CASE);
    while(CASE--)
    {
        scanf("%d%d%d",&r,&c,&num);
        s_edge=1;
        memset(head,0,sizeof(head));
        s=0,t=r+c+1;
        for(i=1;i<=r;i++)
        {
            scanf("%lf",&w);
            addedge(s,i,log(w));
        }
        for(i=1;i<=c;i++)
        {
            scanf("%lf",&w);
            addedge(r+i,t,log(w));
        }
        for(i=1;i<=num;i++)
        {
            scanf("%d%d",&j,&k);
            addedge(j,r+k,double(inf));
        }
        double ans=dinic();
        printf("%.4f\n",exp(ans));//my god，G++输出写成f，毛病真多,以后都用f吧
    }
    return 0;
}