pku3177 （边）双连通分量 求桥

初学双连通分量，刷了第一道题，纪念！

题意：有F个牧场，R条道路，形如<A, B>，表示牧场A到牧场B有一条通路（双向的），给出路保证每两个牧场都有通路，有些cow每次都得走同样的路有A到B。

问：再建最少的路使得任意两个牧场间都有两条不同的路，不同，指的是由A到B的过程中没有经过任一同一段路。

分析：双连通分量，正好保证了切断任何一条边，都不影响分量中任何两个顶点的连通性，也就是说双连通分量的连通度 >= 2，先dfs找到所有的桥，删除后，对图进行缩点，将所有双连通分支都缩成一点，再将这些桥加上，构成的图必然是一棵连通树，就是说，任意两个顶点间有且只有一条路可通。这样只需找到有多少叶子节点，即度为1的顶点，个数设为leaf，则需要建的边最少为(leaf + 1) / 2，特别的如果leaf = 1，就是说整个图是个双连通图，这是就不用加边了，ans = 0 即可。

#include<stdio.h>
#include<string.h>
#define MM 20004// 边最大值
#define NN 5004 // 点最大值
int F, R;
int idx; // 边数
int bCnt;// 连通分支个个数 branch cnt
int time;// 记录时间戳

int low[NN]; // Low(u)为u或u的子树中能通过非父子边追溯到的最早的节点，即DFS序号最小的节点
int dfn[NN]; // 深度优先标号
int bel[NN]; // bel[u] = k，节点u属于分支k
int deg[NN]; // 记录每个分支的度

typedef struct node{
   int v;
   int vis;
   struct node *nxt, *op;        
}NODE;

NODE edg[MM];
NODE *link[NN];
NODE *par[NN]; // 记录到达顶点v的树边

int Min(int a, int b){
    return a < b ? a : b;
}

void Add(int u, int v){ // 加边
     idx++;
     edg[idx].v = v;
     edg[idx].vis = 0;
     edg[idx].nxt = link[u];
	 link[u] = edg + idx;
	 edg[idx].op = edg + idx + 1; // 表示反向边
	 idx++;
     edg[idx].v = u;
     edg[idx].vis = 0;
     edg[idx].nxt = link[v];
	 link[v] = edg + idx;
     edg[idx].op = edg + idx - 1;
}

void dfs(int u){// 建深搜树，类似Tarjan求强连通分支
     int v;
     dfn[u] = low[u] = ++time;
     for (NODE *p = link[u]; p; p = p->nxt){
         if (p->vis) continue;
         p->vis = p->op->vis = 1; // 由于加的是双向边，所以正向边访问后，反向边也要标记
         v = p->v;
         if (!dfn[v]){
            par[v] = p->op; // 由于dfs作用是建一棵深搜树，所以每个点的前边是是唯一的
            dfs(v);
            low[u] = Min(low[u], low[v]);
         }else{
            low[u] = Min(low[u], dfn[v]);
         }
     }
}

void Part(int u){// 将每个点归入对于分支
     bel[u] = bCnt;
     for (NODE *p = link[u]; p; p = p->nxt){
         if (p->vis && !bel[p->v]){
            Part(p->v);
         }
     }
}

void Solve(){
     int i, u, v, x, y, leaf;
     time = 0;
     dfs(1);
     for (v = 2; v <= F; v++){// 删除桥
         u = (par[v])->v;
         if (dfn[u] < low[v]){// 桥边<u, v>的条件
            (par[v])->vis = (par[v])->op->vis = 0; 
         }
     }
     bCnt = 0;
     for (u = 1; u <= F; u++){
         if (!bel[u]){
            bCnt++;
            Part(u);
         }
     }
     memset(deg, 0, sizeof(deg));
     for (u = 1; u <= F; u++){ // 找到每个分支的度，任意两个属于不同分支节点的bel都不一样
         for (NODE *p = link[u]; p; p = p->nxt){
             x = bel[u]; y = bel[p->v];
             if (x != y){
                deg[x]++;
             }
         }
     }
     
     leaf = 0;
     for (i = 1; i <= bCnt; i++){// 记录度为1的分支
         if (deg[i] == 1) leaf++;
     }
     
     if (leaf == 1){
        puts("0");
     }else {
           printf("%d\n", (leaf + 1) / 2);
     }
}

int main()
{
    int a, b;
    scanf("%d%d", &F, &R);
    memset(link, 0, sizeof(link));
	idx = 0;
    while(R--){
       scanf("%d%d", &a, &b);
       Add(a, b);
    }
    Solve();
    return 0;
}

如果有讲的不详细的地方，围观大牛博客：beyond the void

转载于:https://www.cnblogs.com/ylfdrib/archive/2010/08/08/1795048.html