POJ 2553

这道题还是和其他题目相似，没有什么差别和软用

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
#include<vector>
using namespace std;
const int maxn=5000+10;
int n,m;
vector<int> G[maxn];
stack<int> S;
int dfs_clock,scc_cnt;
int pre[maxn],low[maxn],sccno[maxn];
bool in0[maxn],out0[maxn];
void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;
    S.push(u);
    for(int i=0;i<G[u].size();i++)
    {
        int v=G[u][i];
        if(!pre[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!sccno[v])
            low[u]=min(low[u],pre[v]);
    }
    if(low[u]==pre[u])//强连通分量起点
    {
        scc_cnt++;
        while(true)
        {
            int x= S.top(); S.pop();
            sccno[x]=scc_cnt;
            if(x==u) break;
        }
    }
}
void find_scc(int n)
{
    scc_cnt=dfs_clock=0;
    memset(pre,0,sizeof(pre));
    memset(sccno,0,sizeof(sccno));
    for(int i=0;i<n;i++)
        if(!pre[i]) dfs(i);
}

bool on[maxn];

int main()
{
    while(scanf("%d",&n)){
        if(!n){break;}
        scanf("%d",&m);
        for(int i=0;i<n;i++)G[i].clear();
        for(int i=0;i<m;i++){
            int u,v;
            scanf("%d%d",&u,&v);
            G[u-1].push_back(v-1);
        }
        find_scc(n);
        memset(out0+1,true,scc_cnt);
        memset(on,false,n);
        for(int u=0;u<n;u++){
            for(int i=0;i<G[u].size();i++){
                int v=G[u][i];
                if(sccno[u]!=sccno[v]){
                    out0[sccno[u]]=false;
                }
            }
        }
        int cnt=0;
        for(int i=0;i<n;i++){
            if(out0[sccno[i]]){
                on[i]=true;
                cnt++;
            }
        }
        int flag=0;
        for(int i=0;i<n;i++){
            if(on[i]){
                flag++;
                if(flag!=cnt)
                    printf("%d ",i+1);
                else printf("%d\n",i+1);
            }
        }
        if(!cnt)
            printf("\n");
    }
    return 0;
}