ABC132 – F – Small Products（dp）

https://atcoder.jp/contests/abc132/tasks/abc132_f
这道题还是有点想不清楚，以后再来看。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5, mod = 1e9 + 7;

//dp[i][j]：假设第i位的数是x，则dp[i][j]表示所有使得n / x <= sr[j]的x的答案之和
//j的取值的个数是sqrt(n)级别的
int dp[105][66666], sl[66666], sr[66666];
int main()
{
    ios::sync_with_stdio(false); cin.tie(false);
    int n, k;
    cin >> n >> k;

    int cnt = 0, l = 1, r = 1;
    while (l <= n){
        r = n / (n / l);
        sl[++cnt] = l;
        sr[cnt] = r;
        l = r + 1;
    }

    for (int i = 1; i <= cnt; ++i) dp[1][i] = sr[i];
    for (int i = 2; i <= k; ++i){
        for (int j = 1; j <= cnt; ++j){
            //首先dp[i][j]的答案至少包括dp[i][j - 1];
            //然后由于j加了1，这其中包含的数字的个数为(sr[j] - sl[j] + 1)
            //j和cnt - j + 1是对称的位置，sr[cnt - j + 1] * sr[j]刚好不超过n
            //dp[i - 1][cnt - j + 1]是上一次计算好的答案
            //两个相乘就是dp[i][j]的结果了
            dp[i][j] = ((ll)dp[i - 1][cnt - j + 1] * (sr[j] - sl[j] + 1) + dp[i][j - 1]) % mod;
        }
    }

    cout << dp[k][cnt] << endl;
    return 0;
}