【luogu4017】【拓扑排序】最大食物链计数

传送门

题目背景

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

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

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

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

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

输入格式

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

输出格式

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

输入输出样例

输入 #1

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

输出 #1

5

说明/提示

各测试点满足以下约定：

【补充说明】

数据中不会出现环，满足生物学的要求。（感谢 @AKEE ）

---

解题思路

做一个拓扑排序，设一个$s[]$来传递食物链
因为出度为$0$的点肯定是终点，那么最后只要累计$s[$所有出度为$0$的点$]$就好了

---

Code

#include<iostream>
#include<cstdio>
using namespace std;
const int mod=80112002;
struct DT{
	int to,next;
}a[510000];
int n,m,x,y,num,ans,ru[5100],chu[5100],s[5100],head[5100];
int q[5100],h,t;
void add(int x,int y){
	a[++num].to=y,a[num].next=head[x];
	head[x]=num;
}
void dfs(){
	while(h++<t){
    	for(int i=head[q[h]];i;i=a[i].next){
    		ru[a[i].to]--,s[a[i].to]=(s[q[h]]+s[a[i].to])%mod;//把食物链传递下去
    		if(!ru[a[i].to])
    		   q[++t]=a[i].to;
    	}
	}
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d",&y,&x);
		add(x,y);
		ru[y]++,chu[x]++;
	}
	for(int i=1;i<=n;i++)
	    if(!ru[i])
	    	q[++t]=i,s[i]=1;
	dfs();
	for(int i=1;i<=n;i++)
	    if(!chu[i])
	        ans=(s[i]+ans)%mod; //累计
	printf("%d",ans);
}