双连通学习笔记

DFS树博客：链接
双连通是对无向图所说的

边双：对于一个图G，和子图G2，如果去除图中任意一个边子图G2仍然连通的话，G2中的边就是边双联通

1. 关于边双点双桥割点的DFS做法
   跟OIwiki学的：链接

对图跑出一颗DFS树。
如图，图源OIwiki，黑边绿边为树边，红边为非树边。被红边所关照两点的树上路径必然不是桥，除此之外的就都是桥了，所以求桥的方法就是用差分。对于非树边u-v, 假设v为更深的树点，那么开个差分数组a[v]++, a[u]–;在DFS回溯的过程中差分，这样就标记好了。两点是否边双连通，就看树上路径是否有桥

如下图，黑边为树边，红遍为非树边，将非树边两点所对应的树上路径加入一个联通块。这样如果一个点周边的所有的边都属于一个联通块，那么它不是割点，此外就是。点双连通，条件为两点间的路径上的边都为一个联通块。

2. Tarjan做法
A. 判断割点：
我们开两个数组一个为dfn时间戳代表我们dfs访问点的顺序，另一个为low，这个数组不知道该怎么定义，我觉得可以这样说，它所连接的路径直接返回时间戳最小值。对于判割点的算法应该是这样的，对于求scc的按板子定义吧。
这样的话一个点u直接相连的点v如果有一个不能返回到比dfn[u]小的时间戳，那么该点为割点，如果这个点是根节点的话，需要特判，当它的儿子数>=2的时候它也是割点。
例题1：洛谷P3388求哪些是割点
代码：

#include <bits/stdc++.h>
using namespace std;
#define IO ios::sync_with_stdio(false);cin.tie(0)
#define forn(i, n) for(int i = 0; i < n; ++i)
#define for1(i, n) for(int i = 1; i <= n; ++i)

const int maxn = 2e4 + 5;
int ans, ind;
int dfn[maxn], low[maxn];
bool pcut[maxn];
vector<int>e[maxn],res;

inline void tarjan(int u, int pre) {
   
    dfn[u] = low[u] = ++ind;
    int son = 0;
    for(auto v : e[u]) {
   
        if(!dfn[v]) {
   
            ++son;
            tarjan(v, u);
            low[u] = min(low[u], low[v]);
            if(u != pre && !pcut[u] && low[v] >= dfn[u]) {
   
                pcut[u] = 1;
                ++ans;
                res.push_back(u);
            }
        }else if(v != pre) low[u] = min(low[u], dfn[v]);
    }
    if(u == pre && !pcut[u] && son > 1) {
   
        pcut[u] = 1;
        ++ans;
        res.push_back(u);
    } 
}

int main() {
   
    IO;
    int n, m; cin >> n >> m;
    forn(i, m) {
   
        int u, v; cin >> u >> v;
        e[u].push_back(v);
        e[v].push_back(u);
    }
    for1(i, n) if(!dfn[i]) tarjan(i, i);
    cout << ans << '\n';
    sort(res.begin(), res.end());
    for(auto x : res) cout << x << ' ';
    return 0;
}

例题2：POJ2117求一个图删除一个点之后，联通块最多有多少。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<cmath>
#include<sstream>
using namespace std;
#define IO ios::sync_with_stdio(false);cin.tie(0)
#define forn(i, n) for(int i = 0; i < n; ++i)
#define for1(i, n) for(int i = 1; i <= n; ++i)

const int maxn = 1e4 + 5;
int ans, ind, ma;
int dfn[maxn], low[maxn], pcut[maxn];
vector<int>e[maxn];
inline void init(int n) {
   
    ans = 0, ind = 0, ma = 0;
    forn(i, n + 5) {
   
        e[i].clear();
        dfn[i] = pcut[i] = 0;
    }
}

inline void tarjan(int u, int pre) {
   
    dfn[u] = low[u] = ++ind;
    int son = 0;
    for(int i = 0; i < e[u].size(); ++i) {
   
        int v = e[u][i];
        if(!dfn[v]) {
   
            ++son;
            tarjan(v, u);
            low[u] = min(low[u], low[v]);
            if(u != pre && low[v] >= dfn[u]) {
   
                //cerr << "!@ " <<  u << ' '<< v << '\n';
                ++pcut[u];
                ma = max(ma, pcut[u]);
            }
        }else if(v != pre) low[u] = min(low[u], dfn[v]);
    }
    if(u == pre && son > 1) {
   
        //cerr << "!@2 " <<  u  << '\n';
        pcut[u] = son - 1;
        ma = max(ma, pcut[u]);
    }
}

int main() {
   
    IO;
    int n, m;while(cin >> n >> m && (n || m)) {
   
        init(n);
        forn(i, m) {
   
            int u, v; cin >> u >> v;
            ++u, ++v;
            e[u].push_back(v);
            e[v].push_back(u);
        }
        for1(i, n) if(!dfn[i]) tarjan(i, i), ++ans;
        ans += ma