COGS——T 1175. [顾研NOIP] 旅游电车

http://www.cogs.pro/cogs/problem/problem.php?pid=1175

★★☆ 输入文件：buss.in 输出文件：buss.out 简单对比
时间限制：1 s 内存限制：256 MB

【问题描述】

Henryy国正致力于首都的一个旅游电车建设工程。首都有N个旅游景区。Henryy国的电车永远只沿道路规定的方向行驶，为了不使投入使用的电车有可能无法回到它的起始站，Henryy希望知道他的首都的可以在哪些景区设置站点。一个景区可以被设置成车站，当且仅当对于任意一个从该景区出发所能到达的景区，均至少有一条路可回到该景区。你的同事已完成了一份景区之间的道路连通情况的报告。报告中将给出首都的景区数目N、道路总数M以及一些形如“景区A和景区B之间有一条从A到B的单向道路”的信息。现在明确你的任务：根据报告中的信息，

按编号从小到大

列出所有可以被设置成车站的景区。

【输入文件】

输入文件由多份报告组成（这些报告相互无任何联系），每份报告包括：N，M，接下来M对整数Ai、Bi (1<=I<=M)表示Ai和Bi之间有一条单向道路Ai->Bi。仅一个包含整数N=0的报告表示你的工作结束，你的程序不应该对此有任何反应。各整数间用空格或空行分隔。对于任意景区，分别以该景区为起点或终点的道路总数均不超过50。

【输出文件】

对于每份报告，输出一行列表包括：所有能被设置成电车站点的景区编号，各编号之间用一个空格隔开。

【样例输入】

【样例输出】

1 2 3 4 5

1

【数据约定】

对于40%的数据，有N<=200。

对于100%的数据，有N<=5000，M<=50000。

缩点后的新图中，出去为零的强连通里包含的点、

 #include <cstring>

 #include <cstdio>

 const int N(+);

 const int M(+);

 int head[N],sumedge;

 struct Edge

 {

     int v,next;

     Edge(int v=,int next=):v(v),next(next){}

 }edge[M];

 inline void ins(int u,int v)

 {

     edge[++sumedge]=Edge(v,head[u]);

     head[u]=sumedge;

 }

 #define min(a,b) (a<b?a:b)

 int tim,dfn[N],low[N];

 int top,Stack[N],instack[N];

 int col[N],sumcol,cd[N];

 void DFS(int u)

 {

     low[u]=dfn[u]=++tim;

     Stack[++top]=u; instack[u]=;

     for(int v,i=head[u];i;i=edge[i].next)

     {

         v=edge[i].v;

         if(!dfn[v]) DFS(v), low[u]=min(low[u],low[v]);

         else if(instack[v]) low[u]=min(low[u],dfn[v]);

     }

     if(low[u]==dfn[u])

     {

         col[u]=++sumcol;

         for(;Stack[top]!=u;top--)

         {

             instack[Stack[top]]=;

             col[Stack[top]]=sumcol;

         }

         instack[u]=; top--;

     }

 }

 inline void init()

 {

     sumedge=tim=top=sumcol=;

     memset(cd,,sizeof(cd));

     memset(low,,sizeof(low));

     memset(dfn,,sizeof(dfn));

 //    memset(col,0,sizeof(col));

     memset(edge,,sizeof(edge));

     memset(head,,sizeof(head));

 //    memset(Stack,0,sizeof(Stack));

     memset(instack,,sizeof(instack));

 }

 inline void read(int &x)

 {

     x=; register char ch=getchar();

     for(;ch>''||ch<'';) ch=getchar();

     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-'';

 }

 int AC()

 {

     freopen("buss.in","r",stdin);

     freopen("buss.out","w",stdout);

     for(int n,m;;init())

     {

         read(n);

         if(!n) break;

         read(m);

         for(int u,v;m--;ins(u,v))

             read(u),read(v);

         for(int i=;i<=n;i++)

             if(!dfn[i]) DFS(i);

         for(int v,u=;u<=n;u++)

             for(int i=head[u];i;i=edge[i].next)

             {

                 v=edge[i].v;

                 if(col[u]!=col[v]) cd[col[u]]++;

             }

         for(int i=;i<=n;i++)

             if(!cd[col[i]]) printf("%d ",i);

         puts("");

     }

     return ;

 }

 int Hope=AC();

 int main(){;}