CF414B Mashmokh and ACM DP

最新推荐文章于 2022-07-20 23:14:44 发布

原创最新推荐文章于 2022-07-20 23:14:44 发布 · 225 阅读

0 ·

CC 4.0 BY-SA版权

动态规划专栏收录该内容

14 篇文章

订阅专栏

本文介绍了一种使用动态规划解决好数列计数问题的方法。好数列定义为一个数列中每个数都能被前一个数整除。文章详细阐述了状态转移方程的设计与实现，并给出完整代码。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目大意：

如果一个数列中，后一个数都能被前面一个数整除，那么就叫这个数列为好数列。输入n，k，求数列中最大元素为n，数列长度为k的好数列的种数（对1000000007取模）

思路

状态设计： $d p [i] [j]$ 结尾为数字 $i$ ，长度为 $j$

方程： $d p [i] [j] = d p [d] [j - 1]$ ， $d$ 为 $i$ 的约数

初始化： $d p [i] [1] = 1$

递推方式： $d p [d] [j] + = d p [i] [j - 1]$ , $d$ 为 $i$ 的倍数,我们可以递推出i的倍数

答案为： $∑i=1ndp[i][k]\sum_{i=1}^ndp[i][k]$

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int INF = 0x3f3f3f3f;
const double Pi = acos(-1);
namespace {
  template <typename T> inline void read(T &x) {
    x = 0; T f = 1;char s = getchar();
    for(; !isdigit(s); s = getchar()) if(s == '-') f = -1;
    for(;  isdigit(s); s = getchar()) x = (x << 3) + (x << 1) + (s ^ 48);
    x *= f;
  }
}
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define _for(n,m,i) for (register int i = (n); i <  (m); ++i)
#define _rep(n,m,i) for (register int i = (n); i <= (m); ++i)
#define _srep(n,m,i)for (register int i = (n); i >= (m); i--)
#define _sfor(n,m,i)for (register int i = (n); i >  (m); i--)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define lowbit(x) x & (-x)
#define pii pair<int,int>
#define fi first
#define se second
const int N = 2e3+5;
LL dp[N][N];
const LL Mod = 1e9+7;
int main() {
  int n, k; read(n); read(k);
  _rep(1, n, i) dp[i][1] = 1;
  _for(1, k, j) {
    _rep(1, n, i) {
      for(int d = i; d <= n; d += i) {
        dp[d][j+1] += dp[i][j];
        if(dp[d][j+1] >= Mod) dp[d][j+1] -= Mod;
      }
    }
  }
  LL ans = 0;
  _rep(1, n, i) {
    ans += dp[i][k];
    if(ans >= Mod) ans -= Mod;
  }
  cout << ans << endl;
}