POJ-3694 Network

给出一个图,还有q个询问,每个询问给出一条边,问加上这条边之后,图中还有几个桥

给出图后,先用Tarjan算法过一遍,中间对两个数组和一个值进行标记

father数组保存父节点,bri保存这个点到父节点的边是不是桥,cnt保存桥的个数

这样结束后我们得到了一个树形结构,所有的桥都在这个树中作为树边存在

每当添加一条边时,树中就会形成一个环,这个环上的所有桥都不再会是桥

接下来对于每个添加的边(u,v)进行lca算法(最近公共祖先),路径上就可以把环上的桥都找出来,标记并统计个数

详见代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
// Tarjan算法   复杂度O(N+M)
const int MAXN = 100010;//点数
const int MAXM = 400010;//边数
struct Edge
{
	int to,next;
}edge[MAXM];
int head[MAXN],tot;
int Low[MAXN],DFN[MAXN],Stack[MAXN],deep[MAXN],father[MAXN];
int Index,top;
int cnt;
bool Instack[MAXN],bri[MAXN];
void addedge(int u,int v)
{
	edge[tot].to = v;edge[tot].next = head[u];head[u] = tot++;
}
void Tarjan(int u)
{
	int v;
	Low[u] = DFN[u] = ++Index;
	deep[u]=deep[father[u]]+1;
	Stack[top++] = u;
	Instack[u] = true;
	for(int i = head[u];i != -1;i = edge[i].next)
	{
		v = edge[i].to;
		if(v==father[u]) continue;
		if( !DFN[v] )
		{
            father[v]=u;
			Tarjan(v);
			if( Low[u] > Low[v] )Low[u] = Low[v];
			if(Low[v]>DFN[u])
			{
                cnt++;
                bri[v]=true;
			}
		}
		else if(Instack[v] && Low[u] > DFN[v])
		Low[u] = DFN[v];
	}
    top--;
    Instack[u] = false;
}
void solve(int N)
{
	memset(DFN,0,sizeof(DFN));
	memset(deep,0,sizeof(deep));
	memset(Instack,false,sizeof(Instack));
	memset(bri,false,sizeof(bri));
	for(int i=1;i<=N;i++) father[i]=i;
	Index = cnt = top = 0;
	for(int i = 1;i <= N;i++)
	if(!DFN[i]) Tarjan(i);
}
void init()
{
	tot = 0;
	memset(head,-1,sizeof(head));
}
void lca(int u,int v)
{
    while(deep[u]>deep[v])
    {
        if(bri[u]) cnt--,bri[u]=false;
        u=father[u];
    }
    while(deep[v]>deep[u])
    {
        if(bri[v]) cnt--,bri[v]=false;
        v=father[v];
    }
    while(u!=v)
    {
        if(bri[u]) cnt--,bri[u]=false;
        if(bri[v]) cnt--,bri[v]=false;
        u=father[u];
        v=father[v];
    }
}
int main()
{
    //freopen("/home/zlwang/test.txt","r",stdin);
    int n,m,q,a,b;
    int kase=0;
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0) break;
        init();
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a,b);
            addedge(b,a);
        }
        solve(n);
        scanf("%d",&q);
        printf("Case %d:\n",++kase);
        while(q--)
        {
            scanf("%d%d",&a,&b);
            lca(a,b);
            printf("%d\n",cnt);
        }
        printf("\n");
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值