bzoj4562【HAOI2016】食物链

4562: [Haoi2016]食物链

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 216  Solved: 173
[Submit][Status][Discuss]

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

9

DP水题

#include<iostream>#include<cstdio>#include<cstdlib>#include<cmath>#include<cstring>#include<algorithm>#define F(i,j,n) for(int i=j;i<=n;i++)#define D(i,j,n) for(int i=j;i>=n;i--)#define ll long long#define N 100005using namespace std;int n,m,cnt,ans,f[N],head[N],in[N];bool tag[N];struct edge{int next,to;}e[N*2];inline int read(){int x=0,f=1;char ch=getchar();while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}inline void add_edge(int x,int y){e[++cnt]=(edge){head[x],y};head[x]=cnt;}void dfs(int x){for(int i=head[x];i;i=e[i].next){int y=e[i].to;f[y]+=f[x];in[y]--;if (!in[y]) dfs(y);}}int main(){n=read();m=read();F(i,1,m){int x=read(),y=read();add_edge(x,y);in[y]++;}F(i,1,n) tag[i]=!in[i]&&head[i];F(i,1,n) if (tag[i]) f[i]=1,dfs(i);F(i,1,n) if (!head[i]) ans+=f[i];printf("%d\n",ans);return 0;}