最大流—EK算法,流网络,残留网络,定理证明,详细代码

文章目录

    • 零、卡车运输
    • 一、流网络
      • 1.1流网络
      • 1.2流
      • 1.3最大流
      • 1.4残留网络
      • 1.5增广路径
      • 1.6流网络的割
      • 1.7最大流最小割定理
        • 1.7.1证明
      • 1.8Ford-Fulkerson方法
    • 二、Edmonds-Karp算法
      • 2.1定义
      • 2.2EK算法的实现
      • 2.3EK算法详细代码
      • 2.4OJ练习

零、卡车运输

Lucky Puck公司有冰球工厂Vancouver,即源点s,Winnipeg仓库是汇点t,冰球由卡车装载经由中间城市运送,但是每天只能有c(u , v)箱从城市u运送到城市v,在图上表示,每条边<u,v>有上有 f / c数字表示,f为沿这条有向边的运输冰球的箱数,c即为c(u,v),每天都有p箱冰球从Vancouver出发,p箱冰球到达Winnipeg。

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

作为公司老板的Puck,应该如何规划每条道路上的运载量才能使得每天从冰球工厂发出的冰球箱数p最大呢?


一、流网络

1.1流网络

流网络G = (V , E)是一个有向图,其中每条边(u , v)∈E均有一非负容量c(u , v) ≥ 0。如果(u , v) ∉ E,则c(u , v) = 0。

流网络中有两个特别的点:源点s汇点t如例中的工厂和仓库。

1.2流

设f(x , y)是定义在节点二元组(x∈V , y∈V)上的实数函数,且满足:

  1. 容量限制 : f ( x , y ) ≤ c ( x , y ) 容量限制:f(x , y) \le c(x , y) 容量限制:f(x,y)c(x,y)

  2. 反对称性 : f ( x , y ) = − f ( y , x ) 反对称性:f(x , y) = -f(y , x) 反对称性:f(x,y)=f(y,x)

  3. 流守恒性 : ∀ x ≠ s , x ≠ t , ∑ ( u , x ) ∈ E f ( u , x ) = ∑ ( x , v ) ∈ E f ( x , v ) 流守恒性:\forall x \ne s,x \ne t,\sum_{(u,x)\in E}f(u,x)=\sum_{(x,v)\in E}f(x,v) 流守恒性:x=s,x=t,(u,x)Ef(u,x)=(x,v)Ef(x,v)

f(x,y)称为流网络的流函数 , 对于(x , y)∈E , f(x , y)称为边的流量 , c(x , y) - f(x , y)称为边的剩余流量.

流 f 值的定义为
∣ f ∣ = ∑ v ∈ V f ( s , v ) \left |f \right | = \sum_{v\in V}f(s,v) f=vVf(s,v)
亦即 , 从源点s发出的总流。如例中每天从工厂发出的冰球的箱数。

1.3最大流

对于一个给定的流网络 , 合法的流函数 f 有很多. 使得流的值最大的流函数被称为网络的最大流 , 此时的流的值被称为网络的最大流量.

1.4残留网络

假定有一个流网络G = (V , E),其源点为s,汇点为t。设f为G中的一个流,一对顶点u, v∈V。在不超过容量c(u,v)的条件下从u到v之间可以压入的额外网络流量,就是**(u,v) 的残留容量(residual capacity)**,由下式定义:
c f ( u , v ) = c ( u , v ) − f ( u , v ) c_{f}(u,v) = c(u,v) - f(u,v) cf(u,v)=c(u,v)f(u,v)

给定一流网络G = (V , E)和流f,由f压得的G的残留网络是Gf = (V , Ef),其中:
E f = { ( u , v ) ∈ V × V : c f ( u , v ) > 0 } E_{f} = \{ (u,v)\in V\times V:c_{f}(u,v)>0 \} Ef={(u,v)V×V:cf(u,v)>0}
即,残留网络包含了流网络的所有点,和残留容量大于0的有向边。

