网络流（最大流-Dinic算法）
摘自https://www.cnblogs.com/SYCstudio/p/7260613.html
网络流定义
最大流
正如可以通过将道路交通图模型化为有向图来找到从一个城市到另一个城市之间的最短路径，我们也可以将一个有向图看做是一个“流网络”并使用它来回答关于物料流动方面的问题。设想一种物料从产生它的源结点经过一个系统，流向消耗该物料的汇点这样一个过程。源结点以某种稳定的速率生成物料，汇点则以同样的速率消耗物料。从直观上看，物料在系统中任何一个点上的“流量”就是物料移动的速率。这种流网络可以用来建模很多实际问题，包括液体在管道中的流动、装配线上部件的流动、电网中电流的流动和通信网络中信息的流动。我们可以把流网络中每条有向边看做是物料的一个流通通道。每条通道有限定的容量，是物料流经该通道时的最大速率，如一条管道每小时可以流过200加仑的液体。流网络中的结点则是通道的连接点。除了源结点和终结点外，物料在其他结点上只是流过，并不积累或聚集。换句话说，物料进入一个结点速率必须与其离开该结点的速率相等。这个性质称为“流量守恒”，这里的流量守恒与Kirchhoff电流定律等价。在最大流问题中，我们希望在不违反任何容量限制的情况下，计算出从源结点运送物料到汇点的最大速率。这是与流网络有关的所有问题中最简单的问题之一(????).，这个问题可以由高效的算法解决。而且，最大流算法中的一些基本技巧可以用来解决其他网络流问题。————《算法导论》
Dinic算法
定义：
源点：只有流出去的点
汇点：只有流进来的点
流量：一条边上流过的流量
容量：一条边上可供流过的最大流量
残量：一条边上的容量-流量
几个基本性质
基本性质一：
对于任何一条流，总有流量<=容量
这是很显然的
基本性质二
对于任何一个不是源点或汇点的点u，总有
∑p∈Ek[p][u]==∑q∈Ek[u][q]（其中k[i][j]表示i到j的流量）∑p∈Ek[p][u]==∑q∈Ek[u][q]（其中k[i][j]表示i到j的流量）
这个也很显然，即一个点（除源点和汇点）的入流和出流相等
基本性质三
对于任何一条有向边(u,v)，总有
k[u][v]==−k[v][u]k[u][v]==−k[v][u]
这个看起来并不是很好理解，它的意思就是一条边的反边上的流是这条边的流的相反数，可以这么想，就是如果有k[u][v]的流从u流向v，也就相当于有-k[v][u]的流从v流向u。这条性质非常重要。
网络流最大流的求解（Ford-Fulkerson算法）
网络流的所有算法都是基于一种增广路的思想，其基本步骤如下：
1.找到一条从源点到汇点的路径，使得路径上任意一条边的残量>0（注意是大于而不是大于等于，这意味着这条边还可以分配流量），这条路径便称为增广路
2.找到这条路径上最小的F[u][v]（我们设F[u][v]表示u->v这条边上的残量即剩余流量），下面记为flow
3.将这条路径上的每一条有向边u->v的残量减去flow，同时对于起反向边v->u的残量加上flow
4.重复上述过程，直到找不出增广路，此时我们就找到了最大流
这个算法是基于增广路定理(Augmenting Path Theorem): 网络达到最大流当且仅当残留网络中没有增广路
在寻找增广路的过程中，如果我们没有建反边，那么有时候会发现找出来的路径与实际上的最大流不同，我们建这个反边是为了将原先得到的流推回去，从而得到新的流，就是给了我们反悔的机会。
这种朴素的FF算法在某种情况下可能非常低效，比如这个栗子
我们一眼可以看出最大流是999(s->v->t)+999(s->u->t)，但如果程序采取了不恰当的增广策略：s->v->u->t
我们发现中间会加一条u->v的边
而下一次增广时：
若选择了s->u->v->t
然后就变成
如果边上的容量更大，这个算法可能会更加低效，这时引入Dinic算法
Dinic算法
为Dinic算法引入了一个叫做分层图的概念。具体就是对于每一个点，我们根据从源点开始的bfs序列，为每一个点分配一个深度，然后我们进行若干遍dfs寻找增广路，每一次由u推出v必须保证v的深度必须是u的深度+1。
Dinic算法总结
例题：HDU-3549 Flow Problem http://acm.hdu.edu.cn/showproblem.php?pid=3549
1 #include<cstring> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cmath> 5 #include<queue> 6 using namespace std; 7 typedef long long ll; 8 typedef unsigned long long ull; 9 #define INF 0x3f3f3f3f 10 const ll MAXN = 200 + 7; 11 const ll MAXM = 1e3 + 7; 12 const ll MOD = 1e9 + 7; 13 const double pi = acos(-1); 14 int cnt = -1, head[MAXM], dis[MAXN], cur[MAXM]; 15 int n, m; 16 struct edge 17 { 18 int to, value, net; 19 } e[MAXM << 1]; ///共有n*2条边 20 void add(int from, int to, int value) 21 { ///链式前向星 22 cnt++; 23 e[cnt].to = to; 24 e[cnt].value = value; 25 e[cnt].net = head[from]; 26 head[from] = cnt; 27 } 28 int bfs(int st, int ed) 29 { ///建立层次图 30 queue<int> que; 31 memset(dis, -1, sizeof(dis)); 32 dis[st] = 0; 33 que.push(st); 34 while (!que.empty()) 35 { 36 int x = que.front(); 37 que.pop(); 38 for (int i = head[x]; i > -1; i = e[i].net) 39 { 40 int now = e[i].to; 41 if (dis[now] == -1 && e[i].value) 42 { 43 que.push(now); 44 dis[now] = dis[x] + 1; 45 } 46 } 47 } 48 return dis[ed] != -1; 49 } 50 int dfs(int x, int t, int maxflow) 51 { 52 if (x == t) 53 return maxflow; 54 int ans = 0; 55 for (int i = cur[x]; i > -1; i = e[i].net) 56 { ///当前弧优化 57 int now = e[i].to; 58 if (dis[now] != dis[x] + 1 || e[i].value == 0 || ans >= maxflow) 59 continue; 60 cur[x] = i; 61 int f = dfs(now, t, min(e[i].value, maxflow - ans)); 62 e[i].value -= f; 63 e[i ^ 1].value += f; ///反向边加流量 64 ans += f; 65 } 66 if (!ans) 67 dis[x] = -1; ///炸点优化 68 return ans; 69 } 70 int Dinic(int st, int ed) 71 { 72 int ans = 0; 73 while (bfs(st, ed)) 74 { 75 memcpy(cur, head, sizeof(head)); 76 int k; 77 while ((k = dfs(st, ed, INF))) 78 ans += k; 79 } 80 return ans; 81 } 82 void init() 83 { 84 cnt=-1; 85 memset(head,-1,sizeof(head)); 86 } 87 int main() 88 { 89 int t; 90 scanf("%d", &t); 91 for (int i = 1; i <= t; i++) 92 { 93 init(); 94 scanf("%d%d", &n, &m); 95 while (m--) 96 { 97 int u, v, t; 98 scanf("%d%d%d", &u, &v, &t); 99 add(u, v, t); 100 add(v, u, 0); 101 } 102 printf("Case %d: %d\n", i, Dinic(1, n)); 103 } 104 return 0; 105 }
