割边判定法则 搜索桥与求e-DCC(边双连通分量) 缩点

若(x,y)是桥,当且仅当 y是x的子节点,且dfn[x]<low[y]

#include<iostream>
using namespace std;
const int N =10010;
int head[N],nex[N<<1],e[N<<1],tot;
int dfn[N],low[N],n,m,num;
bool bridge[N<<1];

void add(int x,int y){
    e[++tot]=y,nex[tot]=head[x],head[x]=tot;
}

void tarjan(int x,int f){
    dfn[x]=low[x]=++num;//时间戳 并且初始化low
    for(int i=head[x];i;i=nex[i]){
        int y=e[i];
        if(!dfn[y]){//如果此点没被经过,那么向下递归
            tarjan(y,i);
            low[x]=min(low[x],low[y]);//向前迭代回溯
            if(low[y]>dfn[x])//割边判定法则
                bridge[i]=bridge[i^1]=1;
        }
        else if (i!=(f^1))//如果和父节点有重边,可以更新
            low[x]=min(low[x],dfn[y]);
    }
}


int main()
{
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int x,y;
        cin>>x>>y;
        add(x,y),add(y,x);
    }
    for(int i=1;i<=n;i++)//多个图遍历
        if(!dfn[i]) tarjan(i,0);
    for(int i=1;i<tot;i+=2)
        if(bridge[i])
            cout<<e[i^1]<<' '<<e[i]<<endl;
    return 0;
}

求e-DCC

int c[N],dcc;
void dfs(int x){
    c[x]=dcc;
    for(int i=head[x];i;i=nex[i]){
        int y=e[i];
        if(c[y]||bridge[i]) continue;
        dfs(y);
    }
}

int main(){
	for(int i=1;i<=n;i++){
	    if(!c[i]){
	        ++dcc;
	        dfs(i);
	    }
	}
}

缩点


int hc[N],vc[N<<1],nc[N<<1],tc;
void add_c(int x,int y){
    vc[++tc]=y,nc[tc]=hc[x],hc[x]=tc;
}

int main()
{
	tc=1;//缩点
	for(int i=1;i<=tot;i++){
 	   int x=e[i^1],y=e[i];
  	  if(c[x]==c[y]) continue;
  	  add_c(c[x],c[y]);
	}
}

综合

#include<iostream>
using namespace std;
const int N =10010;
int head[N],nex[N<<1],e[N<<1],tot;
int dfn[N],low[N],n,m,num;
bool bridge[N<<1];

void add(int x,int y){
    e[++tot]=y,nex[tot]=head[x],head[x]=tot;
}

void tarjan(int x,int f){
    dfn[x]=low[x]=++num;//时间戳 并且初始化low
    for(int i=head[x];i;i=nex[i]){
        int y=e[i];
        if(!dfn[y]){//如果此点没被经过,那么向下递归
            tarjan(y,i);
            low[x]=min(low[x],low[y]);//向前迭代回溯
            if(low[y]>dfn[x])//割边判定法则
                bridge[i]=bridge[i^1]=1;
        }
        else if (i!=(f^1))//如果和父节点有重边,可以更新
            low[x]=min(low[x],dfn[y]);
    }
}

int c[N],dcc;
void dfs(int x){
    c[x]=dcc;
    for(int i=head[x];i;i=nex[i]){
        int y=e[i];
        if(c[y]||bridge[i]) continue;
        dfs(y);
    }
}

int hc[N],vc[N<<1],nc[N<<1],tc;
void add_c(int x,int y){
    vc[++tc]=y,nc[tc]=hc[x],hc[x]=tc;
}


int main()
{
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int x,y;
        cin>>x>>y;
        add(x,y),add(y,x);
    }
    for(int i=1;i<=n;i++)//多个图遍历
        if(!dfn[i]) tarjan(i,0);
    for(int i=1;i<tot;i+=2)//求桥
        if(bridge[i])
            cout<<e[i^1]<<' '<<e[i]<<endl;

    for(int i=1;i<=n;i++){//求e-DCC
        if(!c[i]){
            ++dcc;
            dfs(i);
        }
    }

    tc=1;//缩点
    for(int i=1;i<=tot;i++){
        int x=e[i^1],y=e[i];
        if(c[x]==c[y]) continue;
        add_c(c[x],c[y]);
    }


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值