【BZOJ4562】【食物链】【dp+拓扑排序】

最新推荐文章于 2026-03-23 23:04:31 发布

原创最新推荐文章于 2026-03-23 23:04:31 发布 · 1k 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

dp 专栏收录该内容

63 篇文章

订阅专栏

这篇博客详细介绍了如何利用动态规划和拓扑排序解决BZOJ4562题目的食物链问题，通过实例展示了输入、输出格式，并给出了样例输入和输出。

Description

如图所示为某生态系统的食物网示意图，据图回答第1小题

现在给你n个物种和m条能量流动关系，求其中的食物链条数。

物种的名称为从1到n编号

M条能量流动关系形如

a1 b1

a2 b2

a3 b3

......

am-1 bm-1

am bm

其中ai bi表示能量从物种ai流向物种bi,注意单独的一种孤立生物不算一条食物链

Input

第一行两个整数n和m,接下来m行每行两个整数ai bi描述m条能量流动关系。

（数据保证输入数据符号生物学特点，且不会有重复的能量流动关系出现）

1<=N<=100000 0<=m<=200000

题目保证答案不会爆 int

Output

一个整数即食物网中的食物链条数

Sample Input

10 16
1 2
1 4
1 10
2 3
2 5
4 3
4 5
4 8
6 5
7 6
7 9
8 5
9 8
10 6
10 7
10 9

Sample Output

题解：按拓扑序统计路径条数即可.

代码：

#include<iostream>
#include<cstdio>
#define N 100010
#define M 200010
using namespace std;
int n,m,x,y,ans,f[N],in[N],out[N],a[N],point[N],next[M],cnt,q[N*10];
struct use{
   int st,en;
}e[M];
void add(int x,int y){
  in[y]++;out[x]++;a[y]++;next[++cnt]=point[x];point[x]=cnt;e[cnt].en=y;
}
void dp(){
  int h(0),t(0);
  for (int i=1;i<=n;i++) if (!in[i]) q[++t]=i,f[i]=1;
  while (h<t){
    int u=q[++h];
    for (int i=point[u];i;i=next[i]){
       f[e[i].en]+=f[u];in[e[i].en]--;
       if (!in[e[i].en]) q[++t]=e[i].en;
    }
  }
  for (int i=1;i<=n;i++) if (!out[i]&&a[i]) ans+=f[i];
  cout<<ans<<endl;  
}
int main(){
   scanf("%d%d",&n,&m);
   for (int i=1;i<=m;i++){
     scanf("%d%d",&x,&y);
     add(x,y);
   }
   return dp(),0;   
}