题目链接:
题目背景
你知道食物链吗?Delia 生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条。于是她来就来求助你,然而你也不会啊!写一个程序来帮帮她吧。
题目描述
给你一个食物网,你要求出这个食物网中最大食物链的数量。
(这里的“最大食物链”,指的是生物学意义上的食物链,即最左端是不会捕食其他生物的生产者,最右端是不会被其他生物捕食的消费者。)
Delia 非常急,所以你只有 1 1 1 秒的时间。
由于这个结果可能过大,你只需要输出总数模上 80112002 80112002 80112002 的结果。
输入格式
第一行,两个正整数 n 、 m n、m n、m,表示生物种类 n n n 和吃与被吃的关系数 m m m。
接下来 m m m 行,每行两个正整数,表示被吃的生物A和吃A的生物B。
输出格式
一行一个整数,为最大食物链数量模上 80112002 80112002 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;
}