HDOJ 3639 Hawk-and-Chicken

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3639

此题可能拿到的时候很容易想到直接深搜可能就能解决，但是实际上不行，如果要深搜的话，那我们必须要满足一个条件，那就是不能成环，但是这个题显然是可能有成环的(如样例2)，那怎么办呢，这时候我们可以去借助强连通分量来解决这个问题，因为如果我们用Tarjan算法把所有的强连通找到后，再讲所有强连通分量缩点处理后，每个点就满足不会成环的性质，我们就可以进行DFS深搜了，当然在这之前还得反向建图求入度(可以不用反向建图求出度)，所以总结此题就是强连通+缩点+反向构图+DFS。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <stack>
#include <vector>
#define MAXN 5010 
using namespace std; 
vector<int> g[MAXN],g2[MAXN];
stack<int> st; 
int n,index,num,dfn[MAXN],low[MAXN],belong[MAXN]; 
bool vis[MAXN],instack[MAXN]; 
int ans[MAXN],in[MAXN],scc[MAXN],sum; 
void tarjan(int u)//ÇóÇ¿Á¬Í¨·ÖÖ§
{ 
    int v; 
    dfn[u] = low[u] = index++;  
    st.push(u);
    instack[u] = true; 
    vis[u] = true; 
    for(int i=0; i<g[u].size(); i++) 
    {
        v = g[u][i];  
        if(!vis[v]) 
        {          
            tarjan(v);       
            low[u] = min(low[u], low[v]);      
        }    
        else if(instack[v]) 
            low[u] = min(low[u], dfn[v]);  
    } 
    if(dfn[u] == low[u]) 
    {         
        do 
        {     
            v = st.top();        
            instack[v] = false; 
            st.pop();
            belong[v]=num;
            scc[num]++;
        }while(v != u); 
        num++;
    }
}
void dfs(int u)//ÉîËÑÇó³ö¿Éµ½´ï¸ÃÁ¬Í¨·ÖÖ§µÄµãÊý
{ 
    vis[u]=true; 
    sum+=scc[u];
    for(int i=0; i<g2[u].size(); i++) 
        if(!vis[g2[u][i]]) 
            dfs(g2[u][i]);
}
int main() 
{
    int T,m,a,b,cas=0; 
    scanf("%d",&T); 
    while(T--) 
    {
        scanf("%d %d",&n,&m); 
        for(int i=0;i<=n;i++)
        { 
            g[i].clear();
            g2[i].clear();
        }
        while(m--) 
        {
            scanf("%d %d",&a,&b); 
            g[a].push_back(b);
        }
        index=num=0; 
        memset(vis,false,sizeof(vis)); 
        memset(instack,false,sizeof(instack)); 
        memset(scc,0,sizeof(scc)); 
        for(int i=0; i<n; i++) 
            if(!vis[i]) 
                tarjan(i);
        memset(in,0,sizeof(in)); 
        for(int i=0; i<n; i++)
            for(int j=0; j<g[i].size(); j++) 
            {
                if(belong[i] != belong[g[i][j]]) 
                {
                    g2[belong[g[i][j]]].push_back(belong[i]);
                    in[belong[i]]++; 
                }
            }
        int maxans=-1; 
        memset(ans,-1,sizeof(ans)); 
        printf("Case %d: ",++cas); 
        for(int i=0; i<num; i++) 
        {
            if(in[i]==0) 
            {
                sum=0; 
                memset(vis,false,sizeof(vis)); 
                dfs(i);
                ans[i]=sum;
                if(sum>maxans) 
                    maxans=sum;
            }
        }
        printf("%d\n",maxans-1); 
        int flag=0; 
        for(int i=0; i<n; i++) 
            if( ans[belong[i]] == maxans) 
            { 
                if(!flag)
                {
                    printf("%d", i); 
                    flag=1; 
                }
                else printf(" %d",i); 
            }
        printf("\n"); 
    }
    return 0; 
}