zoj4011 动态规划

https://vjudge.net/problem/ZOJ-4011

给定一个n，m，要求一个m为长度，并且从第一项起有b[i+1]%b[i]==0的数列，求出组合个数(1<=b[i]<=n)

利用dp来做，数组dp[i][j]表示在长度为i时，以j为结尾的数列，边界条件为dp[1][i]=1,因为当长度为1时，以i为结尾的数列只有一个。转移方程为dp[i][k]=d[i][k]+dp[i-1][j]   (k==a*j)

首先，循环条件为 k=j;k<=n;k+=j 因此k=a*j  ,dp[i][k]就是在长度为i，以k为结尾的数列，那么除了自己的个数以外，还有当长度为i-1,结尾为j的数列的个数dp[i-1][j]，因为当k放入第i个位置的时候，且k=a*j，那么我们就变成了dp[i][k],因此两个加上就是长度为i，结尾为k的数列的个数dp[i][k];

最后求长度为m，以1-n为结尾的个数就直接sum+=dp[m][i](1<=i<=n),记得取模就可以了

#include <iostream>
#include <vector>
#include <set>
#include <cstdio>
#include <cstring>
using namespace std;
const int mod=1e9+7;
long long dp[2009][2009];
//dp[i][j]表示在长度为i，以j结尾的数列的个数
int main()
{
    long long t,n,m,i,j,k;
    scanf("%lld",&t);
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        scanf("%lld%lld",&n,&m);
        for(i=1; i<=n; i++)
        {
            dp[1][i]=1;//长度为1的数列只有一个
        }
        for(i=2; i<=m; i++)
        {
            for(j=1; j<=n; j++)
            {
                for(k=j; k<=n; k+=j)
                {
                    dp[i][k]=(dp[i][k]+dp[i-1][j])%mod;
                }
            }
        }
        long long ans=0;
        for(i=1; i<=n; i++)
        {
            ans+=dp[m][i];
            ans%=mod;
        }
        printf("%lld\n",ans);
    }
}