求点双连通分量

来源:互联网 发布:python vim 编辑:程序博客网 时间:2024/05/16 13:38
 #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);}

0 0
原创粉丝点击