GDUT专题学习2 C - 最大食物链计数

C - 最大食物链计数

Background

你知道食物链吗？Delia 生物考试的时候，数食物链条数的题目全都错了，因为她总是重复数了几条或漏掉了几条。于是她来就来求助你，然而你也不会啊！写一个程序来帮帮她吧。

Description

给你一个食物网，你要求出这个食物网中最大食物链的数量。

（这里的“最大食物链”，指的是生物学意义上的食物链，即最左端是不会捕食其他生物的生产者，最右端是不会被其他生物捕食的消费者。）

Delia 非常急，所以你只有 1 秒的时间。

由于这个结果可能过大，你只需要输出总数模上80112002 的结果。

Input

第一行，两个正整数 n、m，表示生物种类 n 和吃与被吃的关系数 m。

接下来 m 行，每行两个正整数，表示被吃的生物A和吃A的生物B。

Output

一行一个整数，为最大食物链数量模上 80112002 的结果。

Sample 1

Inputcopy	Outputcopy
5 7 1 2 1 3 2 3 3 5 2 5 4 5 3 4	5

Hint

各测试点满足以下约定：

【补充说明】

数据中不会出现环，满足生物学的要求。（感谢 @AKEE ）

思路

        题目大意是找到完全食物链（从生产者到最高级消费者）的数目，可用记忆化搜索解决。首先根据捕食关系确认每个点的入度和出度，入度为0的为生产者（起点），出度为0的为最高级消费者（终点），用记忆化搜索找到从起点到各个终点的路径数目即为完全食物链的数目。

#include<iostream>
using namespace std;
const int mod=80112002;
int n,m,a,b,ans=0;
int flag[5005][5005];//记录捕食关系
int from[5005];//入度
int to[5005];//出度
int index[5005];//标记
int dfs(int x){
	if(!to[x])return 1;
	if(index[x])return index[x];
	index[x]=0;
	for(int i=1;i<=n;i++){
		if(flag[i][x]){
			index[x]=(index[x]+dfs(i))%mod;
		}
	}
	return index[x];
}
int main(){

	cin>>n>>m; 
	while(m--){
		cin>>a>>b;
		flag[a][b]=1;
		from[a]++;
		to[b]++;	
	}
	for(int i=1;i<=n;i++){
		if(!from[i]){//入度为0
			ans=(ans+dfs(i))%mod;
		}
	}
	cout<<ans<<endl;
	return 0;
}