数论分块

数论分块+两道简单的例题

作用：可以在O$\sqrt{n}$的时间复杂度求出以下式子的值：${\sum }_{i=1}^n$$\lfloor$$\frac{n}{i}$$\rfloor$

证明

传送门
题意：输入n与m，输入${\sum }_{i=1}^n$m%i

分析：m%i=m-$\lfloor$$\frac{m}{i}$$\rfloor$ $\times$ i

所以可以化简结果为m$\times$n-${\sum }_{i=1}^n$$\lfloor$$\frac{m}{i}$$\rfloor$ $\times$ i

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long 
int main()
{
    int n,k;cin>>n>>k;
    ll ans=1ll*n*k;
    int l,r;
    for(int l=1;l<=n;l=r+1)
    {
        if(k/l==0) break;
        r=min(n,k/(k/l));
        ans-=1ll*(k/l)*(r-l+1)*(r+l)/2;
    }
    cout<<ans<<endl;
}

传送门

这个题两个版本，数据范围不同，对于数据较小的这个我们可以利用分块，O(N$\times$$\sqrt{n}$)是可以过的

思路：分别开两个数组dp,pre,dp代表的是第i位到1的情况，pre代表的是前缀和

分析：由题意可知，有两种转移方式可以推出转移方程

dp[i]=pre[i-1]+${\sum }_{x=2}^i$dp$\lfloor$$\frac{i}{x}$$\rfloor$

代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=2e5+10;
int f[N],pre[N];
int n;
ll m;
void init(int n)
{
	f[1]=1;
	pre[1]=1;
	for(int i=2;i<=n;i++)
	{
		f[i]=pre[i-1];
		int l,r;
		for(int l=2;l<=i;l=r+1)
		{
		    r=min(i,i/(i/l));
		    f[i]=(f[i]+1ll*f[i/l]*(r-l+1))%m;
		}
		pre[i]=1ll*(f[i]+pre[i-1])%m;
	}
}
int main()
{
	cin>>n>>m;
	init(n);
	cout<<f[n]<<endl;
	return 0;
}

下面讲O(Nlog(N))的做法

方法一是从1开始正向推，现在开始反向推

对于k可以从k+1-n通过方法1转移，前缀和就可以，然后对于方法二，考虑对于k来说，当除数是x时，哪些可以通过除x转移到k，很显然范围是k$\times$x – k$\times$x+x-1

根据这个范围再加上前缀和我们就可以在O(Nlog(N))求出结果

代码：

#include<bits/stdc++.h>
#define ll long long
#define LOCAL
using namespace std;
const int N=4e6+10;
ll dp[N],pre[N];
int n,mod;
int main()
{
    #ifdef ONLINE_JUDGE
    #else
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    #endif
    cin>>n>>mod;
    dp[n]=1,pre[n]=1;
    for(int i=n-1;i>=1;i--)
    {
        dp[i]=pre[i+1];
        for(int x=2;x*i<=n;x++)
        {
            int l=x*i,r=min(n,x*i+x-1);
            dp[i]=(dp[i]+pre[l]-pre[r+1])%mod;
        }
        pre[i]=(dp[i]+pre[i+1])%mod;
    }
    cout<<dp[1]<<endl;
    return 0;
}