P3183 [HAOI2016]食物链

时间:2023-03-09 01:33:39
P3183 [HAOI2016]食物链

题目描述

P3183 [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: 复制
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

拓扑排序+dp
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define inf 2147483647
const ll INF = 0x3f3f3f3f3f3f3f3fll;
#define ri register int
template <class T> inline T min(T a, T b, T c)
{
return min(min(a, b), c);
}
template <class T> inline T max(T a, T b, T c)
{
return max(max(a, b), c);
}
template <class T> inline T min(T a, T b, T c, T d)
{
return min(min(a, b), min(c, d));
}
template <class T> inline T max(T a, T b, T c, T d)
{
return max(max(a, b), max(c, d));
}
#define pi acos(-1)
#define me(x, y) memset(x, y, sizeof(x));
#define For(i, a, b) for (int i = a; i <= b; i++)
#define FFor(i, a, b) for (int i = a; i >= b; i--)
#define mp make_pair
#define pb push_back
const int maxn = ;
#define mod 100003
const int N=; // name*******************************
int n,m;
struct edge
{
int to,next;
} e[];
int in[N],out[N];
int tot=;
int Head[N];
int f[N];
int ans=;
// function******************************
void add(int u,int v)
{
e[++tot].to=v;
e[tot].next=Head[u];
Head[u]=tot;
out[u]++;
in[v]++;
}
void dfs(int u)
{
for(int p=Head[u]; p!=-; p=e[p].next)
{
int v=e[p].to;
f[v]+=f[u];
if(!--in[v])
dfs(v);
}
} //***************************************
int main()
{
//freopen("test.txt", "r", stdin);
cin>>n>>m;
me(Head,-);
For(i,,m)
{
int a,b;
cin>>a>>b;
add(a,b);
}
For(i,,n)
{
if(in[i]==&&out[i]!=&&f[i]==)
{
f[i]=;
dfs(i);
}
}
For(i,,n)
{
if(out[i]==)
ans+=f[i];
}
cout<<ans; return ;
}