给出每条边的下界 求最小流

板题

提供两个板子代码 虽然这个题 第一个比较快

但在loj上https://loj.ac/problem/117

的板题  第一个1700+ms

第二个才600+ms  用输入挂后400+ms

  1. #include <iostream>
  2. #include <cstdio>
  3. #include <sstream>
  4. #include <cstring>
  5. #include <map>
  6. #include <cctype>
  7. #include <set>
  8. #include <vector>
  9. #include <stack>
  10. #include <queue>
  11. #include <algorithm>
  12. #include <cmath>
  13. #include <bitset>
  14. #define rap(i, a, n) for(int i=a; i<=n; i++)
  15. #define rep(i, a, n) for(int i=a; i<n; i++)
  16. #define lap(i, a, n) for(int i=n; i>=a; i--)
  17. #define lep(i, a, n) for(int i=n; i>a; i--)
  18. #define rd(a) scanf("%d", &a)
  19. #define rlld(a) scanf("%lld", &a)
  20. #define rc(a) scanf("%c", &a)
  21. #define rs(a) scanf("%s", a)
  22. #define rb(a) scanf("%lf", &a)
  23. #define rf(a) scanf("%f", &a)
  24. #define pd(a) printf("%d\n", a)
  25. #define plld(a) printf("%lld\n", a)
  26. #define pc(a) printf("%c\n", a)
  27. #define ps(a) printf("%s\n", a)
  28. #define MOD 2018
  29. #define LL long long
  30. #define ULL unsigned long long
  31. #define Pair pair<int, int>
  32. #define mem(a, b) memset(a, b, sizeof(a))
  33. #define _ ios_base::sync_with_stdio(0),cin.tie(0)
  34. //freopen("1.txt", "r", stdin);
  35. using namespace std;
  36. const int maxn = 1e5 + , INF = 0x7fffffff;
  37. int n, m, s, t;
  38. int head[maxn], cur[maxn], vis[maxn], nex[maxn << ], in[maxn];
  39. int d[maxn];
  40. int cnt;
  41. struct node
  42. {
  43. int u, v;
  44. int c;
  45. }Node[maxn << ];
  46.  
  47. void add_(int u, int v, int c)
  48. {
  49. Node[cnt].u = u;
  50. Node[cnt].v = v;
  51. Node[cnt].c = c;
  52. nex[cnt] = head[u];
  53. head[u] = cnt++;
  54. }
  55.  
  56. void add(int u, int v, int c)
  57. {
  58. add_(u, v, c);
  59. add_(v, u, );
  60. }
  61.  
  62. bool bfs()
  63. {
  64. queue<int> Q;
  65. mem(d, );
  66. Q.push(s);
  67. d[s] = ;
  68. while(!Q.empty())
  69. {
  70. int u = Q.front(); Q.pop();
  71. for(int i = head[u]; i != -; i = nex[i])
  72. {
  73. int v = Node[i].v;
  74. if(!d[v] && Node[i].c > )
  75. {
  76. d[v] = d[u] + ;
  77. Q.push(v);
  78. if(v == t) return ;
  79. }
  80. }
  81. }
  82. return d[t] != ;
  83. }
  84.  
  85. int dfs(int u, int cap)
  86. {
  87. int ret = ;
  88. if(u == t || cap == )
  89. return cap;
  90. for(int &i = cur[u]; i != -; i = nex[i])
  91. {
  92. int v = Node[i].v;
  93. if(d[v] == d[u] + && Node[i].c > )
  94. {
  95. int V = dfs(v, min(cap, Node[i].c));
  96. Node[i].c -= V;
  97. Node[i ^ ].c += V;
  98. ret += V;
  99. cap -= V;
  100. if(cap == ) break;
  101. }
  102. }
  103. if(cap > ) d[u] = -;
  104. return ret;
  105. }
  106.  
  107. int Dinic()
  108. {
  109. int ans = ;
  110. while(bfs())
  111. {
  112. memcpy(cur, head, sizeof head);
  113. ans += dfs(s, INF);
  114. }
  115. return ans;
  116. }
  117.  
  118. int tranf(string x)
  119. {
  120. if(x == "+") return ;
  121. if(x == "-") return n + ;
  122. int ans = , len = x.size();
  123. for(int i = ; i < len; i++)
  124. {
  125. ans += (x[i] - '') * pow(, len - i - );
  126. }
  127. return ans + ;
  128. }
  129.  
  130. int main()
  131. {
  132. int w;
  133. string u, v;
  134. int u_, v_;
  135. while(~scanf("%d%d", &n, &m))
  136. {
  137. if(n == && m == ) break;
  138. mem(head, -);
  139. cnt = ;
  140. mem(in, );
  141. int sum = ;
  142. s = , t = n + ;
  143. rap(i, , m)
  144. {
  145. // scanf(" %c %c %d", &u, &v, &w);
  146. cin >> u >> v >> w;
  147. u_ = tranf(u);
  148. v_ = tranf(v);
  149. add(u_, v_, INF - w);
  150. in[v_] += w;
  151. in[u_] -= w;
  152. }
  153. for(int i = ; i <= n + ; i++)
  154. {
  155. if(in[i] > )
  156. {
  157. sum += in[i];
  158. add(s, i, in[i]);
  159. }
  160. else
  161. add(i, t, -in[i]);
  162. }
  163. int ans = Dinic();
  164. add(n + , , INF);
  165. ans += Dinic();
  166. if(sum == ans)
  167. cout << Node[head[(n + )] ^ ].c << endl;
  168. else
  169. cout << "impossible" << endl;
  170.  
  171. }
  172. return ;
  173. }
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <sstream>
  4. #include <cstring>
  5. #include <map>
  6. #include <cctype>
  7. #include <set>
  8. #include <vector>
  9. #include <stack>
  10. #include <queue>
  11. #include <algorithm>
  12. #include <cmath>
  13. #include <bitset>
  14. #define rap(i, a, n) for(int i=a; i<=n; i++)
  15. #define rep(i, a, n) for(int i=a; i<n; i++)
  16. #define lap(i, a, n) for(int i=n; i>=a; i--)
  17. #define lep(i, a, n) for(int i=n; i>a; i--)
  18. #define rd(a) scanf("%d", &a)
  19. #define rlld(a) scanf("%lld", &a)
  20. #define rc(a) scanf("%c", &a)
  21. #define rs(a) scanf("%s", a)
  22. #define rb(a) scanf("%lf", &a)
  23. #define rf(a) scanf("%f", &a)
  24. #define pd(a) printf("%d\n", a)
  25. #define plld(a) printf("%lld\n", a)
  26. #define pc(a) printf("%c\n", a)
  27. #define ps(a) printf("%s\n", a)
  28. #define MOD 2018
  29. #define LL long long
  30. #define ULL unsigned long long
  31. #define Pair pair<int, int>
  32. #define mem(a, b) memset(a, b, sizeof(a))
  33. #define _ ios_base::sync_with_stdio(0),cin.tie(0)
  34. //freopen("1.txt", "r", stdin);
  35. using namespace std;
  36. const int maxn = 1e5 + , INF = 0x7fffffff;
  37. int n, m, s, t;
  38. int head[maxn], cur[maxn], vis[maxn], nex[maxn << ], in[maxn];
  39. int d[maxn];
  40. int cnt;
  41. struct node
  42. {
  43. int u, v, bz;
  44. int c;
  45. }Node[maxn << ];
  46.  
  47. void add_(int u, int v, int c, int bz)
  48. {
  49. Node[cnt].u = u;
  50. Node[cnt].v = v;
  51. Node[cnt].c = c;
  52. Node[cnt].bz = bz;
  53. nex[cnt] = head[u];
  54. head[u] = cnt++;
  55. }
  56.  
  57. void add(int u, int v, int c, int bz)
  58. {
  59. add_(u, v, c, bz);
  60. add_(v, u, , bz);
  61. }
  62.  
  63. bool bfs()
  64. {
  65. queue<int> Q;
  66. mem(d, );
  67. Q.push(s);
  68. d[s] = ;
  69. while(!Q.empty())
  70. {
  71. int u = Q.front(); Q.pop();
  72. for(int i = head[u]; i != -; i = nex[i])
  73. {
  74. int v = Node[i].v;
  75. if(!d[v] && Node[i].c > )
  76. {
  77. d[v] = d[u] + ;
  78. Q.push(v);
  79. if(v == t) return ;
  80. }
  81. }
  82. }
  83. return d[t] != ;
  84. }
  85.  
  86. int dfs(int u, int cap)
  87. {
  88. int ret = ;
  89. if(u == t || cap == )
  90. return cap;
  91. for(int &i = cur[u]; i != -; i = nex[i])
  92. {
  93. int v = Node[i].v;
  94. if(d[v] == d[u] + && Node[i].c > )
  95. {
  96. int V = dfs(v, min(cap, Node[i].c));
  97. Node[i].c -= V;
  98. Node[i ^ ].c += V;
  99. ret += V;
  100. cap -= V;
  101. if(cap == ) break;
  102. }
  103. }
  104. if(cap > ) d[u] = -;
  105. return ret;
  106. }
  107.  
  108. int Dinic()
  109. {
  110. int ans = ;
  111. while(bfs())
  112. {
  113. memcpy(cur, head, sizeof head);
  114. ans += dfs(s, INF);
  115. }
  116. return ans;
  117. }
  118.  
  119. int tranf(string x)
  120. {
  121. if(x == "+") return ;
  122. if(x == "-") return n + ;
  123. int ans = , len = x.size();
  124. for(int i = ; i < len; i++)
  125. {
  126. ans += (x[i] - '') * pow(, len - i - );
  127. }
  128. return ans + ;
  129. }
  130.  
  131. int main()
  132. {
  133. int w;
  134. string u, v;
  135. int u_, v_;
  136. while(~scanf("%d%d", &n, &m))
  137. {
  138. if(n == && m == ) break;
  139. mem(head, -);
  140. cnt = ;
  141. mem(in, );
  142. int sum = ;
  143. s = , t = n + ;
  144. rap(i, , m)
  145. {
  146. // scanf(" %c %c %d", &u, &v, &w);
  147. cin >> u >> v >> w;
  148. u_ = tranf(u);
  149. v_ = tranf(v);
  150. add(u_, v_, INF - w, );
  151. in[v_] += w;
  152. in[u_] -= w;
  153. }
  154. for(int i = ; i <= n + ; i++)
  155. {
  156. if(in[i] > )
  157. {
  158. sum += in[i];
  159. add(s, i, in[i], );
  160. }
  161. else
  162. add(i, t, -in[i], );
  163. }
  164. add(n + , , INF, );
  165. if(sum != Dinic())
  166. ps("impossible");
  167. else
  168. {
  169. sum = Node[head[n + ] ^ ].c;
  170. for(int i = ; i < cnt; i++)
  171. {
  172. if(!Node[i].bz) Node[i].v = ;
  173. }
  174. head[s] = head[t] = -;
  175. s = n + ;
  176. t = ;
  177. cout << sum - Dinic() << endl;
  178. }
  179.  
  180. }
  181. return ;
  182. }

