桥——无向图边的双连通分量

求解加多少边才能将整个图变成边双联通分量(边的双连通分量：极大的不包含桥的连通块)

#include<iostream>
#include<cstring>
#include<algorithm>
#include<stack>
using namespace std;
const int N = 5010,M = 20010;
int h[N],ne[M],e[M],idx;
int dnt[N]; //dnt[u]表示当前到达u
int low[N]; //low[u]表示从u节点出发能遍历到的最小时间戳
int id[N];  //记录节点i所在的双连通分量的编号
int d[N]; //记录节点i的度
bool is_bridge[N]; //记录边是否是桥边
stack<int> sta;
int dcc_cnt; //为双连通分量的编号
int timestamp; //时间戳
int n,m;

void add(int a,int b){
    e[idx] = b;
    ne[idx] = h[a];
    h[a] = idx++;
}

void tarjan(int u,int from){   //from记录的是当前节点由哪条边过来的(防止反向遍历)

    low[u] = dnt[u] = ++timestamp;

    sta.push(u);

    //先将图的所有节点的low和dnt都预处理出来
    for(int i = h[u];i != -1;i = ne[i]){
        int j = e[i];

        if(!dnt[j]){

            tarjan(j,i); 

            low[u] = min(low[u],low[j]);

            //表示j节点永远都走不到u节点(u,j这两个点只能从u走到j)，所以边u-j(w[i])是一条桥边
            if(dnt[u] < low[j]){  
                is_bridge[i] = true;
                is_bridge[i ^ 1] = true; //反向边同样是桥边
                //因为一对正向边和反向边的编号都是成对出现，所以直接^1即可
                // 这里i==idx 如果idx==奇数 则反向边=idx+1 = idx^1
                //            如果idx==偶数 则反向边=idx-1 = idx^1
            } 

        //from ^ 1是边from的反向边，防止再遍历回去
        }else if(i != (from ^ 1)){

            //j节点有可能是之前遍历过的点
            low[u] = min(low[u],dnt[j]);
        }
    }


    //预处理完之后判断哪个节点的low[u] = dnt[u],说明是双联通分量的最高点
    if(low[u] == dnt[u]){

        int y;
        ++dcc_cnt;
        do{
            y = sta.top();
            sta.pop();
            id[y] = dcc_cnt;
        }while(y != u);
    }

}

int main(){
    memset(h,-1,sizeof h);
    cin >> n >> m;
    for(int i = 0;i < m;++i){
        int a,b;
        cin >> a >> b;
        add(a,b); add(b,a);
    }

    for(int i = 1;i <= n;++i){
        if(!dnt[i]) tarjan(i,-1);  //引入from防止反向遍历回遍历过的边
    }

    //遍历所有边，如果边i是桥边,在其所连的出边的点j所在双连通分量的度+1
    for(int i = 0;i < idx;++i){   //包含正向边和反向边，所以度为1的节点一定是叶子节点
        if(is_bridge[i]){
            d[id[e[i]]]++;
        }
    }

    int cnt;
    //枚举所有双连通分量，需要加的边数就是:(度为1的节点个数 + 1 ) / 2
    for(int i = 1;i <= dcc_cnt;++i){
        if(d[i] == 1) cnt++;
    }
    cout << (cnt + 1) / 2 << endl;
    return 0;
}