kuangbin专题九 UVA796 Critical Links(无向图求割边的算法模板题)

题解:
无向图求割边的算法模板,要你输出割边数和割边
从网上学了两个模板,推荐使用第一个。

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
const int MAXN=100005;
struct Edge{
    int to,next;
    bool cut;//是否为桥的标记 
}edge[MAXN];
int head[MAXN],tot;
int low[MAXN],dfn[MAXN],stack[MAXN];
int Index,top;
bool Instack[MAXN];
bool cut[MAXN];
int add_block[MAXN];//删除一个点后增加的连通块,这里还不太懂,等到时候找道题来做一下。 
int bridge;
int n,m;
void addedge(int u,int v)
{
    edge[tot].to=v;
    edge[tot].next=head[u];
    edge[tot].cut=false;
    head[u]=tot++; 
}
void init()
{
    memset(head,-1,sizeof(head));
    memset(dfn,0,sizeof(dfn));
    memset(Instack,false,sizeof(Instack));
    memset(add_block,0,sizeof(add_block));
    memset(cut,false,sizeof(cut));
    tot=0;
    Index=top=0;
    bridge=0;
}
void Tarjan(int u,int pre)
{
    low[u]=dfn[u]=++Index;
    stack[top++]=u;
    Instack[u]=true;
    int son=0;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(v==pre) continue;
        if(!dfn[v])
        {
            son++;
            Tarjan(v,u);
            low[u]=min(low[u],low[v]);
            //桥
            //一条无向边(u,v)是桥,当且仅当(u,v)为树枝边,且满足DFS(u)<Low(v)。          
            if(low[v]>dfn[u])
            {
                bridge++;
                edge[i].cut=true;//这两行的操作相当于可以找到这条边的反向边,但是要从0开始存放,至于为什么,你可以这样想第一条边存放是不是0,1,第二条是不是2,3,你用2^1和3^1就知道怎么回事了。 
                edge[i^1].cut=true;
            }
            //割点
            //一个顶点u是割点,当且仅当满足(1)或(2) (1) u为树根,且u有多于一个子树。
            //(2) u不为树根,且满足存在(u,v)为树枝边(或称父子边,
            //即u为v在搜索树中的父亲),使得DFS(u)<=Low(v)            
            if(u!=pre&&low[v]>=dfn[u])
            {
                cut[u]=true;
                add_block[u]++;
            }
        }
        else if(low[u]>dfn[v])
        low[u]=dfn[v];
    }
    //树根,分支数大于1
    if(u==pre&&son>1) cut[u]=true;
    if(u==pre) add_block[u]=son-1;//这个操作我还不是很理解。 
    Instack[u]=false;
    top--;
}
void solve()
{
    for(int i=1;i<=n;i++)
    if(!dfn[i])
    Tarjan(i,i);
    printf("%d critical links\n",bridge);
    vector< pair<int,int> >ans;//这里的 vector< pair<int,int> >ans;中的 int> >ans必须分开,因为系统会认为是位运算>>ans。所以必须空格弄开。 
    for(int u=1;u<=n;u++)
        for(int i=head[u];i!=-1;i=edge[i].next)
        if(edge[i].cut && edge[i].to>u){//edge[i].to>u可以去重,不懂的话删去试试。 
            ans.push_back(make_pair(u,edge[i].to));
        }
    sort(ans.begin(),ans.end());
    for(int i=0;i<ans.size();i++)
    printf("%d - %d\n",ans[i].first-1,ans[i].second-1);//因为是从1开始的,题目的编号是从0开始的。 
    printf("\n"); 
}
int main()
{
    while(~scanf("%d",&n))
    {
        init();
        int u,v,m;
        for(int i=1;i<=n;i++)
        {
            scanf("%d (%d)",&u,&m);
            u++;//因为题目要求的是从0开始的,做题是从1开始的。所以加1,下面也一样。 
            while(m--)
            {
                scanf("%d",&v);
                v++;
                addedge(u,v);
                addedge(v,u);
            }
        }
        solve(); 
    }
    return 0;
} 
#include<stdio.h>
#include<string.h>
#include<vector>
#include<string.h>
#include<algorithm> 
using namespace std;
#define MAXN 1010
struct node
{
    int x,y;
}a[MAXN];
vector<int>map[MAXN];
int dfn[MAXN];
int low[MAXN];
int fa[MAXN]; 
int Index;
int n;
bool cmp(node c,node d)
{
    if(c.x!=d.x)
    return c.x<d.x;
    else
    return c.y<d.y;
}
void init()
{
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(fa,0,sizeof(fa));
    for(int i=0;i<n;i++) map[i].clear();
    Index=0;    
}
void Tarjan(int u,int pre)
{
    low[u]=dfn[u]=++Index;
    fa[u]=pre;
    for(int i=0;i<map[u].size();i++)
    {
        int v=map[u][i];
        if(!dfn[v])
        {
            Tarjan(v,u);
            low[u]=min(low[u],low[v]);
        }
        else if(pre!=v)
        low[u]=min(low[u],dfn[v]);
    }
}
void slove()
{
    int ans=0;
    for(int i=0;i<n;i++)
    {
        if(!dfn[i])
        Tarjan(i,-1);
    }
    for(int i=0;i<n;i++)
    {
        int v=fa[i];
        if(v!=-1&&low[i]>dfn[v])//父节点为-1说明该点为根节点,与其父节点不可能存在桥
        {                       //dfn[ni]==low[i]说明i还有其他路径到ni,则说明i-ni不属于桥
            a[ans].x=i;
            a[ans].y=v;
            int t;
            if(a[ans].x>a[ans].y)
            {
                t=a[ans].x;
                a[ans].x=a[ans].y;
                a[ans].y=t;
            }
            ans++;
        }
    }
    sort(a,a+ans,cmp);
    printf("%d critical links\n",ans);
    for(int i=0; i<ans; i++)
        printf("%d - %d\n",a[i].x,a[i].y);
    printf("\n");
}
int main()
{
    while(~scanf("%d",&n))
    {
        init();
        int u,v,m;
        for(int i=0;i<n;i++)
        {
            scanf("%d (%d)",&u,&m);
            while(m--)
            {
                scanf("%d",&v);
                map[u].push_back(v);
                map[v].push_back(u);
            }
        }
        slove();
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值