ZJU1311 Network - 无向图的割点

题目大意：

给出一个无向图，求图中割点的个数。（若去掉顶点i及其邻边，导致剩下的图不连通，则i是割点）

分析：

求割点有很成熟的DFS算法，照书写了一个，整理成模板备用。

 

1. /*
2. ZJU1311 Network
3. */
5. #include <stdio.h>
6. #include <memory.h>
8. #define clr(a) memset(a,0,sizeof(a))
9. #define MIN(a,b) ((a)>(b)?(b):(a))
10. #define N 105
12. int n;
13. int a[N][N];
14. int cut[N];
16. char* Next(char str[],char *sp){
17.     while(*sp&&*sp!=' ') sp++;
18.     while(*sp&&*sp==' ') sp++;
19.     return sp;
20. }
22. /**************************************/
24. int ancestor[N];
25. int mark[N];
26. int deep[N];
28. int DFS(int i,int father,int dth,int b[][N],int cut[]){
29.     int j,k,sons=0,count=0;
30.     mark[i]=1;
31.     deep[i]=ancestor[i]=dth;
32.     for(k=b[i][0];k>=1;k--){
33.         j=b[i][k];
34.         if(j!=father && mark[j]==1) ancestor[i]=MIN(ancestor[i],deep[j]);
35.         if(mark[j]==0){
36.             count+=DFS(j,i,dth+1,b,cut);
37.             sons++;
38.             ancestor[i]=MIN(ancestor[i],ancestor[j]);
39.             if((father==-1&&sons>1)||(father!=-1&&ancestor[j]>=deep[i]))
40.                 cut[i]=1;
41.             //if(ancestor[j]>deep[i]) brige[i][j]=1;
42.         }
43.     }
44.     mark[i]=2;
45.     return count+cut[i];
46. }
47. /*
48. 求割点
49.     参数：a[][]邻接矩阵，非0表示连通。n顶点个数。
50.     返回：割点数量，cut[i]==1表示i是割点。
51.     说明：DFS算法，先将邻接矩阵转换为邻接表，可以处理图本身不连通的情况。复杂度O(E)
52.     可稍作修改求图的桥。 
53. */
54. int cutpoints(int a[][N],int n,int cut[]){
55.     int i,j,b[N][N],count=0;
56.     for(i=0;i<n;i++) cut[i]=b[i][0]=0;
57.     memset(mark,0,sizeof(mark));
58.     for(i=0;i<n;i++)for(j=0;j<n;j++) if(a[i][j]) b[i][++b[i][0]]=j;
59.     for(i=0;i<n;i++)if(mark[i]==0) count+=DFS(0,-1,0,b,cut);
60.     return count;
61. }
63. /**************************************/
65. int main()
66. {
67.     int i,j,k;
68.     char str[N*10],*sp;
69.     
70.     while(scanf("%d",&n)!=EOF && n){
71.         //init
72.         clr(a);
73.         //input
74.         while(scanf("%d",&i),i){
75.             gets(str);
76.             for(sp=Next(str,str);sscanf(sp,"%d",&j)!=EOF;sp=Next(str,sp)){
77.                 a[i-1][j-1]=a[j-1][i-1]=1;
78.             }
79.         }
80.         //work
81.         printf("%d/n",cutpoints(a,n,cut));
82.     }
83.     
84.     return 0;
85. }