poj3308 最小割

因为行可以了，那列就不行，所以根据行列建立最小割模型。

然后这题精妙之处在于把乘法取对数后转化为加法，瞬间就简单了。

保证精度，C++AC ,16MS G++WA。

#include<stdio.h>
#include<string.h>
#include<queue>
#include<math.h>
#define maxn 120
#define INF 10000000
using namespace std;
struct node
{
    int to;
    double v;
    int flag;
    int next;
}edge[];
int pre[maxn],index,vis[maxn],n,m,L,st,se;
double col[],row[];
double min(double x,double y)
{
    return x<y?x:y;
}
void add(int x,int y,double z)
{
    edge[index].to=y;
    edge[index].v=z;
    edge[index].flag=index+;
    edge[index].next=pre[x];
    pre[x]=index++;
    edge[index].to=x;
    edge[index].v=;
    edge[index].flag=index-;
    edge[index].next=pre[y];
    pre[y]=index++;
}
double dfs(int u,double low)
{
    int i;
    double used=;
    if(u==se)return low;
    for(i=pre[u];i!=-;i=edge[i].next)
    {
        if(vis[edge[i].to]==vis[u]+&&edge[i].v>)
        {
            double a=dfs(edge[i].to,min(low-used,edge[i].v));
            edge[i].v-=a;
            edge[edge[i].flag].v+=a;
            used+=a;
            if(used-low>1e-)break;
        }
    }
    if(!used)vis[u]=-;
    return used;
}
int BFS()
{
    int i;
    queue<int>q;
    memset(vis,-,sizeof(vis));
    vis[]=;
    q.push();
    while(!q.empty())
    {
        int t=q.front();
        q.pop();
        for(i=pre[t];i!=-;i=edge[i].next)
        {
            if(vis[edge[i].to]<&&edge[i].v>)
            {
                vis[edge[i].to]=vis[t]+;
                q.push(edge[i].to);
            }
        }
    }
    return vis[se]!=-;
}
void init()
{
    int i;
    index=;
    st=,se=;
    memset(pre,-,sizeof(pre));
    scanf("%d%d%d",&n,&m,&L);
    for(i=;i<=n;i++)
    {
        scanf("%lf",&row[i]);
        add(st,i,log(row[i]));
    }
    for(i=;i<=m;i++)
    {
        scanf("%lf",&col[i]);
        add(i+,se,log(col[i]));
    }
    for(i=;i<L;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,+y,INF*1.0);
    }
}
void slove()
{
    double ans=;
    while(BFS())
    {
        while()
        {
            double a=dfs(,INF*1.0);
            if(a<1e-)break;
            ans+=a;
        }
    }
    printf("%.4lf\n",exp(ans));
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        init();
        slove();
    }
}