无向图 求EBC+ 求桥 + 缩点 【模版】

BCC 为 点-双联通 ; EBC 为边-双联通

(和有向图中强联通 的类似)

Eg 条件 给定n个点,m个无向图 ,求EBC +求桥+缩点

代码

#include <cstdio>
#include <cstring>
#include <vector>
#include <stack>
#include <queue>
#include <algorithm>
#define MAXN 200000+10 
#define MAXM 2000000+10 
#define INF 100000000
using namespace std;
struct Edge
{
    int from, to, next;
    bool cut;//是否为桥 
}edge[MAXM];
int head[MAXN], edgenum; 
int dfs_clock;//dfn次数 时间戳的 应用 
int low[MAXN], dfn[MAXN];
int ebcno[MAXN], ebc_cnt;//ebcno[i]表示i属于哪个EBC    ebc_cnt是EBC计数器 
int bridge;//桥数 
int n, m;//点数 边数
stack<int> S;//存储当前ebc的所有点
bool Instack[MAXN];//标记该点是否在栈里面 
vector<int> G[MAXN];//存储缩点后新图 
vector<int> ebc[MAXN];//存储EBC里面的点 
void init()
{
    edgenum = 0;
    memset(head, -1, sizeof(head));
} 
void addEdge(int u, int v)
{
    Edge E = {u, v, head[u], 0};
    edge[edgenum] = E;
    head[u] = edgenum++; 
}
void getMap()//建图 
{
    int a, b;
    while(m--)
    {
        scanf("%d%d", &a, &b);
        addEdge(a, b);
        addEdge(b, a);
    } 
}
void tarjan(int u, int fa)
{
    int v;
    low[u] = dfn[u] = ++dfs_clock;//时间戳 
    S.push(u);
    Instack[u] = true;//进栈
    int have = 1;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        v = edge[i].to;
        if(have && v == fa)//重边有效 无效的话直接判断 v == fa continue; 
        {
            have = 0;
            continue;
        }
        if(!dfn[v])
        {
            tarjan(v, u);
            low[u] = min(low[u], low[v]);
            if(low[v] > dfn[u])//桥 
            {
                bridge++;
                edge[i].cut = edge[i^1].cut = true;
            }
        }
        else if(Instack[v])//已搜索过 且 在栈里面 
        low[u] = min(low[u], dfn[v]);
    } 
    if(low[u] == dfn[u]) 
    {
        ebc_cnt++;
        ebc[ebc_cnt].clear();//清空 开始记录点 
        for(;;)//弹出属于同一个EBC的点 
        {
            v = S.top(); S.pop();
            ebcno[v] = ebc_cnt;//属于当前EBC
            ebc[ebc_cnt].push_back(v); 
            Instack[v] = false;
            if(v == u) break; 
        } 
    }
}
void suodian()//缩点  建新图 
{
    for(int i = 1; i <= ebc_cnt; i++) G[i].clear();//初始化 
    for(int i = 0; i < edgenum; i+=2)//跳过重边 
    {
        int u = ebcno[edge[i].from];
        int v = ebcno[edge[i].to];
        if(u != v)//不在同一个EBC
        G[u].push_back(v), G[v].push_back(u); 
    }
}
void find_cut(int l, int r)
{
    memset(low, 0, sizeof(low));
    memset(dfn, 0, sizeof(dfn));
    memset(Instack, false, sizeof(Instack));
    memset(ebcno, 0, sizeof(ebcno));
    dfs_clock = ebc_cnt = bridge = 0;
    for(int i = l; i <= r; i++)
    if(!dfn[i]) tarjan(i, -1);
} 
int main()
{
    while(scanf("%d%d", &n, &m), n||m)
    {
        init();
        getMap();
        find_cut(1, n);//找桥 以及所有的 EBC 
        suodian();//缩点 
        /*到此已经求出桥数 且 进行过缩点*/ 
    } 
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值