注意Ef中的边既可以是E中的有向边也可以是其反向边,若(u , v)∈E,有f(u , v) < c(u , v),那么根据流网络的性质可知f(v , u) = -f(u,v),那么对应残留容量就是c(v , u) - (-f(u,v)) = c(v , u) + f(u , v) > 0,则其反向边也在残留网络中。

由残留网络可以得出引理

f 为G中的一个流,f‘为Gf中的一个流,那么f + f’仍为流网络G的一个流,其流量为| f + f’ | = | f | + | f‘ |

具体证明可以自己尝试或见《算法导论》

1.5增广路径

已知流网络G = (V , E)和流f,增广路径p残留网络Gf中由源点s到汇点t的一条简单路径

根据残留网络的定义,增广路径上的每条边的剩余容量都大于0,则该路径上的每条边都可以额外容纳一定的流量,这也和我们后续求最大流密切相关。

不难想出,增广路径可以增加的最大流量为该路径上边的最小残留容量

1.6流网络的割

流网络G = (V , E)的割(S , T)将V划分为S和T两部分,使得s∈S,t属于T,通过割的流量为S和T之间边上流量的代数和,但是割的容量仅包含从S到T的边的容量的代数和

如下图,割(S,T)的流量f(S,T) = 12 - 4 + 11 = 19

容量c(S,T) = 12 + 14 = 26

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

我们称容量最小的割为最小割

可以证明f(S , T) = | f | ≤ c(S, T)(证明见《算法导论》)

1.7最大流最小割定理

如果 f是具有源点s和汇点t的流网络G = (V , E)中的一个流,那么下列条件是等价的:

  1. f是G的一个最大流
  2. 残留网络Gf不包含增广路
  3. 存在G的某个割(S , T),有| f | = c(S , T)
1.7.1证明

采用循环证明法,(1) => (2) , (2) => (3) , (3) => (1)

(1) => (2):

很容易证明,采用反证法即可

假设Gf含增广路,那么我们可以在Gf中构造一流f’,| f‘ | = min(cf(u,v) , (u , v) ∈ Ef),那么f + f’仍为流网络的一个流(由1.4中介绍的引理可知),那么|f + f‘| > | f |,那么f就不是最大流,矛盾,则(1) => (2)成立

(2) => (1):

我们只需要在(2)的条件下构造出一个满足(3)的割即可。

选取集合S = {v ∈ V:Gf中从s到v存在一条通路},T = V - S,划分(S , T)为一个割。

对所有<u , v>,u∈S,v∈T,f(u , v) = c(u , v),否则v就属于S。

由此推出| f | = f(S , T) = c(S , T)

(3) => (1):

也很容易证明,由于由于| f | ≤ c(S, T),而此时| f | = c(S , T),故不存在比f更大的流,故f为最大流。

1.8Ford-Fulkerson方法

Ford-Fulkerson方法是最大流的经典求解方法,之所以称之为”方法“而非”算法“,是由于它包含具有不同运行时间的几种实现。

Ford-Fulkerson方法依赖于三种重要思想:残留网络(residual network)、增广路(augmenting path)和割(cut)。这些思想是最大流最小割定理的精髓,这里给出Ford-Fulkerson方法的特定实现。

伪代码如下:

Ford-Fulkerson(G , s , t)
初始化流f = 0
while 流网络中存在增广路pdo 沿着p增加f
return f

二、Edmonds-Karp算法

2.1定义

Ford-Fulkerson方法的第三行寻找p的过程,采用bfs来计算,这种实现方法即Edmonds-Karp算法

算法原理已经不需要介绍了,即最大流最小割定理,不同算法只是Ford-Fulkerson方法的不同实现,我们直接给出EK算法的实现。

