HDOJ3639解题报告【缩点+dfs】

题目地址：

　　http://acm.hdu.edu.cn/showproblem.php?pid=3639

题目概述：

　　题面有点翻译不来……可以用Chrome自带的那个翻译。

大致思路：

　　如果跑bfs或者dfs遇到环的话就会出问题，所以先跑一遍tarjan求强连通分量，把环的情况去掉，在剩下的图上做。

　　这个时候发现样例1里2这个节点指向的两个节点都是答案，所以反向建图，用dfs求出所有入度为0的点的子树大小，最大数即为答案。

　　输出节点的时候直接对0~n-1循环找即可。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
#include <map>
#include <stack>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define sacnf scanf
#define scnaf scanf
#define maxn  5010
#define maxm 26
#define inf 1061109567
#define Eps 0.00001
const double PI=acos(-1.0);
#define mod 7
#define MAXNUM 10000
void Swap(int &a,int &b) {int t=a;a=b;b=t;}
int Abs(int x) {return (x<0)?-x:x;}
typedef long long ll;
typedef unsigned int uint;

vector<int> G[maxn];
vector<int> rG[maxn];
stack<int> S;

struct node
{
    int pos,val;
    bool operator < (const node &a) const
    {
        return val<a.val;
    }
};

int pre[maxn],low[maxn],scc[maxn],scc_val[maxn];
int scc_cnt,dfs_clock;
int n,m;
int ans[maxn],in[maxn],cnt[maxn],vis[maxn];

void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;
    int len=G[u].size();S.push(u);
    for(int i=0;i<len;i++)
    {
        int v=G[u][i];
        if(!pre[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!scc[v]) low[u]=min(low[u],pre[v]);
    }
    if(pre[u]==low[u])
    {
        scc_cnt++;int cnt=0;
        for(;;)
        {
            int t=S.top();S.pop();
            scc[t]=scc_cnt;cnt++;
            if(t==u) break;
        }
        scc_val[scc_cnt]=cnt;
    }
}

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

int Count(int u)
{
    vis[u]=1;
    int len=rG[u].size(),sum=scc_val[u];
    for(int i=0;i<len;i++)
    {
        int v=rG[u][i];
        if(!vis[v]) sum+=Count(v);
    }
    return sum;
}

int main()
{
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    //clock_t st=clock();
    int T;scanf("%d",&T);
    for(int kase=1;kase<=T;kase++)
    {
        int a,b;scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            G[i].clear();
        }
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&a,&b);
            a++;b++;G[a].push_back(b);
        }
        find_scc();

        ///缩点
        for(int i=1;i<=scc_cnt;i++)
        {
            in[i]=0;cnt[i]=-1;
            rG[i].clear();
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=scc_cnt;j++) vis[j]=0;
            int len=G[i].size();
            for(int k=0;k<len;k++)
            {
                int j=G[i][k];
                if(scc[i]!=scc[j]&&!vis[scc[j]])
                {
                    vis[scc[j]]=1;rG[scc[j]].push_back(scc[i]);
                    in[scc[i]]++;
                }
            }
        }

        int ans=-1;bool flag=false;
        for(int i=1;i<=scc_cnt;i++)
        {
            if(!in[i])
            {
                for(int j=1;j<=scc_cnt;j++) vis[j]=0;
                cnt[i]=Count(i);ans=max(ans,cnt[i]);
            }
        }
        printf("Case %d: %d\n",kase,ans-1);
        for(int i=1;i<=n;i++)
        {
            if(ans==cnt[scc[i]])
            {
                if(flag) printf(" ");
                printf("%d",i-1);flag=true;
            }
        }
        printf("\n");
    }
    //clock_t ed=clock();
    //printf("\n\nTime Used : %.5lf Ms.\n",(double)(ed-st)/CLOCKS_PER_SEC);
    return 0;
}