HDU 1827 - Summer Holiday (强连通 + 缩点)

来源:互联网 发布:mac nginx lua 安装 编辑:程序博客网 时间:2024/05/26 09:56

思路

缩点之后有几个入度为0的就要通知几次,然后在那个强连通分量里找最便宜的打过去。

代码

  1. #include <cstdio>
  2. #include <stack>
  3. #include <set>
  4. #include <iostream>
  5. #include <string>
  6. #include <vector>
  7. #include <queue>
  8. #include <functional>
  9. #include <cstring>
  10. #include <algorithm>
  11. #include <cctype>
  12. #include <string>
  13. #include <map>
  14. #include <iomanip>
  15. #include <cmath>
  16. #define LL long long
  17. #define ULL unsigned long long
  18. #define SZ(x) (int)x.size()
  19. #define Lowbit(x) ((x) & (-x))
  20. #define MP(a, b) make_pair(a, b)
  21. #define MS(arr, num) memset(arr, num, sizeof(arr))
  22. #define PB push_back
  23. #define F first
  24. #define S second
  25. #define ROP freopen("input.txt", "r", stdin);
  26. #define MID(a, b) (a + ((b - a) >> 1))
  27. #define LC rt << 1, l, mid
  28. #define RC rt << 1|1, mid + 1, r
  29. #define LRT rt << 1
  30. #define RRT rt << 1|1
  31. #define BitCount(x) __builtin_popcount(x)
  32. #define BitCountll(x) __builtin_popcountll(x)
  33. #define LeftPos(x) 32 - __builtin_clz(x) - 1
  34. #define LeftPosll(x) 64 - __builtin_clzll(x) - 1
  35. const double PI = acos(-1.0);
  36. const int INF = 0x3f3f3f3f;
  37. using namespace std;
  38. const double eps = 1e-6;
  39. const int MAXN = 1100 + 10;
  40. const int MOD = 1000007;
  41. typedef pair<int, int> pii;
  42. typedef vector<int>::iterator viti;
  43. typedef vector<pii>::iterator vitii;
  44. int fee[MAXN];
  45. struct SCC
  46. {
  47. int v, pre[MAXN], lowlink[MAXN], sccno[MAXN], scc_cnt, dfs_clock, in[MAXN], out[MAXN];
  48. vector<int> G[MAXN], scc[MAXN];
  49. stack<int> st;
  50. void init(int v)
  51. {
  52. this->v = v;
  53. for (int i = 1; i <= v; i++) G[i].clear(), scc[i].clear();
  54. MS(sccno, 0); MS(pre, 0); MS(lowlink, 0); MS(in, 0);
  55. scc_cnt = dfs_clock = 0;
  56. }
  57. void DFS(int u)
  58. {
  59. pre[u] = lowlink[u] = ++dfs_clock;
  60. st.push(u);
  61. for (int i = 0; i < SZ(G[u]); i++)
  62. {
  63. int v = G[u][i];
  64. if (!pre[v])
  65. {
  66. DFS(v);
  67. lowlink[u] = min(lowlink[u], lowlink[v]);
  68. }
  69. else if (!sccno[v]) lowlink[u] = min(lowlink[u], pre[v]);
  70. }
  71. if (pre[u] == lowlink[u])
  72. {
  73. ++scc_cnt;
  74. while (1)
  75. {
  76. int x = st.top(); st.pop();
  77. sccno[x] = scc_cnt;
  78. if (x == u) break;
  79. }
  80. }
  81. }
  82. void solve()
  83. {
  84. for (int i = 1; i <= v; i++)
  85. {
  86. scc[sccno[i]].PB(i);
  87. for (int j = 0; j < SZ(G[i]); j++)
  88. {
  89. int x = G[i][j];
  90. if (sccno[i] != sccno[x])
  91. {
  92. out[sccno[i]]++;
  93. in[sccno[x]]++;
  94. }
  95. }
  96. }
  97. int ans = 0, cnt = 0;
  98. for (int i = 1; i <= scc_cnt; i++)
  99. {
  100. int tmp = INF;
  101. if (in[i] == 0)
  102. {
  103. cnt++;
  104. for (int j = 0; j < SZ(scc[i]); j++)
  105. tmp = min(tmp, fee[scc[i][j]]);
  106. ans += tmp;
  107. }
  108. }
  109. printf("%d %d\n", cnt, ans);
  110. }
  111. void find_scc()
  112. {
  113. for (int i = 1; i <= v; i++)
  114. if (!pre[i]) DFS(i);
  115. }
  116. }scc;
  117. int main()
  118. {
  119. //ROP;
  120. int v, e, i, j;
  121. while (~scanf("%d%d", &v, &e))
  122. {
  123. scc.init(v);
  124. for (i = 1; i <= v; i++) scanf("%d", &fee[i]);
  125. for (i = 0; i < e; i++)
  126. {
  127. int a, b;
  128. scanf("%d%d", &a, &b);
  129. scc.G[a].PB(b);
  130. }
  131. scc.find_scc();
  132. scc.solve();
  133. }
  134. return 0;
  135. }
0 0