洛谷 P4017 最大食物链计数

最大食物链计数

题目背景

你知道食物链吗？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

提示

各测试点满足以下约定：

【补充说明】

数据中不会出现环，满足生物学的要求。（感谢 @AKEE ）
思路借鉴：https://www.luogu.com.cn/article/nithzyjc

#include <bits/stdc++.h>
using namespace std;
#define MAX 5005
#define p 80112002 //余数 
vector<int> prey[MAX];//记录每个生物对应的猎物 
int n,m,ru[MAX]={0},chu[MAX]={0};//记录每个点的入度、出度 
queue<int> q;//队列，用于拓扑排序
int f[MAX] = {0}; //到达该点的路径数 
int main()
{
	cin>>n>>m;
	for(int i=0;i<m;i++)
	{
		int a,b;
		cin>>a>>b;
		prey[b].push_back(a);
		ru[a]++;
		chu[b]++;
	}
	for(int i=1;i<=n;i++)//将入度为0的点入队   注意，动物序号从1开始 
	{
//		f[i] = 1;//初始化为1  应该是只有入度为0的点才能初始化为1 
		if(ru[i] == 0)
		{
			f[i] = 1; 
			q.push(i);
		}
	}
	int ans=0;//路径总数 
	while(!q.empty())
	{
		int now = q.front();
		q.pop();
		
		vector<int>::iterator it = prey[now].begin();
		for(;it!=prey[now].end();it++)
		{
			f[*it] +=f[now];//更新路径个数 
			f[*it] %= p;
			ru[*it] --;
			if(ru[*it] == 0)//对于入度为0的点入队
			{
				if(chu[*it] == 0) //如果该点是食物链末端，则计算 
				{
					ans+=f[*it];
					ans%=p;//注意多处优化 
					continue;//不做后续操作 
				}
				q.push(*it); 
			} 
		 } 
	 }
	cout<<ans;
	return 0;
 }