1010 - 网络流 - dining（POJ 3281）

传送门

 

分析

网络流板子很好敲，难就难在建模上

不过ldw老师说建模学好后，对图论很有帮助

所以，还是要好好干网络流啊

这道题建4排点，第一排是食物，第二排和第三排都是牛，第四排是饮料，然后再多个源汇点

为什么第二三排都是牛呢？

因为要保证每头牛都只被匹配过一次

拆点就是一个非常常规的操作了

 

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#define in read()
#define M 4000000
#define inf 10000
#define N 500
using namespace std;
inline int read(){
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>='0'&&ch<='9'){
		res=(res<<3)+(res<<1)+ch-'0';
		ch=getchar();
	}
	return f==1?res:-res;
}
int n,f,d,S=0,T;
int nxt[M],to[M],cap[M],head[N],cur[N],cnt=1,lev[N];
void add(int x,int y,int z){
	nxt[++cnt]=head[x];head[x]=cnt;to[cnt]=y;cap[cnt]=z;
}
bool bfs(){
	memcpy(cur,head,sizeof(head));
	for(int i=S;i<=T;++i) lev[i]=-1;
	queue<int > q;
	q.push(S);lev[S]=0;
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int e=head[u];e;e=nxt[e]){
			int v=to[e];
			if(lev[v]!=-1||cap[e]<=0) continue;
			lev[v]=lev[u]+1;
			if(v==T) return true;
			q.push(v);			
		}
	}
	return false;
}
int dinic(int u,int flow){
	if(u==T) return flow;
	int delta,res=0;
	for(int &e=cur[u];e;e=nxt[e]){
		int v=to[e];
		if(lev[v]>lev[u]&&cap[e]>0){
			delta=dinic(v,min(cap[e],flow-res));
			if(delta){
				res+=delta;cap[e]-=delta;
				cap[e^1]+=delta;if(res==flow) return flow;
			}
		}
	}
	return res;
}
int main(){
	n=in;f=in;d=in;
	T=f+d+2*n+1;
	int i,j,k,a,b;
	for(i=1;i<=f;++i){	add(S,i,1);add(i,S,0);}//稍微画个图就可以了 
	for(i=1;i<=n;++i){ add(f+i,f+i+n,1);add(f+i+n,f+i,0);}
	for(i=1;i<=d;++i){	add(f+2*n+i,T,1);add(T,f+2*n+i,0);}
	for(i=1;i<=n;++i){
		a=in;b=in;
		for(j=1;j<=a;++j) {
			int t=in;
			add(t,i+f,1);add(i+f,t,0);
		}
		for(j=1;j<=b;++j){
			int t=in;
			add(i+f+n,f+2*n+t,1);add(f+2*n+t,i+f+n,0);
		}
	}
	int maxflow=0;
	while(bfs()) maxflow+=dinic(S,inf);
	printf("%d\n",maxflow);
	return 0;
}