hdu - 1269 迷宫城堡 强连通分量 tarjan模板题

有向图强连通分量:在有向图G中,如果两个顶点vi,vj间(vi>vj)有一条从vi到vj的有向路径,同时还有一条从vj到vi的有向路径,则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通,称G是一个强连通图。有向图的极大强连通子图,称为强连通分量(strongly connected components)。

任何一个强连通分量,必定是对原图的深度优先搜索树的子树。那么其实,我们只要确定每个强连通分量的子树的根,然后根据这些根从树的最低层开始,一个一个的拿出强连通分量即可。那么剩下的问题就只剩下如何确定强连通分量的根和如何从最低层开始拿出强连通分量了。
那么如何确定强连通分量的根,在这里我们维护两个数组,一个是indx[1..n],一个是mlik[1..n],其中indx[i]表示顶点i开始访问时间,mlik[i]为与顶点i邻接的顶点未删除顶点j的mlik[j]和mlik[i]的最小值(mlik[i]初始化为indx[i])。这样,在一次深搜的回溯过程中,如果发现mlik[i]==indx[i]那么,当前顶点就是一个强连通分量的根,为什么呢?因为如果它不是强连通分量的根,那么它一定是属于另一个强连通分量,而且它的根是当前顶点的祖宗,那么存在包含当前顶点的到其祖宗的回路,可知mlik[i]一定被更改为一个比indx[i]更小的值。
至于如何拿出强连通分量,如果当前节点为一个强连通分量的根,那么它的强连通分量一定是以该根为根节点的(剩下节点)子树。在深度优先遍历的时候维护一个堆栈,每次访问一个新节点,就压入堆栈。这样,由于当前节点是这个强连通分量中最先被压入堆栈的,那么在当前节点以后压入堆栈的并且仍在堆栈中的节点都属于这个强连通分量。可以用反证法证明这个做法的正确性。假设一个节点在当前节点压入堆栈以后压入并且还存在,同时它不属于该强连通分量,那么它一定属于另一个强连通分量,但当前节点是它的根的祖宗,那么这个强连通分量应该在此之前已经被拿出。

代码中数组dfn为上述indx,low为mlik

题意:中文题

按照题目的描述,相互连接是强连通的,所以判断一下是否是一个强连通图只有一个强连通分量就好了

tarjan算法模板题

链接:hdu 1269

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <set>
#include <map>
#include <stack>
#include <queue>
#define INF 0x3f3f3f3f
#define ll long long
 
using namespace std;
 
const int inf = 0x3f3f3f3f;
const int maxn = 10005;
const int maxm = 100005;
 
struct node {
    int v, next;
}edge[maxm];
 
int sstack[maxn], top;
int belong[maxn], scnt;
int dfn[maxn], low[maxn], cnt;
int head[maxn], tot;
bool instack[maxn];
int n, m;
 
void init() {
    cnt = 0;
    top = tot = scnt = 0; //栈顶指针为0 ,次序计数器,初始化连通分量标号
    memset(head, -1, sizeof(head));
    memset(dfn, 0, sizeof(dfn)); //结点搜索的次序编号数组为0,同时可以当是否访问的数组使用
}
 
void addedge(int u, int v) {
    edge[tot].v = v;
    edge[tot].next = head[u];
    head[u] = tot++;
}
 
void tarjan(int u) { //Tarjan算法求有向图的强连通分量
    dfn[u] = low[u] = ++cnt; //cnt为时间戳
    instack[u] = 1; //标记在栈中
    sstack[top++] = u; //入栈
    for(int e = head[u]; e != -1; e = edge[e].next) { //枚举u的每一条边
        int v = edge[e].v; //u所邻接的边
        if(!dfn[v]) { //未被访问
            tarjan(v); //继续向下找
            if(low[u] > low[v]) low[u] = low[v]; // 更新结点u所能到达的最小次数层
        }
        else if(instack[v] && low[u] > dfn[v]) { //如果v结点在栈内
            low[u] = dfn[v];
        }
    }
    if(low[u] == dfn[u]) { //如果节点u是强连通分量的根
        scnt++; //连通分量标号加1
        while(1) {
            int v = sstack[--top]; //退栈
            instack[v] = 0; //标记不在栈中
            belong[v] = scnt; //出栈结点t属于scnt标号的强连通分量
            if(v == u) break; //直到将u从栈中退出
        }
    }
}
 
int main()
{
    int T, kcase = 0;
    while(~scanf("%d %d", &n, &m) && (n || m)) {
        init();
        int a, b;
        while(m--) {
            scanf("%d %d", &a, &b);
            addedge(a, b);
        }
        for(int i = 1; i <= n; i++) { //枚举每个结点,搜索连通分量
            if(!dfn[i]) { //未被访问
                tarjan(i); //则找i结点的连通分量
            }
        }
        if(scnt == 1) {
            puts("Yes");
        }
        else {
            puts("No");
        }
    }
    return 0;
}

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值