poj1273Drainage Ditches最大流[dinic算法]

这题太经典,可以作为所有最大流算法的入门题...所以昨天搞了好久dinic,愣是没明白建图那一步.后来干脆抛开那些代码,自己根据对dinic的理解写.没想到还能做出来! 好惊喜...

唉,也说说自己的感受吧,网上关于这方面的文章什么的太多了,而且各有各的理解,昨晚好痛苦,看了无数篇博客文章,总算理解到算法最关键的部分了.也给自己留个教训:别总想着模板,好习惯什么的,甚至还担心如果建图甚么的没做好以后学其它算法也会很吃力.其实,既然是算法,那么最关键的是算法思想,是其最精髓的部分,至于怎么实现,完全可以自由发挥,顶多

吸取一下别人或效率或代码量或巧妙一些的实现而已,别瞎折腾,没什么用...........

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <queue>
using namespace std;

const int INF = INT_MAX>>1;//无穷大...
const int MAX_V = 256;
int G[MAX_V][MAX_V]; //邻接矩阵,存图用
int level[MAX_V]; //分层
int n, m;

bool BFS(int s, int t) {
	memset(level, -1, sizeof(level));//初始化为-1,因为分层后要求从0开始,
					//往上一层一层dfs,那层次应该是非负整数,
					//为了防止冲突,就的初始化为负值

	queue<int> Q; int p;

	Q.push(s);
	level[s] = 0;	//起点为0层,单独一层
	
	while (!Q.empty()) {
		p = Q.front(); Q.pop();
		for (int i = 0; i < n; ++i) {
			if (G[p][i] > 0 && level[i] < 0) {//level[i]<0和EK算法的!vis[i]是一样的其实.
				level[i] = level[p] + 1;
				Q.push(i);
				if (i == t) return true;//如果分层分到了终点,换句话说,从起点开始可以一层一层地走到终点了,显然就有路径了.
			}
		}
	}
	return false;
}

int dfs(int s, int t, int flow) {
	if (s == t) return flow;
	int tmp, res = 0;
	for (int i = 0; i < n; ++i) {
		if (G[s][i] > 0 && level[i] == level[s] + 1) {
			tmp = dfs(i, t, min(flow, G[s][i]));//记录残留量
			G[s][i] -= tmp;
			G[i][s] += tmp;//回退残留
			res += tmp; flow -= tmp;//记录.res会记录从当前层开始总的流量变化,作为返回值.而flow表示本层当前流量,每回退一份流量,当前流量自然就要减少一份,从而和上面min(flow,G[s][i])保持同步.其实这里不减,上面改成min(flow-res, G[s][i])也是一样一样的.
		}
	}
	return res;
}

int dinic(int s, int t) {
	int flow = 0;
	while (BFS(s, t)) { flow += dfs(s, t, INF); } //dinic精髓.
	return flow;
}

int main() {
	while (~scanf(" %d %d", &m, &n)) {
		memset(G, 0, sizeof(G));
		int s, t, c;
		/* 建图时我把下标(节点编号统一为0~n-1) */
		for (int i = 0; i < m; ++i) {
			scanf(" %d %d %d", &s, &t, &c);
			G[s-1][t-1] += c;
		}
		printf("%d\n", dinic(0, n-1));
	}
	return 0;
}