最小费用最大流 zkw算法

2024-01-07 10:49
文章标签 算法 最小 最大 费用 zkw

本文主要是介绍最小费用最大流 zkw算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最小费用最大流的算法有很多种,经典的就是EK+SPFA。但是时间效率确实令人不敢恭维~

最近小小的学习了一下传说中的费用流的神算法——zkw算法。(P.s. 以其发明者zkw神牛的名字命名)

 

该算法可以有效地通过修改距离标号的方法避免反复进行最短路运算的部分省去,基本算法类似KM算法的顶标修改。基本思想就是:采用距离标号法,若能增广则无限增广,反之更新距离标号,直到无法更新为止。

本人的嘴比较笨,这里直接引用zkw神牛的原文。

 

【以下为zkw原文:】

最小费用流在 OI 竞赛中应当算是比较偏门的内容,但是 NOI2008 中 employee 的突然出现确实让许多人包括 zkw 自己措手不及。可怜的 zkw 当时想出了最小费用流模型,可是他从来没有实现过,所以不敢写,此题 0 分。zkw 现在对费用流的心得是:虽然理论上难,但是写一个能 AC 题的费用流还算简单。先贴一个我写的 employee 程序:只有不到 70 行,费用流比最大流还好写~

 

 这里使用的是连续最短路算法。最短路算法?为什么程序里没有 SPFA?Dijkstra?且慢,先让我们回顾一下图论中最短路算法中的距离标号。定义 D_i 为点 i 的距离标号,任何一个最短路算法保证,算法结束时对任意指向顶点 i 、从顶点 j 出发的边满足 D_i /le D_j + c_{ij} (条件1),且对于每个 i 存在一个 j 使得等号成立(条件2)。换句话说,任何一个满足以上两个条件的算法都可以叫做最短路,而不仅仅是 SPFA、Dijkstra,算法结束后,恰在最短路上的边满足 D_i = D_j + c_{ij} 

  在最小费用流的计算中,我们每次沿 D_i = D_j + c_{ij} 的路径增广后都不会破坏条件 1,但是可能破坏了条件 2。不满足条件 2 的后果是什么呢?使我们找不到每条边都满足 D_i = D_j + c_{ij} 新的增广路。只好每次增广后使用 Dijkstra,SPFA 等等算法重新计算新的满足条件 2 的距离标号。这无疑是一种浪费。KM 算法中我们可以修改不断修改可行顶标,不断扩大可行子图,这里也同样,我们可以在始终满足条件 1 的距离标号上不断修改,直到可以继续增广(满足条件 2)。 

  回顾一下 KM 算法修改顶标的方法。根据最后一次寻找交错路不成功的 DFS,找到 d = /mathop {/min }/limits_{i /in V,j /notin V} /left/{ { - w_{ij} + A_i + B_j } /right/} ,左边的点增加 d ,右边的点减少 d 。这里也一样,根据最后一次寻找增广路不成功的 DFS,找到 d = /mathop {/min }/limits_{i /in V,j /notin V,u_{ij} > 0} /left/{ {c_{ij} - D_i + D_j } /right/} ,所有访问过的点距离标号增加 d。可以证明,这样不会破坏性质 1,而且至少有一条新的边进入了 D_i = D_j + c_{ij} 的子图。 

  算法的步骤就是初始标号设为 0 ,不断增广,如果不能增广,修改标号继续增广,直到彻底不能增广:源点的标号已经被加到了 + /infty  

  这样我们得到了一个简单的算法,只需要增广,改标号,各自只有 7 行,不需要 BFS,队列,SPFA,编程复杂度很低。由于实际的增广都是沿最短路进行的,所以理论时间复杂度与使用 SPFA 等等方法的连续最短路算法一致,但节省了 SPFA 或者 Dijkstra 的运算时间。实测发现这种算法常数很小,速度较快,employee 这道题所有数据加在一起耗时都在 2s 之内。由于没有找到参考程序,与 SPFA 实现的速度比较大家可以自己评测。 

 

原文地址:http://www.mathlinks.ro/weblog_entry.php?t=256087 

 

下面贴一下本人奇丑无比的代码:

主程序部分:

增广路过程:

顶标调整过程:

这篇关于最小费用最大流 zkw算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/579670

相关文章

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

综合安防管理平台LntonAIServer视频监控汇聚抖动检测算法优势

LntonAIServer视频质量诊断功能中的抖动检测是一个专门针对视频稳定性进行分析的功能。抖动通常是指视频帧之间的不必要运动,这种运动可能是由于摄像机的移动、传输中的错误或编解码问题导致的。抖动检测对于确保视频内容的平滑性和观看体验至关重要。 优势 1. 提高图像质量 - 清晰度提升:减少抖动,提高图像的清晰度和细节表现力,使得监控画面更加真实可信。 - 细节增强:在低光条件下,抖

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

poj 1287 Networking(prim or kruscal最小生成树)

题意给你点与点间距离,求最小生成树。 注意点是,两点之间可能有不同的路,输入的时候选择最小的,和之前有道最短路WA的题目类似。 prim代码: #include<stdio.h>const int MaxN = 51;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int P;int prim(){bool vis[MaxN];

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

秋招最新大模型算法面试,熬夜都要肝完它

💥大家在面试大模型LLM这个板块的时候,不知道面试完会不会复盘、总结,做笔记的习惯,这份大模型算法岗面试八股笔记也帮助不少人拿到过offer ✨对于面试大模型算法工程师会有一定的帮助,都附有完整答案,熬夜也要看完,祝大家一臂之力 这份《大模型算法工程师面试题》已经上传CSDN,还有完整版的大模型 AI 学习资料,朋友们如果需要可以微信扫描下方CSDN官方认证二维码免费领取【保证100%免费