CodeForces - 1561D2 Up the Strip ( dp + 前缀和 )

题目链接：点击进入

题目

题意

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

思路

与此题题意相同，只是 n 变为了4e6，$O(n\sqrt{n})$ 不可过
dp [ i ] : 表示从 n 到 i 的方案数。倒序遍历 i ：
对于减法，依旧是记录前缀和 sum [ i ] ，这个前缀和代表着，i 之前的位置 x 通过减法一步到达 i 的总方案数，$sum={\sum }_{x=i+1}^{n}dp[x]$
对于除法，上一题是用台阶 i ，更新 x 的 dp 值 ( x 是 i 整除到达的位置 )，也就是用 i 的值更新低位置的值。现在尝试用高位置的值更新 i 的值：
对于 z = 2 ，可以通过除法一步到达的高位有 [ 2i , 2i + 1 ]
对于 z = 3 ，可以通过除法一步到达的高位有 [ 3i , 3i + 2 ]
对于 z = 4 ，可以通过除法一步到达的高位有 [ 4i , 4i + 3 ]
…
因此可以枚举 z ,对于 z ，可以通过除法一步到达的高位有 [ zi , zi + z - 1 ]
转换一下，也就是对于位置 i 的 dp 值，可以由位置 [ zi , zi + z - 1 ] 的 dp 值更新，更新的总值就是 ${\sum }_{x=zi}^{zi+z-1}dp[x]$，这个值可以通过前缀和做差求得，即 dp [ i ] += sum [ zi ] - sum [ min ( zi + z , n + 1 ) ] 。

代码

// Problem: Up the Strip
// Contest: Virtual Judge - CodeForces
// URL: https://vjudge.net/problem/CodeForces-1561D2
// 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=4e6+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[maxn];

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