POJ 2553 The Bottom of a Graph 【有向图tarjan 求SCC 处理出度为0的SCC】

【题目】http://poj.org/problem?id=2553

【定义】点v是汇点须满足 --- 对图中任意点u，若v可以到达u则必有u到v的路径；若v不可以到达u，则u到v的路径可有可无。

【题意】在n个点m条边的有向图里面，问有多少个点是汇点。

【分析】首先若SCC里面有一个点不是汇点，那么它们全不是汇点，反之也如此。这也就意味着一个SCC里面的点要么全是，要么全不是。在求出SCC并缩点后，任一个编号为A的SCC若存在指向编号为B的SCC的边，那么它里面所有点必不是汇点（因为编号为B的SCC不可能存在指向编号为A的SCC的边）。若编号为A的SCC没有到达其他SCC的路径，那么该SCC里面所有点必是汇点。因此判断的关键在于SCC的出度是否为0.（来源：https://blog.csdn.net/chenzhenyu123456/article/details/46967485）

【思路】tarjan模板题...T了3个小时....结果M改成50005就过了...

【重点】重点应该是tarjan的邻接表的写法吧...

const int M=50005;
int n,m;
int head[M];
class node
{
public:
    int v,next;
};
node edge[M];
int ei;

void tarjan(int x)
{
    ///省略
    for(int i=head[x]; i!=-1; i=edge[i].next)
    {
        int& t=edge[i].v;///i只是第几条边。用的时候不是拿i用，而是用edge[i].v.
        if(!DFN[t])
        {
            tarjan(t);
            LOW[x]=min(LOW[x],LOW[t]);
        }
        else if(DFN[t]&&inthestack[t])
        {
            LOW[x]=min(LOW[x],DFN[t]);
        }
    }
    ///省略
    return;
}

void addedge(int u,int v)///add的写法也很重要啊
{
    edge[++ei].v=v;
    edge[ei].next=head[u];
    head[u]=ei;
    return;
}

【代码】

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int M=50005;
int n,m;
int head[M];
class node
{
public:
    int v,next;
};
node edge[M];
int ei;
int Stack[M],si;
int DFN[M],LOW[M],tot;
int inthestack[M];
int scc;
int belong[M];

void tarjan(int x)
{
    Stack[++si]=x;
    inthestack[x]=1;
    DFN[x]=LOW[x]=++tot;
    for(int i=head[x]; i!=-1; i=edge[i].next)
    {
        int& t=edge[i].v;
        if(!DFN[t])
        {
            tarjan(t);
            LOW[x]=min(LOW[x],LOW[t]);
        }
        else if(DFN[t]&&inthestack[t])
        {
            LOW[x]=min(LOW[x],DFN[t]);
        }
    }
    if(DFN[x]==LOW[x])
    {
        scc++;
        do
        {
            belong[Stack[si]]=scc;
            inthestack[Stack[si]]=0;
            si--;
        }
        while(x!=Stack[si+1]);
    }
    return;
}

void addedge(int u,int v)
{
    edge[++ei].v=v;
    edge[ei].next=head[u];
    head[u]=ei;
    return;
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(inthestack,0,sizeof(inthestack));
    memset(LOW,0,sizeof(LOW));
    memset(DFN,0,sizeof(DFN));
    scc=tot=ei=si=0;
    return;
}


int main()
{

    while(~scanf("%d",&n)&&n)
    {scanf("%d",&m);
        init();
        for(int i=0; i<m; i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            addedge(a,b);
        }
        for(int i=1; i<=n; i++)
        {
            if(!DFN[i])
                tarjan(i);
        }
        
        ///以下是比较暴力的求法，枚举点看是不是那个强连通分量的儿子，然后放到ans里再sort 141MS
        int in[M],out[M];
        for(int i=1; i<=scc; i++)
        {
            in[i]=out[i]=0;
        }
        for(int i=1; i<=n; i++)
        {
            for(int j=head[i]; j!=-1; j=edge[j].next)
            {
                int &t=edge[j].v;
                if(belong[i]!=belong[t])
                {
                    out[belong[i]]++;
                    in[belong[t]]++;
                }
            }
        }
        
        int ans[M],ai=1;
        int flag=0;
        for(int i=1; i<=scc; i++)
        {
            if(out[i]!=0)
                continue;
            for(int j=1; j<=n; j++)
            {
                if(belong[j]==i)
                {
                    flag=1;
                    ans[ai++]=j;
                }

            }
        }
        if(flag==0)
        {
            printf("\n");
            continue;
        }
        sort(ans+1,ans+ai);
        printf("%d",ans[1]);
        for(int i=2; i<ai; i++)
        {
            printf(" %d",ans[i]);
        }
        printf("\n");
    }

}

int main()
{
    while(~scanf("%d",&n)&&n)
    {
        scanf("%d",&m);
        init();
        for(int i=0; i<m; i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            addedge(a,b);
        }
        for(int i=1; i<=n; i++)
        {
            if(!DFN[i])
                tarjan(i);
        }

        ///比刚才的要骚的写法，（大佬太强了），枚举所有的边，如果跨越了两个强连通，那么出来的那头一定不合题意。
        ///标记这个强连通。最后从1到b扫一遍点如果他的强连通没有被标记就要输出。79MS
        bool book[M]={0};
        for(int i=1; i<=n; i++)
        {
            for(int j=head[i]; j!=-1; j=edge[j].next)
            {
                int &t=edge[j].v;
                if(belong[i]!=belong[t])
                {
                    book[belong[i]]=1;
                }
            }
        }
        int flag=0;
        for(int i=1;i<=n;i++)
        {
            if(book[belong[i]]==1)continue;
            if(flag==0){flag=1;printf("%d",i);}
            else printf(" %d",i);
        }
        printf("\n");
    }

}