網絡流_Edmonds-Karp 增廣路算法

EK增廣路算法的簡單思路:
我們通過bfs去尋找增廣路, 每找到一條,就去更新maxflow, 以及這條路徑上的每條邊;
注意點:使用“成對”存儲的方法。

測試題目 ,洛谷P2470
題目:
題目背景
在農夫約翰的農場上,每逢下雨,貝茜最喜歡的三葉草地就積聚了一潭水。這意味着草地被水淹沒了,並且小草要繼續生長還要花相當長一段時間。因此,農夫約翰修建了一套排水系統來使貝茜的草地免除被大水淹沒的煩惱(不用擔心,雨水會流向附近的一條小溪)。作爲一名一流的技師,農夫約翰已經在每條排水溝的一端安上了控制器,這樣他可以控制流入排水溝的水流量。

題目描述
農夫約翰知道每一條排水溝每分鐘可以流過的水量,和排水系統的準確佈局(起點爲水潭而終點爲小溪的一張網)。需要注意的是,有些時候從一處到另一處不只有一條排水溝。

根據這些信息,計算從水潭排水到小溪的最大流量。對於給出的每條排水溝,雨水只能沿着一個方向流動,注意可能會出現雨水環形流動的情形。

輸入格式
第1行: 兩個用空格分開的整數N (0 <= N <= 200) 和 M (2 <= M <= 200)。N是農夫John已經挖好的排水溝的數量,M是排水溝交叉點的數量。交點1是水潭,交點M是小溪。

第二行到第N+1行: 每行有三個整數,Si, Ei, 和 Ci。Si 和 Ei (1 <= Si, Ei <= M) 指明排水溝兩端的交點,雨水從Si 流向Ei。Ci (0 <= Ci <= 10,000,000)是這條排水溝的最大容量。

輸出格式
輸出一個整數,即排水的最大流量。

輸入輸出樣例
輸入 #1 複製
5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10
輸出 #1 複製
50

AC代碼:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

const int maxn = 205;
const int inf = 0x3f3f3f3f;

int n, m, maxflow, s, t, idx;
int h[maxn], e[maxn * 2], w[maxn * 2], ne[maxn * 2], incf[maxn], pre[maxn];
bool vis[maxn];

inline void add(int a, int b, int c) {
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
	e[idx] = a, w[idx] = 0, ne[idx] = h[b], h[b] = idx ++;
}

inline bool bfs(void) {
	memset(vis, false, sizeof vis);
	queue<int> q;
	q.push(s);
	vis[s] = true;
	incf[s] = inf;
	
	while(q.size()) {
		int x = q.front(); q.pop();
		
		for(int i = h[x]; i != -1; i = ne[i]) {
			if(w[i]) {
				int y = e[i];
				if(!vis[y]) {
					pre[y] = i;
					incf[y] = min(incf[x], w[i]);
					q.push(y);
					vis[y] = true;
					if(y == t) return true;
				}
			}
		}
	}
	
	return  false;
}

inline void update(void) {
	int x = t;
	while(x != s) {
		int i = pre[x];
		w[i] -= incf[t];
		w[i ^ 1] += incf[t];
		x = e[i ^ 1];
	}
	
	maxflow += incf[t];
}

int main(void) {
//	freopen("in.txt", "r", stdin);
	
	scanf("%d%d", &m, &n);
	memset(h, -1, sizeof h);
	s = 1, t = n;
	for(int i = 1; i <= m; i ++) {
		int a, b, c;
		scanf("%d%d%d", &a, &b, &c);
		add(a, b, c);
	}
	
	while(bfs()) update();
	
	printf("%d\n", maxflow);
	
//	fclose(stdin);
	return 0;
}

注意點:
因爲我們在建立邊的時候每次都需要建立2條邊,所以在開數組的時候需要 maxn * 2。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章