[HDU4292] Food 网络流

原题走这里

意思大概就是有多种食品和饮料，每种个有若干个，每个人各有偏好，可能只喜欢某几种食品和饮料

求如何每人分配一件食物和一件饮料，使得各人的偏好都被满足

很容易就能想到网络流

分层图网络流

一开始我们可能会想到建三层图，食物一层，人一层，饮料一层

源点向每一种食物连边，流量为该种食物的个数。类似的，饮料向汇点连边

然后人和他所喜欢的食物和饮料连边，看上去是这样的

然而无法限制每人只选一种食物和饮料

为了限制人这一层点的流量，我们可以将人这一层拆成两层

一层入点，一层出点，入点向出点连一条流量为1的边

直接跑ISAP/Dinic即可，Dinic据说会被卡

然而Dinic还是过了

代码如下：

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
inline int read() {
    int X=0,w=1;char ch=0;
    while(ch>'9'||ch<'0')ch=getchar();
    while(ch<='9'&&ch>='0') X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return X*w;
}
struct edge
{
	int v,p,l;
}e[500010];
int n,f,r,s=1,t,top,head[810],cur[810],d[810],q[1010];
char ch[210];
inline void addedge(int u,int v,int l)
{
	e[++top]=(edge){v,head[u],l};
	head[u]=top;
	e[++top]=(edge){u,head[v],0};
	head[v]=top;
}
bool bfs()
{
	memset(q,0,sizeof(q));
	memset(d,0,sizeof(d));
	int qh=1,qt=0;
	d[s]=1;
	q[++qt]=s;
	while(qh<=qt)
	{
		int u=q[qh++];
		for(register int i=head[u];i;i=e[i].p)
		{
			int v=e[i].v;
			if(e[i].l&&(!d[v]))
			{
				d[v]=d[u]+1;
				q[++qt]=v;
			}
		} 
	}
	return d[t];
}
int dfs(int u,int delta)
{
	if(u==t)
	{
		return delta;
	}
	register int ret=0;
	for(int &i=cur[u];i;i=e[i].p)
	{
		if(e[i].l&&(d[e[i].v]==d[u]+1))
		{
			int dd=dfs(e[i].v,min(e[i].l,delta-ret));
			ret+=dd;
			e[i].l-=dd;
			e[i^1].l+=dd;
			if(ret==delta)break;
		}
	}
	if(ret==0)
	{
		d[u]=-1;
	}
	return ret;
}
int Dinic()
{
	int ret=0; 
	while(bfs())
	{
		memcpy(cur,head,sizeof(head));
		ret+=dfs(s,INF);
	}
	return ret;
}
int main()
{
	while(cin>>n>>f>>r)
	{
		memset(head,0,sizeof(head));
		top=1;
		t=2+f+2*n+r;
		for(register int i=1;i<=f;i++)
		{
			addedge(1,1+i,read());
		}
		for(register int i=1;i<=r;i++)
		{
			addedge(1+f+2*n+i,t,read());
		}
		for(register int i=1;i<=n;i++)
		{
			addedge(1+f+i,1+f+n+i,1);
			scanf("%s",ch+1);
			for(register int j=1;j<=f;j++)
			{
				if(ch[j]=='Y')
				{
					addedge(1+j,1+f+i,1);
				}
			} 
		}
		for(register int i=1;i<=n;i++)
		{
			scanf("%s",ch+1);
			for(register int j=1;j<=r;j++)
			{
				if(ch[j]=='Y')
				{
					addedge(1+f+n+i,1+f+2*n+j,1);
				}
			} 
		}
		cout<<Dinic()<<endl;
	}
}