这题太经典,可以作为所有最大流算法的入门题...所以昨天搞了好久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;
}