poj2553

题目大意：如果v点能够到的点，反过来能够到达v点，则称这个点为sink点，输出所有的sink点

解题思路：求连通分量，然后出度为0的连通分量里面的点就是sink点

一定要注意初始化。

代码：

#include <iostream>
#include <cstdio>
#include <stack>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=5050;

struct Side{
    int to,next;
}side[maxn*maxn];

bool vis[maxn];
int node[maxn],belong[maxn];
int low[maxn],dfn[maxn];
int out[maxn];
stack<int>s;
int n,m;
int t,tp,id;
void add_side(int u,int v){
    side[t].to=v;
    side[t].next=node[u];
    node[u]=t++;
}
void tarjan(int u){
    low[u]=dfn[u]=tp++;
    vis[u]=true;
    s.push(u);
    for(int i=node[u];i!=-1;i=side[i].next){
        int v=side[i].to;
        if(dfn[v]==-1){
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(vis[v]){
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(low[u]==dfn[u]){
        int v;
        do{
            v=s.top();s.pop();
            belong[v]=id;
            vis[v]=false;
        }while(u!=v);
        id++;
    }
}

int main(){
    while(~scanf("%d",&n)&&n){
        scanf("%d",&m);
        memset(node,-1,sizeof(node));
        memset(dfn,-1,sizeof(dfn));
        memset(low,0,sizeof(low));
        memset(vis,false,sizeof(vis));
        memset(out,0,sizeof(out));
        t=0;tp=0;id=0;
        while(m--){
            int a,b;
            scanf("%d%d",&a,&b);
            add_side(a,b);
        }

        for(int i=1;i<=n;i++)
            if(dfn[i]==-1)
                tarjan(i);
        for(int i=1;i<=n;i++){
            for(int j=node[i];j!=-1;j=side[j].next){
                int v=side[j].to;
                if(belong[i]!=belong[v])
                    out[belong[i]]++;
            }
        }
        int flag=1;
        for(int i=1;i<=n;i++){
            if(!out[belong[i]]){
                if(flag){
                    printf("%d",i);
                    flag=0;
                }
                else{
                    printf(" %d",i);
                }
            }
        }
        cout<<endl;
    }
    return 0;
}