“zkw” 费用流算法在哪些图上慢(摘自https://www.cnblogs.com/ECJTUACM-873284962/p/7744943.html

实践中, 上面的这个算法非常奇怪. 在某一些图上, 算法速度非常快,

另一些图上却比纯 SPFA 增广的算法慢. 不少同学经过实测总结的结果是稠密图上比较快,

稀疏图上比较慢, 但也不尽然. 这里我从理论上分析一下, 究竟这个算法用于哪些图可以得到理想的效果.

先分析算法的增广流程. 和 SPFA 直接算法相比, 由于同属于沿最短路增广的算法,

实际进行的增流操作并没有太多的区别, 每次的增流路径也大同小异. 因此不考虑多路增广时,

增广次数应该基本相同. 运行时间上主要的差异应当在于如何寻找增流路径的部分.

那么 zkw 算法的优势在于哪里呢?

与 SPFA 相比, KM 的重标号方式明显在速度上占优, 每次只是一个对边的扫描操作而已.

而 SPFA 需要维护较为复杂的标号和队列操作, 同时为了修正标号, 需要不止一次地访问某些节点, 速度会慢不少.

另外, 在 zkw 算法中, 增广是多路进行的, 同时可能在一次重标号后进行多次增广.

这个特点可以在许多路径都费用相同的时候派上用场, 进一步减少了重标号的时间耗费.

下面想一想 zkw 算法的劣势, 也就是 KM 重标号方式存在的问题.

KM 重标号的主要问题就是, 不保证经过一次重标号之后能够存在增广路.

最差情况下, 一次只能在零权网络中增加一条边而已. 这时算法就会反复重标号,

反复尝试增广而次次不能增广, 陷入弄巧成拙的境地.

接下来要说什么, 大家可能已经猜到了. 对于最终流量较大, 而费用取值范围不大的图,

或者是增广路径比较短的图 (如二分图), zkw 算法都会比较快. 原因是充分发挥优势.

比如流多说明可以同一费用反复增广, 费用窄说明不用改太多距离标号就会有新增广路,

增广路径短可以显著改善最坏情况, 因为即使每次就只增加一条边也可以很快凑成最短路.

如果恰恰相反, 流量不大, 费用不小, 增广路还较长, 就不适合 zkw 算法了.

费用流板题:luoguP3381 【模板】最小费用最大流

zkw版/板:

#include <bits/stdc++.h>
using namespace std;
void read(int &num)
{
char ch; bool flag=0;
while(!isdigit(ch=getchar()))if(ch=='-')flag=!flag;
for(num=0;isdigit(ch);num=num*10+ch-'0',ch=getchar());
if(flag)num=-num;
}
const int MAXN = 5005;
const int MAXM = 50005;
const int inf = 1e9;
int n, m, S, T, Ans;
int cnt = 1, fir[MAXN], nxt[MAXM*2], to[MAXM*2], c[MAXM*2], w[MAXM*2];
int dis[MAXN];
bool vis[MAXN]; void Add(int u, int v, int cc, int wt)
{
to[++cnt] = v; nxt[cnt] = fir[u]; fir[u] = cnt; c[cnt] = cc; w[cnt] = wt;
} int Aug(int u, int flow)
{
if(u == T) { return flow; }
int used = 0, delta = 0;
vis[u] = true;
for(int i = fir[u]; i; i = nxt[i])
if(!vis[to[i]] && c[i] && dis[u] == dis[to[i]] + w[i])
{
delta = Aug(to[i], min(c[i], flow-used));
c[i] -= delta, c[i^1] += delta, used += delta;
Ans += delta * w[i];
if(used == flow) break;
}
return used;
} bool Update()
{
int tmp = inf;
for(int i = 1; i <= n; i++) if(vis[i])
for(int j = fir[i]; j; j = nxt[j])
if(!vis[to[j]] && c[j]) tmp = min(tmp, dis[to[j]]-dis[i]+w[j]);
if(tmp == inf) return false;
for(int i = 1; i <= n; i++) if(vis[i]) dis[i] += tmp;
return true;
} int Cost_flow()
{
int flow = 0, tmp = 0;
do {
do { flow += tmp; memset(vis, 0, sizeof vis); }while(tmp=Aug(S, inf));
}while(Update());
return flow;
} int main ()
{
int u, v, x, y;
read(n), read(m), read(S), read(T);
for(int i = 1; i <= m; i++)
{
read(u), read(v), read(x), read(y);
Add(u, v, x, y);
Add(v, u, 0, -y);
}
int Max_flow = Cost_flow();
printf("%d %d\n", Max_flow, Ans);
}

spfa版:

#include <bits/stdc++.h>
using namespace std;
void read(int &num)
{
char ch; bool flag=0;
while(!isdigit(ch=getchar()))if(ch=='-')flag=!flag;
for(num=0;isdigit(ch);num=num*10+ch-'0',ch=getchar());
if(flag)num=-num;
}
const int MAXN = 5005;
const int MAXM = 50005;
const int inf = 1e9;
int n, m, S, T, Ans;
int cnt = 1, fir[MAXN], nxt[MAXM*2], to[MAXM*2], c[MAXM*2], w[MAXM*2];
int dis[MAXN];
bool vis[MAXN]; void Add(int u, int v, int cc, int wt)
{
to[++cnt] = v; nxt[cnt] = fir[u]; fir[u] = cnt; c[cnt] = cc; w[cnt] = wt;
} bool spfa(int s, int t)
{
memset(vis, 0, sizeof vis);
for(int i = 0; i <= n; i++) dis[i] = inf;
dis[t] = 0; vis[t] = true;
deque<int> Q; Q.push_back(t);
while(!Q.empty())
{
int u = Q.front(); Q.pop_front();
for(int i = fir[u]; i; i = nxt[i])
if(c[i^1] && dis[to[i]] > dis[u] - w[i])
{
dis[to[i]] = dis[u] - w[i];
if(!vis[to[i]])
{
vis[to[i]] = true;
if(!Q.empty() && dis[to[i]] < dis[Q.front()]) Q.push_front(to[i]);
else Q.push_back(to[i]);
}
}
vis[u] = 0;
}
return dis[S] < inf;
} int Aug(int u, int flow)
{
if(u == T) { vis[T] = true; return flow; }
int used = 0, delta = 0;
vis[u] = true;
for(int i = fir[u]; i; i = nxt[i])
if(!vis[to[i]] && c[i] && dis[u] == dis[to[i]] + w[i])
{
delta = Aug(to[i], min(c[i], flow-used));
Ans += delta * w[i], c[i] -= delta, c[i^1] += delta, used += delta;
if(used == flow) break;
}
return used;
} int Cost_flow()
{
int flow = 0;
while(spfa(S, T))
{
vis[T] = 1;
while(vis[T])
{
memset(vis, 0, sizeof vis);
flow += Aug(S, inf);
}
}
return flow;
} int main ()
{
int u, v, x, y;
read(n), read(m), read(S), read(T);
for(int i = 1; i <= m; i++)
{
read(u), read(v), read(x), read(y);
Add(u, v, x, y);
Add(v, u, 0, -y);
}
int Max_flow = Cost_flow();
printf("%d %d\n", Max_flow, Ans);
}

最小费用流模板(zkw与spfa)的更多相关文章

  1. [转]从入门到精通: 最小费用流的“zkw算法”

    >>>> 原文地址:最小费用流的“zkw算法” <<<< 1. 网络流的一些基本概念 很多同学建立过网络流模型做题目, 也学过了各种算法, 但是对于基本 ...

  2. POJ 3068 运送危险化学品 最小费用流 模板题

    "Shortest" pair of paths Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 1215 ...

  3. LG P2285 [模板]负环(spfa判负环)

    题目描述 寻找一个从顶点1所能到达的负环,负环定义为:一个边权之和为负的环. 输入格式 第一行一个正整数T表示数据组数,对于每组数据: 第一行两个正整数N M,表示图有N个顶点,M条边 接下来M行,每 ...

  4. HYSBZ 1061 志愿者招募 【最小费用流】【差分】【最小费用流模板】

    #include<stdio.h> #include<queue> #define MAXN 1003 #define MAXM 10002*4 #define INF 100 ...

  5. 洛谷P3385 [模板]负环 [SPFA]

    题目传送门 题目描述 暴力枚举/SPFA/Bellman-ford/奇怪的贪心/超神搜索 输入输出格式 输入格式: 第一行一个正整数T表示数据组数,对于每组数据: 第一行两个正整数N M,表示图有N个 ...

  6. [模板] zkw线段树

    zkw线段树 code1简单版本 code2差分版本(暂无) code1:(有注释) //By Menteur_Hxy #include<cstdio> #include<iostr ...

  7. 【模板】裸SPFA

    SPFA可以处理带负边权的图,可以判负环,然而SPFA容易被卡,即使加了各种优化. 队列优化的贝尔福德曼:裸SPFA //SPFA #include<bits/stdc++.h> usin ...

  8. 网络流之最大流与最小费用流入门&&模板

    理解处 刷题处 模板处 最大流模板 处理重边的+(优化) #include<bits/stdc++.h> using namespace std; ; const int INF = 0x ...

  9. POJ 2449Remmarguts' Date K短路模板 SPFA+A*

    K短路模板,A*+SPFA求K短路.A*中h的求法为在反图中做SPFA,求出到T点的最短路,极为估价函数h(这里不再是估价,而是准确值),然后跑A*,从S点开始(此时为最短路),然后把与S点能达到的点 ...

随机推荐

  1. Python 判断字符串是否包含中文

    一.摘要 使用 xlrd 模块打开带中文的excel文件时,会报错. FileNotFoundError: [Errno 2] No such file or directory: 'xx.xlsx' ...

  2. SpringBoot指定额外需要扫描的包

    我们都知道,SpringBoot主启动类标注了@SpringBootApplication注解,该注解引入了@ComponentScan注解 所以默认的包扫描规则是,程序会自动扫描主启动类所在包及其子 ...

  3. Consul 的安装与基本使用

    什么是 Consul ​ Consul是一种服务网格解决方案,提供具有服务发现,配置和分段功能的全功能控制平面.这些功能中的每一个都可以根据需要单独使用,也可以一起使用以构建全服务网格.Consul需 ...

  4. 深度学习 Bottleneck layer / Bottleneck feature

    最近在学习deeplearning的时候接触到了bottle-neck layer,好奇它的作用于是便扒了一些论文(论文链接放在文末吧),系统的了解一下bottle-neck feature究竟有什么 ...

  5. ByteBuf源码

    ByteBuf是顶层的抽象类,定义了用于传输数据的ByteBuf需要的方法和属性. AbstractByteBuf 直接继承ByteBuf,一些公共属性和方法的公共逻辑会在这里定义.例如虽然不同性质的 ...

  6. java之hibernate之基于外键的一对一单向关联映射

    这篇讲解基于外键的一对一单向关联映射 1.考察如下信息,人和身份证之间是一个一对一的关系.表的设计 注意:基于外键的一对一关联的表结构和多对一的表结构是一致的,但是,外键是唯一的. 2.类的结构 Pe ...

  7. java之hibernate之配置讲解

    1.映射文件:User.hbm.xml <!-- package 指向class中所有类的包名,可以直接在指定类名时同时指定包名 --> <hibernate-mapping pac ...

  8. MySQL分库备份

    1.需求概述 每天00:00备份MySQL数据库数据: 每一个库生成一个文件,使用gzip压缩,文件名:backup_库名_yyyymmdd.sql.gz,注意yyyymmdd需要是前一天: 备份文件 ...

  9. Switch开关在element-ui表格中点击没有效果解决方法

    <el-table-column label="三审" align="center"> <template slot-scope=" ...

  10. 通过 Web Deploy 发布的配置

    罩着别人的配置弄了一下午,死活认证通不过,后来好不容易试出来,备忘. 服务端:安装IIS,启动管理程序,安装Web Deploy, 建立网站,建立IIS用户,进网站的权限管理里面给IIS用户授权. V ...