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。

【输出文件】

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

【样例输入】

5 6
1 2
2 3
3 4
4 1
2 5
5 2
1 0
0

【样例输出】

1 2 3 4 5
1

【数据约定】

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

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

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

 1 #include <cstring>
 2 #include <cstdio>
 3
 4 const int N(5000+15);
 5 const int M(50000+5);
 6 int head[N],sumedge;
 7 struct Edge
 8 {
 9     int v,next;
10     Edge(int v=0,int next=0):v(v),next(next){}
11 }edge[M];
12 inline void ins(int u,int v)
13 {
14     edge[++sumedge]=Edge(v,head[u]);
15     head[u]=sumedge;
16 }
17
18 #define min(a,b) (a<b?a:b)
19 int tim,dfn[N],low[N];
20 int top,Stack[N],instack[N];
21 int col[N],sumcol,cd[N];
22 void DFS(int u)
23 {
24     low[u]=dfn[u]=++tim;
25     Stack[++top]=u; instack[u]=1;
26     for(int v,i=head[u];i;i=edge[i].next)
27     {
28         v=edge[i].v;
29         if(!dfn[v]) DFS(v), low[u]=min(low[u],low[v]);
30         else if(instack[v]) low[u]=min(low[u],dfn[v]);
31     }
32     if(low[u]==dfn[u])
33     {
34         col[u]=++sumcol;
35         for(;Stack[top]!=u;top--)
36         {
37             instack[Stack[top]]=0;
38             col[Stack[top]]=sumcol;
39         }
40         instack[u]=0; top--;
41     }
42 }
43
44 inline void init()
45 {
46     sumedge=tim=top=sumcol=0;
47     memset(cd,0,sizeof(cd));
48     memset(low,0,sizeof(low));
49     memset(dfn,0,sizeof(dfn));
50 //    memset(col,0,sizeof(col));
51     memset(edge,0,sizeof(edge));
52     memset(head,0,sizeof(head));
53 //    memset(Stack,0,sizeof(Stack));
54     memset(instack,0,sizeof(instack));
55 }
56 inline void read(int &x)
57 {
58     x=0; register char ch=getchar();
59     for(;ch>‘9‘||ch<‘0‘;) ch=getchar();
60     for(;ch>=‘0‘&&ch<=‘9‘;ch=getchar()) x=x*10+ch-‘0‘;
61 }
62
63 int AC()
64 {
65     freopen("buss.in","r",stdin);
66     freopen("buss.out","w",stdout);
67     for(int n,m;1;init())
68     {
69         read(n);
70         if(!n) break;
71         read(m);
72         for(int u,v;m--;ins(u,v))
73             read(u),read(v);
74         for(int i=1;i<=n;i++)
75             if(!dfn[i]) DFS(i);
76         for(int v,u=1;u<=n;u++)
77             for(int i=head[u];i;i=edge[i].next)
78             {
79                 v=edge[i].v;
80                 if(col[u]!=col[v]) cd[col[u]]++;
81             }
82         for(int i=1;i<=n;i++)
83             if(!cd[col[i]]) printf("%d ",i);
84         puts("");
85     }
86     return 0;
87 }
88
89 int Hope=AC();
90 int main(){;}