网络流 poj 3308 最小割

t个样例

n*m的矩阵 L个伞兵

给出每行每列装激光的花费

伞兵的位置

要求杀死所有伞兵 总费用为这些费用的乘积  求花费最小

建图  源点 ->   行   -> 列  -> 汇点

权            花费    inf       花费

花费相乘等于 log(花费)相加  都变成double  最后取e

 #include<stdio.h>
 #include<algorithm>
 #include<string.h>
 #include<queue>
 #include<math.h>

 using namespace std;

 int S,T,cnt;
 #define inf 100000000
 #define MAXN 2000
 #define MAXN1 3500
 int head[MAXN+];
 int vis[MAXN+];

 struct edg
 {
     int to,next;
     double w;
 }x[MAXN1*];

 void add(int u,int v,double w)
 {
     x[cnt].next=head[u];
     x[cnt].to=v;
     x[cnt].w=w;
     head[u]=cnt++;
 }
 int bfs()
 {
     queue<int>q1;
     memset(vis,-,sizeof(vis));
     vis[S]=;
     q1.push(S);

     while(!q1.empty())
     {
         int now=q1.front();
         q1.pop();
         int j;

         for(j=head[now];j!=-;j=x[j].next)
         {
             if(vis[x[j].to]<&&x[j].w>)
             {
                 vis[x[j].to]=vis[now]+;
                 q1.push(x[j].to);
             }
         }
     }
     return vis[T]!=-;

 }
 double dfs(int u,double w)
 {
     double ans=;
     if(u==T)
         return w;
     int j;
     for(j=head[u];j!=-;j=x[j].next)
     {
         if(vis[x[j].to]==vis[u]+&&x[j].w>)
         {
             double b=dfs(x[j].to,min(w-ans,x[j].w));
             x[j].w-=b;
             x[j^].w+=b;
             ans+=b;
         }
     }
     return ans;
 }

 int main()
 {
     int t;
     scanf("%d",&t);

     while(t--)
     {
         int n,m,l,i;
         scanf("%d%d%d",&n,&m,&l);
         S=,T=n+m+;
         cnt=;
         memset(head,-,sizeof(head));
         for(i=;i<=n;i++)
         {
             double a;
             scanf("%lf",&a);
             add(S,i,log(a)),add(i,S,);
         }
         for(i=;i<=m;i++)
         {
             double a;
             scanf("%lf",&a);
             add(i+n,T,log(a)),add(T,i+n,);
         }

         for(i=;i<=l;i++)
         {
             int a,b;
             scanf("%d%d",&a,&b);
             add(a,b+n,inf);
             add(b+n,a,);
         }
         double ans=;
         while(bfs())
              ans+=dfs(S,inf);

         printf("%.4lf\n",exp(ans));
     }
     return ;
 }