Advertisement
Galebickosikasa

Untitled

Jul 13th, 2020
66
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.68 KB | None | 0 0
  1. #pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
  2. #pragma GCC target("sse,sse2,sse3,ssse3,sse4,sse4.1,sse4.2,popcnt,abm,mmx,avx")
  3. // #pragma comment(linker, "/stack:200000000"]
  4.  
  5. #include <iostream>
  6. #include <vector>
  7. #include <cmath>
  8. #include <algorithm>
  9. #include <unordered_set>
  10. #include <unordered_map>
  11. #include <set>
  12. #include <map>
  13. #include <queue>
  14. #include <deque>
  15. #include <bitset>
  16. #include <stack>
  17. #include <random>
  18. #include <fstream>
  19. #include <sstream>
  20. #include <chrono>
  21.  
  22. #define fi first
  23. #define se second
  24. #define pb push_back
  25. #define ll long long
  26. #define ld long double
  27. #define hm unordered_map
  28. #define pii pair<int, int>
  29. #define sz(a) (int)a.size()
  30. #define all(a) a.begin(), a.end()
  31. #define cinv(v) for (auto& x: v) cin >> x
  32. #define fr(i, n) for (int i = 0; i < n; ++i)
  33. #define fl(i, l, n) for (int i = l; i < n; ++i)
  34.  
  35. #define int ll
  36.  
  37. using namespace std;
  38.  
  39. #ifdef __LOCAL
  40. #define dbg(x) cerr << #x << " : " << x << '\n'
  41. const int maxn = 20;
  42. #else
  43. #define dbg(x)
  44. const int maxn = 500 + 20;
  45. #endif
  46.  
  47. //tg: @galebickosikasa
  48.  
  49. const ll inf = (ll) 2e9;
  50. const ld pi = asin (1) * 2;
  51. const ld eps = 1e-8;
  52. const ll mod = (ll)1e9 + 7;
  53. const ll ns = 97;
  54.  
  55. mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
  56.  
  57. struct edge {
  58. int to, c, f;
  59. };
  60.  
  61. vector<edge> edges;
  62. vector<int> g[maxn];
  63. int it[maxn], d[maxn];
  64.  
  65. int bfs (int n, int x) {
  66. for (auto& x: d) x = inf;
  67. d[0] = 0;
  68. queue<int> a;
  69. a.push (0);
  70. while (!a.empty ()) {
  71. int v = a.front ();
  72. a.pop ();
  73. for (auto& u: g[v]) {
  74. auto& e = edges[u];
  75. if (d[e.to] == inf && e.c - e.f >= x) {
  76. d[e.to] = d[v] + 1;
  77. a.push (e.to);
  78. }
  79. }
  80. }
  81. return d[n - 1] < inf;
  82. }
  83.  
  84. int dfs (int v, int t, int x, int c = inf, int p = -1) {
  85. if (v == t) return c;
  86. while (it[v] < sz (g[v])) {
  87. auto& e = edges[g[v][it[v]]];
  88. if (e.to == p || e.c - e.f < x || d[v] + 1 != d[e.to]) {
  89. ++it[v];
  90. continue;
  91. }
  92. int delta = dfs (e.to, t, x, min (c, e.c - e.f), v);
  93. if (delta > 0) {
  94. auto& _e = edges[g[v][it[v]]^1];
  95. e.f += delta;
  96. _e.f -= delta;
  97. return delta;
  98. }
  99. ++it[v];
  100. }
  101. return 0;
  102. }
  103.  
  104. signed main () {
  105. ios_base::sync_with_stdio(false);
  106. cin.tie(nullptr);
  107. cout.tie(nullptr);
  108. int n, m;
  109. cin >> n >> m;
  110. fr (i, m) {
  111. int a, b, c;
  112. cin >> a >> b >> c;
  113. --a, --b;
  114. g[a].pb (sz (edges));
  115. edges.pb ({b, c, 0});
  116. g[b].pb (sz (edges));
  117. edges.pb ({a, 0, 0});
  118. }
  119. for (int x = 30; x >= 0; --x) {
  120. while (bfs (n, (1<<x))) {
  121. while (dfs (0, n - 1, (1<<x))) {}
  122. for (auto& y: it) y = 0;
  123. }
  124. }
  125. int f = 0;
  126. for (auto& x: g[0]) {
  127. auto& e = edges[x];
  128. f += e.f;
  129. }
  130. cout << f << '\n';
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement