热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

增广路算法(最大流问题)

Edmonds-Karp算法:计算机科学中,Edmonds–Karp算法通过实现Ford–Fulkerson算法来计算网络中的最大流,其时间复杂度为O(VE2).该算法由

Edmonds-Karp算法:

计算机科学中, Edmonds–Karp算法通过实现Ford–Fulkerson算法来计算网络中的最大流,其时间复杂度为O(V E2). 该算法由Yefim (Chaim) Dinic 在1970年最先提出并由Jack Edmonds和Richard Karp 在1972年独立发表。

视频:Edmonds-Karp算法视频

最大流问题的目标:把最多的物品从s运送到t,其它点都是中转站。

算法思路:从0流开始不断增加流量,保持每次增加流量后都满足容量限制·斜对称性·流量平衡3个条件。

增广路:残量网络中任何一条从s到t的有向道路都对应一条原图中的增广路。

增广:只要求出增广路中所有残量的最小值d,把对应的所有边都加上d。

最大流判断条件:当残量网络中不存在增广路,则当前流就是最大流。

 

POJ 1459

模板题!

 1 #include 
2 #include
3 #include
4 #include
5 #include
6 #include
7 #include <string>
8 #include
9 #include
10 #include
11 #include
12 #include <set>
13
14 #define c_false ios_base::sync_with_stdio(false); cin.tie(0)
15 #define INF 0x3f3f3f3f
16 #define INFL 0x3f3f3f3f3f3f3f3f
17 #define zero_(x,y) memset(x , y , sizeof(x))
18 #define zero(x) memset(x , 0 , sizeof(x))
19 #define MAX(x) memset(x , 0x3f ,sizeof(x))
20 #define swa(x,y) {LL s;s=x;x=y;y=s;}
21 using namespace std ;
22 #define N 505
23
24 const double PI = acos(-1.0);
25 typedef long long LL ;
26
27 int n, np, nc, m, ans;
28 int Map[N][N], pre[N], que[N];
29 bool vis[N];
30
31 bool bfs(){
32 int head, tail;
33 zero(vis);
34 head = tail = 1;
35 que[tail++] = n;
36 vis[n] = true;
37 while(tail > head){
38 int u = que[head ++];
39 for(int i = 0; i <= n+1; i++){
40 if(!vis[i] && Map[u][i]){
41 pre[i] = u;
42 if(i == n+1) return true;
43 que[tail ++] = i;
44 vis[i] = true;
45 }
46 }
47 }
48 return false;
49 }
50
51 void End(){
52 int i, sum = INF;
53 for(i = n + 1; i != n; i = pre[i])
54 sum = min(sum, Map[pre[i]][i]);
55 for(i = n + 1; i != n; i = pre[i]){
56 Map[pre[i]][i] -= sum;
57 Map[i][pre[i]] += sum;
58 }
59 ans += sum;
60 }
61
62 int main(){
63 //freopen("in.txt","r",stdin);
64 //freopen("out.txt","w",stdout);
65 //ios_base::sync_with_stdio(false); cin.tie(0);
66 int u, v, w;
67 while(~scanf("%d%d%d%d", &n, &np, &nc, &m)){
68 zero(Map);
69 while(m--){
70 while(getchar() != '(');
71 scanf("%d,%d)%d", &u, &v, &w);
72 Map[u][v] += w;
73 }
74 while(np--){
75 while(getchar() != '(');
76 scanf("%d)%d", &u, &w);
77 Map[n][u] = w;
78 }
79 while(nc--){
80 while(getchar() != '(');
81 scanf("%d)%d", &u, &w);
82 Map[u][n+1] = w;
83 }
84 ans = 0;
85 while(bfs()) End();
86 printf("%d\n", ans);
87 }
88 return 0;
89 }
View Code

 

POJ 3041

题意:给定N*N矩阵,其中M个点,每次删除整行或者整列,求最小删除数;

思路:将行看成是x集合的点,列看成是y集合的点,M个点为对应边,就变成了二分图的最大匹配问题,

   这样就方便的转化成了最大流问题,只要在图中加入一个源点,和一个汇点即可;

   用0-n-1存储行,n-(2*n-1)存储列,2*n存储源点,2n+1存储汇点;

 

 1 #include 
2 #include
3 #include
4 #include
5 #include
6 #include
7 #include <string>
8 #include
9 #include
10 #include
11 #include
12 #include <set>
13
14 #define c_false ios_base::sync_with_stdio(false); cin.tie(0)
15 #define INF 0x3f3f3f3f
16 #define INFL 0x3f3f3f3f3f3f3f3f
17 #define zero_(x,y) memset(x , y , sizeof(x))
18 #define zero(x) memset(x , 0 , sizeof(x))
19 #define MAX(x) memset(x , 0x3f ,sizeof(x))
20 #define swa(x,y) {LL s;s=x;x=y;y=s;}
21 using namespace std ;
22 #define N 1005
23
24 const double PI = acos(-1.0);
25 typedef long long LL ;
26
27 int n, m, ans;
28 int Map[N][N], pre[N], que[N];
29 bool vis[N];
30
31 bool bfs(){
32 int head, tail;
33 zero(vis);
34 head = tail = 1;
35 que[tail++] = 2*n;
36 vis[2*n] = true;
37 while(tail > head){
38 int u = que[head ++];
39 for(int i = 0; i <= 2*n+1; i++){
40 if(!vis[i] && Map[u][i]){
41 pre[i] = u;
42 if(i == 2*n+1) return true;
43 que[tail ++] = i;
44 vis[i] = true;
45 }
46 }
47 }
48 return false;
49 }
50
51 void End(){
52 int i, sum = INF;
53 for(i = 2*n + 1; i != 2*n; i = pre[i])
54 sum = min(sum, Map[pre[i]][i]);
55 for(i = 2*n + 1; i != 2*n; i = pre[i]){
56 Map[pre[i]][i] -= sum;
57 Map[i][pre[i]] += sum;
58 }
59 ans += sum;
60 }
61
62 int main(){
63 //freopen("in.txt","r",stdin);
64 //freopen("out.txt","w",stdout);
65 //ios_base::sync_with_stdio(false); cin.tie(0);
66 int u, v;
67 while(~scanf("%d%d", &n, &m)){
68 zero(Map);
69 for(int i = 0; i ){
70 scanf("%d%d", &u, &v);
71 Map[u-1][v-1+n] = 1;
72 Map[2*n][u-1] = 1;
73 Map[n+v-1][(2*n)+1] = 1;
74 }
75 ans = 0;
76 while(bfs()) End();
77 printf("%d\n", ans);
78 }
79 return 0;
80 }
View Code

 

 


