POJ 1459 Power Network (最大流, 网络流, EdmondsKarp算法求解最大流)

题意:  给定np个发电站, nc个消费者, m条从u流到v的边, 其权值为边的最大容量, 求解, 消费者总的最大可用电量

求解:  增加两个点:源点以及汇点,  源点指向所有的发电站(即给发电站提供电量， 但是每个发电站最大接收量已被限制), 所有的消费者指向汇点(即消费者所能用的电量全部汇聚到该汇点上,每个消费者能用的电量有最大限制), 构建网络流模型， 用EdmondsKarp算法求解

code:

/*
   s(u)  u被提供的能量数
   p(u)  0<=p(u)<=Pmax(u) u产生的能量数
   c(u)  u消耗的能量数  (0<=c(u)<=minn(s(u), Cmax(u)))
   d(u)  u向外传递的能量数 = s(u)+p(u)-c(u)
   每个发电站         c(u) = 0;
   每个消费者(用电者)  p(u)=0;
   每个调度员   c(u)=0 && p(u)=0;
   u到v至多一条电力通路
   l(u, v) u传送到v的总能量数
   Con = EuC(u), 消耗总和
   求Con最大值
 */
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int maxn = 205;
const int inf = 0x3f3f3f3f3f;
int n, np, nc, m, s, e;//s->源点,  e->汇点
//n:nodes
//np: power stations 发电站数
//nc: consumers  消费者数
//m: 电力通路数
int edge[maxn][maxn];// u-v之间电力通路传送最大值
int flow[maxn], vis[maxn];
int pre[maxn];
//找出所有增广路, 每找出一条就在各个边上去掉汇点上所获得的能量
int bfs(){
	int i, v;
	queue<int>que;
	memset(flow, -1, sizeof(flow));
	memset(vis, 0, sizeof(vis));
	while(!que.empty()) que.pop();
	que.push(s);
	vis[s] = 1;
	flow[s] = inf;
	pre[s] = s;
	while(!que.empty()){
		v = que.front();
		que.pop();
		for(i=s; i<=e; i++){
			if(vis[i]&&edge[v][i]){
				vis[i] = 1;
				flow[i] = min(flow[v], edge[v][i]);
				que.push(i);
				pre[i] = v;
			}
		}
	}
	if(flow[e]!=-1){
		return flow[e];
	} else {
		return 0;
	}
}
int EdmondsKarp(){
	int i, MIN;
	int maxflow = 0;
	while(MIN=bfs()){
		for(i=e; i!=s; i=pre[i]){
			edge[pre[i]][i] -= MIN;//正向边
			edge[i][pre[i]] += MIN;//反向边
		}
		maxflow += MIN;
	}
	return maxflow;
}
int main()
{
	int i, j, k;
	char c;
	while(~scanf("%d %d %d %d%*c", &n, &np, &nc, &m)){
		s=0, e=n+1;
		memset(edge, 0, sizeof(edge));
		int x, y, z;
		for(i=1; i<=m; i++){
			while((c=getchar())!='(');
			scanf("%d,%d)%d", &x, &y, &z);
			if(x==y) continue;
			edge[x+1][y+1]+=z;
		}
		for(i=1; i<=np; i++){
			while((c=getchar())!='(');
			scanf("%d)%d", &x, &z);
			edge[s][x+1] = z;//源点与每个发电站建立联系
		}
		for(i=1; i<=nc; i++){
			while((c=getchar())!='(');
			scanf("%d)%d", &x, &z);
			edge[x+1][e] = z;//汇点与每个消费者建立联系
		}
		printf("%d\n", EdmondsKarp());
	}
	return 0;
}