您的位置:首页 > 其它

poj 3281 最大流 拆点构图,关键在于构图

2012-07-28 09:38 302 查看
/*

最大流,模版题目,要拆点

题目大意:给出N个牛喜欢的食物和饮料列表,每种食物和饮料只能用一种,求能满足多少只牛?

思路:把牛拆成两个点,连边e(i,n+i)=1

食物和i连i边e(food,i)=1,n+i和饮料连边e(n+i,drink)=1

食物和s连边e(s,i)=1, 饮料和t连边(i,)=1

*/

#include <iostream>

#include <cstdio>

#include <cstring>

#include <queue>

#include <algorithm>

using namespace std;

#define M 40010

int gap[M],dis[M],pre[M],cur[M];

int NE,NV;

int head[M];

struct Node{

int c,pos,next;

}E[999999];

#define FF(i,NV) for(int i=0;i<NV;i++)

int sap(int s,int t) {

//memset(pre,-1,sizeof(pre));

memset(dis,0,sizeof(int)*(NV+1));

memset(gap,0,sizeof(int)*(NV+1));

FF(i,NV) cur[i] = head[i];

int u = pre[s] = s,maxflow = 0,aug = -1;

gap[0] = NV;

while(dis[s] < NV) {

loop:for(int &i = cur[u]; i != -1; i = E[i].next) {

int v = E[i].pos;

if(E[i].c && dis[u] == dis[v] + 1) {

aug=min(aug,E[i].c);

pre[v] = u;

u = v;

if(v == t) {

maxflow += aug;

for(u = pre[u];v != s;v = u,u = pre[u]) {

E[cur[u]].c -= aug;

E[cur[u]^1].c += aug;

}

aug = 999999;

}

goto loop;

}

}

int mindis = NV;

for(int i = head[u]; i != -1 ; i = E[i].next) {

int v = E[i].pos;

if(E[i].c && mindis > dis[v]) {

cur[u] = i;

mindis = dis[v];

}

}

if( (--gap[dis[u]]) == 0)break;

gap[ dis[u] = mindis+1 ] ++;

u = pre[u];

}

return maxflow;

}

void addEdge(int u,int v,int c ) {

E[NE].c = c;E[NE].pos = v;

E[NE].next = head[u];head[u] = NE++;

E[NE].c = 0;E[NE].pos = u;

E[NE].next = head[v];head[v] = NE++;

}

int N, F, D,idx;

int main()

{

while (scanf("%d %d %d", &N, &F, &D) != EOF)

{

NV=F+2*N+D+2;

FF(i,NV) head[i] = -1;NE = 0;

memset(head, -1, sizeof(head));

idx = 0;

int a, b, f, d;

for (int i = 1; i <= N; ++i)

{

scanf("%d %d", &a, &b);

while (a--)

{

scanf("%d", &f);

addEdge(f, F + i, 1);

}

while (b--)

{

scanf("%d", &d);

addEdge(F + N + i, F + N + N + d, 1);

}

}

int source = 0, sink = N + N + F + D + 1;

//int n = sink + 1;

for (int i = 1; i <= F; ++i)

addEdge(source, i, 1);

for (int i = 1; i <= N; ++i)

addEdge(F + i, F + N + i, 1);

for (int i = 1; i <= D; ++i)

addEdge(F + N + N + i, sink, 1);

printf("%d\n", sap(source, sink));

}

return 0;

}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: