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

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 <iostream>
 2 #include <algorithm>
 3 #include <stdlib.h>
 4 #include <time.h>
 5 #include <cmath>
 6 #include <cstdio>
 7 #include <string>
 8 #include <cstring>
 9 #include <vector>
10 #include <queue>
11 #include <stack>
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 <iostream>
 2 #include <algorithm>
 3 #include <stdlib.h>
 4 #include <time.h>
 5 #include <cmath>
 6 #include <cstdio>
 7 #include <string>
 8 #include <cstring>
 9 #include <vector>
10 #include <queue>
11 #include <stack>
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 < m; 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

 

 

转载于:https://www.cnblogs.com/yoyo-sincerely/p/5399446.html

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值