node_sltfl.cpp

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=10007;
int iscut[maxn];  //割点标志
int bccno[maxn];  //i点属于哪个点双连通分块
int pre[maxn], dfs_clock, bcc_cnt;   //bcc_cnt代表点双连通分块数
vector<int> G[maxn], bcc[maxn];  //bcc[i]代表第i个点双连通分块，注意i从1开始编号
//割点可以属于任何一个点双连通分块
struct Edge{
	int u,v;
};
stack<Edge> S;
void init(int n){
	for(int i=0;i<=n;i++)
		G[i].clear();
}
void addEdge(int u,int v){
	 G[u].push_back(v);
}
int dfs(int u,int fa){
	int lowu=pre[u]=++dfs_clock;
	int child=0;
	for(int i=0;i<G[u].size();i++){
		 int v=G[u][i];
		 Edge e=(Edge){u,v};
		 if(!pre[v]){
			  S.push(e);
			  child++;
			  int lowv=dfs(v,u);
			  lowu=min(lowu,lowv);
			  if(lowv>=pre[u]){
					iscut[u]=true;
					bcc_cnt++; bcc[bcc_cnt].clear();
					for(;;){
						 Edge x=S.top(); S.pop();
						 if(bccno[x.u]!=bcc_cnt){ bcc[bcc_cnt].push_back(x.u);bccno[x.u]=bcc_cnt; }
						 if(bccno[x.v]!=bcc_cnt){ bcc[bcc_cnt].push_back(x.v);bccno[x.v]=bcc_cnt; }
						 if(x.u==u&&x.v==v) break;
					}
			  }
		 }
		 else if(pre[v]<pre[u] && v!=fa){
			  S.push(e);
			  lowu=min(lowu,pre[v]);
		 }
	}
	if(fa<0&&child==1) iscut[u]=0;
	return lowu;
}
void find_bcc(int n){    //寻找点双连通分块
	 memset(pre,0,sizeof(pre));
	 memset(iscut,0,sizeof(iscut));
	 memset(bccno,0,sizeof(bccno));
	 dfs_clock=bcc_cnt=0;
	 for(int i=1;i<=n;i++)
		if(!pre[i]) dfs(i,-1);
}