洛谷P4017 最大食物链计数

题目链接：

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

提示

各测试点满足以下约定：

思路：

标签：拓扑排序+dp。

转移方程：
以生产者为食物链结尾的食物链条数为1
以 某消费者 为 食物链结尾 的食物链条数 = 他所能捕食的猎物为结尾的食物链数 之和。

最终答案为 终极消费者为食物链结尾 的 食物链条数 之和。

注意：在计算a为食物链结尾的食物链条数时，需要确保a的猎物的食物链条数都已计算完毕，所以要用到拓扑排序，用deg_in记录每个动物的入度，只有入度为0时，将他入队，计算食物链数目。

AC代码

#include<iostream>
#include<queue>
using namespace std;

const int MAXN = 5005;
const int MAXM = 500005;
const int MOD = 80112002;

bool food[MAXN];
long long dp[MAXN]={0};
int deg_in[MAXN]={0};
bool map[MAXN][MAXN];
queue<int>q;
long long ans=0;

int main()
{
    int n,m;
    int a,b,head;
    cin>>n>>m;

    for(int i=1;i<=m;++i)
    {
        cin>>a>>b;
        map[a][b]=true;
        food[a]=true;
        deg_in[b]++;
    }
    for(int i=1;i<=n;++i)
    {
        if(deg_in[i]==0) 
        {
            dp[i]=1;
            q.push(i);
        }
    }

    while(!q.empty())
    {
        head=q.front();
        q.pop();
        for(int i=1;i<=n;++i)
        {
            if(map[head][i])
            {
                deg_in[i]--;
                dp[i]=(dp[i]+dp[head])%MOD;
                if(deg_in[i]==0) q.push(i);
            }
        }
    }

    for(int i=1;i<=n;++i)
    {
        if(food[i]==false) ans=(dp[i]+ans)%MOD;
    }

    cout<<ans<<endl;

    return 0;
}