poj 3308 Paratroopers

题意：在一个矩阵中某些坐标点存在你的敌人，现有一些炮弹，他们可以摆在任意一行或者任意一列，并且可以消灭掉一行或者一列上的敌人。要求消灭掉所都敌人需要的最小花费。

最小点覆盖：
建立超级源点和每行连接，权值为放在该行炮弹的花费。
建立超级汇点和每列连接，权值为放在该列炮弹的花费。
存在敌人的行和列之间连一条边，权值为inf。
跑一遍最小割。

#include<cmath>
#include<stdio.h>
#include<string.h>
#include<iostream>
#define inf 10000000
using namespace std;
struct node{
int u,v,next;
double c;
}edge[1200];
int head[1200],pre[1200],cur[1200],dis[1200],gap[1200];
double flow[1200],a[1200],b[1200];
int T,m,n,l,x,y,e,start,end;
double ans,tot;

void add_adge(int u,int v,double c)
{
edge[e].u=u;  edge[e].v=v;  edge[e].c=c;  edge[e].next=head[u]; head[u]=e++;
edge[e].u=v;  edge[e].v=u;  edge[e].c=0;  edge[e].next=head[v]; head[v]=e++;
}

double sap()
{
double flow=0,aug=inf;
bool flag; int u;
for(int i=0; i<=n+m+1; i++)
{
cur[i]=head[i];
gap[i]=dis[i]=0;
}
gap[start]=n+m+2;
u=pre[start]=start;
while(dis[start]<n+m+2)
{
flag=0;
for(int &j=cur[u]; j!=-1; j=edge[j].next)
{
int v=edge[j].v;
if(edge[j].c>0&&dis[u]==dis[v]+1)
{
flag=1;
if(edge[j].c<aug) aug=edge[j].c;
pre[v]=u;
u=v;
if(u==n+m+1)
{
flow+=aug;
while(u!=0)
{
u=pre[u];
edge[cur[u]].c-=aug;
edge[cur[u]^1].c+=aug;
}
aug=inf;
}
break;
}
}
if(flag) continue;
int mindis=n+m+2;
for(int j=head[u]; j!=-1; j=edge[j].next)
{
int v=edge[j].v;
if(edge[j].c>0&&dis[v]<mindis)
{
mindis=dis[v];
cur[u]=j;
}
}
if((--gap[dis[u]])==0)
break;
gap[dis[u]=mindis+1]++;
u=pre[u];
}
return flow;
}

int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d%d%d",&m,&n,&l);
e=0;start=0,end=n+m+1;
memset(head,-1,sizeof(head));
for(int i=1;i<=m;i++)
{
scanf("%lf",&a[i]);
add_adge(0,i,log(a[i]));
}
for(int i=1;i<=n;i++)
{
scanf("%lf",&b[i]);
add_adge(m+i,m+n+1,log(b[i]));
}
for(int i=1;i<=l;i++)
{
scanf("%d%d",&x,&y);
add_adge(x,m+y,inf);
}
printf("%.4lf\n",exp(sap()));
}
return 0;
}