題目鏈接：http://codeforces.com/problemset/problem/414/B

?

dp[i][j]表示長度為i、最后一個數字為j的合法序列總數

dp[1][1...2000]都是1

后面用dp[i-1][j] 去更新 dp[i][j*t] (1 <= j*t <= 2000) 即用因子去更新它的倍數

表面上看是2000^3的復雜度會爆 其實不用算那么多次

最外層循環是2000

分析第二層和第三層 需要算 2000/1 + 2000/2 + 2000/3 + 2000/4 + ... + 2000/2000 次

即2000 * (1/1 + 1/2 + 1/3 + ... + 1/2000)

后面的括號是一個【調和級數】 利用高等數學知識可知 它是有上限的 為自然對數e

所以本算法的時間復雜度實際上僅為 e*2000*2000?

綽綽有余

?

因為下標從1開始計數

然后一開始fill的時候沒注意看wa了幾炮T^T

?

#include <cstdio>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <stack>
#include <set>
#include <queue>
#include <vector>using namespace std;const int maxn = 2010;
const int M = 1000000007;int dp[maxn][maxn];int main()
{//freopen("in.txt", "r", stdin);int n, k;scanf("%d%d", &n, &k);fill(dp[1] + 1, dp[1] + 2001, 1);for(int i = 2; i <= k; i++){for(int j = 1; j <= 2000; j++){for(int t = 1; j*t <= 2000; t++){dp[i][j*t] = (dp[i][j*t] + dp[i-1][j]) % M;}}}int ans = 0;for(int i = 1; i <= n; i++)ans = (ans + dp[k][i]) % M;printf("%d\n", ans);return 0;
}

?

轉載于:https://www.cnblogs.com/dishu/p/4295089.html

總結

以上是生活随笔為你收集整理的cf 414B Mashmokh and ACM 动态规划的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。