POJ3281 Dining（网络流）

【题目大意】

给你N头牛，F种食物，D种饮料，每种食物和饮料只能给一头牛，每头牛都有食物和饮料的喜好，只有食物和饮料都是一头牛喜欢的，我们称这头牛就是满足的，问最多有多少头牛可以满足。

【输入格式】

第一行三个整数：N,F,D。

接下来N行，每行前两个数为fi,di，表示喜好数量，接下来fi个数表示该牛喜欢的食物，接下来di个数表示该牛喜欢的饮料。

【输出格式】

一个整数表示最多有多少头牛满足。

【样例输入】

4 3 3

2 2 1 2 3 1

2 2 2 3 1 2

2 2 1 3 1 2

2 1 1 3 3

【样例输出】

3

【备注】

N,F,D<=100.

【题目分析】

一道很好的网络流板题，关键在于建图。

首先，如果只有食物，那么这个题就是一道二分图最大匹配问题，现在加了饮料，因为是牛与食物和饮料一一对应，所以我们网络流建图的顺序一定为：食物（饮料）----牛----饮料（食物）。

但这么做有一个问题：一头牛可能会匹配多个食物或饮料，为了解决这个问题，我们将一头牛拆成两个点，这也是网络流的一种经典做法，每头牛向自己连一条流量为1的边，这样就能限制一头牛最多只能匹配一个食物和一个饮料。

最后建完图就跑一个最大流就行了。

【代码】

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
const int MAXN=410;
const int MAXM=1e5+10;
const int INF=0x3f3f3f3f;

int n,m,cnt,s,t;
int head[MAXN],cur[MAXN],depth[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];
queue<int> q;

void Add(int x,int y,int z)
{
	nxt[cnt]=head[x];
	head[x]=cnt;
	w[cnt]=z;
	to[cnt]=y;
	cnt++;
}

void add(int x,int y,int z)
{
	Add(x,y,z);
	Add(y,x,0);
}

bool bfs()
{
	while(!q.empty())
	  q.pop();
	memset(depth,0,sizeof(depth));
	depth[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=head[u];i!=-1;i=nxt[i])
		{
			int v=to[i];
			if(w[i]&&depth[v]==0)
			{
				depth[v]=depth[u]+1;
				q.push(v);
			}
		}
	}
	if(depth[t]==0)
	  return false;
	return true;
}

int dfs(int u,int dist)
{
	if(u==t)
	  return dist;
	for(int &i=cur[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(w[i]&&depth[v]==depth[u]+1)
		{
			int di=dfs(v,min(dist,w[i]));
			if(di>0)
			{
				w[i]-=di;
				w[i^1]+=di;
				return di;
			}
		}
	}
	return 0;
}

int dinic()
{
	int ans=0;
	while(bfs())
	{
		for(int i=s;i<=t;++i)
		  cur[i]=head[i];
		while(int d=dfs(s,INF))
		  ans+=d;
	}
	return ans;
}

int main()
{
	memset(head,-1,sizeof(head));
	int f,d;
	scanf("%d%d%d",&n,&f,&d);
	s=0,t=n*2+f+d+1;
	for(int i=1;i<=f;++i)
	  add(s,i,1);
	for(int i=1;i<=d;++i)
	  add(i+n*2+f,t,1);
	for(int i=1;i<=n;++i)
	  add(f+i,f+i+n,1);
	for(int i=1;i<=n;++i)
	{
		int fc,dc;
		scanf("%d%d",&fc,&dc);
		for(int j=1;j<=fc;++j)
		{
			int x;
			scanf("%d",&x);
			add(x,i+f,1);
		}
		for(int j=1;j<=dc;++j)
		{
			int x;
			scanf("%d",&x);
			add(i+n+f,x+n*2+f,1);
		}
	}
	printf("%d",dinic());
	return 0;
}