题目描述
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;
}