拓扑排序介绍及其应用

拓扑排序简介

拓扑排序是一种图论排序方式，要了解拓扑排序，首先要了解两个图论概念：入度和出度。在有向无环图（$DAG$）中，一个点是多少弧的弧头称为入度，是多少弧的弧尾是出度。
而作拓扑排序时，只需做以下步骤：

1. 将入度为0的点存在一个队列内
2. 每次取出一个点，将其能达到的点的入度减$1$
3. 将上一步操作中入度降为$0$的点放在队尾
4. 删除第$2$步取出的那个点

最后，被删除的点的先后顺序就是这个图的拓扑序。
例如下面这个有向无环图的拓扑序便是：$4\to 3\to 6\to 5\to 1\to 2$

下面附上一般拓扑排序的部分代码

void topo(){//拓扑排序 
	for(int i=1;i<=n;i++)
	{
		if(in[i]==0)
		{
			q.push(i);
		}
	}
	while(!q.empty())
	{
		int p=q.front(); q.pop();
		for(int i=head[p];i>0;i=nxt[i])
		{
			int go=to[i];
			in[go]--;
			if(in[go]==0)
			q.push(go);
		}
	}
}

其中$in[i]$用来存储入度，$go[i]$用来存储出度。

应用

我们看这样一道例题：

最大食物链计数

题目背景

你知道食物链吗？Delia 生物考试的时候，数食物链条数的题目全都错了，因为她总是重复数了几条或漏掉了几条。于是她来就来求助你，然而你也不会啊！写一个程序来帮帮她吧。

题目描述

给你一个食物网，你要求出这个食物网中最大食物链的数量。

（这里的“最大食物链”，指的是生物学意义上的食物链，即最左端是不会捕食其他生物的生产者，最右端是不会被其他生物捕食的消费者。）

Delia 非常急，所以你只有 $1$ 秒的时间。

由于这个结果可能过大，你只需要输出总数模上 $80112002$ 的结果。

输入格式

第一行，两个正整数 $n、m$，表示生物种类 $n$ 和吃与被吃的关系数 $m$。

接下来 $m$ 行，每行两个正整数，表示被吃的生物A和吃A的生物B。

输出格式

一行一个整数，为最大食物链数量模上 $80112002$ 的结果。

样例 #1

样例输入 #1

5 7
1 2
1 3
2 3
3 5
2 5
4 5
3 4

样例输出 #1

5

提示

各测试点满足以下约定：

首先不难看出食物网是一张有向无环图（$DAG$),可以考虑通过$dp$来求解，我们由相对若的生物指向相对强的生物来建边，如何建边可以考虑用邻接矩阵或者是我之前文章中提到的链式前向星算法引用链接
而状态转移方程并不难想，用数组$dp[i]$存储以$i$为终点的食物链条数，如果以$i$为弧尾的弧的弧头是$j$,那么我们有（在本题中注意取模）
$$dp[j]+=dp[i]$$
而最终的答案$ans=\sum dp[i](whilego[i]=0）$
最后附上本题解答

#include<bits/stdc++.h>
using namespace std;

queue<int> q;

int cnt=0,n,m;
int nxt[500005],to[500005],head[5005],dp[5005],in[5005];
bool has_son[5005]; 
const int mod=80112002;
long long ans=0;

void add(int u,int v){//链式前向星 
	nxt[++cnt]=head[u];
	head[u]=cnt;
	to[cnt]=v;
	has_son[u]=true;
}

void topo(){//拓扑排序 
	for(int i=1;i<=n;i++)
	{
		if(in[i]==0)
		{
			q.push(i);
			dp[i]=1;
		}
	}
	while(!q.empty())
	{
		int p=q.front(); q.pop();
		for(int i=head[p];i>0;i=nxt[i])
		{
			int go=to[i];
			dp[go]=(dp[go]+dp[p])%mod;
			in[go]--;
			if(in[go]==0)
			q.push(go);
		}
	}
}

int read(){
	int x=0,f=1;
    char c=getchar();
    if(c<'0'||c>'9'){
    	if(c=='-') f=-1;
    	c=getchar();
	}
	while(c>='0'&&c<='9'){
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}


int main(){
      n=read(),m=read(); 
	for(int i=1;i<=m;i++){
		int u=read(),v=read();
		add(u,v);
		in[v]++;
	}
	topo();
	for(int i=1;i<=n;i++){
		if(!has_son[i]){
			ans=(ans+dp[i])%mod;
		}
	}
	cout<<ans;
	return 0;
}