题目背景
你知道食物链吗?Delia 生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条。于是她来就来求助你,然而你也不会啊!写一个程序来帮帮她吧。
题目描述
给你一个食物网,你要求出这个食物网中最大食物链的数量。
(这里的“最大食物链”,指的是生物学意义上的食物链,即最左端是不会捕食其他生物的生产者,最右端是不会被其他生物捕食的消费者。)
Delia 非常急,所以你只有 111 秒的时间。
由于这个结果可能过大,你只需要输出总数模上 801120028011200280112002 的结果。
输入格式
第一行,两个正整数 n、mn、mn、m,表示生物种类 nnn 和吃与被吃的关系数 mmm。
接下来 mmm 行,每行两个正整数,表示被吃的生物A和吃A的生物B。
输出格式
一行一个整数,为最大食物链数量模上 801120028011200280112002 的结果。
样例 #1
样例输入 #1
5 71 21 32 33 52 54 53 4
样例输出 #1
5
提示
各测试点满足以下约定:
【补充说明】
数据中不会出现环,满足生物学的要求。(感谢 @AKEE )
代码:
#include<iostream>#include<queue>#include<cstring>using namespace std;const int N = 1e6 + 10;int e[N], ne[N], h[N], idx = 0;int din[N],dout[N];int f[N];const int mod = 80112002;void add(int a, int b){e[idx] = b;ne[idx] = h[a];h[a] = idx++;}int n, m;bool st[N];int seq[N];void top(){queue<int>q;int cnt = 0;for (int i = 1; i <= n; i++){if (din[i] == 0){f[i] = 1;q.push(i);}if (dout[i] == 0){seq[cnt++] = i;}}while (!q.empty()){int t = q.front();q.pop();if (st[t]){continue;}st[t] = true;for (int i = h[t]; i != -1; i = ne[i]){int j = e[i];f[j] = (f[t] + f[j]) % mod;din[j]--;if (din[j] == 0){q.push(j);}}}int res = 0;for (int i = 0; i < cnt; i++){res += f[seq[i]];res %= mod;}cout << res << endl;}int main(){memset(h, -1, sizeof h);cin >> n >> m;for (int i = 0; i < m; i++){int a, b;cin >> a >> b;add(a, b);din[b]++;dout[a]++;}top();}