hdu3639 Hawk-and-Chicken（强联通分量+反建图dfs）

Hawk-and-Chicken

题目大意：给一个有向图，如果u是v的粉丝，v是w的粉丝，那么

u是w的粉丝，求最大粉丝的数量和点。

可以想到最大粉丝点肯定是存在在出度为零的强联通分量中，因为如果不是这样肯定有另一个他通往的点比他的粉丝多，而求最大粉丝点这里，用到了反建图，因为方便dfs来搜索每个点，搜索过程中，每个点首先有他自己强联通分量内的n-1个点以及他到达的强联通分量的n个点，然后搜索即可。

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <stack>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=5e3+5;
vector<int>a[maxn];
vector<int>p[maxn];
int t,n,m,u,v,tmp,cnt,tol,vis[maxn],dfn[maxn],low[maxn],id[maxn],out[maxn],num[maxn],ans[maxn];
stack<int>s;
void tarjan(int x)//targan算法缩点，把所有点缩成强联通分量，以强联通分量建图
{
    dfn[x]=low[x]=(++tol);
    vis[x]=1;
    s.push(x);
    for(int i=0; i<a[x].size(); i++)
    {
        if(!dfn[a[x][i]])tarjan(a[x][i]),low[x]=min(low[x],low[a[x][i]]);
        else if(vis[a[x][i]])low[x]=min(low[x],dfn[a[x][i]]);
    }
    if(low[x]==dfn[x])
    {
        cnt++;
        while(1)
        {
            int ss=s.top();
            s.pop();
            id[ss]=cnt,vis[ss]=0;
            if(ss==x)break;
        }
    }
}
void init()//初始化
{
    memset(vis,0,sizeof(vis));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(out,0,sizeof(out));
    memset(num,0,sizeof(num));
    memset(ans,0,sizeof(ans));
}
void dfs(int x)//搜索建完后的图，tmp代表粉丝数
{
    for(int i=0; i<p[x].size(); i++)
    {
        int y=p[x][i];
        if(!vis[y])tmp+=num[y],vis[y]=1;
        dfs(y);
    }
}
int main()
{
    ios::sync_with_stdio(false);
    while(cin>>t)
    {
        int ca=0;
        while(t--)
        {
            cin>>n>>m;
            init();
            for(int i=1; i<=m; i++)
            {
                cin>>u>>v;
                u++,v++;//嫌麻烦直接加了个1
                a[u].push_back(v);
            }
            cnt=tol=0;
            for(int i=1; i<=n; i++)
                if(!dfn[i])
                    tarjan(i);
            for(int i=1; i<=n; i++)num[id[i]]++;
            for(int i=1; i<=n; i++)
                for(int j=0; j<a[i].size(); j++)
                    if(id[i]!=id[a[i][j]])p[id[a[i][j]]].push_back(id[i]),out[id[i]]++;//对每个点所在的强联通分量反向建图，出度加一
            int maxn=0;
            for(int i=1; i<=cnt; i++)
            {
                if(out[i]==0)//对出度为0的分量搜索
                {
                    memset(vis,0,sizeof(vis));
                    tmp=num[i]-1;
                    dfs(i);
                    ans[i]=tmp;
                    if(maxn<tmp)maxn=tmp;
                }
                else maxn=max(maxn,num[i]-1),ans[i]=num[i]-1;
            }
            cout<<"Case "<<++ca<<": "<<maxn<<endl;
            int fl=0;
            for(int i=1; i<=n; i++)
            {
                if(ans[id[i]]==maxn&&fl==0)cout<<i-1,fl=1;
                else if(ans[id[i]]==maxn&&fl==1)cout<<" "<<i-1;//记得处理，否则pe
            }
            cout<<endl;
            for(int i=1; i<=n; i++)a[i].clear(),p[i].clear();
        }
    }
    return 0;
}