hdu2460Network

链接:http://acm.hdu.edu.cn/showproblem.php?pid=2460

题意:给出一个n个点m条无向边,然后逐步加入一些边,求每次加入当前边后图内还剩多少条割边。

分析:先用tarjan缩点,将图缩成一颗树,然后在树上加个边权每次缩到lca就好了。详见代码。

代码:

#include<map>
#include<set>
#include<stack>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1e5+10;
const int M=2e5+10;
const int mod=1000000007;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int INF=1000000010;
const ll MAX=1ll<<55;
const double eps=1e-5;
const double inf=~0u>>1;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned int uint;
typedef unsigned long long ull;
struct Edge {
    int u,v,id,nex;
    Edge() {}
    Edge(int u,int v,int id,int nex) {
        this->u=u;this->v=v;this->id=id;this->nex=nex;
    }
}edge[M<<1];
stack<int>S;
int tot,first[N];
int sccn,dfs_clock,pre[N],dfn[N],sccno[N],vis[M];
void dfs_scc(int u) {
    pre[u]=dfn[u]=++dfs_clock;S.push(u);
    for (int i=first[u];~i;i=edge[i].nex)
    if (!vis[edge[i].id]) {
        vis[edge[i].id]=1;
        if (!pre[edge[i].v]) {
            dfs_scc(edge[i].v);
            dfn[u]=min(dfn[u],dfn[edge[i].v]);
        } else if (!sccno[edge[i].v]) dfn[u]=min(dfn[u],dfn[edge[i].v]);
    }
    if (pre[u]==dfn[u]) {
        sccn++;
        while (1) {
            int x=S.top();
            sccno[x]=sccn;S.pop();
            if (x==u) break ;
        }
    }
}
void find_scc(int n,int m) {
    dfs_clock=sccn=0;
    for (int i=1;i<=m;i++) vis[i]=0;
    for (int i=1;i<=n;i++) pre[i]=sccno[i]=0;
    for (int i=1;i<=n;i++)
    if (!pre[i]) dfs_scc(i);
}
int dep[N],fa[N];
vector<int>V[N];
void dfs(int u,int v) {
    dep[u]=dep[v]+1;fa[u]=v;
    for (int i=0;i<V[u].size();i++)
    if (V[u][i]!=v) dfs(V[u][i],u);
}
int get(int u,int v) {
    int a=u,b=v,ret=0;
    while (a!=b)
    if (dep[a]>dep[b]) ret+=vis[a],a=fa[a];
    else ret+=vis[b],b=fa[b];
    a=u;b=v;
    while (a!=b)
    if (dep[a]>dep[b]) vis[a]=0,a=fa[a];
    else vis[b]=0,b=fa[b];
    return ret;
}
int main()
{
    int u,v,i,n,m,q,ca=0,ans;
    while (scanf("%d%d", &n, &m)&&(n+m)) {
        for (tot=0,i=1;i<=n;i++) first[i]=-1;
        for (i=1;i<=m;i++) {
            scanf("%d%d", &u, &v);
            edge[tot]=Edge(u,v,i,first[u]);first[u]=tot++;
            edge[tot]=Edge(v,u,i,first[v]);first[v]=tot++;
        }
        find_scc(n,m);
        for (i=1;i<=sccn;i++) V[i].clear();
        for (i=0;i<tot;i+=2)
        if (sccno[edge[i].u]!=sccno[edge[i].v]) {
            V[sccno[edge[i].u]].push_back(sccno[edge[i].v]);
            V[sccno[edge[i].v]].push_back(sccno[edge[i].u]);
        }
        dfs(1,0);
        scanf("%d", &q);ans=sccn-1;
        for (i=2;i<=sccn;i++) vis[i]=1;
        printf("Case %d:\n", ++ca);
        while (q--) {
            scanf("%d%d", &u, &v);
            u=sccno[u];v=sccno[v];
            ans-=get(u,v);
            printf("%d\n", ans);
        }
        printf("\n");
    }
    return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值