浅谈网络最大流

网络最大流是一种比较常见的高级知识点吧。
别的概念我就先不解释了，直接简单的说一下大概是什么。
其实可以这么说：已知一个图，每一条边都有一个流量限制，换句话讲就是最多只能有那么多的流量经过，请问从源点开始最多能有多少单位流量流到汇点（其实就是指定点）。

就比如说下图最大流是2。

下面说一说我们怎么实现最大流。

Edmonds-Karp

其实本质很简单：就是找源点到汇点的路（路径上每一个边的流量都要大于0），然后减去相应的流量。

就这么简单？？？

想多了，有些情况下单单这样做是错的，就像这样：

像这样最大流就只有1了，但正确答案是2，怎么办呢？

建反向边

对跑过的边建立反向边，流量为当前这一条路流的流量，像这样：

然后再跑一次就可以啦。
可以证明，这样做是对的。

End？

等等……如果是这样，求点1到点5的最大流。

运气不好的话会跑1998次。
像这样：

1->4->3->5
1->2->3->4->5
1->4->3->5
1->2->3->4->5
1->4->3->5
1->2->3->4->5
......

其实这个问题很好解决，每一次只要找最短的路径就可以了，其实用BFS来寻找路径就可以啦。

哦对，该算法的时间复杂度是$O(n{m}^2)$，不过通常情况下跑不满，所以通常来说不会超时，嗯嗯，的确结束了。

当然，这只是网络流的一个基本算法，挺好实现的，代码就不给了。

Dinic

可以说是网络流中最常用的一种求最大流的算法。

分为两部：分层和增广（增广就是寻找路径）

分层

这一个还是比较容易的，就是算出每一个点到源点的最短距离（经过的边的流量必须大于0），就像这样：

然后按照距离用DFS来增广，这样子可以保证走的是最短的增广路

然后发现没有增广路了，所以就结束了，最后的最大流就是1998。
说一下关于反向边的问题：
因为我们可以让每一对边的编号长这样：
(2,3) (4,5) (6,7)……
这样子我们找反向边的时候直接异或1就可以了。

以LOJ101最大流为例

#include <cstdio>
#include <algorithm>
#define to k[i].ar
#define ll long long
#define maxint (ll)0x3f3f3f3f3f3f3
#define N 10010
using namespace std;
struct node {
    int next, ar, st;
    long long flow;
} k[N * 100];
int first[N], dis[N], len = 1;//编号从2开始
int head, tail, d[N], x, mx;
void add(int a, int b, int t) {
    len++;
    k[len].ar = b;
    k[len].next = first[a];
    first[a] = len;
    k[len].flow = t;
    k[len].st = a;
	//建图的时候记得先建好反向边
    len++;
    k[len].ar = a;
    k[len].next = first[b];
    first[b] = len;
    k[len].st = b;
}
int a, b, t;
//BFS给图分层
bool bfs() {
    head = 0, tail = 1;
    d[0] = a;
    for (register int i = 1; i <= mx; i++) dis[i] = 0;
    dis[a] = 1;
    while (head < tail) {
        t = d[head];
        head++;
        for (register int i = first[t]; i; i = k[i].next) {
            if (dis[to] == 0 && k[i].flow > 0) {
                dis[to] = dis[t] + 1;
                if (to == b)
                    return true;
                d[tail] = to;
                tail++;
            }
        }
    }
    //没有找到增广路，Byebye
    return false;
}
//DFS增广
ll dfs(int xx, ll flow) {
	//到终点了的话就返回当前流量
    if (xx == b)
        return flow;
    ll h, s = 0;
    for (register int i = first[xx]; i; i = k[i].next) {
        if (flow == 0) {
            dis[xx] = -1;
            return s;
        }
        if (dis[to] == dis[xx] + 1 && k[i].flow > 0) {
            h = dfs(to, min(flow, k[i].flow));
            s += h;
            flow -= h;
            k[i].flow -= h;
            k[i ^ 1].flow += h;
        }
    }
    if (s == 0)
        dis[xx] = -1;
    return s;
}
int n, m, u, v, w;
ll ans;
//Dinic主体
void dinic() {
    while (bfs()) {
        ans += dfs(a, maxint);
    }
}
int main() {
    scanf("%d%d%d%d", &n, &m, &a, &b);
    mx = n;
    for (int i = 1; i <= m; i++) {
        scanf("%d%d%d", &u, &v, &w);
        add(u, v, w);
    }
    dinic();
    printf("%lld", ans);
}

然后最大流就讲完了