Advertisement
Guest User

Untitled

a guest
Jul 24th, 2014
414
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.16 KB | None | 0 0
  1. #include <iostream>
  2. #include <cmath>
  3. #include <algorithm>
  4. #include <cstdio>
  5. #include <string>
  6. #include <vector>
  7. #include <set>
  8. #include <cstring>
  9. #include <cstdlib>
  10. #include <map>
  11. #include <ctime>
  12. #include <queue>
  13. #include <bitset>
  14.  
  15. using namespace std;
  16.  
  17. #define f first
  18. #define s second
  19. #define pb push_back
  20. #define mp make_pair
  21. #define forit(s) for(__typeof(s.begin()) it = s.begin(); it != s.end(); it ++)
  22. #define forn(i, l, r) for(int i = l; i < r; i ++)
  23. #define pii pair <int, int>
  24. #define vi vector <int>
  25. #define N 1005
  26. #define y1 stupid_cmath
  27. #define sz(a) (int)a.size()
  28. #define ll long long
  29. #define all(a) a.begin(), a.end()
  30.  
  31. const int inf = (int)1e9;
  32. const int mod = (int)1e9 + 7;
  33. const double pi = acos(-1.0);
  34. const double eps = 1e-9;
  35.  
  36. int n, m, col[N];
  37. vi g[N], gr[N], order;
  38. pii e[10 * N];
  39. int used[N], mt[N], c, sz, ver = 1;
  40. int was[N][N];
  41. vector <bitset <N> > b;
  42.  
  43. void clear() {
  44. b.clear();
  45. for(int i = 0; i < N; i++) {
  46. g[i].clear();
  47. gr[i].clear();
  48. }
  49. order.clear();
  50. memset(col, 0, sizeof col);
  51. sz = 0;
  52. c = 1;
  53. memset(mt, -1, sizeof mt);
  54. memset(used, 0, sizeof used);
  55. }
  56.  
  57. void dfs1(int v) {
  58. used[v] = 1;
  59. forit(gr[v]) {
  60. if(!used[*it]) dfs1(*it);
  61. }
  62. order.pb(v);
  63. }
  64.  
  65. void dfs2(int v) {
  66. used[v] = 1;
  67. col[v] = sz;
  68. forit(g[v]) {
  69. if(!used[*it]) dfs2(*it);
  70. }
  71. }
  72.  
  73. void dfs3(int v) {
  74. used[v] = 1;
  75. forit(gr[v]) {
  76. int to = *it;
  77. if(used[to]) continue;
  78. dfs3(to);
  79. }
  80. order.pb(v);
  81. }
  82.  
  83. void read() {
  84. scanf("%d%d", &n, &m);
  85. for(int i = 0, x, y; i < m; i++) {
  86. scanf("%d%d", &x, &y);
  87. x--; y--;
  88. e[i] = (mp(x, y));
  89. g[x].pb(y); gr[y].pb(x);
  90. }
  91. for(int i = 0; i < n; i++)
  92. if(!used[i]) dfs1(i);
  93. reverse(all(order));
  94. memset(used, 0, sizeof used);
  95. for(int i = 0; i < n; i++) {
  96. int v = order[i];
  97. if(!used[v]) {
  98. dfs2(v);
  99. sz++;
  100. }
  101. }
  102.  
  103. for(int i = 0; i < N; i++) {gr[i].clear(); g[i].clear();}
  104.  
  105. for(int i = 0; i < m; i++) {
  106. if(col[e[i].f] != col[e[i].s]) {
  107. int x = col[e[i].f], y = col[e[i].s];
  108. if(was[x][y] == ver) continue;
  109. was[x][y] = ver;
  110. gr[x].pb(y);
  111. }
  112. }
  113.  
  114.  
  115.  
  116.  
  117. order.clear();
  118. memset(used, 0, sizeof used);
  119. for(int i = 0; i < sz; i++) {
  120. if(!used[i]) dfs3(i);
  121. }
  122. b.assign(sz, bitset <N>(0));
  123. for(int i = 0; i < sz; i++) {
  124. int v = order[i];
  125. b[v][v] = 1;
  126. forit(gr[v]) {
  127. b[v] |= b[*it];
  128. }
  129. }
  130. for(int i = 0; i < sz; i++) {
  131. for(int j = 0; j < sz; j++) if(b[i][j] == 1 && i != j) g[i].pb(j);
  132. }
  133. }
  134.  
  135. bool dfs(int v) {
  136. if(used[v] == c) return false;
  137. used[v] = c;
  138. forit(g[v]) {
  139. int i = *it;
  140. if(mt[i] == -1) {
  141. mt[i] = v;
  142. return true;
  143. }
  144. }
  145. forit(g[v]) {
  146. int i = *it;
  147. if((mt[i] == -1 || dfs(mt[i]))) {
  148. mt[i] = v;
  149. return true;
  150. }
  151. }
  152. return false;
  153. }
  154.  
  155. void solve(int num) {
  156. printf("Case %d: ", num);
  157.  
  158. clear();
  159. read();
  160.  
  161. int ans = 0;
  162. memset(used, 0, sizeof used);
  163. for(int i = 0; i < sz; i++) {
  164. ans += dfs(i);
  165. c++;
  166. }
  167.  
  168. printf("%d\n", sz - ans);
  169. ver++;
  170. }
  171.  
  172. int main () {
  173.  
  174. #ifdef LOCAL
  175. freopen("a.in", "r", stdin);
  176. freopen("a.out", "w", stdout);
  177. #endif
  178.  
  179. int T;
  180. scanf("%d", &T);
  181. for(int i = 1; i <= T; i++) solve(i);
  182.  
  183. return 0;
  184. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement