网络流简述
网络流主要可以拿来解决一些跟有向关系相关的问题,例如液体在管道中的流动、货物的运载、网络中的信息波动等。
简单介绍一下它:在一个有向图上选择一个源点s、一个汇点t。源点只流出,汇点只流进。同时,一条边经过的流量记为,也有允许通过的最大流量称为容量,记为。(若该边不存在,则)。除了源点和汇点以外的每个店入流和出流都相等。一条边上的剩余流量(没有用完的)称为残量,即 容量-流量。
因此网络流模型可以形象地描述为:在每一条边都不超过容量限制的前提下,“流”从源点源源不断地产生,最终全部归于汇点。
基本性质
- (容量限制)
- 对于任何一个不是源点或汇点的点,总有
因为入流和出流相等(流量平衡) - 对于任何一条有向边,总有 (斜对称性)
最大流
先讲一下比较简单的最大流问题。
直接从字面意思即可理解,从源点流到汇点的流量最大就是最大流。
算法思想:从零流(所有边的流量均为0)开始不断增加流量,保持每次增加流量后都满足以上性质(增加流量,也就是减去容量时要相应的给反向边加上等量的容量,便于反悔)。计算每条边上的残量,得到残量网络,再继续在残量网络中尝试增加流量。
增广路算法基于:残量网络中任何一条从到的有向道路都对应一条原图中的增广路——只要求出该道路中所有容量的最小值,把对应的所有边上的流量减去,在答案里加上它即可,这个过程被称为增广。
显然只要残量网络中存在增广路,流量就可以增大;反之,如果不存在增广路,流量就已经最大。(增广路定理)
顺便也讲一下最大流最小割定理。
最小割
- 割:对于一个网络流图,其割的定义为一种点的划分方式:将所有的点划分为和两个集合,其中源点,汇点。
- 割的容量:定义割的容量表示所有从到的边的容量之和,即。当然也可以用表示割的容量。
- 最小割:使得容量最小的割。也可以理解为使得和不联通所需要删去边权最小的割。
最大流最小割
- 定理:
- 证明:对于任意一个可行流和任意割,有:。而当达到最大流时,残量网络中不存在从到的增广路,所以的出边都是满流,上式等号成立。同时,上式的另一表达形式为。又因等号可以取到,则,证毕。
问题模型
一般在最小割的问题中,割掉一条边表示选择某个条件,而边权表示的是选择这个条件的价值。至于该价值的贡献是好是坏,则根据你如何使用这个价值而定。
举个例子,对于二分图,有点集。表示u、v中至少选择一个,保证无孤立点,求最小点权覆盖集。
设左边为点集,右边为点集,如图所示连边。割掉红色边表示选择中某个点,其点权设为;割掉蓝色边表示选择该边两端的点,其点权设为;割掉橙色边表示选择中某个点,边权设为。这样一来,若一条边不从连通到,则表示至少有一个点被选中,求得的最小割即为最小点权之和,覆盖集可通过边的使用情况求得。
事实上,我们都清楚,为了使点权之和最小化,蓝色的边是不会被割掉的。而仍然要设置它的原因,一是为了结合问题背景来建模,二是让它承担通道的角色,保证网络流图的性质。问题建模是非常值得琢磨的。
如果转换一下,边的限制变成u、v中至多选一个,求最大点权独立集。其本质相同,因为至多选一个等价于至少删去一个。由此可以看出,最大点权独立集为最小点权覆盖集的补集,其解法自然也就不言而喻了。
这里只是提供了一个参考的思路,具体的问题还要具体分析。
最小割的边数
将边权设为1重新求一遍最小割即可。
最小割集的求解
由于使用网络流解决此类问题效率较为低下,普遍使用算法进行求解,有兴趣可自行查阅。
一些小细节
-
存图的时候要从偶数开始存,同时存正向边和反向边(这样就可以保证正向边编号全为偶数,反向边编号为 i^1(奇偶性相反))
-
反向边怎么用?因为找到的增广路不一定是最优的,反边给你“反悔”的机会。如果一条边边权为0,那么往回走的时候并不会对流量有所影响(走不回去)。所以一开始反边的边权应该存0。当正向边减去流过这条边的增广路上容量最小值d(此时这个值已经被加入到了答案)的时候,反向边应该加上d(因为对于源点和汇点来说中间流量的这些变化都是无差别的,为了保证反向正向相加得原边权,也就是不改变原本的条件就得这么做)
Edmonds-Karp(EK) 便是不断用来寻找增广路,直到图中不存在增广路的算法。
但是如果一条一条地找出增广路,万一有一些极(毒)端(瘤)数据(比如几条相邻的边容量相差特别大),这个时间复杂度就是无法承受的。Dinic的高效之处在于它能够同时找出几条增广路.
关于最大流,我还没有讲完!当然,实际上在大部分情况下以上两个算法已经够用(我认为)。可以先跳过剩下有关最大流的算法。
(这里应该有ISAP和HLPP)
最大流/最小割 练习题
费用流
假如流经一条边有对应流量的花费,那么问题就可以有更多的变式了。
- 定义一条边的费用表示单位流量流经所需花费的费用。即,当边的流量为时,需要花费的费用。
类似的,我们先从最小费用最大流引入。即在最大化流量的基础上使得总花费最小。
这当然和上面的最大流相关。之前寻找增广路的方式是,对于来说,就是在边权为1的图上找到的一条最短路。那么,将 更换为寻找最短路的算法,边权设为,会产生什么样的效果呢?
这样的方式,既保证了最短路的性质,也保证了增广路的性质。由于每单位流到汇点的流量都要花费途经的边权之和,则最短路的性质使得了这些每次从最短路流过的流量是最划算的。在这样的情境下,反向边相当于退钱,所以要将边权设为。
而 是同时找到多条增广路的算法,只需将 更改为最短路算法的同时,限制流量只能由当前点流向到汇点的最短路上的点即可。
因为有负权边的存在,所以不能直接采用 ,应该采用 或经由 处理的 。同时,由于向下走的条件发生改变,有可能会往回走,所以需要标记当前链上的点,防止陷入无限循环。
Reference
- EK不够快? 再学个Dinic吧 by 钱逸凡
- Dinic算法 by SYCstudio
- 从入门到精通:最小费用流的“zkw”算法 by zkw
- 基于 Capacity Scaling 的弱多项式复杂度最小费用流算法 by ouuan
- OI Wiki
- 《算法竞赛入门经典第二版》
- 《算法竞赛进阶指南》
- 《算法导论》