PKU 2553 The Bottom of a Graph - 极大强连通分量

PKU2553 ZJU1979 The Bottom of a Graph

题目描述：

给出一个有向图，定义：若节点v所有能到达的点{wi}，都能反过来到达v，那么称节点v是sink。题目要求所有的sink点。

分析：

很类似这样一个题，足球队队员之间某些人可以联系到另一些人，问教练至少要通知多少个队员，可以让所有队员都收到训练通知。

这个经典题目在吴文虎的图论书里面有讲解，就是求有向图的最小点基。方法如下：

首先将有向图划分成若干个极大强连通分量。那么对于每个强连通分量Ci来说，只需要通知一个人。而这些强连通分量之间要么相互独立，要么存在单向边从连通分量Ci指向Cj，构成森林。对于每棵树来说，只需要通知根结点Croot当中的任意一人，整棵树都可以被通知到。那么，只需要统计入度为0的连通分量的个数即可。

 

现在再回到这个题目，类似地，题目要求的实际上就是所有在出度为0的连通分量当中的节点。按序号输出即可。

此题的关键在求有向图的强连通分量。后面的工作都很简单。

另外，很高兴的是在此题ZJU刷到了第一：）~

Submit Time                  Language    Run Time(ms)       Run Memory(KB)         User Name
2008-11-05 21:53:04      C++           30                            628                               Tiaotiao

1. /*
2. PKU2553 The Bottom of a Graph
3. */
5. #include <stdio.h>
6. #include <memory.h>
8. #define clr(a) memset(a,0,sizeof(a))
9. #define N 5005
10. #define M 50000
12. typedef struct StrNode{
13.     int j;
14.     struct StrNode* next;
15. }Node;
17. Node mem[M];
18. int memp;
20. void addEdge(Node* e[],int i,int j){
21.     Node* p = &mem[memp++];
22.     p->j=j;
23.     p->next = e[i];
24.     e[i]=p;
25. }
27. int g_DFS_First;
28. void DFS_conn(Node* e[],int i,int mark[],int f[],int* nf){
29.     int j; Node* p;
30.     if(mark[i]) return;
31.     
32.     mark[i]=1;
33.     if(!g_DFS_First) f[i]=*nf;  //反向搜索，获取连通分量编号
35.     for(p=e[i];p!=NULL;p=p->next){
36.         j=p->j;
37.         DFS_conn(e,j,mark,f,nf);
38.     }
39.     if(g_DFS_First) f[(*nf)++]=i;   //正向搜索，获取时间戳
40. }
42. /*
43. 有向图极大强连通分量
44. 参数：
45.     邻接表e[]，节点数n。返回极大强连通分支的个数ncon。 
46.     返回con[i]表示节点i所属强连通分量的编号，0 ~ ncon-1。
47. */
48. int Connection(Node* e[],int n,int con[]){
49.     int i,j,k,mark[N],ncon;
50.     int time[N],ntime;  //time[i]表示时间戳为i的节点 
51.     Node *p,*re[N]; //反向边
52.     
53.     //构造反向边邻接表
54.     clr(re);
55.     for(i=0;i<n;i++){
56.         for(p=e[i];p!=NULL;p=p->next)
57.             addEdge(re,p->j,i);
58.     }
59.     
60.     //正向DFS，获得时间戳 
61.     g_DFS_First = 1;
62.     clr(mark); clr(time); ntime=0;
63.     for(i=0;i<n;i++){
64.         if(!mark[i])
65.             DFS_conn(e,i,mark,time,&ntime);
66.     }
67.     
68.     //反向DFS，获得强连通分量 
69.     g_DFS_First = 0;
70.     clr(mark); clr(con); ncon=0;
71.     for(i=n-1;i>=0;i--){
72.         if(!mark[time[i]]){
73.             DFS_conn(re,time[i],mark,con,&ncon);
74.             ncon++;
75.         }
76.     }
77. }
80. //vars
81. Node* e[N];
82. int con[N],ncon;
83. int n,m;
85. int d[N];
87. int main()
88. {
89.     int i,j,k,x,y;
90.     Node *p;
91.     
92.     while(scanf("%d",&n)!=EOF && n){
93.         //init
94.         memp=0; clr(e);
95.         //input
96.         scanf("%d",&m);
97.         for(k=0;k<m;k++){
98.             scanf("%d%d",&i,&j); i--; j--;
99.             addEdge(e,i,j);
100.         }
101.         //Connection
102.         ncon=Connection(e,n,con);
103.         //work d[]
104.         clr(d);
105.         for(i=0;i<n;i++){
106.             x=con[i];
107.             for(p=e[i];p!=NULL;p=p->next){
108.                 y=con[p->j];
109.                 if(x!=y) d[x]++;
110.             }
111.         }
112.         //output
113.         j=0;
114.         for(i=0;i<n;i++){
115.             k=con[i];
116.             if(d[k]==0){
117.                 if(j) printf(" "); j=1;
118.                 printf("%d",i+1);
119.             }
120.         }
121.         puts("");
122.     }
123.     
124.     return 0;
125. }