计算点-双连通分量算法

//割顶的bccno无意义#include<cstdio>#include<cstring>#include<stack>#include<vector>using namespace std;#define maxn 1100struct Edge{int u,v;};int pre[maxn],iscut[maxn],bccno[maxn],dfs_clock,bcc_cnt;stack<Edge> S;vector<int> G[maxn],bcc[maxn];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;e.u=u;e.v=v;if(!pre[v])//没有访问过v{S.push(e);child++;int lowv=dfs(v,u);lowu=min(lowu,lowv);//用后代的low函数更新自己if(lowv>=pre[u]){iscut[u]=true;bcc_cnt++;bcc[bcc_cnt].clear();//bcc从1开始编号while(1){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,(n+10)*sizeof(pre[0]));memset(iscut,0,(n+10)*sizeof(iscut[0]));memset(bccno,0,(n+10)*sizeof(bccno[0]));dfs_clock=bcc_cnt=0;for(int i=0;i<n;i++){if(!pre[i])dfs(i,-1);}}