最小费用最大流的实现之Dijkstra+Ford-fulksonff算法

 

文章写的很详细，对于学习最小费用最大流的求解算法有很大的启发作用

原帖地址附上http://www.cppblog.com/guojingjia2006/archive/2009/11/12/57905.html还要感谢博主啊

最小费用最大流 修改的dijkstra  +  Ford - Fulksonff算法
修改的dijkstra其实和Johnson算法的思想是一致的。
 
一个求最小费用最大流的朴素算法是这样的：
1  求最小费用增广路
2  判断是否存在增广路，否的话算法终止。
3  增加增广路上边的流量
4  在增广路上添加必要的逆向负权边
5   goto   1
 
因为负权边的存在，求最小费用增广路就不可以用dijkstra算法。当然，我们可以用bellman - ford算法，可是这样的话求一次最短路的时间代价就是O(e * n)，e是边数，n是顶点数。代价大了点，如果能用dijkstra算法就好了。利用Johnson算法的思想，这是可以做到的。
 
第一次求最短路可以用dijkstra算法（如果一开始就有负权边，那就用bellman - ford算法，这没关系），求出源点到所有点的距离，嗯，我说的距离是指路径上边的费用之和的最小值。注意，要求出到所有点的距离，而不是求出到汇点的距离就完事了。
假设有一条边u -> v，源点到u的距离是d[u]，到v的距离是d[v]，边的费用（权值）是w(u，v)。很显然，d[u] + w(u,v) >= d[v]，不然的话，你会发现一条更好的路径从源点到v。问题是，什么时候取等呢？当u -> v在v的最优路径上，范围说小一点，当u -> v在从源点到汇点的最优路径，即最小费用增广路上。
好的，如果u -> v被你增载了，你要开始添负权边v -> u了，权值取负，就是 - w(u,v)。负权就是讨厌，是正的就好了，dijkstra算法就可以再用了。怎么办呢，把负权边加个权值，让它非负。要加多少呢，d[v] - d[u]。当然不能只加一条边，对所有边，无论原有的还是新添的，按这个规则加，构造一个新的图：
                对边a -> b，新的边权w ' (a,b)=w(a,b)+d[a]-d[b]
现在来看看你的杰作：
                对原来的边u -> v, w ' (u,v)=w(u,v)+d[u]-d[v]: 记得么d[u]+w(u,v)>=d[v], 所以 w ' (u,v)  >=   0
                对新加的负权边v -> u, w ' (v,u)=w(v,u)+d[v]-d[u]=-w(u,v)+d[v]-d[u]: 记得么d[u]+w(u,v)==d[v]，这里可是取等号的，所以w ' (v,u)  ==   0
哈哈，这下所有边又是非负的了。
可是，问题是，为啥不每个边加个足够大的正数，这样不是所有边也都是正的了么。仔细想想，边权为啥要为正，不就是为了求源点到汇点的最短路方便么，可是，都加大正数的话，你求出的最短路和原来图的最短路能一致么，不能，为啥，画个三角形，自己想想。可是，我的方法就能一致么，能。我证明给你看。
 
假设从源点s到汇点t有一条路径s -> a -> b -> c -> d. -> t,在原图中的路径长度为
                 w(s,a) + w(a,b) + w(b,c) + + w(x,t)
在新图中的路径为
                 w ' (s,a)+w ' (a,b) + w ' (b,c)+w ' (x,t)
展开来就是
                 w(s,a) + d[a] - d[s] + w(a,b) + d[b] - d[a] + w(c,d) + d[d] - d[b] + . + w(x,t) + d[t] - d[x]
消阿消，d[a]和 - d[a]，d[b]和 - d[b]d[x]和 - d[x]，剩下什么呢：
                 w(s,a) + w(a,b) + w(b,c) + + w(x,t) + d[t] - d[s]
噢，不就比原图中多d[t] - d[s]么（其实d[s] == 0 ）。这可是对所有s到t的路径都成立的，既然所有路径，在新图中的权值都比在原图中的权值多了d[t]，那么，新图的最短路，也就对应原图的最短路，只不过路径长度多了d[t]，这不仅对t成立，对所有节点u都成立，只不过新图中到u的最短路长度比原图多了d[u]。
 
好，用dijkstra算法，第二次求出最短路。然后求出新的d’[u]，然后添加新的边，然后准备第三次的dijkstra算法。。。为什么第二次可以这样做，第三次还可以这样做，第三次的原图可能有很多负权边啊？我可没说过w(u,v) >= 0这样的限制，所以，即使原图有负权边还是可以这样做的。
 
好了，第一次dijkstra算法（或者bellman - ford算法，如果有负权边的话，只用一次，不会成为瓶颈的），然后每次求最小增广路用一次修改的dijkstra算法。这个算法求最小费用最大流复杂度是O(m * n * n), m是最大流量，或者是求增广路次数的上界。最后，如果用这个算法来求最优匹配问题，复杂度是O(n ^ 3 )的。