POJ3281，Dining（二分匹配）

不难看出，这道题属于网络流的二分匹配应用，不过，在建模的时候要小心注意一点。此题当中有两个二分图，分别是1 ~ n->1 ~ f, 1 ~n -> 1 ~ d，建模要把两个二分图给联系，同时应建模过程中要注意，保持每个food和drink只能选一次以及每只cow只能挑选一次的特点。

本人所建立的网络模型是，引入源点与汇点s=0,t=f+2n+d+1，以及四个点集V1,V2,V3,V4，V1编号从1 ~ f，表示food，V2编号从f+1 ~ f+n，表示奶牛，V3编号从f+n+1 ~ f+2n，V4编号从f+2n+1 ~ f+2n+d。初始化时，s与V1中各点均连接一条容量为1的边，V2中编号为f+i的点与V3中编号为f+n+i的点连接一条容量为1的边，V4中各点与汇点t连接一条容量为1的边，这么做就可以保持上述两个特点（在AC之前本人WA了很多次，就是建模时忘了去保持这些特点）。接下来按照数据的输入连接V1与V2，V3与V4之间的边即可。

代码如下：

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<iostream>
using namespace std;
const int maxn=1e3+5;
const int maxm=1e6+5;
const int inf=0x7f7f7f7f;
int head[maxn],cur[maxn],cnt;
int n,f,d;
struct edge
{
	int from,to,val,next;
 }e[maxm];
 
void addedge(int u,int v,int w)
{
	e[cnt].from=u; e[cnt].to=v; e[cnt].val=w; e[cnt].next=head[u]; head[u]=cnt++;
	e[cnt].from=v; e[cnt].to=u; e[cnt].val=0; e[cnt].next=head[v]; head[v]=cnt++;		//注意，反向边的容量为0
}

void init()		//0,1~f, f+1~f+n, f+n+1~f+2n,f+2n+1~f+2n+d,f+2n+d+1; 
{
	memset(head,0,sizeof(head));
	cnt=2;
	for(int i=1;i<=f;i++)
		addedge(0,i,1);
	for(int i=1;i<=n;i++)
		addedge(f+i,f+n+i,1);
	for(int i=1;i<=d;i++)
		addedge(f+2*n+i,f+2*n+d+1,1);
	for(int i=1;i<=n;i++)
	{
		int a,b,c;
		scanf("%d%d",&a,&b);
		for(int j=1;j<=a;j++)
		{
			scanf("%d",&c);
			addedge(c,f+i,1);
		}
		for(int j=1;j<=b;j++)
		{
			scanf("%d",&c);
			addedge(f+n+i,f+2*n+c,1);
		}
	}
}

struct Dinic
{
	int s,t,dis[maxn];
	
	bool bfs()
	{
		memset(dis,0,sizeof(dis));
		queue<int> Q;
		Q.push(s);
		dis[s]=1;
		while(!Q.empty())
		{
			int fro=Q.front(); Q.pop();
			for(int i=head[fro]; i ;i=e[i].next)
			{
				if(e[i].val&&!dis[e[i].to])
				{
					Q.push(e[i].to);
					dis[e[i].to]=dis[fro]+1;
				}
			}
		}
		return dis[t]!=0; 
	}
	
	int dfs(int pos,int flow)
	{
		if(pos==t)	return flow;
		int ret=0;
		for(int i=cur[pos]; flow&&i ;i=e[i].next)
		{
			cur[pos]=i;
			if(e[i].val&&dis[e[i].to]==dis[pos]+1)
			{
				int temp=dfs(e[i].to,min(e[i].val,flow));
				if(temp>0)
				{
					e[i].val-=temp;
					e[i^1].val+=temp;
					flow-=temp;
					ret+=temp;
				}
				else dis[e[i].to]=-1;
			}
		}
		return ret;
	}
	
	int max_flow()
	{
		int ans=0;
		while(bfs())
		{
			for(int i=0;i<maxn;i++)	cur[i]=head[i];
			ans+=dfs(s,inf);
		}
		return ans;
	}
}dinic;

int main()
{
	scanf("%d%d%d",&n,&f,&d);
	init();
	dinic.s=0; dinic.t=f+2*n+d+1;
	printf("%d\n",dinic.max_flow());
	return 0;
}