先总结一下,杜教筛的的精髓之处我认为在于通过两个积性函数做狄利克雷卷积以后就可以对其进行整除分块了,又因为一般用到杜教筛的题目数据量都特别大,是o(n)时间都跑不过来的数据,所以肯定不能预处理。但是这样的题样例数量不会太大,你只能每一次都计算结果,不能与处理出来结果,所以你需要用到记忆化搜索或hashmap去存整除分块的过程中出现的s(n)。以后用到s(n)可以o(sqrt(n))查询调用。
记忆化搜索实现:
我们需要先预处理出来一部分莫比乌斯前缀和,因为你预处理出来的越多,后面记忆化的时候需要查询是否已经记忆化搜索出来的部分越少,同时需要你去记忆化搜索去递归调用的次数也少,大大优化了时间复杂度,就拿本题为例,我先预处理出来1e5的莫比乌斯函数前缀和就T了,而我预处理出来5*1e6的莫比乌斯函数前缀和就1s过题。
我们都知道整除分块的时间复杂度为O(sqrt(n)),因为n / x 一共只有 2*sqrt(n)组解。
令s(n)为莫比乌斯函数的前n项和
代码如下:
#include<bits/stdc++.h>
const int maxn = 1e6+5;
const int mod = 857777;
typedef long long ll;
using namespace std;
int mu[maxn], prime[maxn],nxt[mod],had[mod],tot;
ll tmp[mod],val[mod];
bool isprime[maxn];
void init(){//线性筛预处理出来1e6范围内的莫比乌斯函数和,优化使搜索加速
int cnt=0;
mu[1]=1;
for(int i=2;i<maxn;i++)
{
if(!isprime[i]){
prime[cnt++]=i;
mu[i]=-1;
}
for(int j=0;j<cnt&&i*prime[j]<maxn;j++)
{
isprime[i*prime[j]]=1;
if(i%prime[j])
mu[i*prime[j]]=-mu[i];
else
{
mu[i*prime[j]]=0;
break;
}
}
}
for(int i=2;i<maxn;i++)
mu[i]=mu[i-1]+mu[i];
}
void add(ll x,ll y,ll z){//记录下来记忆化过程中搜索到的结果
val[++tot] = y;//记录被计算的n值
nxt[tot] = had[x];//记忆化中被记录x的下一个位置(也就是没有0)
had[x] = tot;//记录下上一次被搜索到的值的位置
tmp[tot] = z;//记录结果
}
ll cal(ll x){//记忆化搜索,i从2到n,很多n/i是相等的。
//所以我们在这里用到分块优化,大佬算出杜教筛时间复杂度大致达到o(n^(3/4))
if(x<maxn) return mu[x];
int t = x%mod,ret=1;
for(int i=had[t];i;i=nxt[i])
if(val[i]==x)
return tmp[i];
ll l=2,r;
while(l <= x)//计算每一个s(x)的值
{
r = x / (x / l);
ret -= 1ll*cal(x/l)*(r-l+1);
l = r + 1;
}
add(t,x,ret);
return ret;
}
int main() {
init();
ll l,r;
tot=0;
scanf("%lld %lld",&l,&r);
printf("%lld\n",cal(r)-cal(l-1));
return 0;
}