最大食物链计数——拓扑排序

你知道食物链吗？Delia 生物考试的时候，数食物链条数的题目全都错了，因为她总是重复数了几条或漏掉了几条。于是她来就来求助你，然而你也不会啊！写一个程序来帮帮她吧。
给你一个食物网，你要求出这个食物网中最大食物链的数量。（这里的“最大食物链”，指的是生物学意义上的食物链，即最左端是不会捕食其他生物的生产者，最右端是不会被其他生物捕食的消费者。）
Delia 非常急，所以你只有 1 秒的时间。
由于这个结果可能过大，你只需要输出总数模上 80112002 的结果。数据中不会出现环，满足生物学的要求。

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

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

样例输入 
5 7
1 2
1 3
2 3
3 5
2 5
4 5
3 4
样例输出
5

解析：

对于通过最佳捕食者的路径，取决于通过图中的前节点的路径的数量，将通过这些节点的路径的数量相加就是通过最佳捕食者的路径数量；以此类推一直到通过最佳生产者的路径，为1.

#include <bits/stdc++.h>
using namespace std;
#define int long long 
#define ios ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
typedef pair<int,int> PII;
const int N=1e6+10;
vector <int> g[N];
int n,m;
int d[N];   //记录每个点的入度，入度为零的点为 “最佳生产者”
int s[N];   //记录每个点的出度，出度为零的点为 “最佳捕食者”
int w[N];
queue <int> q;
int mod=80112002;
signed main()
{
    ios;
    cin>>n>>m;
    while (m--)
    {
        int a,b;
        cin>>a>>b;
        g[a].push_back(b);
        d[b]++;
        s[a]++;
    }
    for (int i=1;i<=n;i++) 
    {
        if (!d[i])
        {
            q.push(i);
            w[i]=1;
        }
    }
    while (q.size())
    {
        int t=q.front();
        q.pop();
        for (auto x:g[t])
        {
            w[x]=(w[x]+w[t])%mod;
            d[x]--;
            if (!d[x]) q.push(x);
        }
    }
    int ans=0;
    for (int i=1;i<=n;i++) 
    {
        if (!s[i]) ans=(ans+w[i])%mod;
    }
    cout<<ans;
    return 0;
}