2.2EK算法的实现

  • 建图,建立<u ,v , w>的有向边,初始容量为w,同时建立<u , v , 0>的反向边,初始容量为0,最大流maxflow = 0
  • 路径数组pre记录增广路上点的前驱边,数组incf存储每条边的剩余容量,标记数组vis用于记录点是否访问
  • bfs找增广路
    • 初始,源点s进入队列q,vis[s] = 1
    • 从队头取出u,遍历u发出的边i,incf[v] = min(w(i) , incf[u]),v入队,pre[v] = i
    • 如果u = t,说明找到增广路,return true
    • 否则队空,return false
  • 找到增广路,maxflow += incf[t] , 从t开始向前更新剩余容量,w(pre[x]) = incf[x]
  • 没找到增广路,说明最大流计算完毕,返回maxflow

2.3EK算法详细代码

采用链式前向星存图,每条边和其反向边的编号间的关系为异或1关系

关于链式前向星详见:一种实用的边的存储结构–链式前向星-CSDN博客

#define N 205
#define M 5005
const int MOD = 10000007;
const int inf = 0x3f3f3f3f3f3f3f3f;struct edge
{int v, w, nxt;
} edges[M << 1];
int head[N], idx = 0;
inline void addedge(int u, int v, int w)
{edges[idx] = {v, w, head[u]};head[u] = idx++;
}int n, m, s, t, incf[N], pre[N];
bool vis[N];bool bfs()
{memset(vis, 0, sizeof(vis));queue<int> q;q.emplace(s), vis[s] = true, incf[s] = inf;while (q.size()){int u = q.front();q.pop();for (int i = head[u]; ~i; i = edges[i].nxt){int v = edges[i].v;if (!vis[v] && edges[i].w){vis[v] = 1;incf[v] = min(incf[u], edges[i].w);pre[v] = i, q.emplace(v);if (v == t)return true;}}}return false;
}int EK()
{int maxflow = 0;while (bfs())//找增广路{int x = t;while (x != s)//更新剩余容量{int i = pre[x];edges[i].w -= incf[t];edges[i ^ 1].w += incf[t];x = edges[i ^ 1].v;}maxflow += incf[t];}return maxflow;
}//mainint a, b, c;memset(head, -1, sizeof(head));cin >> n >> m >> s >> t;for (int i = 0; i < m; i++)cin >> a >> b >> c, addedge(a, b, c), addedge(b, a, 0);cout << EK();

2.4OJ练习

P3376 【模板】网络最大流 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

