网络流初步 增广路算法求最大流 hdoj3549

内容详见刘汝佳的《算法竞赛入门经典》，这里主要结合题目解释算法。

Edmond-Karp算法思想：

1.从起点出发利用bfs寻找一条路到终点，之后从终点回溯每条路，为每条路增加终点增量；

2.重复1直到终点增量为0。

Edmond-Karp算法用bfs代替了Ford-Fulkerson算法中的dfs，时间复杂度为O(nm^2)（n是端点数，m是边数）。

hdoj 3549： http://acm.hdu.edu.cn/showproblem.php?pid=3549

AC代码：

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int M = 20;
const int INF = 0xfffffff;
int d[M][M], p[M], add[M], f[M][M];//d是通道容量上限，f是通道使用容量，add是本次通过i点的容量
queue<int> q;
main()
{
    int t, n, m;
    scanf("%d", &t);
    for(int j = 1; j <= t; j++)
    {
        int ans = 0;
        while(!q.empty())
            q.pop();
        memset(d, 0, sizeof(d));
        memset(f, 0, sizeof(f));
        scanf("%d %d", &n, &m);
        while(m--)
        {
            int a, b, c;
            scanf("%d %d %d", &a, &b, &c);
            d[a][b] += c;//注意这题可能有重边的情况出现
        }
        for(;;)
        {
            memset(add, 0, sizeof(add));
            add[1] = INF; q.push(1);
            while(!q.empty())
            {
                int u = q.front(); q.pop();
                for(int i = 1; i <= n; i++)
                    if( !add[i] && d[u][i] > f[u][i])
                    {
                        p[i] = u; q.push(i);
                        add[i] = min(add[u], d[u][i] - f[u][i]);
                    }
            }
            if(add[n] == 0)//终点增量为0则退出
                break;
            for(int i = n; i != 1; i = p[i])//回溯并更新正向容量与反向容量
            {
                f[p[i]][i] += add[n];
                f[i][p[i]] -= add[n];
            }
            ans += add[n];
        }
        printf("Case %d: %d\n", j, ans);
    }
}

关于这个算法有个【反向容量】，在《算法竞赛》中并没有对它的工作原理作出具体解释，这里记一下：

当出现“走错”的路线时，可以通过反向容量将流量“向回推“。

如图，按照最短路就会走1-2-4-6，结果为10，但是很显然这个图的路线应该为1-2-5-6和1-3-4-6，结果为20。

当走过1-2-4-6后，1-2、 2-4、 4-6会形成权值为10的反向边，在1-3-4-2-5-6中会将2-6的反向边"推回去"，这样反向边就只存在于1-3-4-6和1-2-5-6，这样根据残余网络就可以得出最大流为20。