8.3重要城市——图的割点

#include<iostream>
using namespace std;
int n,m,e[9][9],root;
int num[9],low[9],flag[9],index;//index用来进行时间戳的递增
//求两个数中较小一个数的函数
int min(int a,int b){
	return a<b?a:b;
} 
//割点算法的核心
void dfs(int cur,int father){//需要传入两个参数,当前顶点编号和父顶点的编号
	int child=0,i,j;//child用来记录在生成树中当前顶点cur的儿子个数
	
	index++;//时间戳加1//时间戳:指图中的某结点是第几个被访问到的 
	num[cur]=index;//当前顶点cur的时间戳
	low[cur]=index;//当前顶点cur能够访问到最早顶点的时间戳,刚开始当然是自己啦
	for(i=1;i<=n;i++){//枚举与当前顶点cur有边相连的顶点i
		if(e[cur][i]==1){
			if(num[i]==0){//如果顶点i的时间戳不为0,说明顶点i还没有被访问过
				child++;
				dfs(i,cur);//继续往下深度优先遍历
				//更新当前顶点cur能否访问到最早顶点的时间戳
				low[cur]=min(low[cur],low[i]);
				//如果当前顶点不是根结点并且满足low[i]>=num[cur],则当前顶点为割点
				if(cur!=root&&low[i]>=num[cur]){
					flag[cur]=1;
				} 
				//如果当前顶点是根结点,在生成树中根结点必须要有两个儿子,那么这个根结点才是割点
				if(cur==root&&child==2){
					flag[cur]=1;
				}
			}
			else if(i!=father){//否则如果顶点i曾经被访问过,并且这个顶点不是当前顶点cur的父亲,则
			//需要更新当前结点cur能否访问到最早顶点的时间戳
				low[cur]=min(low[cur],num[i]); 
			}
		} 
	} 
} 

int main(){
	int i,j,x,y;
		
	//读入n和m,n表示顶点个数,m表示边的条数
	cout<<"请分别输入原图顶点个数n,边的条数m："<<endl;
	cin>>n>>m;

	//初始化
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++){
			e[i][j]=0;
		}
	} 
 	
	//读入边 
	cout<<"请依次输入m条边，每条边起点为b1，终点为b2，权重为z3。即b1->b2(z3)："<<endl; 
	for(i=1;i<=m;i++){
		//读入每一条边 
		cin>>x>>y;
		//这里是无向图,所以需要将边反向再存储一遍 
		e[x][y]=1;
		e[y][x]=1;
	}
	root=1;
	dfs(1,root);//从1号顶点开始进行深度优先遍历
	
	for(i=1;i<=n;i++){
		if(flag[i]==1){
			cout<<i<<" ";
		}
	} 
	
	getchar();
	return 0;
}

/*
总结:
1、在一个无向连通图中,如果删除某个顶点后,图不再连通,称这样的顶点为割点(或称为割顶) 
2、思路:如果k点事割点,那么剩下的没有被访问过的点中至少会有一个点在不经过k点的情况下,是无论如何再也回不到已访问过的点 
3、核心:对于某个顶点u,如果存在至少一个顶点v(即顶点u的儿子),使得low[v]>=num[u],即不能回到祖先,那么u点为割点 
4、num[]数组用来记录每一个顶点的时间戳;low[]数组用来记录每个顶点在不经过父顶点时,能够回到的最小"时间戳"
5、上面的代码使用的邻接矩阵来存储图,这显然是不对的,因为这样无论如何时间复杂度都会在O(N^2),因为边的处理就需要N^2的时间;
这里这样写是为了突出割点算法部分,实际应用中需要改为使用邻接表来存储,这样整个算法的时间复杂度是O(N+M)
*/