牛客：小易喜欢的数列

最新推荐文章于 2022-07-07 10:30:14 发布

ShellDawn

最新推荐文章于 2022-07-07 10:30:14 发布

阅读量265

点赞数

分类专栏： NewCoder OJ 文章标签：牛客

本文链接：https://blog.csdn.net/ShellDawn/article/details/93377473

版权

NewCoder OJ 专栏收录该内容

37 篇文章 0 订阅

订阅专栏

动态规划入门题：

#include<bits/stdc++.h>

using namespace std;

#define maxn 100005
#define LL long long
#define MM(x) memset(x,0,sizeof(x));
#define mod 1000000007

LL N,M;
LL dp[12][maxn];// dp[i][j] i位末尾为j的序列的数量

int main(){
    cin>>N>>M;
    MM(dp);
    for(int i=1;i<=M;i++) dp[1][i] = 1;
    for(int i=2;i<=N;i++){
        LL Sum = 0;
        for(int j=1;j<=M;j++) Sum += dp[i-1][j];
        Sum %= mod;
        for(int j=1;j<=M;j++){
            dp[i][j] = Sum;
            // 只有a为b的2倍以上倍数才非法
            for(int k=2;j*k<=M;k++){
                dp[i][j] = (dp[i][j] - dp[i-1][j*k] + mod)%mod;
            }
        }
    }
    LL ans = 0;
    for(int i=1;i<=M;i++){
        ans += dp[N][i];
        ans %= mod;
    }
    cout<<ans<<endl;
    return 0;
}