推荐阅读
  • 题目描述Takuru是一名情报强者,所以他想利用他强大的情报搜集能力来当中间商赚差价。Takuru的计划是让Hinae帮他去市场上买一个商品,然后再以另一个价格卖掉它。Takur ... [详细]
  • 如何自行分析定位SAP BSP错误
    The“BSPtag”Imentionedintheblogtitlemeansforexamplethetagchtmlb:configCelleratorbelowwhichi ... [详细]
  • 在Android开发中,使用Picasso库可以实现对网络图片的等比例缩放。本文介绍了使用Picasso库进行图片缩放的方法,并提供了具体的代码实现。通过获取图片的宽高,计算目标宽度和高度,并创建新图实现等比例缩放。 ... [详细]
  • HDU 2372 El Dorado(DP)的最长上升子序列长度求解方法
    本文介绍了解决HDU 2372 El Dorado问题的一种动态规划方法,通过循环k的方式求解最长上升子序列的长度。具体实现过程包括初始化dp数组、读取数列、计算最长上升子序列长度等步骤。 ... [详细]
  • 本文主要解析了Open judge C16H问题中涉及到的Magical Balls的快速幂和逆元算法,并给出了问题的解析和解决方法。详细介绍了问题的背景和规则,并给出了相应的算法解析和实现步骤。通过本文的解析,读者可以更好地理解和解决Open judge C16H问题中的Magical Balls部分。 ... [详细]
  • 796.[APIO2012]派遣在一个忍者的帮派里,一些忍者们被选中派遣给顾客,然后依据自己的工作获取报偿。在这个帮派里,有一名忍者被称之为Master。 ... [详细]
  • 向QTextEdit拖放文件的方法及实现步骤
    本文介绍了在使用QTextEdit时如何实现拖放文件的功能,包括相关的方法和实现步骤。通过重写dragEnterEvent和dropEvent函数,并结合QMimeData和QUrl等类,可以轻松实现向QTextEdit拖放文件的功能。详细的代码实现和说明可以参考本文提供的示例代码。 ... [详细]
  • 本文讨论了使用差分约束系统求解House Man跳跃问题的思路与方法。给定一组不同高度,要求从最低点跳跃到最高点,每次跳跃的距离不超过D,并且不能改变给定的顺序。通过建立差分约束系统,将问题转化为图的建立和查询距离的问题。文章详细介绍了建立约束条件的方法,并使用SPFA算法判环并输出结果。同时还讨论了建边方向和跳跃顺序的关系。 ... [详细]
  • 本文介绍了一种划分和计数油田地块的方法。根据给定的条件,通过遍历和DFS算法,将符合条件的地块标记为不符合条件的地块,并进行计数。同时,还介绍了如何判断点是否在给定范围内的方法。 ... [详细]
  • 《数据结构》学习笔记3——串匹配算法性能评估
    本文主要讨论串匹配算法的性能评估,包括模式匹配、字符种类数量、算法复杂度等内容。通过借助C++中的头文件和库,可以实现对串的匹配操作。其中蛮力算法的复杂度为O(m*n),通过随机取出长度为m的子串作为模式P,在文本T中进行匹配,统计平均复杂度。对于成功和失败的匹配分别进行测试,分析其平均复杂度。详情请参考相关学习资源。 ... [详细]
  • 动态规划算法的基本步骤及最长递增子序列问题详解
    本文详细介绍了动态规划算法的基本步骤,包括划分阶段、选择状态、决策和状态转移方程,并以最长递增子序列问题为例进行了详细解析。动态规划算法的有效性依赖于问题本身所具有的最优子结构性质和子问题重叠性质。通过将子问题的解保存在一个表中,在以后尽可能多地利用这些子问题的解,从而提高算法的效率。 ... [详细]
  • CF:3D City Model(小思维)问题解析和代码实现
    本文通过解析CF:3D City Model问题,介绍了问题的背景和要求,并给出了相应的代码实现。该问题涉及到在一个矩形的网格上建造城市的情景,每个网格单元可以作为建筑的基础,建筑由多个立方体叠加而成。文章详细讲解了问题的解决思路,并给出了相应的代码实现供读者参考。 ... [详细]
  • Linux环境变量函数getenv、putenv、setenv和unsetenv详解
    本文详细解释了Linux中的环境变量函数getenv、putenv、setenv和unsetenv的用法和功能。通过使用这些函数,可以获取、设置和删除环境变量的值。同时给出了相应的函数原型、参数说明和返回值。通过示例代码演示了如何使用getenv函数获取环境变量的值,并打印出来。 ... [详细]
  • 本文介绍了一个题目的解法,通过二分答案来解决问题,但困难在于如何进行检查。文章提供了一种逃逸方式,通过移动最慢的宿管来锁门时跑到更居中的位置,从而使所有合格的寝室都居中。文章还提到可以分开判断两边的情况,并使用前缀和的方式来求出在任意时刻能够到达宿管即将锁门的寝室的人数。最后,文章提到可以改成O(n)的直接枚举来解决问题。 ... [详细]
  • 3.223.28周学习总结中的贪心作业收获及困惑
    本文是对3.223.28周学习总结中的贪心作业进行总结,作者在解题过程中参考了他人的代码,但前提是要先理解题目并有解题思路。作者分享了自己在贪心作业中的收获,同时提到了一道让他困惑的题目,即input details部分引发的疑惑。 ... [详细]
author-avatar
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有