拓扑排序
灵魂画图
文字
类似于宽度优先搜索的思想,区别在于,当从u找到v时,并不马上把v放入队列,而是找到v的时候就需要把v的“入度” 减去一个1,当入度为零时则入队并输出。
例如②中我们从2开始搜索(从任意一个入度为0的节点开始即可),找到了4,但是此时4应该还需要3完成之后才能点亮4,所以4不入队,但他的入度会减一,从原本的2变成了1,然后再从1开始搜索,搜索到3,点亮,搜索到4,点亮(因为此时的入度已经为零了),最后搜索到5。
对于一个有向无环图(DAG)而言,一定存在至少一个拓扑序,并且入度为零的节点被称为起点,出度为零的节点被称为终点。
题目 模板
描述
给定一个n个点m条边的有向图,点的编号是1到n,图中可能存在重边和自环。
请输出任意一个该有向图的拓扑序列,如果拓扑序列不存在,则输出-1。
若一个由图中所有点构成的序列A满足:对于图中的每条边(x, y),x在A中都出现在y之前,则称A是该图的一个拓扑序列。
输入格式
第一行包含两个整数n和m
接下来m行,每行包含两个整数x和y,表示存在一条从点x到点y的有向边(x, y)。
输出格式
共一行,如果存在拓扑序列,则输出拓扑序列。
否则输出-1。
数据范围
1≤n,m≤105
输入样例:
3 3
1 2
2 3
1 3
输出样例:
1 2 3
参考代码
#include <bits/stdc++.h>
using namespace std;
const int N = 100010;
int e[N], h[N], ne[N], idx;
int din[N], st[N];
void add(int a, int b) {
e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}
int n, m;
vector<int> top_sort() {
vector<int> res;
queue<int> q;
for (int i = 1; i <= n; i++) {
if (din[i] == 0) {
// 入度为0 加入队列
q.push(i);
st[i] = 1;
}
}
while (q.size()) {
auto d = q.front();
res.push_back(d);
q.pop();
for (int i = h[d]; ~i; i = ne[i]) {
int j = e[i];
if (!st[j]) {
din[j]--;
if (din[j] == 0) {
st[j] = 1;
q.push(j);
}
}
}
}
if (res.size() == n) return res;
res.resize(1);
res[0] = -1;
return res;
}
int main() {
memset(h, -1, sizeof h);
ios::sync_with_stdio(false);
cin >> n >> m;
while (m--) {
int a, b;
cin >> a >> b;
add(a, b);
din[b] ++;
}
auto re = top_sort();
for (auto x : re) {
cout << x << " ";
}
return 0;
}
注意
需要注意的是,如果图中成环则无拓扑序。
当然,成环可以看作一个强连通分量,在进行缩点之后也可以跑拓扑排序。