HDU 3639 SCC Hawk-and-Chicken

求SCC缩点，统计出每个SCC中的点的个数。

然后统计能到达u的最多的点的个数，可以反向建图，再dfs一遍统计出来。

最后说一下，有必要开一个标记数组，因为测试数据中有重边，结果无限WA。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #include <vector>
 #include <stack>
 using namespace std;

 const int maxn =  + ;
 int n, m;
 vector<int> G[maxn];

 stack<int> S;
 int pre[maxn], lowlink[maxn], sccno[maxn], sz[maxn], dfs_clock, scc_cnt;

 void dfs(int u)
 {
     pre[u] = lowlink[u] = ++dfs_clock;
     S.push(u);

     for(int i = ; i < G[u].size(); i++)
     {
         int v = G[u][i];
         if(!pre[v])
         {
             dfs(v);
             lowlink[u] = min(lowlink[u], lowlink[v]);
         }
         else if(!sccno[v]) lowlink[u] = min(lowlink[u], pre[v]);
     }

     if(pre[u] == lowlink[u])
     {
         scc_cnt++;
         for(;;)
         {
             int x = S.top(); S.pop();
             sccno[x] = scc_cnt;
             sz[scc_cnt]++;
             if(x == u) break;
         }
     }
 }

 void find_scc()
 {
     dfs_clock = scc_cnt = ;
     memset(pre, , sizeof(pre));
     memset(sccno, , sizeof(sccno));
     memset(sz, , sizeof(sz));
     for(int i = ; i < n; i++) if(!pre[i]) dfs(i);
 }

 int support[maxn];
 int indeg[maxn];
 bool vis[maxn];
 vector<int> G2[maxn];

 int dfs2(int u)
 {
     int ans = sz[u];
     vis[u] = true;
     for(int i = ; i < G2[u].size(); i++)
     {
         int v = G2[u][i];
         if(vis[v]) continue;
         ans += dfs2(v);
     }
     return ans;
 }

 int main()
 {
     int T; scanf("%d", &T);
     for(int kase = ; kase <= T; kase++)
     {
         scanf("%d%d", &n, &m);
         for(int i = ; i < n; i++) G[i].clear();
         while(m--)
         {
             int u, v; scanf("%d%d", &u, &v);
             G[u].push_back(v);
         }

         find_scc();

         memset(indeg, , sizeof(indeg));
         for(int i = ; i <= scc_cnt; i++) G2[i].clear();
         for(int i = ; i < n; i++)
             for(int j = ; j < G[i].size(); j++)
             {
                 int u = sccno[i], v = sccno[G[i][j]];
                 if(u == v) continue;
                 indeg[u]++;
                 G2[v].push_back(u);
             }

         int ans = ;
         for(int i = ; i <= scc_cnt; i++) if(!indeg[i])
         {
             memset(vis, false, sizeof(vis));
             support[i] = dfs2(i);
             ans = max(ans, support[i]);
         }

         printf("Case %d: %d\n", kase, ans - );
         bool flag = false;
         for(int i = ; i < n; i++)
         {
             int u = sccno[i];
             if(!indeg[u] && support[u] == ans)
             {
                 if(flag) printf(" ");
                 printf("%d", i);
                 flag = true;
             }
         }
         puts("");
     }

     return ;
 }

代码君