链接: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;
}