【图论】割点,桥,边双联通

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<string>
using namespace std;
const int maxn=20010;
struct Edge{
    int from,to;
//    不带边权可删除var
    int var;
    bool is_bridge;
};
vector<int>G[maxn];
vector<Edge>edges;
vector<int>Set;
int dfs_clock;//初始化为0
int pre[maxn];//初始化为-1
int low[maxn];
bool cut_point[maxn];
int color[maxn];
bool vis[maxn];
int res;

void add(int a,int b,int c){
    Edge e;
    e.from=a;e.to=b;e.var=c;e.is_bridge=false;
    edges.push_back(e);
    G[a].push_back(edges.size()-1);
    e.from=b;e.to=a;
    edges.push_back(e);
    G[b].push_back(edges.size()-1);
}


// u为当前点
// fa为父亲节点
// EdgeNo为到当前点的边的序号
int dfs(int u,int fa,int EdgeNO){
    int lowu=pre[u]=++dfs_clock;
    int child=0;
    for(int i=0;i<G[u].size();i++){
        int v=edges[G[u][i]].to;
        if(pre[v]==-1){
            child++;
            int lowv=dfs(v,u,G[u][i]);
            lowu=min(lowv,lowu);
            if(lowv>=pre[u]){
                cut_point[u]=true;
                if(lowv>pre[u]){
                    res++;
                    edges[G[u][i]].is_bridge=true;
                    edges[G[u][i]^1].is_bridge=true;
                }
            }
        }
        else if(pre[v]<pre[u]&&G[u][i]!=(EdgeNO^1))
            lowu=min(pre[v],lowu);
    }
    if(fa<0&&child==1)cut_point[u]=0;
    low[u]=lowu;
    return lowu;
}

int Find_Color(int u,int fa){
    int minC=u;
    Set.push_back(u);
    for(int i=0;i<G[u].size();i++){
        int v=edges[ G[u][i] ].to;
        if(edges[G[u][i]].is_bridge)continue;
        if(vis[v])continue;
        vis[v]=true;
        int col=Find_Color(v,u);
        minC=min(minC,col);
    }
    return minC;
}
void Paint(int n){
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++){
        if(vis[i])continue;
        vis[i]=true;
        int col=Find_Color(i,-1);
        for(int i=0;i<Set.size();i++)
            color[Set[i]]=col;
        Set.clear();
    }
}

void init(){
    res=1;
    dfs_clock=0;
    for(int i=0;i<maxn;i++)
        G[i].clear();
    edges.clear();
    memset(pre,-1,sizeof(pre));
    memset(cut_point,0,sizeof(cut_point));
    memset(color,-1,sizeof(color));
}

/**************************************
用add 完成建图之后
初始化
调用dfs()函数找出割点和桥
调用Paint(n)进行边双联通着色
结点序号从1开始
输出:
cut_point[i]       判断i是否为割点
edges[j].is_bridge 判断j边是否是桥(反向边也有标记)
color[i]           判断i属于哪个集合
***************************************/

int main(){
//    freopen("C:\\Users\\金柯\\Desktop\\input.txt","r",stdin);
    int n,m;
    while(cin>>n>>m){
        init();
        for(int i=1;i<=m;i++){
            int a,b;
            cin>>a>>b;
            add(a,b,0);
        }
        dfs(1,-1,-1);
        Paint(n);
        cout<<res<<endl;
        for(int i=1;i<=n;i++){
            cout<<color[i]<<" ";
        }
        cout<<endl;
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值