• Tarjan总结


    在图论中,连通图基于连通的概念。在一个无向图G中,若从顶点到顶点有路径相连(当然从到也一定有路径),则称和是连通的。如果G是有向图,那么连接和的路径中所有的边都必须同向。如果图中任意两点都是连通的,那么图被称作连通图。图的连通性是图的基本性质。
    有向图的所有的有向边替换为无向边,所得到的图称为原图的基图。如果一个有向图的基图是连通图,则有向图是弱连通图,弱有向图的任意两点都可以相互到达那么称这个有向图为强连通图,如果一个有向图的子图是强连通图,那么这个子图称为该有向图的强连通分量。有向图中一个单个的点也是一个强连通分量
    Tarjan在无向图里可以求割点和割桥。在有向图里可以求强连通分量。
     
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    #include<cstdio>
    #include<cstring>
    #include<iostream>
    #include<cstring>
    #include<vector>
    #include<stack>
    #include<algorithm>
    using namespace std;
     
    #define N 10003
    int dfn[N],low[N],ins[N],Time,num;//ins是否在栈里
    vector<int>gra[N];
    stack<int>sta;
     
    void Tarjan(int s)
    {
        dfn[s] = low[s] = ++Time;
        sta.push(s);
        ins[s] = 1;
        for(int i=0;i<gra[s].size();i++)
        {
            int k = gra[s][i];
            if(dfn[k] == 0){
                Tarjan(k);
                low[s] = min(low[s] ,low[k]);
            }
            if(dfn[k] != 0 && ins[k] == 1){
                low[s] = min(low[s] ,dfn[k]);//low[s] = min(low[s] ,low[k]);好像也是对的
            }
        }
        if(dfn[s] == low[s])
        {
            num++;
            while(!sta.empty())//一次性弹出来的所有点属于一个强连通分量
            {
                int temp = sta.top();
                sta.pop();
                ins[temp] = 0;
                if(temp == s) break;
            }
        }
    }
     
    int main()
    {
        int n,m;
        while(scanf("%d%d",&n,&m)&&(n + m))
        {
            memset(dfn,0,sizeof(dfn));
            memset(ins,0,sizeof(ins));
            memset(low,0,sizeof(low));
            while(!sta.empty()) sta.pop();
            for(int i=1;i<=n;i++) gra[i].clear();
            int a,b;
            while(m--)
            {
                scanf("%d%d",&a,&b);
                gra[a].push_back(b);
            }
            num = Time = 0;
            for(int i=1;i<=n;i++)
            {
                if(dfn[i]==0) Tarjan(i);
            }
    //num就是强连通分量的个数
        }
    }
     
    1.HDU1269 
    裸Tarjan求强连通分量个数,跑N遍Tarjan得到的NUM就是结果
     
    2.POJ2186 
    Tarjan后重建图判断出度为0的节点个数,若为1则OK。
       for(int i=1;i<=n;i++)
        {
            for(int j=0;j<Gra[i].size();j++){
                int k = Gra[i][j];
                if(belong[i] != belong[k]){
                    outDegree[belong[i]]++;
                }
            }
        }
    其实也算不上重建图....只是求了一下各个节点的出度而已......
     
    3.POJ2762 
        Tarjan + 拓扑排序
     这个题是真的烦人,都给我WA哭了有木有,从中午2点错到晚上10点,中间气的我玩了几把游戏,真的气。
     还有一个就是有人说缩晚点以后的DAG是一条链,也有好多人是这么写的,都AC了,说是dfs点的数目等于强连通块的数目就可以。
    但是感觉明显不对啊,整个就是一个树,你dfs的点数肯定等于强连通块的数目啊,不懂他们怎么A的。
    感觉和2186差不多,但是最多只有1000个顶点。
     我怎么这么愚蠢,居然想到要暴力!!!我太蠢啦!!!
    学长一语点醒

     
    我这都没想到-。-||
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    #include<cstdio>
    #include<cstring>
    #include<iostream>
    #include<string>
    #include<vector>
    #include<stack>
    #include<set>
    #include<algorithm>
    using namespace std;
    #define N 1002
    vector<int>Gra[N];
    stack<int>Sta;
    int map[N][N];
    int dfn[N],low[N],inStack[N],belong[N],Time,cnt;
    int inDegree[N];
     
    void init()
    {
        Time = cnt = 0;
        memset(dfn,0,sizeof(dfn));
        memset(low,0,sizeof(dfn));
        memset(inStack,0,sizeof(inStack));
        memset(inDegree,0,sizeof(inDegree));
        memset(belong,0,sizeof(belong));
        for(int i=0;i<N;i++) Gra[i].clear();
        memset(map,0,sizeof(map));
        while(!Sta.empty()) Sta.pop();
    }
     
    void Tarjan(int s)
    {
        dfn[s] = low[s] = ++Time;
        inStack[s] = 1;
        Sta.push(s);
        for(int i=0;i<Gra[s].size();i++)
        {
            int j = Gra[s][i];
            if(dfn[j] == 0){
                Tarjan(j);
                low[s] = min(low[s], low[j]);
            }
            else if(inStack[j] == 1){
                low[s] = min(low[s], dfn[j]);
            }
        }
        if(dfn[s] == low[s])
        {
            cnt ++;
            while(!Sta.empty()){
                int temp = Sta.top(); Sta.pop();
                inStack[temp] = 0;
                belong[temp] = cnt;
                if(temp == s) break;
            }
        }
        return;
    }
     
    void tsort()
    {
        for(int k=0;k<cnt;k++){
            int fuck = 0,pos;
            for(int i=1;i<=cnt;i++)
            {
                if(inDegree[i] == 0)
                {
                    fuck ++;
                    pos = i;
                }
            }
            if(fuck > 1){
                printf("No ");
                return ;
            }
            inDegree[pos ] = -1;
            for(int i=1;i<=cnt;i++)
            {
                if(map[pos][i] == 1)
                    inDegree[i]--;
            }
        }
        printf("Yes ");
    }
     
    int main()
    {
        int noc;
        cin>>noc;
        while(noc--)
        {
            init();
            int n,m,x,y;
            scanf("%d%d",&n,&m);
            for(int i=0;i<m;i++)
            {
                scanf("%d%d",&x,&y);
                Gra[x].push_back(y);
            }
            for(int i=1;i<=n;i++) if(dfn[i] == 0) Tarjan(i);
            if(cnt == 1) {
                printf("Yes ");
                continue;
            }
            for(int i=1;i<=n;i++)
            {
                for(int j=0;j<Gra[i].size();j++)
                {
                    int k = Gra[i][j];
                    if(belong[i]!=belong[k]){
                        if(map[belong[i]][belong[k]] == 0){
                            map[belong[i]][belong[k]] = 1;
                            inDegree[belong[k]]++;
                        }
                    }
                }
            }
            for(int i=1;i<=cnt;i++) printf("%d %d ",i,inDegree[i]);
            tsort();
        }
    }

     
    4.HYSBZ 1179
        这个就是tarjan + bfs就可以了,我感觉我的代码命名都很清晰,一看就明白了。
        还有就是网上好多人都是 tarjan + spfa,但我感觉spfa没有必要,因为缩点以后的图一定是一个DAG图。他的搜索树一定是一个没有向后边的的树,所以一遍vfs     就够了。
     
    5.最后说一下就是,这个缩点的时候,如果只用入度出度是否为0就能得到答案的话直接
    for(int i=1;i<=n;i++)
        {
            for(int j=0;j<Gra[i].size();j++){
                int k = Gra[i][j];
                if(belong[i] != belong[k]){
                    outDegree[belong[i]]++;//inDegree[belong[i]]++;
                }
            }
        }
    就可以了。
    但是假如要对入度出度进行操作就不可以了,因为强连通块1是由1,2,3节点构成,强连通块2是由4,5构成
    现在缩点后块1 ----> 块2。其中节点2 ---> 4,2 ---> 5,3 ---> 5。这样下来outdegree[1] == 3;
    indegree[2] == 3;这其实就错了,如果拓扑排序是不出结果的。比如上面的poj2762.
     
    那么这个类型的就告一段落了,感觉掌握的还不错,因为本身Tarjan比较好理解(当然一开始看的时候各种网上找讲解
    也是看得我一脸懵逼,一开始看还是有点不好想的)但是只要理解了以后代码,就很好打了。
  • 相关阅读:
    Codeforces Round #336 B
    Codeforces Round #336 A
    hiho一下157
    Codeforces Round #420 E
    Codeforces Round #420 C
    Codeforces Round #420 B
    Codeforces Round #420 A
    Codeforces Round #418 C
    Codeforces Round #418 B
    CodeForces 811D Vladik and Favorite Game bfs,模拟
  • 原文地址:https://www.cnblogs.com/liwenchi/p/7259276.html
Copyright © 2020-2023  润新知