动态规划dp P2193 HXY和序列洛谷题解 dp入门

题目描述
HXY突发奇想，她想要找到一个正整数序列，满足序列中所有的数不超过n，序列长度为p，且除了第一个数外，所有的数都能被前一个数整除（即是前一个数的倍数）。很快她找到了一个这样的序列。可是她觉得还不够，想要知道这样的序列有多少个，可她被惊人的数据范围吓怕了。现在她找到了你，请你来帮助她解决这个问题。（因为结果可能会很大，请输出对1e9+7取模后的值）

输入格式
第一行，两个正整数n，p。

输出格式
仅一行，一个正整数，表示满足条件序列的个数对1e9+7取模后的值。

输入输出样例
输入 #1复制
3 2
输出 #1复制
5
输入 #2复制
6 4
输出 #2复制
39
说明/提示
数据范围：

对于10%的数据，p=1；

对于30%的数据，1<=n，p<=10；

对于60%的数据，1<=n，p<=500;

对于100%的数据，1<=n，p<=2000。

我们设f[i][j]为长度为i ，j开头的定义
知道f[i][j] = f[i - 1][j] + f[i -1][ j 2] + f[i-1][j3]…

在i == 1时 不管j为多少他的序列都只能为
for (int i = 1; i <= n; i++ ) f[1][i] = 1;
这位初始化

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
const int MOD = 1e9 + 7;
int n, p;
int f[2005][2005];
int main() {
	cin >> n >> p;
	//i长度j为第一个数
	for (int i = 1; i <= n; i++) f[1][i] = 1;
	for (int i = 2; i <= p; i++) // 序列的长度
		for (int j = 1; j <= n; j++) //序列的开头
			for (int k = j; k <= n; k += j) {
				f[i][j] += f[i - 1][k];
				f[i][j] %= MOD;
			}

	int ans = 0;
	//这里我们需要把所有长度为p的序列的数量累加起来
	for (int i = 1; i <= n; i++){
	//f[p][i]长度为p，i开头的序列的数量
		ans = (ans + f[p][i]) % MOD;
	}	
	cout << ans;
	return 0;
}