B. Mashmokh and ACM（cf）dp

原题链接：Problem - 414B - Codeforces

题目大意：给出两个数n,k，可以取1~n中的数，组成长度为k的数组，数组满足a1 <= a2 <= ..<= ak,且相邻两数满足前面这个数可以除尽后面那个数ai +1 % ai == 0.问你有多少种情况。

一开始我想的是完全背包，因为1~n中的数可以用不止一次，但是限制是个数，我又觉得无从下手了。参考了其他大佬的代码，最后总结出来：

可以令数组dp【i】【z】//i表示当前数组长度，z表示以z值作为最后一个数的情况。z是j的倍数，所以有：

dp【i】【z】 += dp【i - 1】【j】

有时候是以右端点为第一维，有时候是以数组长度为第一维，视题目而定。

​​​​​​​注意初始化！

AC代码：

#include<bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ios ios::sync_with_stdio(false);cin.tie(0); cout.tie(0);
typedef pair<int, int> PII;
const double pi = acos(-1.0);
#define rep(i, n) for (int i = 1; i <= (n); ++i)
#define rrep(i, n) for (int i = n; i >= (1); --i)
typedef long long ll;
#define sqar(x) ((x)*(x))

const int N = 2010;
long long M = 1e9 + 7;
int a[N];
int dp[N][N]; //dp[i][z] ： 长度为i的序列中最后一个数为下标为z的情况数量
//状态转移方程：dp[i][z] = (dp[i][z] + dp[i -  1][j]) % M (z是j的倍数)
int main(){
  int n, k;
  cin >> n >> k;
  rep(i, n) dp[1][i] = 1; //初始化(长度为1且以i为最后一位元素的数列只有一种)
  rep(i, k)
   for(int j = 1; j <= n; j++)
    for(int z = j; z <= n; z += j){
      dp[i][z] = (dp[i][z] + dp[i - 1][j]) % M;
    }
  ll res = 0;
  rep(i, n) res = (res + dp[k][i]) % M;
  cout << res;
  return 0;
}