题目背景
你知道食物链吗?
D
e
l
i
a
Delia
Delia 生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条。于是她来就来求助你,然而你也不会啊!写一个程序来帮帮她吧。
题目描述
给你一个食物网,你要求出这个食物网中最大食物链的数量。
(这里的“最大食物链”,指的是生物学意义上的食物链,即最左端是不会捕食其他生物的生产者,最右端是不会被其他生物捕食的消费者。)
D
e
l
i
a
Delia
Delia 非常急,所以你只有
1
1
1 秒的时间。
由于这个结果可能过大,你只需要输出总数模上
80112002
80112002
80112002 的结果。
输入格式
第一行,两个正整数
n
、
m
n、m
n、m,表示生物种类
n
n
n 和吃与被吃的关系数 mm。
接下来
m
m
m 行,每行两个正整数,表示被吃的生物
A
A
A和吃
A
A
A的生物
B
B
B。
输出格式
一行一个整数,为最大食物链数量模上
80112002
80112002
80112002 的结果。
输入输出样例
输入
5 7
1 2
1 3
2 3
3 5
2 5
4 5
3 4
输出
5
说明/提示
各测试点满足以下约定:
【补充说明】
数据中不会出现环,满足生物学的要求。(感谢 @AKEE )
解题思路
做一个拓扑排序,设一个
f
[
]
f[]
f[]来传递食物链
因为出度为
0
0
0的点肯定是终点,那么最后只要累计所有出度为
0
0
0的点的
f
[
]
f[]
f[]就好了
代码
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int mod=80112002;
int n,m,h,head[5100],t,sum,ans,xx,yy,k,ru[5100],chu[5100],f[5100],q[5100];
struct c{
int x,next;
}a[501000];
void add(int x,int y)
{
++k;
a[k].x=y;
a[k].next=head[x];
head[x]=k;
}
void bfs(){
while(h<t){
h++;
for(int i=head[q[h]];i;i=a[i].next)
{
ru[a[i].x]--;
f[a[i].x]=(f[q[h]]+f[a[i].x])%mod;//把食物链传递下去
if(!ru[a[i].x])
q[++t]=a[i].x;
}
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d",&yy,&xx);
add(xx,yy);
ru[yy]++;
chu[xx]++;
}
for(int i=1;i<=n;i++)
if(!ru[i])
{
q[++t]=i;
f[i]=1;
}
bfs();
for(int i=1;i<=n;i++)
if(!chu[i])
ans=(f[i]+ans)%mod;//累加
printf("%d",ans);
}