D19【模板】vDCC 缩点 Tarjan 算法

D19【模板】vDCC 缩点 Tarjan 算法_哔哩哔哩_bilibili

 

参考：D16【模板】割点 Tarjan 算法 - 董晓 - 博客园

在一张连通的无向图中，对于两个点 x 和 y，如果无论删去哪一个点（不删 x和y）都不能使它们不连通，就说 x 和 y 点双连通

对于一个无向图中的 极大 点双连通的子图，称这个子图为一个 点双连通分量（vDCC）

图中有两个割点 1,5，对割点要裂点，分割成 4 个的 vDCC: {1,2,3,4}, {5,6}, {5,7,8}, {1,5}

用一个栈存点，若遍历回到 x 时，割点判定法则 low[y]≥dfn[x] 成立，则从栈中弹出节点，直到 y 被弹出

那么弹出的节点和 x 一起构成一个 vDCC

三类缩点归纳：

 

P8435 【模板】点双连通分量 - 洛谷

对于一个 n 个节点 m 条无向边的图，请输出其点双连通分量的个数，并且输出每个点双连通分量

#include<bits/stdc++.h>
using namespace std;

const int N=500010;
int n,m;
int dfn[N],low[N],stk[N],top,cnt,root,cut[N];
vector<int> e[N],dcc[N];

void tarjan(int x){
  dfn[x]=low[x]=++dfn[0]; stk[++top]=x;
  
  if(x==root&&!e[x].size()){ //孤立点
    dcc[++cnt].push_back(x);
    return;
  }
  int son=0;
  for(int y:e[x]){
    if(!dfn[y]){ //若y未访问
      tarjan(y);
      low[x]=min(low[x],low[y]);
      
      if(low[y]>=dfn[x]){ //x是割顶
        son++;
        if(x!=root||son>1)cut[x]=1; //记录割点
        
        // 每次回到割顶时记录vDCC
        ++cnt;
        for(int z=-1;z!=y;){ //让x留在栈中
          z=stk[top--];
          dcc[cnt].push_back(z);
        }
        dcc[cnt].push_back(x); //记录vDCC
      }
    }
    else //若y已访问
      low[x]=min(low[x],dfn[y]);
  }
}
int main(){
  ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
  cin>>n>>m;
  for(int a,b;m--;){
    cin>>a>>b;
    if(a==b) continue; //忽略自环
    e[a].push_back(b),
    e[b].push_back(a);
  }
  
  for(root=1;root<=n;root++)if(!dfn[root])tarjan(root);
  
  cout<<cnt<<"\n";
  for(int i=1;i<=cnt;i++){
    cout<<dcc[i].size()<<" ";
    for(int j:dcc[i])cout<<j<<" ";
    cout<<"\n";
  }
}

 

vDCC 缩点

vDCC：不删除边，分裂割点，构成的连通块

把 vDCC 缩成点，把缩点和割点对应连边，构成树

1. cut[x] 记录割点，dcc[] 记录 vDCC，id[] 割点的编号

2. 每次回溯：判割点，记录 vDCC

3. 特判孤立点，忽略自环

#include<bits/stdc++.h>
using namespace std;

const int N=10010;
int n,m,a,b;
vector<int> e[N],ne[N],dcc[N];
int dfn[N],low[N],tot,stk[N],top,cut[N],root,cnt;
int id[N];

void tarjan(int x){
  dfn[x]=low[x]=++tot;
  stk[++top]=x;
  if(x==root&&!e[x].size()){ //孤立点
    dcc[++cnt].push_back(x);
    return;
  }
  int son=0;
  for(int y:e[x]){
    if(!dfn[y]){ //若y未访问
      tarjan(y);
      low[x]=min(low[x],low[y]); 
      if(low[y]>=dfn[x]){
        son++;
        if(x!=root||son>1)cut[x]=1; //割点
        
        ++cnt;
        while(1){
          int z=stk[top--];
          dcc[cnt].push_back(z);
          if(z==y) break; //让x留在栈中
        }
        dcc[cnt].push_back(x); //vDCC
      }
    }
    else //若y已访问
      low[x]=min(low[x],dfn[y]);
  }
}
int main(){
  cin>>n>>m;
  while(m--){
    cin>>a>>b;
    e[a].push_back(b),
    e[b].push_back(a);
  }
  for(root=1;root<=n;root++)if(!dfn[root])tarjan(root);
      
  //给每个割点一个新编号（cnt+1开始）
  int num=cnt;
  for(int i=1;i<=n;i++)if(cut[i])id[i]=++num;
  //建树，从每个缩点向对应割点连边
  for(int i=1;i<=cnt;i++){
    for(int j=0;j<dcc[i].size();j++){
      int x=dcc[i][j];
      if(cut[x]){
        ne[i].push_back(id[x]),
        ne[id[x]].push_back(i);        
      }
    }
  }
}