poj2553 求汇点

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

其实就是求联通分量罢了。

想要互相到达，就是说明v这点所在的联通分量没有出度，如果有的话，v点一定可以顺着那条出度的边指向某个点G，但是G是没法到达v的。如果可以到达，那么又构成联通了，G会归与这个联通图的。所以说，找所有出度为0的连通图所包含的点即可

注意点1：排序输出

注意点2：不要用数组存图，内存会超

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#include<stack>
using namespace std;
const int MAXN=5100;
int n;
int DFN[MAXN];
int LOW[MAXN];
int vis[MAXN];
int belong[MAXN];//belong[i]表示i属于缩点后的哪个节点
int cnt;
int out[MAXN];
int tot;
struct Edge
{
      int v;
      int next;
}edge[MAXN*MAXN];
int edgecount;
int head[MAXN];
void Init()
{
      edgecount=0;
      memset(head,-1,sizeof(head));
}
void Add_edge(int u,int v)
{
      edge[++edgecount].v=v;
      edge[edgecount].next=head[u];
      head[u]=edgecount;
}
stack<int > St;
void Tarjan(int u)//从节点x开始搜索
{
     DFN[u]=LOW[u]=++tot;
     vis[u]=1;//为1表示在队列里面
     St.push(u);
     for(int k=head[u];k!=-1;k=edge[k].next)
     {


           int v=edge[k].v;
           if(!DFN[v])//还未访问过
           {
                 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 x;
           ++cnt;
           while(1)
           {
                 x=St.top();
                 St.pop();
                 vis[x]=0;
                 belong[x]=cnt;
                 if(x==u)break;
           }
     }
}
void Solve()
{
    tot=0;
    cnt=0;//缩点后的点数
    memset(DFN,0,sizeof(DFN));
    memset(LOW,0,sizeof(LOW));
    memset(vis,0,sizeof(vis));


    memset(out,0,sizeof(out));
    while(!St.empty()) St.pop();
    for(int i=1;i<=n;i++)
    {
          if(DFN[i]==0)Tarjan(i);
    }
    for(int u=1;u<=n;u++)
    {
          for(int k=head[u];k!=-1;k=edge[k].next)
          {
                int v=edge[k].v;
                if(belong[u]!=belong[v])
                {
                      out[belong[u]]++;


                }
          }
    }
    int ans1=0;
    int ans2=0;
    int ff=0;
    int sum[MAXN];
    int numm=0;
    for(int i=1;i<=cnt;i++)
    {
          if(out[i]==0)
          {
              for(int j=1;j<=n;j++)
              {
                  if(belong[j]==i)
                  {
                     sum[numm++]=j;
                  }
              }
          }
    }
    sort(sum,sum+numm);
    cout<<sum[0];
    for(int i=1;i<numm;i++)
    {
        cout<<" "<<sum[i];
    }


    cout<<endl;


}
int main()
{
    int v;
    while(scanf("%d",&n)!=EOF&&n)
    {
          Init();
          int mm,aa,bb;
          cin>>mm;
          for(int u=1;u<=mm;u++)
          {
              cin>>aa>>bb;
              Add_edge(aa,bb);
          }
          Solve();
    }
    return 0;
}