二维dp 之 完美序列

来源：中国石油大学训练平台

4874: 完美序列

时间限制: 1 Sec  内存限制: 128 MB
提交: 28  解决: 17
[ 提交][ 状态][ 讨论版]

题目描述

已知一个长度为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取余后的结果。

样例输入

33 26 42 1

样例输出

5392

【分析】：

根据题意推得，可以打出一张表格，输入n，l直接询问。

用dp[i][j]表示长度为i时，以数字j结尾的，序列数量

状态方程：dp[i][j] = sum ( dp[i-1][k] ) ，sum()的意思是累加所有的dp[i-1][k]，其中k是j的因子。

初态：dp[1][j]=1，长度为1时的数量

根据初态，递推打出2000*2000的表格。

当输入n，l时，只需把第n行前l项累加即可

【代码】：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll dp[2020][2020];
int T,n,m;
int main()
{
    for(int i=1;i<=2002;i++)dp[1][i]=1;
    for(int i=2;i<=2002;i++)//长度
    {
        for(int j=1;j<=2002;j++)
            for(int k=1;k*j<=2002;k++)
            {
                dp[i][j*k]+=dp[i-1][j];
                dp[i][j*k]%=mod;
            }
    }
    for(int i=1;i<=2002;i++)
    {
        for(int j=2;j<=2002;j++)
        {
            dp[i][j]+=dp[i][j-1];
            dp[i][j]%=mod;
        }
    }
    /*for(int i=1;i<=8;i++)
    {
        printf("i = %d  ",i);
        for(int j=1;j<=8;j++)
            printf("%6lld",dp[i][j]);
        puts("");
    }*/
    cin>>T;
    while(T--)
    {
        cin>>n>>m;
        cout<<dp[m][n]<<endl;
    }
}
 
/**************************************************************
    Problem: 4874
    User: summer17083
    Language: C++
    Result: 正确
    Time:932 ms
    Memory:33576 kb
****************************************************************/