hdu3639 强连通

题意：有 n 个人，m 组支持关系，已知支持关系可以传递，比如 A 支持 B，则所有支持 A 的人也同时支持 B，问哪些人获得的支持数最多，最多获得多少支持（自己不能获得自己的支持）。

首先，如果一些人他们互相支持，那么他们的支持数肯定都是一样的，所有支持他们其中一个人的也同时支持他们所有人，所以对于这些人我们可以强连通缩点。然后就获得了一个有向无环图，每个强连通分量有各自的人数。由于支持的人数可以传递，所以某个人支持了别人，那么他获得的支持数就一定不如他支持的人多，那么其实我们需要获得的就只是那些没有支持过别人的人，所以我们可以通过反向建图，那些没有支持别人的人就是入度为 0 的点，从每个这样的点开始DFS下去统计一共有多少个人是他能够到达的，然后再加上他自己的人数 - 1（不能获得自己的支持），这样我们再从中找出最大值以及这些人就行了。

 #include<stdio.h>
 #include<string.h>
 #include<stack>
 #include<queue>
 using namespace std;

 const int maxn=;
 const int maxm=;

 int head[][maxn],point[][maxm],nxt[][maxm],size[];
 int n,t,scccnt,maxx;
 int stx[maxn],low[maxn],scc[maxn];
 int dp[maxn],id[maxn],vis[maxn];
 stack<int>S;

 int max(int a,int b){return a>b?a:b;}

 void init(){
     memset(head,-,sizeof(head));
     size[]=size[]=;
     memset(dp,,sizeof(dp));
     memset(id,,sizeof(id));
     maxx=;
 }

 void add(int a,int b,int c=){
     if(c){
         for(int i=head[][a];~i;i=nxt[][i])if(point[][i]==b)return;
     }
     point[c][size[c]]=b;
     nxt[c][size[c]]=head[c][a];
     head[c][a]=size[c]++;
 }

 void dfs(int s){
     stx[s]=low[s]=++t;
     S.push(s);
     for(int i=head[][s];~i;i=nxt[][i]){
         int j=point[][i];
         if(!stx[j]){
             dfs(j);
             low[s]=min(low[s],low[j]);
         }
         else if(!scc[j]){
             low[s]=min(low[s],stx[j]);
         }
     }
     if(low[s]==stx[s]){
         scccnt++;
         while(){
             int u=S.top();S.pop();
             scc[u]=scccnt;
             dp[scccnt]++;
             if(s==u)break;
         }
     }
 }

 void setscc(){
     memset(stx,,sizeof(stx));
     memset(scc,,sizeof(scc));
     t=scccnt=;
     for(int i=;i<n;++i)if(!stx[i])dfs(i);
     for(int i=;i<n;++i){
         for(int j=head[][i];~j;j=nxt[][j]){
             int k=point[][j];
             if(scc[i]!=scc[k]){
                 add(scc[i],scc[k],);
                 id[scc[k]]++;
             }
         }
     }
 }

 int Dp(int s){
     int ans=;
     vis[s]=;
     for(int i=head[][s];~i;i=nxt[][i]){
         if(!vis[point[][i]])ans+=Dp(point[][i]);
     }
     return ans+dp[s];
 }

 int main(){
     int T;
     scanf("%d",&T);
     for(int q=;q<=T;q++){
         int m;
         scanf("%d%d",&n,&m);
         init();
         while(m--){
             int a,b;
             scanf("%d%d",&a,&b);
             add(b,a);
         }
         setscc();
         for(int i=;i<=scccnt;++i)if(!id[i]){
             memset(vis,,sizeof(vis));
             dp[i]=Dp(i)-;
             if(dp[i]>maxx)maxx=dp[i];
         }
         int ans=;
         printf("Case %d: %d\n",q,maxx);
         int cnt=;
         for(int i=;i<n;++i){
             if(!id[scc[i]]&&dp[scc[i]]==maxx){
                 if(cnt++)printf(" ");
                 printf("%d",i);
             }
         }
         printf("\n");
     }
     return ;
 }