板子题,直接跑板子即可

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <queue>
#include <unordered_set>
using namespace std;
#define sc scanf
#define int long long
#define N 205
#define M 5005
const int MOD = 10000007;
const int inf = 0x3f3f3f3f3f3f3f3f;struct edge
{int v, w, nxt;
} edges[M << 1];
int head[N], idx = 0;
inline void addedge(int u, int v, int w)
{edges[idx] = {v, w, head[u]};head[u] = idx++;
}int n, m, s, t, incf[N], pre[N];
bool vis[N];bool bfs()
{memset(vis, 0, sizeof(vis));queue<int> q;q.emplace(s), vis[s] = true, incf[s] = inf;while (q.size()){int u = q.front();q.pop();for (int i = head[u]; ~i; i = edges[i].nxt){int v = edges[i].v;if (!vis[v] && edges[i].w){vis[v] = 1;incf[v] = min(incf[u], edges[i].w);pre[v] = i, q.emplace(v);if (v == t)return true;}}}return false;
}int EK()
{int maxflow = 0;while (bfs()){int x = t;while (x != s){int i = pre[x];edges[i].w -= incf[t];edges[i ^ 1].w += incf[t];x = edges[i ^ 1].v;}maxflow += incf[t];}return maxflow;
}void solve()
{int a, b, c;memset(head, -1, sizeof(head));cin >> n >> m >> s >> t;for (int i = 0; i < m; i++)cin >> a >> b >> c, addedge(a, b, c), addedge(b, a, 0);cout << EK();
}signed main()
{ios::sync_with_stdio(false);cin.tie(0), cout.tie(0);freopen("in.txt", "r", stdin);int _ = 1;// cin >> _;while (_--)solve();
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/636609.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Unity导出Android项目踩坑记录

导出的时候需要注意以下地方的配置&#xff1a; 1、buildSetting-> 设置ExportProject 2、buildsetting ->playerSetting ->设置IL2CPP 3、设置ndk edit->preferences->external tools->ndk 如果unity的ndk版本和android项目里的ndk版本不一致会报错&…

【Qt开发】初识Qt

文章目录 1. Qt的背景1.1 Qt是什么1.2 Qt的发展史1.3 Qt支持的平台 2. Qt开发环境的搭建2.1 Qt SDK下载2.2 Qt SDK的安装 3. 一个简单的Qt模板程序的创建4. Qt模板程序的代码讲解4.1 main.cpp4.2 widget.h4.3 widget.cpp4.4 widget.ui4.5 test_1_18.pro4.6 一些中间文件 5. Qt在…

keil软件仿真

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 文章目录 前言一、pandas是什么&#xff1f;二、使用步骤 1.引入库2.读入数据总结 前言 例如&#xff1a;随着人工智能的不断发展&#xff0c;机器学习这门技术也越来越重要…

vue基于Spring Boot框架的甘肃敦煌文化旅游管理系统

本敦煌文化旅游管理系统是为了提高用户查阅信息的效率和管理人员管理信息的工作效率&#xff0c;可以快速存储大量数据&#xff0c;还有信息检索功能&#xff0c;这大大的满足了用户和管理员这两者的需求。操作简单易懂&#xff0c;合理分析各个模块的功能&#xff0c;尽可能优…

(蓝桥杯每日一题)平方末尾及补充(常用的字符串函数功能)

能够表示为某个整数的平方的数字称为“平方数 虽然无法立即说出某个数是平方数&#xff0c;但经常可以断定某个数不是平方数。因为平方数的末位只可能是:0,1,4,5,6,9 这 6 个数字中的某个。所以&#xff0c;4325435332 必然不是平方数。 如果给你一个 2 位或 2 位以上的数字&am…

六、标准对话框、多应用窗体

一、标准对话框 Qt提供了一些常用的标准对话框&#xff0c;如打开文件对话框、选择颜色对话框、信息提示和确认选择对话框、标准输入对话框等。1、预定义标准对话框 &#xff08;1&#xff09;QFileDialog 文件对话框 QString getOpenFileName() 打开一个文件QstringList ge…

You need to add dependency of ‘poi-ooxml‘ to your project, and version >= 4.1.2

原因 由于在依赖中引用了多个版本的 hutool,导致在最终打包时使用的版本不是由在开发时所引用的版本 cn.hutool.core.exceptions.DependencyException: You need to add dependency of poi-ooxml to your project, and version > 4.1.2at cn.hutool.poi.excel.ExcelUtil.get…

MyBatis-Plus 日常操作

本文主要介绍 mybatis-plus 日常操作。 一、快速开始 本文基于 springboot、maven、jdk1.8、mysql 环境。 新建如下数据库&#xff1a; 建议大家选择 utf8mb4 这种字符集&#xff0c;做过微信的同学应该会知道&#xff0c;微信用户名称的表情&#xff0c;是需要这种字符集才…

基于python旅游推荐系统 协同过滤算法 爬虫 Echarts可视化 Django框架(源码)✅

毕业设计&#xff1a;2023-2024年计算机专业毕业设计选题汇总&#xff08;建议收藏&#xff09; 毕业设计&#xff1a;2023-2024年最新最全计算机专业毕设选题推荐汇总 &#x1f345;感兴趣的可以先收藏起来&#xff0c;点赞、关注不迷路&#xff0c;大家在毕设选题&#xff…

Android14之DefaultKeyedVector实现(一百八十二)

简介&#xff1a; CSDN博客专家&#xff0c;专注Android/Linux系统&#xff0c;分享多mic语音方案、音视频、编解码等技术&#xff0c;与大家一起成长&#xff01; 优质专栏&#xff1a;Audio工程师进阶系列【原创干货持续更新中……】&#x1f680; 优质专栏&#xff1a;多媒…

计算机组成原理 第一弹

ps&#xff1a;本文章的图片来源都是来自于湖科大教书匠高老师的视频&#xff0c;声明&#xff1a;仅供自己复习&#xff0c;里面加上了自己的理解 这里附上视频链接地址&#xff1a;1-2 计算机的发展_哔哩哔哩_bilibili ​​ 目录 &#x1f680;计算机系统 &#x1f680;计…

基于SpringBoot的欢乐校园管理系统

文章目录 项目介绍主要功能截图&#xff1a;部分代码展示设计总结项目获取方式 &#x1f345; 作者主页&#xff1a;超级无敌暴龙战士塔塔开 &#x1f345; 简介&#xff1a;Java领域优质创作者&#x1f3c6;、 简历模板、学习资料、面试题库【关注我&#xff0c;都给你】 &…

SpringMVC搭建环境

idea创建java项目后添加webapp怎么配置 1.首先在main下创建一个普通文件webapp 2. 3.选中你的项目&#xff0c;添加Web 4.修改这两处的路径&#xff0c;修改为你webapp所在的路径 先修改左下角的路径&#xff0c;然后再添加web.xml. 然后再修改右上角的地址&#xff0c;注…

基于SpringBoot Vue家政服务预约平台系统

大家好✌&#xff01;我是Dwzun。很高兴你能来阅读我&#xff0c;我会陆续更新Java后端、前端、数据库、项目案例等相关知识点总结&#xff0c;还为大家分享优质的实战项目&#xff0c;本人在Java项目开发领域有多年的经验&#xff0c;陆续会更新更多优质的Java实战项目&#x…

在vue中使用echarts渲染地图,geo点击某个区域可高亮,取消

一、安装echarts npm install echarts --save二、main.js引入注册 import Vue from "vue";import * as echarts from "echarts";Vue.prototype.$echarts echarts;三、vue文件中使用echarts <template><div class"page-warp"><…

Vue-30、Vue非单文件组件。

非单文件组件&#xff1a; 一个组件包含n个组件 <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><title>非单文件组件</title><script type"text/javascript" src"https://cdn.jsde…

BossCMS RCE

简介 BossCMS是一款基于自主研发PHP框架MySQL架构的内容管理系统&#xff0c;能够满足各类网站开发建设的需求。系统开源、安全、稳定、简洁、易开发、专注为中小型企业及政企单位、个人站长、广大开发者、建站公司提供一套简单好用的网站内容管理系统解决方案。官网提供了大量…

SystemVerilog中数组内置函数sum()的一个注意点

Systemverilog内置了数组求和运算方法(sum())&#xff0c;将数组的所有元素累加起来&#xff0c;返回一个最终值。在使用时要注意数组类型的位宽&#xff0c;通常情况下&#xff0c;如果你将一组单bit的值加起来&#xff0c;Systemverilog会使用足够的精度来确保不丢失任何bit的…

INS-06003错误处理

在麒麟V10操作系统上安装Oracle RAC 19C&#xff0c;安装GI的建立互信步骤中&#xff0c;遇到INS-06003错误&#xff1a; [INS-06003] Failed to setup password SSH connectivity with following node(s) 查看详细信息&#xff1a; PRVG-11001: PRCZ-2136: PRCZ-2006: 此时在操…

实验四 SQL语言

&#x1f57a;作者&#xff1a; 主页 我的专栏C语言从0到1探秘C数据结构从0到1探秘Linux &#x1f618;欢迎关注&#xff1a;&#x1f44d;点赞&#x1f64c;收藏✍️留言 &#x1f3c7;码字不易&#xff0c;你的&#x1f44d;点赞&#x1f64c;收藏❤️关注对我真的很重要&…