poj3281--Dining（最大流（EK算法））

题目来源：http://poj.org/problem?id=3281

题意

给出n头牛，d杯饮料（都比我过得好。。。T^T），以及f份食物，每头牛有各自喜欢的食物和饮料，那么每头牛吃饱喝好才算可以，所以问给出的东西能令多少牛吃饱。。。（吃完就没啦。。。）

思路

或许这道题可以弄成三分图匹配（胡诌的。。）。。
但是这是一道最大流模板题，本菜采用的最糟糕的EK算法，主要的一点就是如何去构图，，，
对于样例，构建如下图：

为啥cow要到cow，为了防止cow被重复使用，用边去判断是否使用过。。

代码

#include<cmath>#include<stack>#include<queue>#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int maxn=100+10;int n,f,d,s,t;int F[maxn][maxn],D[maxn][maxn];int c[maxn*5][maxn*5],vis[maxn*5],pre[maxn*5];void init(){    scanf("%d%d%d",&n,&f,&d);    memset(F,0,sizeof(F));    memset(D,0,sizeof(D));    for(int i=1; i<=n; i++)    {        int ff,dd;        scanf("%d%d",&ff,&dd);        while(ff--)        {            int x;            scanf("%d",&x);            F[i][x]=1;        }        while(dd--)        {            int y;            scanf("%d",&y);            D[i][y]=1;        }    }}bool bfs(){    memset(vis,0,sizeof(vis));    queue<int> Q;    while(!Q.empty()) Q.pop();    Q.push(s);    vis[s]=1;    while(!Q.empty())    {        int u=Q.front();        Q.pop();        for(int v=s;v<=t;v++)        {            if(!vis[v]&&c[u][v]>0)            {                pre[v]=u;                if(v==t)                    return 1;                Q.push(v);                vis[v]=1;            }        }    }    return 0;}void print_EK(){    s=0,t=f+d+2*n+1;    int maxflow=0;    while(bfs())    {        int tt=t;        while(tt!=s)        {            c[pre[tt]][tt]--;            c[tt][pre[tt]]++;            tt=pre[tt];        }        maxflow++;    }    printf("%d\n",maxflow);}void build_graph()//一共f+2*n+d+2个点，0~f+2*n+d+1{    memset(c,0,sizeof(c));    for(int i=1; i<=f; i++) //s->food        c[0][i]=1;    for(int i=f+2*n+1; i<=f+2*n+d; i++) //drink->t        c[i][f+2*n+d+1]=1;    for(int i=1; i<=n; i++)    {        for(int j=1; j<=f; j++) //food-->cow            if(F[i][j])                c[j][f+i]=1;            c[f+i][f+n+i]=1;//cow-->cow        for(int j=1; j<=d; j++) //cow-->drink            if(D[i][j])                c[f+n+i][f+2*n+j]=1;    }}int main(){    init();    build_graph();//建图。。    print_EK();    return 0;}//网络流第一道题，，，患得患失。。好菜，，，