CodeForces - 1561D1 Up the Strip (simplified version) ( dp + 整除分块 + 前缀和 )

题目链接：点击进入

题目

题意

n 个台阶，问从 n 到 1 有多少种方案，答案 % m。
若是位于 x ，则有两种操作：
1、通过减法到达 [ 1 , x - 1 ]
2、通过除法到达 $\lfloor \frac{x}{z}\rfloor$ ( 下取整 , z 取值 [ 2 , x ] )

思路

dp [ i ] : 表示从 n 到 i 的方案数。倒序遍历 i ：
对于减法，记录一个前缀和 sum ，这个前缀和代表着，i 之前的位置 x 通过减法一步到达 i 的总方案数，$sum={\sum }_{x=i+1}^{n}dp[x]$
对于除法，对于位置 i ，它状态转移到：$\lfloor \frac{i}{2}\rfloor$,$\lfloor \frac{i}{3}\rfloor$…$\lfloor \frac{i}{i}\rfloor$。
整除分块 ？！
也就是对于 i ，我可以$\sqrt{n}$更新 dp [ x ]（ x 是 i 整除到达的位置 ）
时间复杂度$O(n\sqrt{n})$

代码

// Problem: Up the Strip (simplified version)
// Contest: Virtual Judge - CodeForces
// URL: https://vjudge.net/problem/CodeForces-1561D1
// Memory Limit: 131 MB
// Time Limit: 6000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define base 233 
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=1e6+10;
const int mod=1e9+7;
const double eps=1e-9;
const int N=5e3+10;
/*--------------------------------------------*/
inline int read()
{
    int k = 0, f = 1 ;
    char c = getchar() ;
    while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
    while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
    return k * f ;
}
/*--------------------------------------------*/

int n,m,dp[maxn],sum;

signed main() 
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);cout.tie(0);
	cin>>n>>m;
	dp[n]=1;
	for(int i=n;i>=1;i--)
	{
		dp[i]=(dp[i]+sum)%m;
		for(int l=2,r;l<=i;l=r+1)
		{
			r=i/(i/l);
			dp[i/l]=(dp[i/l]+(r-l+1)*dp[i])%m;
		}
		sum=(sum+dp[i])%m;
	}
	cout<<dp[1]<<endl;
    return 0;
}