题目:HDU-4292
题意:有N个人,F种食品,D种饮料,接下来一行表示F种食品的数量,再下一行D种饮料的数量,接下来的N行F列表示N个人对F种食品的可以接受情况:Y为可接受,N为不可接受。再接下来N行D列表示N个人对D种饮料可以接受的情况,Y与N和上面的代表相同。问最多能满足多少人,一人任意一种可接受的一个食品和一杯饮料。
题解:源点vs与每种食品建边,边权为当前食品的数量;每种饮料与汇点相连,边权为饮料的数量;N 个人点拆成一条有向边(两点先命名为左点与右点),边权为1 。每个人: 左点与能接受的食品种类点相连,边权为1;右点与能接受的饮料品种点相连,边权为1。最后跑一次最大流,最大流值就是最多能满足的人数。给每个点编号时注意一下。
代码:
#include<bits/stdc++.h>
#define N 100005
#define INF 0x7fffffff
using namespace std;
int n,F,D,s,t,c;
int head[N],cur[N],deep[N];
char ss[N];
struct ljh
{
int to,next,w;
}e[N*2];
inline void add(int x,int y,int z)
{
e[c].next=head[x];
e[c].w=z;
e[c].to=y;
head[x]=c++;
}
bool bfs(int s,int t)
{
memset(deep,0,sizeof(deep));
queue<int>q;
deep[s]=1;
q.push(s);
while(!q.empty())
{
int x=q.front();
q.pop();
for(int i=head[x];i!=-1;i=e[i].next)
{
int nex=e[i].to;
if(e[i].w>0&&deep[nex]==0)
{
deep[nex]=deep[x]+1;
q.push(nex);
}
}
}
return deep[t];
}
int dfs(int x,int t,int maxflow)
{
if(x==t)return maxflow;
int ans=0;
for(int i=cur[x];i!=-1;i=e[i].next)
{
int nex=e[i].to;
if(deep[nex]!=deep[x]+1||ans>=maxflow||e[i].w<=0)continue;
cur[x]=i;
int k=dfs(nex,t,min(e[i].w,maxflow-ans));
e[i].w-=k;
e[i^1].w+=k;
ans+=k;
}
if(ans==0)deep[x]=-2;
return ans;
}
int Dinic(int s,int t)
{
int ans=0;
while(bfs(s,t))
{
memcpy(cur,head,sizeof(head));
ans+=dfs(s,t,INF);
}
return ans;
}
int main()
{
while(~scanf("%d%d%d",&n,&F,&D))
{
c=0;
memset(head,-1,sizeof(head));
s=0;t=F+D+2*n+1;
for(int i=1;i<=F;i++)
{
int z;
scanf("%d",&z);
add(s,i,z);
add(i,s,0);
}
for(int i=1;i<=D;i++)
{
int z;
scanf("%d",&z);
add(i+F+2*n,t,z);
add(t,i+F+2*n,0);
}
for(int i=1;i<=n;i++)
{
add(i+F,i+n+F,1);
add(i+n+F,i+F,0);
scanf("%s",ss);
for(int j=1;j<=F;j++)
if(ss[j-1]=='Y')
{
add(j,i+F,1);
add(i+F,j,0);
}
}
for(int i=1;i<=n;i++)
{
scanf("%s",ss);
for(int j=1;j<=D;j++)
if(ss[j-1]=='Y')
{
add(i+F+n,j+F+2*n,1);
add(j+F+2*n,i+F+n,0);
}
}
printf("%d\n",Dinic(s,t));
}
}