完美序列

題目描述

已知一個長度為l的序列：b1，b2，b3，…，bl (1<=b1<=b2<=b3<=…<=bl<=n)。若這個序列滿足每個元素是它後續元素的因子，換句話說就是對於任意的i (2<=i<=l)都滿足bi%bi-1=0 (其中“%”代表求餘)，則稱這個序列是完美的。你的任務是對於給定的n和l，計算出一共有多少序列是完美序列。由於答案很大，所有輸出答案對1000000007取餘後的結果。

輸入

輸入的第一行為一個正整數T (T<=1000)，代表一共有T組測試資料。

每組測試資料包含兩個正整數n，l (1<=n, l<=2000)，分別代表序列中元素大小的最大值和序列的長度。

輸出

對於每組測試資料，輸出一行包含一個整數，代表答案對1000000007取餘後的結果。

樣例輸入

3 2
6 4
2 1

樣例輸出

5

39

2

思路:二維dp[i][j],i表示序列的長度，j表示序列最後一個數的值

首先打表:

i到maxn所有的dp[1][i]賦初始值為1

接著三個for迴圈，i,j,k,分別代表序列長度，序列最後一個數的值，滿足ai%ai-1=0序列需要的值(k+=j)

狀態轉移方程：dp[i][k]=(d[i-1][j]%mod+dp[i][k]%mod)%mod

最後結果把從1到最大值val,長度為len的dp[len][i]方案數累加即可

AC Code:

#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<map>
#include<vector>
#include<stack>
#include<queue>
#include<set>
#define Max 0x3f3f3f3f
#define Min 0xc0c0c0c0
#define mst(a) memset(a,0,sizeof(a))
#define f(i,a,b) for(int i=a;i<b;i++)
using namespace std;
typedef long long ll;
const int maxn=2005;
const int mod=1e9+7;
int dp[maxn][maxn];
void init(){
    mst(dp);
    for(int i=1;i<maxn;i++){
        dp[1][i]=1;
    }
    for(int i=2;i<maxn;i++){
        for(int j=1;j<maxn;j++){
            for(int k=j;k<maxn;k+=j){
                dp[i][k]=(dp[i-1][j]%mod+dp[i][k]%mod)%mod; //狀態轉移方程
            }
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    int T,val,len;
    init();
    cin>>T;
    while(T--){
        cin>>val>>len;
        ll ans=0;
        for(int i=1;i<=val;i++){
            ans=(ans+dp[len][i])%mod;
        }
        cout<<ans<<endl;
    }
    return 0;
}