POJ1459多源最大流

普通最大流详解

多源最大流的原理和普通最大流相同，只要建一个虚点S和一个虚点T，点S与所有的源点相连，点T与所有的汇点相连，即可。

题目传送门

本题的大意就是有np个源点，每个源点有限定的流出量；有nc个汇点，每个汇点有限定的流入量，求流入量的最大值。

（对于汇点u,如果流入量为in,流出量为out，则实际流入量为in-out）

我们只要将S与每一个源点连上等于源点流出量的边，T与每一个汇点连上等于限定流入量的边。做一遍DINIC即可。

#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
int S,T,h,t,ans,cnt;
int n,np,nc,m;
int dis[115],head[115],nxt[100005],q[500005];
struct node{
	int to,val;
}L[100005];
void add(int x,int y,int c){
	L[cnt].to=y;L[cnt].val=c; 
	nxt[cnt]=head[x];head[x]=cnt;cnt++;
	L[cnt].to=x;L[cnt].val=0;
	nxt[cnt]=head[y];head[y]=cnt;cnt++;
}
int BFS(){	
	memset(dis,0,sizeof(dis));
	dis[S]=1;h=0;t=1;q[t]=S;
	while(h<t){
		int front=q[++h];
		for(int i=head[front];i!=-1;i=nxt[i]){
			int to=L[i].to;
			if(!dis[to]&&L[i].val){
				q[++t]=to;
				dis[to]=dis[front]+1;
			}
		}
	}
	return dis[T];
}
int DFS(int now,int x){
	if(now==T) return x;
	int res=0;
	for(int i=head[now];i!=-1&&x;i=nxt[i]){
		int to=L[i].to;
		if(dis[to]==dis[now]+1&&L[i].val){
			int fd=DFS(to,min(x,L[i].val));
			L[i].val-=fd;L[i^1].val+=fd;
			res+=fd;x-=fd;
		} 
	}
	if(!res) dis[now]=0;
	return res;
}
int main()
{
	while(~scanf("%d%d%d%d",&n,&np,&nc,&m)){
		memset(head,-1,sizeof(head));
		S=0;T=n+1;cnt=0;ans=0;
		for(int i=1;i<=m;i++){
			int x,y,c;char t;
			cin>>t>>x>>t>>y>>t>>c;
			add(x+1,y+1,c);
		}
		for(int i=1;i<=np;i++){
			int x,c;char t;
			cin>>t>>x>>t>>c;
			add(S,x+1,c); 
		}
		for(int i=1;i<=nc;i++){
			int x,c;char t;
			cin>>t>>x>>t>>c; 
			add(x+1,T,c);
		}
		while(BFS()){
			ans+=DFS(S,2e9);
		}
		printf("%d\n",ans);
	}
	return 0;
}