poj3281--Dining（最大流（EK算法））

题目来源：http://poj.org/problem?id=3281

题意

给出n头牛，d杯饮料（都比我过得好。。。T^T），以及f份食物，每头牛有各自喜欢的食物和饮料，那么每头牛吃饱喝好才算可以，所以问给出的东西能令多少牛吃饱。。。（吃完就没啦。。。）

思路

或许这道题可以弄成三分图匹配（胡诌的。。）。。

但是这是一道最大流模板题，本菜采用的最糟糕的EK算法，主要的一点就是如何去构图，，，

对于样例，构建如下图：



为啥cow要到cow，为了防止cow被重复使用，用边去判断是否使用过。。

代码

#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100+10;
int n,f,d,s,t;
int F[maxn][maxn],D[maxn][maxn];
int c[maxn*5][maxn*5],vis[maxn*5],pre[maxn*5];
void init()
{
scanf("%d%d%d",&n,&f,&d);
memset(F,0
4000
,sizeof(F));
memset(D,0,sizeof(D));
for(int i=1; i<=n; i++)
{
int ff,dd;
scanf("%d%d",&ff,&dd);
while(ff--)
{
int x;
scanf("%d",&x);
F[i][x]=1;
}
while(dd--)
{
int y;
scanf("%d",&y);
D[i][y]=1;
}
}
}
bool bfs()
{
memset(vis,0,sizeof(vis));
queue<int> Q;
while(!Q.empty()) Q.pop();
Q.push(s);
vis[s]=1;
while(!Q.empty())
{
int u=Q.front();
Q.pop();
for(int v=s;v<=t;v++)
{
if(!vis[v]&&c[u][v]>0)
{
pre[v]=u;
if(v==t)
return 1;
Q.push(v);
vis[v]=1;
}
}
}
return 0;
}

void print_EK()
{
s=0,t=f+d+2*n+1;
int maxflow=0;
while(bfs())
{
int tt=t;
while(tt!=s)
{
c[pre[tt]][tt]--;
c[tt][pre[tt]]++;
tt=pre[tt];
}
maxflow++;
}
printf("%d\n",maxflow);

}
void build_graph()//一共f+2*n+d+2个点，0~f+2*n+d+1
{
memset(c,0,sizeof(c));
for(int i=1; i<=f; i++) //s->food
c[0][i]=1;
for(int i=f+2*n+1; i<=f+2*n+d; i++) //drink->t
c[i][f+2*n+d+1]=1;
for(int i=1; i<=n; i++)
{
for(int j=1; j<=f; j++) //food-->cow
if(F[i][j])
c[j][f+i]=1;
c[f+i][f+n+i]=1;//cow-->cow
for(int j=1; j<=d; j++) //cow-->drink
if(D[i][j])
c[f+n+i][f+2*n+j]=1;
}
}
int main()
{
init();
build_graph();//建图。。
print_EK();
return 0;
}
//网络流第一道题，，，患得患失。。好菜，，，