P3183 [HAOI2016]食物链-bfs+图论+dp

[HAOI2016]食物链

题目描述

如图所示为某生态系统的食物网示意图，据图回答第1小题现在给你n个物种和m条能量流动关系，求其中的食物链条数。物种的名称为从1到n编号M条能量流动关系形如a1 b1a2 b2a3 b3…am-1 bm-1am bm其中ai bi表示能量从物种ai流向物种bi,注意单独的一种孤立生物不算一条食物链

输入格式

第一行两个整数n和m,接下来m行每行两个整数ai bi描述m条能量流动关系。（数据保证输入数据符号生物学特点，且不会有重复的能量流动关系出现）1<=N<=100000 0<=m<=200000题目保证答案不会爆 int

输出格式

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

样例 #1

样例输入 #1

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

样例输出 #1

9

#include <bits/stdc++.h>
using namespace std;
#define de(x) cout<<x<<" ";
#define Pu puts("");
#define sf(x) scanf("%d",&x);
typedef long long ll;
const int N=1e6+10,mod=100003;
int n,m,ans;
int in[N],out[N];//记录入度
int var[N],hed[N],nxt[N],tot;
int cnt[N];//记录当前点所含有的食物链条数
void _add(int x,int y){
    var[++tot]=y;
    nxt[tot]=hed[x];
    hed[x]=tot;
}
void fun(){
    queue<int>q;
    for(int i=1;i<=n;i++){
        if(in[i]==0){
            cnt[i]=1;//加入入度为0的点，把其初始值赋值为1
            q.push(i);
        } 
    }
    int w,t;
    while(!q.empty()){
        w=q.front();q.pop();
        for(int i=hed[w];i;i=nxt[i]){
            t=var[i];

            if(out[t]==0){
                ans+=cnt[w];//当下一个节点是食物链的末端时
                //加上这个末端的前一个cnt[w]（表示是到这这个点时含有的链数）
            }else{//如果下一个点t不是食物链的末端，则点t汇总它之前所有的链数
                cnt[t]+=cnt[w];
            }
            in[t]--;
            if(in[t]==0){
                q.push(t);
            }
        }
    }
}
int main(){
    cin>>n>>m;
    int x,y;
    for(int i=1;i<=m;i++){
        sf(x)sf(y)
        _add(x,y);
        in[y]++;
        out[x]++;
    }
    fun();
    de(ans)
    return 0;
}