【算法模板】拓扑排序

可以使用邻接矩阵或邻接表

步骤：

• 初始时，先第一遍扫描入度为零的点，存入队列中

• 以该点为父点，通过邻接矩阵遍历他的所有子点

• 子点的入度减一，并对该点进行处理（如求最长路，可以通过父点更新他到起点的长度，或如求从起点到该点的路数，即可以累加上父点的路数）

• 当子点的入度为零时，将他加入到队列中

• 循环直到队列为空

模板题1：P4017 最大食物链计数

#include<bits/stdc++.h>
using namespace std;
const int N=5000005,M=80112002;
queue<int> q;
int n,m,head[N],cnt,d[N],f[N],D[N];
struct Edge{
	int nxt,to;
}g[N*2];
void add(int from,int to){
	g[++cnt].nxt=head[from];
	g[cnt].to=to;
	head[from]=cnt;
}
void init(){
	cin>>n>>m;
	int x,y;
	for(int i=1;i<=m;i++){
		cin>>x>>y; add(x,y); d[y]++; D[x]++;
	}
}
void work(){
	for(int i=1;i<=n;i++){
		if(d[i]==0){
			q.push(i); f[i]=1;
		}
	}
	int ans=0;
	while(!q.empty()){
		int x=q.front(); q.pop();
		for(int i=head[x];i;i=g[i].nxt){
			int v=g[i].to;
			d[v]--; f[v]=(f[v]+f[x])%M;
			if(d[v]==0) q.push(v);
		}
	}
	for(int i=1;i<=n;i++){
		if(!D[i]) ans=(ans+f[i])%M;
	}
	cout<<ans<<endl;
}
int main(){
	init();
	work();
	return 0;	
}