Dining POJ - 3281

【题目大意】
有 F 种食物和 D 种饮料，每种食物或饮料只能供一头牛享用，且每头牛只享用一 种食物和一种饮料。现在有 N 头牛，每头牛都有自己喜欢的食物种类列表和饮 料种类列表，问最多能使几头牛同时享用到自己喜欢的食物和饮料。（ 1 <= F <= 100, 1 <= D <= 100, 1 <= N <= 100）

三分图，可以用网络流写 ，还要保证每个牛只有一个还需要拆点。

此题的建模方法比较有开创性。以往一般都是左边一个点集表示供应并与源相连， 右边一个点集表示需求并与汇相连。现在不同了，供应有两种资源，需求仍只有 一个群体，怎么办？其实只要仔细思考一下最大流的建模原理，此题的构图也不 是那么难想。最大流的正确性依赖于它的每一条 s-t 流都与一种实际方案一一对 应。那么此题也需要用 s-t 流将一头牛和它喜欢的食物和饮料“串”起来，而食 物和饮料之间没有直接的关系，自然就想到把牛放在中间，两边是食物和饮料， 由 s, t 将它们串起来构成一种分配方案。至此建模的方法也就很明显了：每种食 物 i 作为一个点并连边(s, i, 1)，每种饮料 j 作为一个点并连边(j, t, 1)，将每头牛 k 拆成两个点 k’, k’’并连边(k’, k’’, 1), (i, k’, 1), (k’’, j, 1)，其中 i, j 均是牛 k 喜欢的食物 或饮料。求一次最大流即为结果。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1000;
const int maxm=2*maxn*maxn;
const int inf=0x3f3f3f3f;
int head[maxn],cnt=0;
struct edge
{
    int v,nxt,w;
}edge[maxm*3+100];
void add_edge(int u,int v,int w)
{
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].nxt=head[u];
    head[u]=cnt++;
    edge[cnt].v=u;
    edge[cnt].w=0;
    edge[cnt].nxt=head[v];
    head[v]=cnt++;
}
int numh[maxn],h[maxn],curedge[maxn],pre[maxn];
int sap(int s,int t,int n)
{
    //printf("111\n");
    memset(numh,0,sizeof(numh));
    memset(h,0,sizeof(h));
    memset(pre,-1,sizeof(pre));
    int cur_flow,flow_ans=0,u,tmp,neck,i;
    for(i=1;i<=n;i++)
    {
        curedge[i]=head[i];
    }
    numh[0]=n;
    u=s;
    while(h[s]<n)
    {
        if(u==t)
        {
            cur_flow=inf;
            for(i=s;i!=t;i=edge[curedge[i]].v)
            {
                if(cur_flow>edge[curedge[i]].w)
                {
                    neck=i;
                    cur_flow=edge[curedge[i]].w;
                }
            }
            for(i=s;i!=t;i=edge[curedge[i]].v)
            {
                tmp=curedge[i];
                edge[tmp].w-=cur_flow;
                edge[tmp^1].w+=cur_flow;
            }
            flow_ans+=cur_flow;
            u=neck;
        }
        for(i=curedge[u];i!=-1;i=edge[i].nxt)
        {
            if(edge[i].w&&h[u]==h[edge[i].v]+1)
            {
                break;
            }
        }
        if(i!=-1)
        {
            curedge[u]=i;
            pre[edge[i].v]=u;
            u=edge[i].v;
        }
        else
        {
            if(0==--numh[h[u]])
                break;
            curedge[u]=head[u];
            for(tmp=n,i=head[u];i!=-1;i=edge[i].nxt)
            {
                if(edge[i].w)
                {
                    tmp=min(tmp,h[edge[i].v]);
                }
            }
            h[u]=tmp+1;
            ++numh[h[u]];
            if(u!=s)
                u=pre[u];
        }
    }
    return flow_ans;
}
int main ()
{
    int f,d,n;
    while(~scanf("%d%d%d",&n,&f,&d))
    {
        int xx,yy;
        memset(head,-1,sizeof(head));
        cnt=0;
        for(int i=1;i<=f;i++)
        {
            add_edge(f+2*n+1+d,i,1);
        }
        for(int i=1;i<=n;i++)
        {
            //printf("%d\n",i);
            scanf("%d%d",&xx,&yy);
            for(int j=1;j<=xx;j++)
            {
                int zz;
                scanf("%d",&zz);
                add_edge(zz,f+i,1);
            }
            add_edge(f+i,f+n+i,1);
            for(int j=1;j<=yy;j++)
            {
                //printf("%d %d\n",j,yy);
                int zz;
                scanf("%d",&zz);
                add_edge(f+n+i,f+2*n+zz,1);
            }
            //printf("111\n");
        }
        //printf("111\n");
        for(int i=1;i<=d;i++)
        {
            add_edge(f+2*n+i,f+2*n+d+2,1);
        }
        int res=sap(f+2*n+d+1,f+2*n+d+2,f+2*n+d+2);
        printf("%d\n",res);     
    }
}