杜教筛

杜教筛

     

             

于是μ的前缀和就可以求了， 其他就是乘了一些奇怪的数也差不多是这样化

不用哈希，因为我们预处理了至少√n个，所以我们要算的所有的数都是小于√n的，

当我们在计算n的答案时，在计算过程中所要访问的所有数x，它的n/x的值都是唯一的，

所以我们就直接用n/x作为它的哈希值， 计算下一个答案时再清空即可

 

【51nod 1244】莫比乌斯函数之和

求：

思路：这个直接求出1~b的减去1~(a-1)的即可

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=2000500,M=1333331;
LL m,m1,n;
bool prz[N];
int pr[N/5],mu[N],Hx[M+1][2],co;
int HX(LL q){
    int i=q%M;
    while(Hx[i][0]&&Hx[i][0]!=q)i=(i+1)%M;
    return i;
}
LL Gmu(LL q){
    if(q<=n)return mu[q];
    int t=HX(q);
    if(Hx[t][0])return Hx[t][1];
    Hx[t][0]=q;
    LL ans=0,i=2;
    while(i<=q){
        LL ni=q/(q/i);
        ans+=(ni-i+1)*Gmu(q/i);
        i=ni+1;
    }
    return Hx[t][1]=1-ans;
}
int main(){
    scanf("%lld%lld",&m1,&m);
    n=sqrt(m)*20;
    mu[1]=1;
    for(int i=2;i<=n;i++){
        if(!prz[i])pr[++pr[0]]=i,mu[i]=-1;
        for(int j=1;j<=pr[0];j++){
            LL t=pr[j]*i;
            if(t>n)break;
            prz[t]=1;
            if(i%pr[j]==0)break;
            mu[t]=-mu[i];
        }
    }
    LL q=0;
    for(int i=2;i<=n;i++)mu[i]+=mu[i-1];
    printf("%lld\n",Gmu(m)-Gmu(m1-1));
    return 0;
}