无向图的连通分量(BCC)模版

点双连通:如果任意两点之间至少存在两条”点不重复”的路径,则说这个图是点双连通的,这个要求等价于任意两条边都在同一个简单环中,即内部无割点

边双连通:如果任意两点至少存在两条“边不重复“的路径,我们说这个图是边-双连通的,这个要求低一点,只需要每条边都至少在一个简单环中,即所有的边都不是桥

对于一张无向图,每条边恰好属于一个双连通分量,但不同的双连通分量之间可能会有公共点,但是不同的双连通分量之间的公共点最多只有一个,且这个点一定是割顶,另一方面,任意割顶都是至少两个不同双连通分量的公共点

点-双连通

#include <cstdio>
#include <cstring>
#include <vector>
#include <stack>
#include <algorithm>
using namespace std;
const int MAXNODE = 10005;
const int MAXEDGE = 100010;
const int INF = 0x3f3f3f3f;

struct Edge{
    int v, id, next;
    bool bridge;
    Edge() {}
    Edge(int v, int id, int next): v(v), id(id), next(next), bridge(false){}
}E[MAXEDGE * 2], cut[MAXEDGE];

struct Node {
    int u, v;
    Node() {}
    Node(int u, int v): u(u), v(v) {}
};

stack<Node> Stack;
//bcc表示的是一个bcc里面的点
vector<int> bcc[MAXNODE];
//pre纪录的是时间戳,lowlink纪录的是该点及其该子孙节点所能返回的最早时间戳是多少,bccno纪录的是该点当前是属于哪个bcc的
int head[MAXNODE], pre[MAXNODE], lowlink[MAXNODE], bccno[MAXNODE];
int n, m, tot, bcc_cnt, dfs_clock, cut_cnt; 
bool iscut[MAXNODE];

//双向边的添加
void AddEdge(int u, int v, int id) {
    E[tot] = Edge(v, id, head[u]);
    head[u] = tot++;
    E[tot] = Edge(u, id, head[v]);
    head[v] = tot++;
}

void init() {
    memset(head, -1, sizeof(head));
    tot = 0;
}

