网络最大流是一种比较常见的高级知识点吧。
别的概念我就先不解释了,直接简单的说一下大概是什么。
其实可以这么说:已知一个图,每一条边都有一个流量限制,换句话讲就是最多只能有那么多的流量经过,请问从源点开始最多能有多少单位流量流到汇点(其实就是指定点)。
就比如说下图最大流是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 ) O(nm^2) O(nm2),不过通常情况下跑不满,所以通常来说不会超时,嗯嗯,的确结束了。
当然,这只是网络流的一个基本算法,挺好实现的,代码就不给了。
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);
}
然后最大流就讲完了