求最大食物链

#include<stdio.h>
#include<iostream>
#include<math.h>
#include<algorithm>
#include<queue>
#include<vector>
#define MOD 80112002
using namespace std;
int n,m,ans;
vector<int> p[5005];
queue<int> q;
int f[500005],ind[500005],outd[500005];//记录结果、入度、出度
int main(){
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int x,y;
        cin>>x>>y;
        outd[x]++; ind[y]++; //出度入度改变
        p[x].push_back(y); //邻接表存储图
    }
    memset(f,0,sizeof(f));
    for(int i=1;i<=n;i++){  //首先将所有入度为0的加入队列，代表食物链底层
        if(ind[i] == 0){
            q.push(i);
            f[i]=1;
        }
    }
    while(!q.empty()){
        int x = q.front();
        for(int i=0,sz=p[x].size();i<sz;i++){
            int y = p[x][i];
            f[y] = (f[x] + f[y]) % MOD;
            ind[y]--;
            if(ind[y]==0)  //如果入度为0了也加入队列
                q.push(y);
        }
        q.pop();
    }
    for(int i=1;i<=n;i++) //统计所有出度为0的点即为结果
        if(outd[i] == 0) ans = (ans+f[i]) % MOD;
    cout<<ans<<endl;
    return 0;
}