The Bottom of a Graph
Time Limit: 3000MS | Memory Limit: 65536K | |
Total Submissions: 7103 | Accepted: 2907 |
Description
We will use the following (standard) definitions from graph theory. Let
V be a nonempty and finite set, its elements being called vertices (or nodes). Let
E be a subset of the Cartesian product
V×V, its elements being called edges. Then
G=(V,E) is called a directed graph.
Let n be a positive integer, and let p=(e1,...,en) be a sequence of length n of edges ei∈E such that ei=(vi,vi+1) for a sequence of vertices (v1,...,vn+1). Then p is called a path from vertex v1 to vertex vn+1 in G and we say that vn+1 is reachable from v1, writing (v1→vn+1).
Here are some new definitions. A node v in a graph G=(V,E) is called a sink, if for every node w in G that is reachable from v, v is also reachable from w. The bottom of a graph is the subset of all nodes that are sinks, i.e., bottom(G)={v∈V|∀w∈V:(v→w)⇒(w→v)}. You have to calculate the bottom of certain graphs.
Let n be a positive integer, and let p=(e1,...,en) be a sequence of length n of edges ei∈E such that ei=(vi,vi+1) for a sequence of vertices (v1,...,vn+1). Then p is called a path from vertex v1 to vertex vn+1 in G and we say that vn+1 is reachable from v1, writing (v1→vn+1).
Here are some new definitions. A node v in a graph G=(V,E) is called a sink, if for every node w in G that is reachable from v, v is also reachable from w. The bottom of a graph is the subset of all nodes that are sinks, i.e., bottom(G)={v∈V|∀w∈V:(v→w)⇒(w→v)}. You have to calculate the bottom of certain graphs.
Input
The input contains several test cases, each of which corresponds to a directed graph
G. Each test case starts with an integer number
v, denoting the number of vertices of
G=(V,E), where the vertices will be identified by the integer numbers in the set
V={1,...,v}. You may assume that
1<=v<=5000. That is followed by a non-negative integer
e and, thereafter,
e pairs of vertex identifiers
v1,w1,...,ve,we with the meaning that
(vi,wi)∈E. There are no edges other than specified by these pairs. The last test case is followed by a zero.
Output
For each test case output the bottom of the specified graph on a single line. To this end, print the numbers of all nodes that are sinks in sorted order separated by a single space character. If the bottom is empty, print an empty line.
Sample Input
3 3 1 3 2 3 3 1 2 1 1 2 0
Sample Output
1 3 2
求出度为零的强连通分量的点。
#include<iostream> #include<string.h> #include<stdio.h> #include<vector> using namespace std; #define max_n 5005 #define max_e 250002 #define inf 99999999 int stack[max_n],top;//栈 int isInStack[max_n];//是否在栈内 int low[max_n],dfn[max_n],tim;//点的low,dfn值;time从1开始 int node_id; int head[max_n],s_edge;//邻接表头 s_edge从1开始 int gro_id[max_n]; int n,m; int val[max_n]; int out[max_n]; struct Node { int to; int next; } edge[max_e]; void init() { top=0; node_id=0; memset(isInStack,0,sizeof(isInStack)); memset(low,0,sizeof(low)); memset(dfn,0,sizeof(dfn)); tim=0; memset(out,0,sizeof(out)); memset(head,0,sizeof(head)); s_edge=0; memset(edge,0,sizeof(edge)); } void addedge(int u,int v) { s_edge++; edge[s_edge].to=v; edge[s_edge].next=head[u]; head[u]=s_edge; } int min(int a,int b) { if(a<b)return a; else return b; } void tarjan(int u) { //low值为u或u的子树能够追溯到得最早的栈中节点的次序号 stack[top++]=u; isInStack[u]=1; dfn[u]=++tim; //记录点u出现的记录,并放在栈中 low[u]=tim; int e,v; for(e=head[u]; e; e=edge[e].next) //如果是叶子节点,head[u]=0,edge[e].next=0; { v=edge[e].to; if(!dfn[v]) { tarjan(v); low[u]=min(low[u],low[v]); } else if(isInStack[v]) low[u]=min(low[u],dfn[v]); } int j; if(dfn[u]==low[u]) { node_id++; while(j=stack[--top]) { isInStack[j]=0; gro_id[j]=node_id; if(j==u)break; } } } void find() { for(int i = 1 ; i <=n ; ++i) { if(!dfn[i]) { tarjan(i); } } } vector<int> vec[max_n]; int main() { int a,b; while(scanf("%d",&n)!=EOF) { if(n==0)break; cin>>m; init(); for(int i=1; i<=n; i++) { vec[i].clear(); } for(int i = 0 ; i <m ; ++i) { scanf("%d%d",&a,&b); vec[a].push_back(b); addedge(a,b); } find(); int sum=0; for(int i=1; i<=n; i++) { for(int j=0; j<vec[i].size(); j++)//求入度为零的点 { if(gro_id[i]!=gro_id[vec[i][j]]) { if(out[gro_id[i]]==0) sum++; out[gro_id[i]]=1; } } } if(sum!=node_id) for(int i=1; i<=n; i++) { if(out[gro_id[i]]==0) { cout<<i<<" "; } } cout<<endl; } return 0; }