拓扑排序

一、什么是拓扑排序

在图论中，拓扑排序（Topological Sorting）是一个有向无环图（DAG, Directed Acyclic Graph）的所有顶点的线性序列。且该序列必须满足下面两个条件：

1. 每个顶点出现且只出现一次。
2. 若存在一条从顶点 A 到顶点 B 的路径，那么在序列中顶点 A 出现在顶点 B 的前面。

有向无环图（DAG）才有拓扑排序，非DAG图没有拓扑排序一说。

例如，下面这个图：

 
它是一个 DAG 图，那么如何写出它的拓扑排序呢？这里说一种比较常用的方法：

1. 从 DAG 图中选择一个 没有前驱（即入度为0）的顶点并输出。
2. 从图中删除该顶点和所有以它为起点的有向边。
3. 重复 1 和 2 直到当前的 DAG 图为空或当前图中不存在无前驱的顶点为止。后一种情况说明有向图中必然存在环

于是，得到拓扑排序后的结果是 { 1, 2, 4, 3, 5 }。

通常，一个有向无环图可以有一个或多个拓扑排序序列。

二、拓扑排序的应用

拓扑排序通常用来“排序”具有依赖关系的任务。

比如，如果用一个DAG图来表示一个工程，其中每个顶点表示工程中的一个任务，用有向边<A,B><A,B>表示在做任务 B 之前必须先完成任务 A。故在这个工程中，任意两个任务要么具有确定的先后关系，要么是没有关系，绝对不存在互相矛盾的关系（即环路）。

三、拓扑排序的实现

根据上面讲的方法，我们关键是要维护一个入度为0的顶点的集合。
下面我们来看一下C++代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int MAX_N = 100010;
vector <int> vec[MAX_N];
int top[MAX_N],in[MAX_N];
int n, m,cnt = 0;;
void topsort(){
	queue<int> q;
	for(int i = 1; i <= n; i++){
		if(!in[i]){
			q.push(i);
		}
	}
	while(!q.empty()){
		int now = q.front();
		q.pop();
		top[++cnt] = now; 
		int sz = vec[now].size();
		for(int i = 0; i < sz; i++){
			int u = vec[now][i];
			in[u]--;
			if(!in[u]) q.push(u);
		}
	}
}
int main(){
	cin >> n >> m;
	for(int i = 0; i < m; i++){
		int a,b;
		scanf("%d%d",&a,&b);
		vec[a].push_back(b);
		in[b]++;
	}
	topsort();
	for(int i = 1; i <= cnt; i++){
		cout << top[i] << " ";
	}
	return 0;
} 

拓扑排序还可以用于判有向图是否有环，我们只需要在上面的代码上判断下：cnt < n，如果成立，说明有环，否则不成环。

下面来看一道例题，更深入的理解下拓扑序的应用：

旅行计划

题目要求出以一个城市为终点最多游览城市个数，而且使得线路上除了第一个城市，每个城市都在路线前一个城市东面，

根据这句，我们知道这里有一个前驱后继关系，所以可以根据关系建有向图，因为拓扑序表示的就是先后关系，所以不难想到用拓扑序求到每个点的最远距离（因为整个图都遍历后每个点的所有前驱一定都找到了），我们只需要按拓扑序中点的关系计算下距离（每个点的所有前驱个数+1即是距离）。

下面看下代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int MAX_N = 100010;
vector <int> vec[MAX_N];
int dis[MAX_N],in[MAX_N];
int n, m;
void topsort(){
	queue<int> q;
	for(int i = 1; i <= n; i++){
		if(!in[i]){
			q.push(i);
			dis[i] = 1;//第一个点没有其他点可以到达 
		}
	}
	while(!q.empty()){
		int now = q.front();
		q.pop();
		int sz = vec[now].size();
		for(int i = 0; i < sz; i++){
			int u = vec[now][i];
			in[u]--;
			dis[u] = dis[now] + 1;//递推式，now的后继u由now更新 
			if(!in[u]) q.push(u);
		}
	}
}
int main(){
	cin >> n >> m;
	for(int i = 0; i < m; i++){
		int a,b;
		scanf("%d%d",&a,&b);
		vec[a].push_back(b);
		in[b]++;
	}
	topsort();
	for(int i = 1; i <= n; i++){
		printf("%d\n",dis[i]);
	}
	return 0;
}