http://poj.org/problem?id=3281
题意:牛只能吃一种食物和喝一个饮料,牛喜欢多种食物和饮料,现在每种食物饮料都只有一个,问你一共有多少牛可以吃上喝上
题解:难在建图
我想的是源点->食物->牛->饮料->汇点,容量都为1 ,但这样搞会出现一头牛吃了好几种食物(会有多种食物指向一头牛)再流向汇点所以不行,
这时候拆点!!(见到网络流题一定要想拆点),把牛点拆成牛->牛(容量为1)这样从牛往下面流的时候,因为容量限制1,你只能流1了,这样就解决了上面的问题,最后跑到汇点的值就是该网络的最大流了(多少牛吃上喝上)。
#include<cstdio>
#include<cstring>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<stdlib.h>
#include<queue>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxn=1000000;
typedef long long LL;
const int N = 20005;
const int INF = 0x3f3f3f3f;
bool vis[N];
struct Edge{
int to, cap, flow, next;
}edge[N*50];
int n, m, cnt;//n是点 m是边 cnt是加边操作后的边
int head[N];//邻接表
int dis[N];//分层 等级
int cur[N];//弧优化
void add(int u, int v, int w){
edge[cnt] = (struct Edge){v, w, 0, head[u]};
head[u] = cnt++;
edge[cnt] = (struct Edge){u, 0, 0, head[v]};
head[v] = cnt++;
}
bool bfs(int start, int endd){//分层
memset(dis, -1, sizeof(dis));
memset(vis, false, sizeof(vis));
queue<int>que;
dis[start] = 0;
vis[start] = true;
que.push(start);
while(!que.empty()){
int u = que.front();
que.pop();
for(int i = head[u]; i != -1; i = edge[i].next){
Edge E = edge[i];
if(!vis[E.to] && E.flow<E.cap){
dis[E.to] = dis[u]+1;
vis[E.to] = true;
if(E.to == endd) return true;
que.push(E.to);
}
}
}
return false;
}
int dfs(int x, int res, int endd){ //增广
if(x == endd || res == 0) return res;
int flow = 0, f;
for(int& i = cur[x]; i != -1; i = edge[i].next){
Edge E = edge[i];
if(dis[E.to] == dis[x]+1){
f = dfs(E.to, min(res, E.cap-E.flow), endd);
if(f>0){
edge[i].flow += f;
edge[i^1].flow -= f;
flow += f;
res -= f;
if(res == 0) break;
}
}
}
return flow;
}
int max_flow(int start, int endd){
int flow = 0;
while(bfs(start, endd)){
memcpy(cur, head, sizeof(head));//初始化弧优化数组
flow += dfs(start, INF, endd);
}
return flow;
}
void init(){//初始化
cnt = 0;
memset(head, -1, sizeof(head));
}
int main(){
init();
int n,f,d;
cin>>n>>f>>d;
int sp=0;
int tp=2*n+f+d+1;
for(int i=1;i<=f;i++){
add(sp,i,1);
}
for(int i=1;i<=n;i++){
int x,y;
cin>>x>>y;
for(int j=1;j<=x;j++){
int xx;
cin>>xx;
add(xx,f+i*2-1,1); //食物和牛
}
add(f+i*2-1,f+i*2,1);//牛和牛 拆点
for(int k=1;k<=y;k++){
int yy;
cin>>yy;
add(f+i*2,f+2*n+yy,1);//牛和饮料
}
}
for(int i=1;i<=d;i++){
add(f+2*n+i,tp,1);//饮料喝汇点
}
int ans=max_flow(sp,tp);
cout<<ans<<endl;
}