博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
POJ 3204 Ikki's Story I - Road Reconstruction
阅读量:5113 次
发布时间:2019-06-13

本文共 2789 字,大约阅读时间需要 9 分钟。

题意:给出一个流网络,问有仅将一条边的容量增大可以使得最大流增大的边有多少条。

首先,肯定将是满流的边的容量增大才可能使得最大流增大。但是这些满流的边不一定能是最小割,即使它出现在最小割中,把它的容量增大也不一定能使最大流增大。那么考虑一下求完最大流后不能增广的原因,或者说,怎么样才能在残量网络中增广呢?从S出发,dfs出所有满流的正向边,给dfs经过的点标上f1。再从T出发,逆向dfs出所有满流的边,给dfs经过的点标上f2。那么如果有一条边是这样的:满流且有f1标记和f2标记,那么增加它的容量,残量网络就又可以增广了,最大流会增加,因此这样的边就是我们需要寻找的边。

1 #include 
2 #include
3 #include
4 #include
5 #define INF 1<<30 6 #define maxn 510 7 #define maxm 100000 8 using namespace std; 9 10 int u[maxm],v[maxm],next[maxm],w[maxm]; 11 int first[maxn],d[maxn],work[maxn],q[maxn]; 12 int e,S,T; 13 bool f1[maxn],f2[maxn]; 14 int n,m; 15 void init(){ 16 memset(first,-1,sizeof(first)); 17 e = 0; 18 } 19 20 void add_edge(int a,int b,int c){ 21 u[e] = a;v[e] = b;next[e] = first[a];w[e] = c;first[a] = e++; 22 u[e] = b;v[e] = a;next[e] = first[b];w[e] = 0;first[b] = e++; 23 } 24 25 int bfs(){ 26 int rear = 0; 27 memset(d,-1,sizeof(d)); 28 d[S] = 0;q[rear++] = S; 29 for(int i = 0;i < rear;i++){ 30 for(int j = first[q[i]];j != -1;j = next[j]) 31 if(w[j] && d[v[j]] == -1){ 32 d[v[j]] = d[q[i]] + 1; 33 q[rear++] = v[j]; 34 if(v[j] == T) return 1; 35 } 36 } 37 return 0; 38 } 39 40 int dfs(int cur,int a){ 41 if(cur == T) return a; 42 for(int &i = work[cur];i != -1;i = next[i]) 43 if(w[i] && d[v[i]] == d[cur] + 1) 44 if(int t = dfs(v[i],min(a,w[i]))){ 45 w[i] -= t;w[i^1] += t; 46 return t; 47 } 48 return 0; 49 } 50 51 int dinic(){ 52 int ans = 0; 53 while(bfs()){ 54 memcpy(work,first,sizeof(first)); 55 while(int t = dfs(S,INF)) ans += t; 56 } 57 return ans; 58 } 59 60 void dfs1(int u){ 61 f1[u] = 1; 62 for(int i = first[u];i != -1;i = next[i]){ 63 if(!f1[v[i]] && w[i]) dfs1(v[i]); 64 } 65 } 66 67 void dfs2(int u){ 68 f2[u] = 1; 69 for(int i = first[u];i != -1;i = next[i]){ 70 if(!f2[v[i]] && w[i^1]) dfs2(v[i]); 71 } 72 } 73 74 int solve(){ 75 memset(f1,0,sizeof(f1)); 76 memset(f2,0,sizeof(f2)); 77 dfs1(S); 78 dfs2(T); 79 int ret = 0; 80 for(int i = 0;i < e;i += 2){ 81 if(f1[u[i]] && f2[v[i]] && !w[i]) 82 ret++; 83 } 84 return ret; 85 } 86 87 int main() 88 { 89 init(); 90 int n,m; 91 scanf("%d%d",&n,&m); 92 S = 0,T = n-1; 93 for(int i = 0;i < m;i++){ 94 int a,b,c; 95 scanf("%d%d%d",&a,&b,&c); 96 add_edge(a,b,c); 97 } 98 dinic(); 99 printf("%d\n",solve());100 return 0;101 }
View Code

 

转载于:https://www.cnblogs.com/zhexipinnong/p/3417156.html

你可能感兴趣的文章
js 字符串转化成数字
查看>>
duobango-tinySAK,20121214
查看>>
突破自己的知识体系
查看>>
TensorFlowSharp
查看>>
Egret Wing3 FTP使用方法
查看>>
使用PCF8563实现硬件看门狗
查看>>
vue devtools is not detected
查看>>
LeetCode "Restore IP Addresses"
查看>>
java连接mysql
查看>>
eclipse makefile
查看>>
Git fetch和git pull的区别
查看>>
leetcode - Flatten Binary Tree to Linked List
查看>>
LINUX 日志服务器的搭建
查看>>
MySQL中concat函数
查看>>
页面自动适应大小&&获取页面的大小
查看>>
jquery获取table,遍历输出tr中各个td的内容(转载)
查看>>
Leetcode No.85 *
查看>>
HTTP协议学习笔记-1
查看>>
面向对象
查看>>
jsf初学解决faces 中文输入乱码问题
查看>>