//点双连通
void dfs(int u, int fa) {
    pre[u] = lowlink[u] = ++dfs_clock;
    int child = 0;//纪录当前节点有多少个子节点
    for (int i = head[u]; ~i; i = E[i].next) {
        int v = E[i].v;
        if (!pre[v]) {
            Stack.push(Node(u, v));
            child++;
            dfs(v, u);
            lowlink[u] = min(lowlink[u], lowlink[v]);//更新
            //子节点最多返回到该点
            if (lowlink[v] >= pre[u]) {
                //该边为桥
                if (lowlink[v] > pre[u]) {
                    E[i].bridge = E[i ^ 1].bridge = true;
                    cut[cut_cnt++] = E[i];
                }
                iscut[u] = true;
                bcc_cnt++; bcc[bcc_cnt].clear();
                while (1) {
                    Node x = Stack.top(); Stack.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 (v != fa && pre[v] < pre[u]) {//反向边
            Stack.push(Node(u, v));
            lowlink[u] = min(lowlink[u], pre[v]);
        }
    }
    //u是根结点,且只有一个孩子,那就不是割点了
    if (fa < 0 && child == 1) iscut[u] = 0;
}

void find_bcc() {
    memset(pre, 0, sizeof(pre));
    memset(iscut, 0, sizeof(iscut));
    memset(bccno, 0, sizeof(bccno));
    dfs_clock = bcc_cnt = cut_cnt = 0;
    for (int i = 0; i < n; i++) 
        if (!pre[i]) dfs(i, -1);
}

int main() {
    init();
    return 0;
}

边双连通

学长的:

#include <cstdio>  
#include <cstring>  
#include <algorithm>  
using namespace std;  

const int N = 10005;  
const int M = 20005;  

struct Edge{  
    int u, v, id, next;  
    bool iscut;  
    Edge() {}  
    Edge(int u, int v, int id, int next): u(u), v(v), id(id), next(next), iscut(false) {}  
}E[M * 2], cut[M];  

int n, m, tot, bcc_cnt, dfs_clock, cut_cnt;
int head[N], val[N], pre[N], lowlink[N], bccno[N]; 

void init() {  
    memset(head, -1, sizeof(head));
    tot = 0;
}  

void AddEdge(int u, int v, int id) {  
    E[tot] = Edge(u, v, id, head[u]);
    head[u] = tot++;
    E[tot] = Edge(v, u, id, head[v]);
    head[v] = tot++;
}  

void dfs_cut(int u, int fa) {  
    pre[u] = lowlink[u] = ++dfs_clock;  
    for (int i = head[u]; ~i; i = E[i].next) {
        int v = E[i].v;
        if (E[i].id == fa) continue;
        if (!pre[v]) {  
            dfs_cut(v, E[i].id);  
            lowlink[u] = min(lowlink[u], lowlink[v]);  
            if (lowlink[v] > pre[u]) {  
                E[i].iscut = E[i^1].iscut = true;//标记割边  
                cut[cut_cnt++] = E[i];//存放割边  
            }  
        } else lowlink[u] = min(lowlink[u], pre[v]);  
    }  
}  

void find_cut() {  
    dfs_clock = cut_cnt = 0;  
    memset(pre, 0, sizeof(pre));  
    for (int i = 0; i < n; i++)  
        if (!pre[i]) dfs_cut(i, -1);  
}  

void dfs_bcc(int u) {  
    pre[u] = 1;  
    bccno[u] = bcc_cnt;  
    for (int i = head[u]; ~i; i = E[i].next) {  
        if (E[i].iscut) continue;  
        int v = E[i].v;  
        if (pre[v]) continue;  
        dfs_bcc(v);  
    }  
}  

void find_bcc() {  
    bcc_cnt = 0;  
    memset(pre, 0, sizeof(pre));  
    for (int i = 0; i < n; i++) {  
        if (!pre[i]) {  
            dfs_bcc(i);  
            bcc_cnt++;  
        }  
    }  
}  

int main() {  
    return 0;  
}  

另一份

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXNODE = 10005;
const int MAXEDGE = 100010;
const int INF = 0x3f3f3f3f;

struct Edge{
    int v, id, next;
    bool bridge;
    Edge() {}
    Edge(int v, int id, int next): v(v), id(id), next(next), bridge(false){}
}E[MAXEDGE * 2], cut[MAXEDGE];

//pre纪录的是时间戳,lowlink纪录的是该点及其该子孙节点所能返回的最早时间戳是多少,bccno纪录的是该点当前是属于哪个bcc的
int head[MAXNODE], pre[MAXNODE], lowlink[MAXNODE], bccno[MAXNODE], Stack[MAXNODE], belong[MAXNODE];
int n, m, tot, bcc_cnt, dfs_clock, top; 

void AddEdge(int u, int v, int id) {
    E[tot] = Edge(v, id, head[u]);
    head[u] = tot++;
    E[tot] = Edge(u, id, head[v]);
    head[v] = tot++;
}

void init() {
    memset(head, -1, sizeof(head));
    tot = 0;
}

//边双连通
void dfs(int u, int fa) {
    pre[u] = lowlink[u] = ++dfs_clock;
    Stack[++top] = u;

    for (int i = head[u]; ~i; i = E[i].next) {
        int v = E[i].v;
        if (!pre[v]) {
            dfs(v, u);
            lowlink[u] = min(lowlink[u], lowlink[v]);
            if (lowlink[v] > pre[u]) { 
                E[i].bridge = E[i ^ 1].bridge = true;
                bcc_cnt++;
                while (1) {
                    int x = Stack[top--];
                    belong[x] = bcc_cnt;
                    if (x == v) break;
                }
            }
        }
        else if (pre[v] < pre[u] && v != fa) lowlink[u] = min(lowlink[u], pre[v]);
    }
}

void find_bcc() {
    memset(pre, 0, sizeof(pre));
    memset(bccno, 0, sizeof(bccno));
    dfs_clock = bcc_cnt = 0;
    for (int i = 0; i < n; i++) 
        if (!pre[i]) dfs(i, -1);
}

int main() {
    init();
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值