Tarjan(边的双连通分量) - Redundant Paths - POJ 3177

Tarjan(边的双连通分量) - Redundant Paths - POJ 3177

为了从F个草场中的一个走到另一个，奶牛们有时不得不路过一些她们讨厌的可怕的树。

奶牛们已经厌倦了被迫走某一条路，所以她们想建一些新路，使每一对草场之间都会至少有两条相互分离的路径，这样她们就有多一些选择。

每对草场之间已经有至少一条路径。

给出所有R条双向路的描述，每条路连接了两个不同的草场，请计算最少的新建道路的数量,路径由若干道路首尾相连而成。

两条路径相互分离，是指两条路径没有一条重合的道路。

但是，两条分离的路径上可以有一些相同的草场。

对于同一对草场之间，可能已经有两条不同的道路，你也可以在它们之间再建一条道路，作为另一条不同的道路。

输入格式

第1行输入F和R。

接下来R行，每行输入两个整数，表示两个草场，它们之间有一条道路。

输出格式

输出一个整数，表示最少的需要新建的道路数。

数据范围

$1\le F\le 5000,F-1\le R\le 10000$

输入样例：

7 7
1 2
2 3
3 4
2 5
4 5
5 6
5 7

输出样例：

2

---

分析：

$边的双连通分量：任意两点间存在至少两条不相交的路径。即不存在桥。$

$题意：在一个无向图中至少添加几条边，可以使得整个图是一个边的双连通图。$

$tarjan算法先缩点，那么在每个连通分量中，要使得各点之间连通，至少要存在一个环。$

$有环就满足边双连通分量的定义。$

$那么缩点后得到的新图就是一棵树。$

$树中边均为原图中的桥。$

如何判断桥：

$条件：dfn[u]<low[j]$

$即遍历到u的时间戳要严格小于u的孩子j向下遍历的最小时间戳。$

$也就是从j向下遍历，无论如何都不会返回到u。如下图：$

$大圆表示连通分量，小圆表示点。$

$缩点后的树中，根到叶子节点是不存在两条不相交的路径的。$

$因此我们需要在叶子节点与其父节点之间新建一条边。$

结论：$\lfloor \frac{叶子的数量}{2}\rfloor$

$缩点后，我们统计度为1的连通分量的数量即可。$

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;

const int N=5010, M=20010;

int n,m;
int e[M],ne[M],h[N],idx;
int stk[N],top;
bool is_bridge[M];
int dfn[N],low[N],timestamp;
int id[N],dcc_cnt;
int d[N];

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

void tarjan(int u,int from)
{
    dfn[u]=low[u]=++timestamp;
    stk[++top]=u;
    for(int i=h[u];~i;i=ne[i])
    {
        int j=e[i];
        if(!dfn[j])
        {
            tarjan(j,i);
            low[u]=min(low[u],low[j]);
            if(dfn[u]<low[j])
                is_bridge[i]=is_bridge[i^1]=true;
        }
        else if(i!=(from ^ 1)) 
            low[u]=min(low[u],dfn[j]);
    }
    
    if(dfn[u]==low[u])
    {
        ++dcc_cnt;
        int y;
        do
        {
            y=stk[top--];
            id[y]=dcc_cnt;
        }while(y!=u);
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    
    memset(h,-1,sizeof h);
    int a,b;
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&a,&b);
        add(a,b),add(b,a);
    }
    
    tarjan(1,-1);
    
    for(int i=0;i<idx;i++)
        if(is_bridge[i])
            d[id[e[i]]]++;
            
    int cnt=0;
    for(int i=1;i<=dcc_cnt;i++)
        if(d[i]==1)
            cnt++;
    
    printf("%d\n",(cnt+1)/2);
    
    return 0;
}