P4017 最大食物链计数
入口
题目背景
你知道食物链吗?Delia 生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条。于是她来就来求助你,然而你也不会啊!写一个程序来帮帮她吧。
题目描述
给你一个食物网,你要求出这个食物网中最大食物链的数量。
(这里的“最大食物链”,指的是生物学意义上的食物链,即最左端是不会捕食其他生物的生产者,最右端是不会被其他生物捕食的消费者。)
Delia 非常急,所以你只有 11 秒的时间。
由于这个结果可能过大,你只需要输出总数模上 80112002的结果。
输入格式
第一行,两个正整数 n、mn、m,表示生物种类 nn 和吃与被吃的关系数 mm。
接下来 mm 行,每行两个正整数,表示被吃的生物A和吃A的生物B。
输出格式
一行一个整数,为最大食物链数量模上 8011200280112002 的结果。
输入
5 7
1 2
1 3
2 3
3 5
2 5
4 5
3 4
输出
5
n<=5000,m<=500000;
这题就是拓扑排序的变形题吧,入度为0就把点入队,遍历这个点的所有出边,当出度为0时就保存答案。
#include <bits/stdc++.h>
using namespace std;
const int mod = 80112002;
typedef long long ll;
int in[5005], out[5005];
vector<int> adj[1005];
int cnt[5005];
int main()
{
int n, m, a, b;
cin >> n >> m;
for (int i = 0; i < m; i++)
{
cin >> a >> b;
adj[a].push_back(b);
out[a]++;
in[b]++;
}
queue<int> q;
int ans = 0;
for (int i = 1; i <= n; i++)
{
if (in[i] == 0)
{
q.push(i);
cnt[i] = 1;
}
}
while (!q.empty())
{
int u = q.front();
q.pop();
for (int i = 0; i < adj[u].size(); i++)
{
int v = adj[u][i];
cnt[v] += cnt[u];
cnt[v] %= mod;
in[v]--;
if (in[v] == 0)
{
if (out[v] == 0)
{
ans += cnt[v];
ans %= mod;
}
q.push(v);
}
}
}
cout << ans << endl;
return 0;
}