网络流初步——最大流&费用流

网络流初步——最大流&最小割&费用流

一、基础定义

网络 ：一个网络G=(V,E)是一张有向图；
容量 ：图中每条边(u,v)∈E都有一个给定的权值c(u,v)被称为容量
源/汇点：两个指定的特殊节点,分别设为S/T
流函数/流量/剩余容量 ：设f(x,y)中x、y均为图上节点，且满足：

• f(x,y)<=c(x,y)
• f(x,y)=-f(y,x)
• 对于除源/汇点外的任意节点，∑f(u,x)=∑f(x,v)( (u,x),(x,v)∈E )

则称f(x,y)为网络的 流函数 。对于(x,y)∈E，f(x,y)被称为边的 流量，c(x,y)-f(x,y)被称为边的 剩余容量 。网络中由源点出发的每条路径的流量总和被称为被称为 整个网络的流量

流函数性质： 上述三条性质，分别被称为 容量限制 ，斜对称 与 流量守恒 。尤其是流量守恒定律，它告诉我们网络中除源汇点之外，所有节点都不存储流，流出总量等于流入总量，就像“流”从源点源源不断地产生，在不超过容量限制的前提下不断地流经整个网络，最终全部进入汇点。

注意：每一条边的反向边都有一个负的流量

---

二、最大流

1.定义： 使得某一张网络的流量最大的流函数被称为最大流，此时网络的流量被称为 最大流量
2.算法： 种类很多，下面主要介绍Edmonds_Karp增广路算法 与 Dinic算法

Edmonds_Karp增广路算法

主要思想： 若存在一条路径，其最小剩余容量不为0,则称该路径为增广路。EK算法就是通过BFS不断找到增广路，直到网络中不再有增广路。显然，此时的流函数就是最大流。
流程：
在每轮寻找增广路时，EK算法仅考虑剩余容量不为0的边，找到任意一条由S到T的路径，该条路径上的最小剩余容量（设为minf)就是网络可以增加的流量
需要注意的是，当一条边的流量增加时，根据斜对称性质，其反向边流量f(y,x)<0，此时必有f(y,x)<c(y,x)。故EK算法还需要考虑每条边的反向边。
具体实现时，可以利用“成对存储”的技巧，边权的定义可以直接规定为剩余容量。
复杂度： O(nm² ),实际上远远达不到这个上界。
一般能处理10³ ~10⁴ 级别的网络

Dinic算法

在任意时刻，网络中所有点以及剩余容量大于0的边构成的子图被称为 残量网络。
S到任意点的最短距离被称为这个节点的层次，在残量网络中，满足d[y]=d[x]+1的边(x,y)构成的子图被称为分层图。分层图显然是DAG。
Dinic算法不断重复以下步骤，直到在残量网络中S不能到达T:

• 1.在残量网络上BFS求出节点的层次，构造分层图；
• 2.在分层图上dfs寻找增广路，在回溯时实时更新剩余容量

一点补充

你可能会有这样的疑问：找到一条增广路就更新，怎么就最优了呢？难道不会出现一条增广路，使得选择这条增广路比不选择这条增广路更劣吗？
事实上，在《信息学奥赛一本通》系列的《金牌导航》中就有关于这个问题的解释。答案是：“见到一条增广路就增广”的策略不能保证每次都是最优选择，但是其中隐含的退流思想使得在算法执行完毕时得到的最终结果是最优的。
什么叫 退流思想 ？事实上，我们可以将上面两个算法视为 有反悔机制的贪心，而退流思想扮演着 “反悔机制” 的角色。在一条边(u,v)被增广后，斜对称性质保证了其反向边(v,u)一定可以被增广，而如果其反向边(v,u)也被增广，则等效于先前(u,v)增加的一部分流量被”退回去”了。于是这个算法的正确性得到了保证。如果需要，在《金牌导航》中你会找到更为详细的介绍。

---

代码

洛谷P3376

#include<bits/stdc++.h>
using namespace std;
#define il inline
#define re register
typedef long long LL;
il int read()
{
   
	int s=0,w=1;char c=getchar();
	while(c<'0'||c>'9'){
    if(c=='-') w=-1;c=getchar();}
	while(c>='0'&&c<='9'){
    s=(s<<1)+(s<<3)+c-'0';c=getchar();}
	return s*w;
}
il LL read_ll()
{
   
	LL s=0,w=1;char c=getchar();
	while(c<'0'||c>'9'){
    if(c=='-') w=-1;c=getchar();}
	while(c>='0'&&c<='9'){
    s=(s<<1)+(s<<3)+c-'0';c=getchar();}
	return s*w;
}
const int N=210;
const int M=1e4+10;
const LL MAX=