洛谷-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 ）

//应注意题目对最大食物链的定义：即起点是入度为零的点，终点是出度为零的点
#include<iostream>
#include<vector>
#include<queue>
using namespace std;
const int N = 5*1e5+10;
const int mod = 80112002;    //题目要求模数

vector<int> e[N];            //二维vector数组存图
int indegree[N];             //入度
int outdegree[N];            //出度
int n,m;
int dis[N];                  //食物链长度,dis[i]代表入度为0的点到第i个点的长度

void topsort(){              //拓扑排序
    queue<int>q;             //队列存入度为0的点
    
    for(int i = 1;i <= n;i++){    //搜索所有点
        if(indegree[i] == 0){     //找到入度为0的点
            q.push(i);            //压入队列
            dis[i] = 1;           //入度为0的点的距离设为1
        }
    }
    
    while(q.size()){              //当还有入度为0的点
        int t = q.front();        //取出队头
        q.pop();
        
        for(int i = 0;i < e[t].size();i++){    //从队头扩展点
            int j = e[t][i];                   //取出终点
            indegree[j]--;                     //删掉一个点就要让这个点的终点的入度减1
            dis[j]=(dis[t]+dis[j])%mod;        //加上上一个点的距离，记得mod模数，不然爆炸
            
            if(!indegree[j])q.push(j);         //如果有入读为零的点，那么就是下一次遍历的起点，压入队列
        }
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin >> n >> m;
    
    for(int i = 0;i < m;i++){        //存图
        int a,b;
        cin >> a >> b;
        e[a].push_back(b);
        indegree[b]++;               //记录入度
        outdegree[a]++;              //记录出度
    }
    
    int ans = 0;
    
    topsort();
    
    for(int i = 1;i <= n;i++){       //遍历所有点
        if(outdegree[i] == 0){       //如果这个点出度为0
            ans = (ans+dis[i])%mod;  //答案加上这个点所对应的食物链的条数，即以这个点为终点的最大食物链的数量
        }
    }
    
    printf("%d",ans);
    return 0;
}