序列

最新推荐文章于 2024-03-24 23:25:39 发布

原创最新推荐文章于 2024-03-24 23:25:39 发布 · 267 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#c++

all 同时被 2 个专栏收录

51 篇文章

订阅专栏

2 篇文章

订阅专栏

给定一个整数N和序列长度k，一个'好序列'定义为每个后续元素是前一个元素的倍数。本题要求计算满足条件的'好序列'数量，并对10^9 + 7取模。采用动态规划（DP）思路，状态转移方程为f[i+1][jn]+=f[i][j]，其中f[i][j]表示第i个数选择为j的方案数。边界条件是f[1][i]为1。

题目描述

1138.序列
时间限制:1000MS 内存限制:128000KB

题目描述
一个长度为k的整数序列b1，b2，…，bk（1≤b1≤b2≤…≤bk≤N）称为“好序列”当且仅当后一个数是前一个数的倍数，即bi+1是bi的倍数对任意的i（1≤i≤k-1）成立。
给定N和k，请算出有多少个长度为k的“好序列”，答案对1000000007取模。

输入
输入共1行，包含2个用空格隔开的整数N和k。
输出
输出共1行，包含一个整数，表示长度为k的“好序列”的个数对1000000007取模后的结果。

输入样例
3 2
输出样例
5

说明

【输入输出样例说明】
“好序列”为：[1,1]，[1,2]，[1,3]，[2,2]，[3,3]。

【数据说明】
对于40%的数据，1≤N≤30，1≤k≤10。
对于100%的数据，1≤N≤2000，1≤k≤2000。

思路：DP

设f[i][j]为第i个数放j的方案总数
顺推：
若f[i][j]已经放置好，则f[i+1][jn]的方案数就+=f[i][j] //n为倍数
方程：f[i+1][jn]+=f[i][j]
边界：f[1][i]为1

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
#include <cmath>
using namespace std;

const int N=1e6+100,M=1000000007;
int a[2010][2010];
int num,ans,n,m;

int main ()
{
	scanf ("%d%d",&m,&n);
	
	for (int i=1;i<=m;i++) a[1][i]=1;    //初始化，指若第一个数放i，则必然只有一个数
	
	for (int i=1;i<n;i++)          //放第几个数
	  for (int j=1;j<=m;j++)           //第i个数放j
	    for (int k=1;k*j<=m;k++)            //j的倍数
	      a[i+1][j*k]=(a[i+1][j*k]+a[i][j])%M;        //方程
	      
	for (int i=1;i<=m;i++) ans+=a[n][i],ans%=M;       //答案为最后一个数放i的所有方案的和
	
	printf ("%d",ans%M);
	
	return 0;
}