【无向图 联通分量】

版权声明:本博客主要用于记录思考过程,方便复习,欢迎留言交流讨论~ https://blog.csdn.net/Floraqiu/article/details/82023985

强连通分量
1

2

无向图 双联通分量
1

2 例题

这里写图片描述

关于tarjan算法,一直有一个很大的争议,就是low[u]=min(low[u],dfn[v]);

这句话,如果改成low[u]=min(low[u],low[v])就会wa掉,但是在求强连通分量时却没有问题

根据许多大佬的观点,我想提出自己的一点看法,在求强连通分量时,如果v已经在栈中,那么说明u,v一定在同一个强连通分量中,所以到最后low[u]=low[v]是必然的,提前更新也不会有问题,但是在求割点时,low的定义有了小小的变化,不再是最早能追溯到的祖先,(因为是个无向图)没有意义,应该是最早能绕到的割点,为什么用绕到,是因为是无向边,所以有另一条路可以走,如果把dfn[v]改掉就会上翻过头,可能翻进另一个环中,所以wa掉,仅是本人的一些个人看法,不知道讲的对不对,请各位指教

做了的例题

POJ 1523
luogu P3388 【模板】割点(割顶)
POJ 3694(就是WA???先放着)

无向图 割顶 桥

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <list>
#define clr(x, y)  memset(x, y, sizeof(x))
#define MOD 1000000000+7
#define INF 0x3f3f3f3f

using namespace std;

const int maxn = 1000;
int n, m;
vector<int> G[maxn];
int pre[maxn], low[maxn];//pre是dfs时间序列, low是该节点及其子孙能够连回最早的祖先的dfs值
int dfs_clock; //时间戳
int iscut[maxn];//标记是否为割顶

int dfs(int u, int fa)//u在dfs树中的父节点是fa
{
    int lowu = pre[u] = ++dfs_clock;
    int child = 0;//子节点数目
    for(int i = 0; i < G[u].size(); ++i)
    {
        int v = G[u][i];
        if(!pre[v]) // 没有访问过v
        {
            ++child;
            int lowv = dfs(v, u);
            lowu = min(lowu, lowv);
            if(lowv >= pre[u]) iscut[u] = 1; // u是割顶
            if(lowv > pre[u]) cout << "桥 : " << u << " -> " << v << endl; // u -> v是桥
        }
        else if(pre[v] < pre[u] && v != fa)//用反向边更新u的low函数
            lowu = min(lowu, pre[v]);
    }
    if(fa < 0 && child == 1) iscut[u] = 0;//对根节点的处理
    low[u] = lowu;
    return lowu;
}

int main()
{
    while(~scanf("%d%d", &n, &m))
    {
        clr(pre, 0);
        clr(cut, 0);
        for(int i = 1; i <= n; ++i)
            G[i].clear();
        int u, v;
        for(int i = 1; i <= m; ++i)
        {
            scanf("%d%d", &u, &v);
            G[u].push_back(v);
            G[v].push_back(u);
        }
        dfs(1, -1);//对于根节点的fa,要初始化为-1
        for(int i = 1; i <= n; ++i)//打印割顶
            if(iscut[i]) cout<<i<<endl;
    }
    return 0;
}

点双联通分量

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <vector>
#include <queue>
#include <set>
#include <stack>
#include <map>
#include <list>
#define clr(x, y)  memset(x, y, sizeof(x))
#define MOD 1000000000+7
#define INF 0x3f3f3f3f

using namespace std;

const int maxn = 1000;
int pre[maxn], bccno[maxn], dfs_clock, bcc_cnt;
bool iscut[maxn];
vector<int> G[maxn], bcc[maxn];
struct Edge
{
    int u, v;
    Edge(int uu, int vv): u(uu), v(vv) {};
};

stack<Edge> S;

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];
        if(!pre[v])
        {
            ++child;
            S.push(Edge(u, v));
            int lowv = dfs(v, u);
            lowu = min(lowu, lowv);
            if(lowv >= pre[u])//割顶,有BCC
            {
                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;
                    }
                }
            }
        }
        else if(pre[v] < pre[u] && v != fa)
        {
            S.push(Edge(u, v));
            lowu = min(lowu, pre[v]);
        }
    }
    if(fa < 0 && child == 1)
        iscut[u] = false;
    return lowu;
}

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

void print_bcc()
{
    for(int i = 1; i <= bcc_cnt; ++i)
    {
        for(int j = 0; j < bcc[i].size(); ++j)
            printf("%d ", bcc[i][j]);
        printf("\n");
    }
}

边双联通分量

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <vector>
#include <queue>
#include <set>
#include <stack>
#include <map>
#include <list>
#define clr(x, y)  memset(x, y, sizeof(x))
#define MOD 1000000000+7
#define INF 0x3f3f3f3f

using namespace std;

const int maxn = 1000;
int pre[maxn], bccno[maxn], head[maxn];
int dfs_clock, ebc_cnt, tot;
bool isbridge[maxn];//标记编号为i的边是否为桥
vector<int> ebc[maxn];

struct Edge
{
    int no, v, next;
}edge[maxn];

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

void addEdge(int no, int u, int v)
{
    edge[tot].no = no;//记录边的编号,和tot不一样
    edge[tot].v = v;
    edge[tot].next = head[u];
    head[u] = tot++;
}

void dfs_bridge(int u, int fa)//找出所有桥
{
    int lowu = pre[u] = ++dfs_clock;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(!pre[v])
        {
            int lowv = dfs_bridge(v, u);
            lowu = min(lowu, lowv);
            if(lowv > pre[u])
                isbridge[edge[i].no] = true;
        }
        else if(pre[v] < pre[u] && v != fa)
            lowu = min(lowu, pre[v]);
    }
    return lowu;
}

void dfs_countbridge(int u, int fa)
{
    ebc[ebc_cnt].push_back(u);
    pre[u] = ++dfs_clock;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(!isbridge[edge[i].no] && !pre[v])//除去桥和已经访问过的点
            dfs_countbridge(v, u);
    }
}

void find_ebc()
{
    clr(pre, 0);
    clr(isbridge, false);
    dfs_clock = ebc_cnt = 0;
    for(int i = 1; i <= n; ++i)
        if(!pre[i]) dfs_bridge(i, -1);
    clr(pre, 0);
    dfs_clock = 0;
    for(int i = 1; i <= n; ++i)
    {
        if(!pre[i])
        {
            ++ebc_cnt;
            ebc[ebc_cnt].clear();
            dfs_countbridge(i, -1);
        }
    }
}

void print_ebc()
{
    for(int i = 1; i <= ebc_cnt; ++i)
    {
        for(int j = 0; j < ebc[i].size(); ++j)
            printf("%d ", ebc[i][j]);
        printf("\n");
    }
}

int main()
{
    for(int i = 1; i <= m; ++i)
    {
        scanf("%d%d", &u, &v);
        addEdge(i, u, v);
        addEdge(i, v, u);
    }
}

猜你喜欢

转载自blog.csdn.net/Floraqiu/article/details/82023985