我爱崔老师系列之 最大流 拆点 EK poj3281 dining

题目连接：http://poj.org/submit?problem_id=3281

有人所用邻接矩阵会超时。。。然后。。。哥笑了~(*^__^*) 嘻嘻……

其实这题我觉得这不超时

这题一开始崔老师说把牛放中间我的意思是把人放在最前面，因为如果放在中间1头牛3物三水的话答案就会错误，但是如果把人放在前面感觉又会不大对= =。。然后崔老师就神一般的回答说拆点！顿时天亮了= =。。。其实我真心不知道拆点是神马。。。比赛有一道拆点我枚举匹配一下就过了。。虽然也是拆点的思想= =。。。

总之拆点就是为了先吃每个点的出度唯一（个人理解）而做的重复两次操作。博客的上一个acm computer其实就是拆点。。。但是没人点出来。。。还是崔老师是文化人毕竟出去学习过啊。。。在此给崔神跪了~。。。

代码

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <queue>
#define maxn 100000
using namespace std;
struct node
{
    int fod[105],drk[105],f,d;
}cow[105];
int cap[405][405];
int ek(int n)
{
    int f = 0,i;
    int flow[405][405] = {0};
    memset(flow,0,sizeof(flow));
    int temp[405],pre[405];
    memset(pre,0,sizeof(pre));

    queue<int>q;
    for(;;)
    {
        q.push(0);
        memset(temp,0,sizeof(temp));
        temp[0] = maxn;
        q.push(0);
        while(!q.empty())
        {
            int u = q.front();
            q.pop();
            for(i = 0; i <= n;i++)
            {
                if(!temp[i] && cap[u][i]>flow[u][i])
                {
                    q.push(i);
                    pre[i]=  u;
                    temp[i] = min(temp[u],cap[u][i]-flow[u][i]);
                }
            }
        }
        if(temp[n] == 0)
        break;

        for(i = n;i != 0;i = pre[i])
        {
            flow[pre[i]][i] += temp[n];
            flow[i][pre[i]] -= temp[n];
        }
        f += temp[n];
        //printf("%d****\n",temp[n]);
    }
    return f;
}
int main()
{
    int m,n,t,i,j,k;
    int u;
    int val;
    while(~scanf("%d %d %d",&m,&n,&t))
    {
        memset(cap,0,sizeof(cap));
        for(i = 1;i <= m;i++)
        {
            scanf("%d %d",&cow[i].f,&cow[i].d);
            for(j = 1;j <= cow[i].f;j++)
            {

                scanf("%d",&val);
                cap[val][n+i] = 1;
            }
            for(j = 1;j <= cow[i].d;j++)
            {
                scanf("%d",&val);
                cap[n+m+i][n+m+m+val] = 1;
            }
        }
        for(i = 1;i <= m;i++)
        cap[n+i][n+m+i] = 1;
        for(i = 1;i <= n;i++)
            cap[0][i] = 1;
        for(i = 1;i <= t;i++)
        cap[n+m+m+i][n+2*m+t+1] = 1;
                int res = ek(2*m+t+n+1);
        cout<<res<<endl;
    }
    return 0;
}