poj3281 dining
考点
- 网络流最大流
题意
一些人去吃饭,要求一个人只能拿一份菜和一份肉,每份菜和肉都是独一无二的,每个人有不同的口味,求最多有多少人能吃菜和吃肉。
思路
首先规定食物左边,牛在中间,饮料右边
乍一看很像二分图匹配,但是出现了三个不同的部分,所以舍去
所以我们可以设dp[i][j]为中间i号位选取第j个左边的最大值,不难得出dp不行
网络流看起来是一个不错的替代品但是如果按照原样来画点会发现特殊情况(本图省略超级源点和汇点 ):
我们模拟一下程序,不难发现一头牛竟然吃了两次!(吃的关系用一个直线表示):
那么如何解决呢?
这里我们可以把牛给拆成两个点,来确保一头牛只能吃一次
这个方法对四分图等也都有效
额外
源码中有两行SuperStart=(n+a+b)*2+1;SuperEnd=(n+a+b+2)*2+2;
因为调试的时候输出了4(样例,下同),输出forline(i,11)11为原SuperStart(n+a+b+1)
时发现直接连接到了12//原SuperEnd
,调大了就好了
这是我写过最长的题解博客了吧2333
代码
#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int maxnode=30100;
const int maxline=61100;
struct nn
{
int nxt;
int to;
int value;
}edge[maxnode];
int head[maxnode];
int will,cnt;
#define forline(i,x) for(int i=head[x];i!=-1;i=edge[i].nxt)
int n,a,b;
int kn1[101];
int kn2[101];
int ka[101];
int kb[101];
int SuperStart,SuperEnd;
int CodeName;
queue<int> team;
int dis[maxnode];
bool Dinic_BFS(int startnode)
{
for(int i=0;i<maxnode;++i) dis[i]=-1;
team.push(startnode);
dis[startnode]=0;
int teamhead;
while(!team.empty())
{
teamhead=team.front();
team.pop();
forline(i,teamhead)
{
will=edge[i].to;
if(dis[will]==-1&&edge[i].value>0)
{
team.push(will);
dis[will]=dis[teamhead]+1;
}
}
}
if(dis[SuperEnd]==-1) return 0;
return 1;
}
int Dinic_ZG(int nownode,int wave)
{
if(nownode==SuperEnd)
{
return wave;
}
int tmp;
forline(i,nownode)
{
will=edge[i].to;
if(dis[will]==dis[nownode]+1&&edge[i].value>0)
{
if((tmp=Dinic_ZG(will,1)))
{
edge[i].value-=tmp;
edge[i^1].value+=tmp;
return tmp;
}
}
}
return 0;
}
void addEdge(int x,int y,int val)
{
edge[cnt].to=y;
edge[cnt].nxt=head[x];
edge[cnt].value=val;
head[x]=cnt;
cnt++;
}
int main()
{
cnt=0;
int d;
int ans=0;
for(int i=0;i<maxnode;++i) head[i]=-1;
scanf("%d%d%d",&n,&a,&b);
SuperStart=(n+a+b)*2+1;
SuperEnd=(n+a+b+2)*2+2;
for(int i=1;i<=a;++i)
{
ka[i]=++CodeName;
addEdge(SuperStart,ka[i],1);
addEdge(ka[i],SuperStart,0);
}
for(int i=1;i<=b;++i)
{
kb[i]=++CodeName;
addEdge(kb[i],SuperEnd,1);
addEdge(SuperEnd,kb[i],0);
}
for(int i=1;i<=n;++i)
{
kn1[i]=++CodeName;
kn2[i]=++CodeName;
addEdge(kn1[i],kn2[i],1);
addEdge(kn2[i],kn1[i],0);
}
int x,y;
for(int i=1;i<=n;++i)
{
scanf("%d%d",&x,&y);
for(int j=1;j<=x;++j)
{
scanf("%d",&d);
addEdge(ka[d],kn1[i],1);
addEdge(kn1[i],ka[d],0);
}
for(int j=1;j<=y;++j)
{
scanf("%d",&d);
addEdge(kn2[i],kb[d],1);
addEdge(kb[d],kn2[i],0);
}
}
int tmp;
while(Dinic_BFS(SuperStart))
{
while((tmp=Dinic_ZG(SuperStart,1)))
{
ans+=tmp;
}
}
printf("%d\n",ans);
}