Crazy Circuits HDU - 3157(有源汇有上下界最小流)的更多相关文章

  1. loj #117. 有源汇有上下界最小流

    题目链接 有源汇有上下界最小流,->上下界网络流 注意细节,边数组也要算上后加到SS,TT边. #include<cstdio> #include<algorithm> ...

  2. LOJ.117.[模板]有源汇有上下界最小流(Dinic)

    题目链接 有源汇有上下界最小流 Sol1. 首先和无源汇网络流一样建图,求SS->TT最大流: 然后连边(T->S,[0,INF]),再求一遍SS->TT最大流,答案为新添加边的流量 ...

  3. sgu 176 Flow construction(有源汇的上下界最小流)

    [题目链接] http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=11025 [模型] 有源汇点的上下界最小流.即既满足上下界又满足 ...

  4. Flow construction SGU - 176 有源汇有上下界最小流 二分法和回流法

    /** 题目:Flow construction SGU - 176 链接:https://vjudge.net/problem/SGU-176 题意: 有源汇有上下界的最小流. 给定n个点,m个管道 ...

  5. 【 POJ - 3801】Crazy Circuits(有源汇、上下界最小流)

    Description You’ve just built a circuit board for your new robot, and now you need to power it. Your ...

  6. sgu176 Flow Construction【有源汇有上下界最小流】

    同样是模板题. 首先将有源汇转换为无源汇,假设原来的源汇为st,我们加入的源汇为ST,那么我们应该从t到s连一条流量为+∞的边,使原来的st满足收支平衡,退化为普通节点. 分离必要边和其他边,从S到T ...

  7. bzoj 2502 清理雪道(有源汇的上下界最小流)

    [题意] 有一个DAG,要求每条边必须经过一次,求最少经过次数. [思路] 有上下界的最小流.  边的下界为1,上界为无穷.构造可行流模型,先不加ts边跑一遍最大流,然后加上t->s的inf边跑 ...

  8. LOJ117 有源汇有上下界最小流(上下界网络流)

    跑出可行流后从原来的汇点向原来的源点跑最大流,原图最小流=inf-maxflow.显然超源超汇的相关边对其也没有影响.原图最小流=可行流-原图新增流量,因为t向s流量增加相当于s向t流量减少.但为什么 ...

  9. bzoj2502: 清理雪道(有源汇有上下界最小流)

    传送门 别说话,自己看,我不会->这里 我这里用的建图方法是先跑一次最大流,连上$(t,s,inf)$之后再跑一遍,然后答案就是之前连的那条边的反向边的流量 据说还有种方法是连上$(t,s,in ...

随机推荐

  1. Ubuntu: 软件库(software repositories)

    Linux 生态下,几乎每个发行版都有自己的软件库(software repositories),Ubuntu 当然也不例外.Ubuntu 提供了四个不同的软件库,分别是 main.restricte ...

  2. Php7 开发笔记

    Ubuntu环境安装 http://www.jianshu.com/p/1d312d9f1be1 sudo apt-get install python-software-properties sof ...

  3. H5 21-属性选择器下

    21-属性选择器下 --> <!DOCTYPE html> <html lang="en"> <head> <meta charse ...

  4. K 班1-7,alpha,beta 作业成绩汇总

    K 班1-7,alpha,beta 作业成绩汇总 千帆竞发 详细得分 短学号 名 1 2 3 4 5 6 7 alpha beta TOTAL 505 基智 4.55 1 -2 0 0 -10 4.3 ...

  5. Python_动态参数、名称空间、作用域、作用域链、加载顺序、函数的嵌套、global、nonlocal

    1.动态参数 当实参数量与形参数量相等时,参数传递正常. def func1(a, b, c): pass func1(1, 2, 3) 当实参数量与形参数量不相等时,则会报错. def func1( ...

  6. PhpStorm的注册激活方法

    首先,需要修改本地的hosts文件(路径一般为C:\Windows\System32\drivers\etc\hosts),添加下面这行代码. 0.0.0.0 account.jetbrains.co ...

  7. java list 去重

    Guava ImmutableSet源码去重处理,Hash处理 - 汪小哥 - CSDN博客 https://blog.csdn.net/u012881904/article/details/6895 ...

  8. oracle事务的四个特性(ACID)

    事务产生的背景 当在PL/SQL中同时操作多个SQL语句,比如通过DML语句添加.修改或删除数据时,如何确保数据库数据不会因为意外而倒置错误数据是一个非常重要的问题. 以仓库发料系统为例,如果某一张领 ...

  9. Farm Irrigation

    题目:Farm Irrigation 题目链接:http://210.34.193.66:8080/vj/Problem.jsp?pid=1494 题目思路:并查集 #include<stdio ...

  10. Day 3-4 函数进阶

    1.名称空间 定义:Python使用叫做命名空间的东西来记录变量的轨迹.命名空间是一个 字典(dictionary) ,它的键就是变量名,它的值就是那些变量的值.是存放变量和值的内存地址的绑定关系的空 ...