网络流（EK）

网络流Edmonds_Karp (模板)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#define mem(a) memset(a,0,sizeof(a))
#define inf 999999999
using namespace std;
int c[210][210],f[210][210],p[210];
//f[u][v]为<u,v>流量，c[u][v]为<u,v>容量，p[i]记录i的前驱 
int s,t;//s为源点，t为汇点 
int res[210];//res[i]表示源点s到节点i的路径上的最小残留量 
int sum;//记录最大流 
int m,n;
void init(){
	s=1;t=m;
	sum=0;
	mem(c);
	mem(f);
	mem(p);
}
void Edmonds_Karp(){
	queue <int> q;
	int u;
	while(1){
		memset(res,0,sizeof(res));//每找一次，初始化一次 
		res[s]=inf; 
		q.push(s);
		while(!q.empty()){
			u=q.front();
			q.pop();
			for(int i=1;i<=m;i++){
				if(!res[i] && f[u][i]<c[u][i]){
					p[i]=u;
					q.push(i);
					res[i]=min(res[u],c[u][i]-f[u][i]);//u-i路径上的最小残量
				}
			}
		}
		if(!res[t])break;//找不到增广路,则当前流已经是最大流 
		sum+=res[t];
		for(int i=t;i!=s;i=p[i]){//从汇点顺着这条增广路往回走
			f[p[i]][i]+=res[t];//更新正向流量
			f[i][p[i]]-=res[t];//更新反向流量
		}
	}
	printf("%d\n",sum);
	return ;
}
int main(){
	while(scanf("%d%d",&n,&m)!=EOF){
		init();
		while(n--){
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			c[x][y]+=z;//可能出现相同边 
		}
		Edmonds_Karp();
	}
	return 0;
}