本菜鸟看过一句话:动态规划的本质就算在一个有向无环图上寻找最优解。
本菜鸟很认可这句话,因为写题目的时候的的确确有这种感觉,但本菜鸟对这句话的运用只限于一些简单的题目。
这道题就可以非常简洁明了地诠释这句话,废话不多说,下面就开始题解;
P4017 最大食物链计数 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)
这道题我们可以先设置一个结构体,用来描述一个动物;
struct p{
vector<int> ne;//存储下一个动物的下标
int in=0,out=0,num=0;//in为该点的入度,out为该点的出度,num为从所以捕食者到该点的路径数量
};
动态规划是运筹学的一个分支,为多阶段决策问题,一个阶段的决策往往依赖其他决策的信息,而这个问题中在某一个点的决策就依赖于其父节点,我们要求出其父节点的num,然后有状态转移方程来确定该点的决策,即该点的num。
状态方程如下。
即所有父节点的num相加,那么为了正确地确定该点num,我们得计算出该点所有父节点的num,那么我们怎么确认一个点的父节点的num全部计算完呢?我们可以计算出一个点后,将其所以的子节点的入度-1,当某一个子节点为0时,说明该点的所有父节点num计算完成。
这一个过程有没有很熟悉?对!这就是拓扑排序的过程。
那么接下来就是代码实现了。
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
const int N=1000010;
struct p{
vector<int> ne;
int in=0; int num=0; int out=0;
}st[N];
queue<int> que;//队列存储可更新的点
int n,m;//n为食物数,m为关系数
int ans;//存储答案
void get(){
while(que.size()){
int t=que.front();
que.pop();
for(int i=0;i<st[t].ne.size();i++){
st[st[t].ne[i]].in--;//减去一个in
st[st[t].ne[i]].num+=st[t].num%80112002;//加上父节点的num
st[st[t].ne[i]].num%=80112002;
if(!st[st[t].ne[i]].in)//如果入度为0,则加入队列
que.push(st[t].ne[i]);
}
}
}
int main(){
cin>>n>>m;
for(int i=0;i<m;i++){//接收数据
int x,y; scanf("%d%d",&x,&y);
st[y].ne.push_back(x);
st[x].in++; st[y].out++;
}
for(int i=1;i<=n;i++){//先找到初始时可更新的点
if(!st[i].in)
que.push(i),st[i].num=1;
}
get();
for(int i=1;i<=n;i++){//寻找出度为0的点,即生产者,将其num叠加即可
if(!st[i].out)
ans+=st[i].num%80112002;
ans%=80112002;
}
cout<<ans;
}