poj 3281 Dining(最大流)

题意：有N头牛，F种食物，D种饮料，每头牛有自己喜欢的食物和饮料，每种食物/饮料只能被一头牛享用，每头牛只能享用一种饮料和一种食物，求最后能有几头牛享用到自己的那份套餐（套餐是饮料+食物）。

解法：很容易想到建图方法，食物放左边，饮料放右边，牛放中间，然后按照给定的关系连边，但这里需要对牛拆点，拆点之间流量1.控制牛只选一份套餐；

#include<cstdio>
#include<cstring>
#include<map>
#include<vector>
#include<cmath>
#include<cstdlib>
#include<stack>
#include<queue>
#include <iomanip>
#include<iostream>
#include<algorithm>
using namespace std ;
const int N=500 ;
const int M=5000 ;
const int inf=1<<30 ;
struct node
{
	int u,v,c,next;
}edge[M];
int pre[N],cur[N],gap[N],dis[N],head[N];
int top ;


int add(int u ,int v,int c)
{
	edge[top].u=u;
	edge[top].v=v;
	edge[top].c=c;
	edge[top].next=head[u];
	head[u]=top++;
	edge[top].u=v;
	edge[top].v=u;
	edge[top].c=0;
	edge[top].next=head[v];
	head[v]=top++;
}

int sap(int s,int t,int nv)
{
	memset(dis,0,sizeof(dis));
	memset(gap,0,sizeof(gap));
	for(int i = 0 ; i < nv ; i++)   cur[i]=head[i];
	int minflow=inf ,flow = 0 ,u,v;
	pre[s]=u=s;
	gap[s]=nv;
	while(dis[s]<nv)
	{
		    loop :
		    	 for(int &j = cur[u]; j!=-1;j=edge[j].next)
		    	 {
		    	 	   v=edge[j].v ;
		    	 	   if(edge[j].c > 0 && dis[u] == dis[v]+1)
		    	 	   {
		    	 	   	       pre[v]=u;
		    	 	   	       u=v;
		    	 	   	       minflow = min(minflow,edge[j].c);
		    	 	   	       if(v==t)
		    	 	   	       {
		    	 	   	       	    for(u=pre[v];v!=s;v=u,u=pre[u])
		    	 	   	       	    {
		    	 	   	       	    	edge[cur[u]].c -= minflow ;
		    	 	   	       	    	edge[cur[u]^1].c += minflow ;
		    	 	   	       	    }
		    	 	   	       	    flow += minflow;
		    	 	   	       	    minflow = inf ;
		    	 	   	       }
		    	 	   	     goto loop ;
		    	 	   }
		    	 }
		       int mindis=nv ;
		       for(int i = head[u] ; i!=-1 ; i=edge[i].next)
		       {
		       	    v=edge[i].v;
		       	    if(edge[i].c>0 && dis[v] < mindis)
		       	    {
		       	    	     cur[u]=i;
		       	    	     mindis=dis[v];
		       	    }
		       }
		       if(--gap[ dis[u] ]==0) break;
		       gap[ dis[u]=mindis+1 ]++;
		       u=pre[u];
	}  
	 return flow;
}

int main()
{
	int n,f,d;
	while(~scanf("%d%d%d",&n,&f,&d))
	{
		   top=0;
		   memset(head,-1,sizeof(head));
		   int s=0,t=2*n+f+d+1;
		   for(int i = 1 ; i <= f; i++)
		       add(s,i,1);
		   for(int i = 1 ; i <= d ; i++)
		       add(2*n+f+i,t,1);
		   for(int i = 1 ; i <= n ; i++)	       
		       add(f+i,f+n+i,1);
		   for(int i = 1 ; i <= n ; i++)
		   {
		   	      int F,D,K;
				  scanf("%d%d",&F,&D);
		   	      while(F--)
		   	      {
		   	      	    scanf("%d",&K);
		   	      	    add(K,f+i,1);
		   	      }
		   	      while(D--)
		   	      {
		   	      	    scanf("%d",&K);
		   	      	    add(f+n+i,2*n+f+K,1);
		   	      }
		   }   
		   int ans = sap(s,t,t+1);
		   printf("%d\n",ans);
	}
	return 0;
}