图论-zkw费用流
这是一个求最小费用最大流的算法,因为发明者是神仙zkw,所以叫zkw费用流(就是zkw线段树那个zkw)。有些时候比EK快,有些时候慢一些,没有比普通费用流算法更难,所以学zkw费用流之前,不需要先掌握普通费用流。
前置知识: 网络最大流 \color{8888cc}\texttt{网络最大流} 网络最大流。
在学了网络最大流后,如果在没条边上加个限制,就是 c o s t cost cost,表示这条边上每走过 1 1 1 流量就要付费 c o s t cost cost,求在最大流的情况下,要交的最少费用。
如上费用流图,最大流为 1 1 1,在达到最大流的同时,最小费用为 3 3 3。最小费用路径: s → 1 → 2 → t s\to1\to2\to t s→1→2→t, 1 × ( 2 + 0 + 1 ) = 3 1\times (2+0+1)=3 1×(2+0+1)=3。
像这种图,普通费用流算法更快,因为它一次只增广一条路径。但一些大的费用流图往往最大流要经过好多路径,那么zkw费用流快得多。那么为什么还要有EK的存在呢?因为有些毒瘤题专门卡zkw费用流(路径长,最大流路径单一)。
接下来开始讲zkw费用流的实现。 首先凡是网络流题,都是要建双向边的,要不然就“走不了回头路”了。对于费用流,走回头路可以把钱要回来,所以建边如下:
class Graph{
public:
int top,to[M<<1],fw[M<<1],ct[M<<1];
//top表示边数,因为后面要把互为反边的边通过^1获得,所以正边的编号要为偶数,然后负边的编号要为正边+1
//to表示这条边通往的节点,fw表示边的流量,ct表示边的费用
vector<int> g[V];//比较奇怪的建边方式
Graph(){
top=1;}
void add(int x,int y,int f,int c){
g[x].push_back(++top);
to[top]=y,fw[top]=f,ct[top]=c;
}
void Add(int x,int y,int f,int c){
add(x,y,f,c),add(y,x,0,-c);
//如同最大流,建反边,钱可以走反边要回来,所以反边费用为-c
}
};
然后开始做最小费用最大流。先看看求出答案的框架:
void mcmf(){
while(spfa()){
vis[t