HDU - 3639-D - Hawk-and-Chicken（强连通分量+缩点+重新建图）

题意：有n个孩子，m对关系，关系a b代表孩子a支持孩子b。求被支持数最多的孩子并输出

思路：

*很明显，由于边是单向的，对于在同一个强连通分量里点，它们的支持数是一样的。所以可以先求强连通分量然后缩点

*缩点后重点是怎么求每个点的支持数，这个是比较麻烦的。重点关注出度为0的点，若出度不为0，则它指向的那个点肯定比当前点更优。

*所以问题归结为求出度为0的点的支持数即可。重新构图，建立反向边（方便求出出度为0的点的支持数）

附：本题可以学习的东西感觉还是挺多的。对于缩点后的图，求出度入度的时候，其实只有出度和入度为0才是真正可以判断的（因为同一个分量里面可能会重复指向另一个分量，这样此分量的出度和另一个分量的入度都会增加，并且缩点后来看，其实都是同一条边，即重边）。所以只有出入度为0是比较”有意义“的。

另外，建立反向边这个思路也值得学习，此题建立后跑一遍dfs边可以求出某个点的支持数（有多少个点是指向这个点的，包括传递的点）。重构的图是有重边的（即一个分量可能有多条指向另一个分量的边，这点可以从代码中看出），但是不影响dfs，因为只要把跑过的子节点标记一下，下次就不会选它了，其他重边也就没用了

代码：

//#include<bits/stdc++.h>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<stack>
#include<vector>
#define rep(i,e) for(int i=0;i<(e);++i)
#define rep1(i,e) for(int i=1;i<=(e);++i)
#define repx(i,x,e) for(int i=(x);i<=(e);++i)
#define pii pair<int,int>
#define X first
#define Y second
#define PB push_back
#define MP make_pair
#define mset(var,val) memset(var,val,sizeof(var))
#define scd(a) scanf("%d",&a)
#define scdd(a,b) scanf("%d%d",&a,&b)
#define scddd(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define IOS ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

#ifdef LOCAL
template<typename T>
void dbg(T t){
    cout<<t<<" "<<endl;
}
template<typename T, typename... Args>
void dbg(T t, Args... args){
    cout<<t<<" ";dbg(args...);
}


#else
#define dbg(...)
#endif // local
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3fll;
const int mod = 1e9+7;
const int N = 500+10;
typedef long long ll;
const int maxn = 5e3+10;
const int maxm = 3e4+10;


struct node{
    int v,next;
}edge[maxm],nedge[maxm];
int cnt,ncnt;
int head[maxn];
int nhead[maxn];
int low[maxn],dfn[maxn],dfn_cnt,sccno[maxn],scc_cnt;
stack<int>S;
int outdeg[maxn];
int num[maxn];
int res[maxn];
bool vis[maxn];
void add(int u,int v)
{
    edge[++cnt].v = v;
    edge[cnt].next = head[u];
    head[u] = cnt;
}

void nadd(int u,int v)
{
    nedge[++ncnt].v = v;
    nedge[ncnt].next = nhead[u];
    nhead[u] = ncnt;
}

void tarjan(int u)
{
    low[u] = dfn[u] = ++dfn_cnt;
    S.push(u);
    for(int i = head[u]; ~i; i = edge[i].next)
    {
        int v = edge[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[u] = min(low[u],low[v]);
        }else if(!sccno[v])
        {
            low[u] = min(low[u],dfn[v]);
        }
    }
    if(dfn[u] == low[u])
    {
        ++scc_cnt;
        for(;;){
            int x = S.top();S.pop();
            sccno[x] = scc_cnt;
            if(x == u) break;
        }
    }
}
int sum;
void dfs(int u)
{
    vis[u] = true;
    sum += num[u];
    for(int i = nhead[u]; ~i; i = nedge[i].next)
    {
        int v = nedge[i].v;
        if(!vis[v])
            dfs(v);
    }
}

void work()
{
    int t,n,m,x,y,k=1;
    cin>>t;
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i = 0; i <= n; i++)
        {
            low[i] = dfn[i] = sccno[i] = outdeg[i] = num[i] = res[i] = 0;
            head[i] = nhead[i] = -1;
        }
        scc_cnt = dfn_cnt = cnt = ncnt = 0;
        for(int i = 0; i < m; i++)
        {
            scanf("%d%d",&x,&y);
            add(x,y);
        }

        for(int i = 0; i < n; i++)
        {
            if(!dfn[i])
                tarjan(i);
        }
        //dbg(scc_cnt);
        for(int i = 0; i < n; i++)
        {
            num[sccno[i]]++;
        }
        for(int u = 0; u < n; u++)
        {
            for(int i = head[u]; ~i; i = edge[i].next)
            {
                int v = edge[i].v;
                if(sccno[v] != sccno[u])
                {
                    nadd(sccno[v],sccno[u]);
                    outdeg[sccno[u]]++;
                }
            }
        }

        int mx = 0;
        for(int i = 1; i <= scc_cnt; i++)
        {
            if(!outdeg[i])
            {
                memset(vis,0,sizeof(vis));
                sum = 0;
                dfs(i);
                sum--;
                res[i] = sum;
                mx = max(mx,sum);
            }
        }
        printf("Case %d: %d\n",k++,mx);
        bool tag = true;
        for(int i = 0; i < n; i++)
        {
            if(res[sccno[i]] == mx)
            {
                if(tag)
                    printf("%d",i),tag = false;
                else printf(" %d",i);
            }
        }
        printf("\n");
    }
}

int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    //IOS;
    work();
    return 0;
}