求点双连通分量

 #include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#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);
}