这道题最开始看见是在书上的二分图匹配上面看见的,(我是先学二分图在学网络流),当时一直不知道怎么建这个二部图,直到学了最大流之后才猛然想起这道题;
简单说下建图思路,主要是牛这个地方要拆点,如果不拆的话以下这组数据就不对了:
1 2 2
2 2 1 2 1 2
答案很明显是1,但是如果不拆点的话答案就是2,主要是如果不拆点的话就会出现一只牛可能同时对应多种饮料或者食物,就是会导致增广路多余正确值;
应为做了POJ 2391发现边表无论是时间上还是空间上都要优于二维数组,所以果断这道题目就用边表写了,但是WA死我了。。。。。
至今没想明白怎么回事,沮丧之际,顺手写了一个二维数组的。。。63MS水过。。。。。。不明白为什么。。。。
先是我用边表写的代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<queue>
#define MAXN 105
#define INF 1000000
using namespace std;
struct EDGE
{
int to,next,cost;
}edge[1600005];
int head[4*MAXN];
int leve[4*MAXN];
int judge[4*MAXN][4*MAXN];
int s,t,num_edge;
int N,F,D;
void add_edge(int x,int y)
{
edge[num_edge].to =y;
edge[num_edge].cost=1;
edge[num_edge].next=head[x];
head[x]=num_edge++;
edge[num_edge].to=x;
edge[num_edge].cost=0;
edge[num_edge].next=head[y];
head[y]=num_edge++;
return;
}
int BFS(){
int k;
queue<int> q;
memset(leve,-1,sizeof(leve));
leve[s]=0;
q.push(s);
while(!q.empty()){
k=q.front();
q.pop();
for(int i=head[k];i!=-1;i=edge[i].next)
{
if(leve[edge[i].to]==-1 && edge[i].cost)
{
leve[edge[i].to] = leve[k] + 1;
q.push(edge[i].to);
}
}
if(k==t) return 1;
}
return 0;
}
int DFS(int v,int sum)
{
if(v==t) return sum;
int temp=sum;
for(int i=head[v];i!=-1;i=edge[i].next)
{
int x=edge[i].to;
if(leve[x]==leve[v]+1&&edge[i].cost!=0)
{
int t=DFS(x,min(sum,edge[i].cost));
edge[i].cost-=t;
edge[i+1].cost+=t;
sum-=t;
}
}
return temp-sum;
}
int main()
{
int num_food,num_drink;
int x,ans;
while(scanf("%d%d%d",&N,&F,&D)!=EOF)
{
s=0,t=N+N+F+D+1;
ans=0;
num_edge=0;
memset(head,-1,sizeof(head));
memset(judge,0,sizeof(judge));
memset(edge,0,sizeof(edge));
for(int i=1;i<=N;i++)
{
if(!judge[F+i][F+N+i])
{
add_edge(F+i,F+N+i);
judge[F+i][F+N+i]=1;
}
}
for(int i=1;i<=N;i++)
{
scanf("%d%d",&num_food,&num_drink);
for(int j=0;j<num_food;j++)
{
scanf("%d",&x);
if(!judge[s][x])
{
add_edge(s,x);
judge[s][x]=1;
}
if(!judge[x][F+i])
{
add_edge(x,F+i);
judge[x][F+i]=1;
}
}
for(int j=0;j<num_drink;j++)
{
scanf("%d",&x);
if(!judge[F+N+i][F+N+N+x])
{
add_edge(F+N+i,F+N+N+x);
judge[F+N+i][F+N+N+x]=1;
}
if(!judge[F+N+N+x][t])
{
add_edge(F+N+N+x,t);
judge[F+N+N+x][t]=1;
}
}
}
while(BFS())
{
ans+=DFS(s,INF);
}
printf("%d\n",ans);
}
return 0;
}
说真的,死活看不出哪里有问题。。。。。
以下是AC代码:
#include<stdio.h>
#include<queue>
#include<iostream>
#include<string.h>
#include<algorithm>
#define MAXN 405
#define INF 10000000
using namespace std;
int leve[MAXN];
int map_flow[MAXN][MAXN];
int N,F,D;
int n,s,t;
bool BFS()
{
memset(leve,0,sizeof(leve));
leve[s]=1;
queue<int> Q;
Q.push(s);
while(!Q.empty())
{
int k=Q.front();
Q.pop();
for(int i=0;i<=t;i++)
{
if(!leve[i]&&map_flow[k][i])
{
leve[i]=leve[k]+1;
Q.push(i);
}
}
if(k==t) return true;
}
return false;
}
int DFS(int v,int sum)
{
if(v==t) return sum;
int res=sum;
for(int i=0;i<=t;i++)
{
if(leve[i]==(leve[v]+1)&&map_flow[v][i])
{
int t=DFS(i,min(map_flow[v][i],sum));
map_flow[v][i]-=t;
map_flow[i][v]+=t;
sum-=t;
}
}
return res-sum;
}
int main()
{
int ans,num_food,num_drink;
int x;
while(scanf("%d%d%d",&N,&F,&D)!=EOF)
{
s=0;
t=N+N+F+D+1;
ans=0;
memset(map_flow,0,sizeof(map_flow));
for(int i=1;i<=F;i++)
{
map_flow[s][i]=1;
}
for(int i=1;i<=N;i++)
{
map_flow[F+i][F+N+i]=1;
}
for(int i=1;i<=D;i++)
{
map_flow[F+N+N+i][t]=1;
}
for(int i=1;i<=N;i++)
{
scanf("%d%d",&num_food,&num_drink);
for(int j=0;j<num_food;j++)
{
scanf("%d",&x);
map_flow[x][F+i]=1;
}
for(int j=0;j<num_drink;j++)
{
scanf("%d",&x);
map_flow[F+N+i][F+N+N+x]=1;
}
n=t+1;
}
while(BFS()) ans+=DFS(s,INF);
printf("%d\n",ans);
}
return 0;
}