bzoj4562 [Haoi2016]食物链

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

百度上搜索：bzoj 记忆化搜索，弹出来有人发过这题的题解，于是去bzoj上刷这题…讲道理，这是我在bzoj上做过的除A+B problem最水的一题…这哪是记忆化搜索啊，就是个纯搜索，10mins 一次A…
每次输入两个点要记录下出度入度，最后 ans 只加入读为 0 的点。
代码：

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
bool judge[100004];
int dp[100001],tot,head[100002],tov[200001],nex[200001],ru[100005],chu[100002];
void add(int a,int b)
{
    tov[++tot]=b;
    nex[tot]=head[a];
    head[a]=tot;
}
void dfs(int k)
{
    int t=head[k];
    if(!tov[t])
    {
        judge[k]=true;
        dp[k]=1;
    }
    int v=tov[t];
    while(v)
    {
        if(!dp[tov[t]])
            dfs(tov[t]);
        dp[k]+=dp[tov[t]];
        t=nex[t];
        v=tov[t];
    }
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        ru[b]++;
        chu[a]++;
        add(a,b);
    }
    for(int i=1;i<=n;i++)
        if(!dp[i])
            dfs(i);
    int ans=0;
    for(int i=1;i<=n;i++)
        if(!judge[i]&&ru[i]==0)
            ans+=dp[i];
    printf("%d",ans);
}