Codeforces Round #740 (Div. 2）D1. Up the Strip (simplified version) DP，整除分块优化

题意：你现在的值为n，你可以进行两种操作
1.n - x [ 1 <= x <= n - 1 ]
2.n / x [ 2 <= x <= n ]
问你从n变到1的方案数对m取模
数据范围 2 ≤ n ≤ 2e5 1e8 < m < 1e9，m为素数

题解：设f[i]表示把i变到1的所有方案数,我们可以很容易地写出来O(n^2)的状态转移表达式。然后去优化即可，一个是前缀和优化，一个是整数分块优化。

AC代码：

#include <bits/stdc++.h>
#define int long long
#define x first
#define y second
using namespace std;
typedef pair<int,int> PII;
const int N=1e6+5;
int n,p;
int f[N];
main(){
	cin>>n>>p;
	f[1]=1;
	
	int sum=f[1];
	for(int i=2;i<=n;i++){
		f[i]=(f[i]+sum)%p;
		
		int cnt=0,x=f[i];
		for(int l=1,r;l<=i;l=r+1){//整数分块板子 
			r=i/(i/l);
			cnt=(cnt+(r-l+1)*f[i/l]%p)%p;
		}
		f[i]=(f[i]+cnt-x+p)%p;//这里减去x是因为原来式子是从l=2开始的累加i/l的和，而整数分块是从1开始的，所以要减去 
		sum=(sum+f[i])%p;
	}
	
	cout<<f[n]<<endl;
}