【图论专题】拓扑排序

最新推荐文章于 2021-08-30 16:36:33 发布

繁凡さん

最新推荐文章于 2021-08-30 16:36:33 发布

阅读量429

点赞数 4

分类专栏： # 第三章图论图论 - 拓扑排序 AcWing课程/洛谷题单刷题总结文章标签：拓扑排序

本文链接：https://blog.csdn.net/weixin_45697774/article/details/107707818

版权

本文介绍了拓扑排序的概念及其在解决差分约束问题中的应用。通过拓扑排序，可以消除循环依赖，实现DP递推求解最长路。文章列举了AcWing平台上的几个实例，如家谱树、奖金问题、可达性统计和车站分级，详细阐述了如何利用拓扑排序优化算法，避免超时，并给出了建图优化技巧。

摘要由CSDN通过智能技术生成

拓扑排序
给定一张有向无环图，若一个序列A满足图中的任意一条边(x,y)x都在y的前面呢么序列A就是图的拓扑排序

实际上拓扑排序就是满足所有的边x指向y，x一定在y的前面。这样按照拓扑排序递推，就可以满足每一个状态都没有循环依赖 $- >$ 没有后效性 $- >$ 可以满足DP递推

题目	算法
AcWing 1191. 家谱树	toposort模板题
AcWing 1192. 奖金	差分约束+拓扑序 = 递推求最长路
AcWing 164. 可达性统计	拓扑图递推求答案
AcWing 456. 车站分级	拓扑图递推求解差分约束系统

A、AcWing 1191. 家谱树（toposort模板题）

在这里插入图片描述

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 110, M = 50007;

int ver[M], nex[M], edge[M], head[N], tot;
int n, m;
int d[N];//入度
bool vis[N];
int q[M];
//int ans[M];
int cnt;

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

void toposort(){
   
    int hh = 0, tt = -1;
    for(int i = 1;i <= n;++i)
        if(!d[i]){
   
            q[++ tt] = i;
            if(tt == M)tt = 0;
        }
    
    while(hh <= tt){
   
        int x = q[hh ++ ];
        
        //ans[++ cnt] = x;
        
        if(hh == M )hh = 0;
        for(int i = head[x];~i;i = nex[i]){
   
            int y = ver[i];
            if(-- d[y] == 0)
            q[ ++ tt] = y;
        }
    }
}

int main(){
   
    scanf("%d",&n);
    memset(head, -1, sizeof head);
    
    for(int i = 1;i <= n;++i)
    {
   
        int son;
        while(~scanf("%d",&son) && son){
   
            add(i, son);
        }
    }
    
    toposort();
    
    for(int i = 0;i < n;++i)
    printf("%d ",q[i]);
}

B、AcWing 1192. 奖金（差分约束+拓扑序 = 递推求最长路）

在这里插入图片描述
实际上是一道差分约束的题，要求最小值，也就是取所有下界的最大值，用最长路

如果用差分约束spfa的时间复杂度为 $O (k m)$ ~ $O (n m)$ ，可能会被卡

因为本题中所有的权值都是正值，因此我们可以不用差分约束就求这个最长路。
我们直接拓扑排序得到一张有向无环图DAG $= >$ 每一个状态都没有循环依赖 $= >$ 没有后效性 $= >$ 就可以用DP递推求最长路。

然后别忘了条件 a >= b + 1 等价于b向a连一条权值为1的边。

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 50007, M = 500007;

int n, m;
int dist[N];//longest way
int ver[M], edge[M], nex