poj3694Network(lca+边双连通缩点)

POJ3694Network
这题意思就是给你一个计算机网络,所有计算机互相直接或间接相连,然后给你一些询问,询问每次添加一条边后还有多少桥边。
先把边双连通的部分缩点,然后就变成了一颗树
每一次询问就是删去LCA路径上的桥边,因为树上连一条边会形成一个环,那么环上的边就都不是桥
这题数据比较水,似乎是不考虑重边也能过,而且做法是O(qn)居然也能过.

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn=100010;
const int maxm=400010;
struct edge{
    int u,v,nxt;
    int bri;
}e[maxm];
int h[maxn];
int low[maxn],dfn[maxn],myst[maxn];
int vis[maxn],belong[maxn];
int f[maxn],dep[maxn],a[maxn];
int siz,idx,st,be,ans,bri;
void init(){
    memset(h,-1,sizeof(h));
    memset(e,-1,sizeof(e));
    memset(low,0,sizeof(low));
    memset(dfn,0,sizeof(dfn));
    memset(vis,0,sizeof vis);
    memset(a,0,sizeof(a));
    siz=idx=st=be=bri=0;
}
void add(int u,int v){
    e[siz].nxt=h[u];
    e[siz].u=u;
    e[siz].v=v;
    h[u]=siz++;
}
void tarjan(int u,int fa){
    int v;
    vis[u]++;
    myst[++st]=u;
    dfn[u]=low[u]=++idx;
    int k=0;
    for(int i=h[u];~i;i=e[i].nxt){
        v=e[i].v;
        if(v==fa&&k==0){k++;continue;}
        if(!vis[v]){
            tarjan(v,u);
            if(low[v]<low[u])low[u]=low[v];
            if(low[v]>dfn[u]){
                bri++;
                e[i].bri=1;
                e[i^1].bri=1;
            }
        }
        else if((vis[v]==1)&&(low[u]>dfn[v])){
            low[u]=dfn[v];
        }
    }
    if(low[u]==dfn[u]){
        be++;
        do{
            v=myst[st--];
            vis[v]++;
            belong[v]=be;
        }while(v!=u);
    }
}
vector<int>mp[100010];
void dfs(int rt,int fa){
    f[rt]=fa;
    dep[rt]=dep[fa]+1;
    //printf("be:%d fa:%d dep:%d\n",rt,fa,dep[rt]);
    int sz=mp[rt].size();
    for(int i=0;i<sz;i++){
        int v=mp[rt][i];
        if(v==fa)continue;
        dfs(v,rt);
    }
}
void lca(int u,int v){
    if(dep[u]<dep[v])swap(u,v);
    while(dep[u]>dep[v]){
        if(!a[u]){a[u]=1;ans--;}
        u=f[u];
    }
    while(u!=v){
        if(!a[u]){a[u]=1;ans--;}u=f[u];
        if(!a[v]){a[v]=1;ans--;}v=f[v];
    }
}
void gao(int m){
    for(int i=1;i<=be;i++){
        mp[i].clear();
    }
    for(int i=0;i<2*m;i++){
        int u=belong[e[i].u];
        int v=belong[e[i].v];
        if(u!=v)mp[u].push_back(v);
    }
    dep[1]=0;
    dfs(1,1);
}
int main(){
    int n,m;
    int cas=0;
    while(scanf("%d%d",&n,&m)&&(n||m)){
        int q;
        init();
        for(int i=1;i<=m;i++){
            int u,v;
            scanf("%d%d",&u,&v);
            add(u,v);
            add(v,u);
        }
        tarjan(1,0);
    /*  for(int i=1;i<=n;i++){
            printf("%d %d\n",i,belong[i]);
        }*/
        ans=be-1;
        gao(m);
/*      for(int i=1;i<=n;i++){
            printf("i:%d be:%d\n",i,belong[i]);
        }
        for(int i=1;i<=be;i++){
            printf("be: %d fa: %d\n",i,f[i]);
        }*/
        scanf("%d",&q);
        printf("Case %d:\n",++cas);
        for(int i=1;i<=q;i++){
            int u,v;
            scanf("%d%d",&u,&v);
            lca(belong[u],belong[v]);
            printf("%d\n",ans);
        }
        puts("");
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值