P4017 最大食物链计数

题目背景

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

题目描述

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

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

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

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

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

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

输出格式

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

思路

生物试卷好评
一道非初中生物知识水平做不出的破题
topsort直接过，没啥好说的
code：

#include<iostream>
#include<cstdio>
#include<map>
using namespace std;
int n,m,head[5001],to[500001],nxt[500001],tot=1;
int head2[5001],a[5001],s=0;
void add(int x,int y)
{
	nxt[tot]=head[x];
	to[tot]=y;
	head[x]=tot;
	head2[y]=1;
	tot++;
}
void dfs(int x)
{
	if (a[x]) return;
	if (head[x]==0) a[x]=1;
	for (int i=head[x];i;i=nxt[i])
	{
		dfs(to[i]);//记忆化搜索完成topsort
		a[x]=(a[x]+a[to[i]])%80112002;
	}
	return;
}
int main()
{
	cin>>n>>m;
	for (int i=1;i<=m;i++)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		add(y,x);
	}
	for (int i=1;i<=n;i++)
	{
		dfs(i);
		if (head2[i]==0) s=(s+a[i])%80112002;
	}
	cout<<s;
    return 0;
}