POJ - 3281 Dining （最大流+拆点+建图）

链接：https://cn.vjudge.net/problem/POJ-3281

题意：n头牛，F种食物，D种饮料。每头牛都有若干个喜好的食物和饮料。既给牛喜好的食物又给牛喜好的饮料才会是牛满意，问最多可以使多少头牛满意。

思路：如果只有食物或饮料，那么这个问题就成了简单的二分图最大匹配问题。但此题既有食物又有饮料，有两个限制。那么我们就可以把牛拆成两个点，内部之间连一条边，容量为1；然后一个点和喜爱的食物连边，容量为1；另一个点和喜爱的饮料连边，容量也为1。（如果不拆点，一个牛会匹配到多个组合，也就是牛、饮料、食物都会被重复利用。）跑一遍最大流即可。

盗的图：

PS:太菜了。。。。。。。。。。。。。。。。。。。。。。。想了一个错误的思路（就是把每种食物和饮料的搭配就看做一个点，然后源点和每种食物连边容量为1，然后每种食物和每种搭配连边容量为1。最后，每头牛和他喜好的所有搭配连边容量为1，每头牛和汇点连边容量为1。），现在不知道为啥错，望有大佬指导。

 

#include <cstdio>
#include <cstring>
#include <queue>
#define ll long long
using namespace std;
const int N = 1e3+10;
const int M = 1e5+10;
int head[N],cnt,deep[N],cur[N];
const int inf = 0x3f3f3f3f;
int n,F,D,s,t;
struct node
{
	int to,ca,next;
}g[M];
void Init()
{
	cnt=0;
	s=0;
	t=2*n+F+D+1;
	for(int i=s;i<=t;i++)
		head[i]=-1;
	//memset(head,-1,sizeof(head));
	return ;
}
void add(int u,int v,int ca)
{
	g[cnt].to=v;
	g[cnt].ca=ca;
	g[cnt].next=head[u];
	head[u]=cnt++;
	return ;
}
bool bfs()
{
	//memset(deep,0,sizeof(deep));
	for(int i=s;i<=t;i++)
		deep[i]=0;
	queue<int> q;
	int u,v;
	q.push(s);
	deep[s]=1;
	while(!q.empty())
	{
		u=q.front();
		q.pop();
		if(u==t) return 1;
		for(int i=head[u];i!=-1;i=g[i].next)
		{
			v=g[i].to;
			if(g[i].ca>0&&!deep[v])
			{
				deep[v]=deep[u]+1;
				q.push(v);
			}
		}
	}
	return deep[t]!=0;
}
int dfs(int u,int flow)
{
	if(u==t||!flow) return flow;
	int v,ans=0,nowflow;
	for(int& i=cur[u];i!=-1;i=g[i].next)
	{
		v=g[i].to;
		if(g[i].ca>0&&deep[v]==deep[u]+1)
		{
			nowflow=dfs(v,min(flow,g[i].ca));
			if(nowflow)
			{
				ans+=nowflow;
				flow-=nowflow;
				g[i].ca-=nowflow;
				g[i^1].ca+=nowflow;
				if(!flow) break;
			}
		}
	}
	if(!ans) deep[u]=0;
	return ans;
}
int Dinic()
{
	int maxflow=0,flow;
	while(bfs())
	{
		//memcpy(cur,head,sizeof(head));
		for(int i=s;i<=t;i++)
			cur[i]=head[i];
		while(flow=dfs(s,inf))
			maxflow+=flow;
	}
	return maxflow;
	
}
int main(void)
{
	
	scanf("%d%d%d",&n,&F,&D);
	Init();
	for(int i=1;i<=F;i++)
	{
		add(s,2*n+i,1);
		add(2*n+i,s,0);
	}	
	for(int i=1;i<=n;i++)
	{
		int x,y,u;
		scanf("%d%d",&x,&y);		
		for(int j=1;j<=x;j++)
		{
			scanf("%d",&u);
			add(2*n+u,i,1);
			add(i,2*n+u,0);	
		}
		for(int j=1;j<=y;j++)
		{
			scanf("%d",&u);
			add(n+i,2*n+F+u,1);
			add(2*n+F+u,n+i,0);	
		}			
	}
	for(int i=1;i<=n;i++)
	{
		add(i,n+i,1);
		add(n+i,i,0);		
	}
	for(int i=1;i<=D;i++)
	{
		add(2*n+F+i,t,1);
		add(t,2*n+F+i,0);		
	}
	printf("%d\n",Dinic());
	
	
	
	return 0;
}