51Nod_P1244 莫比乌斯函数之和(数论+杜教筛+哈希)

51Nod传送门

基准时间限制：3 秒 空间限制：131072 KB 分值: 640 难度：8级算法题

莫比乌斯函数，由德国数学家和天文学家莫比乌斯提出。梅滕斯(Mertens)首先使用μ(n)（miu(n)）作为莫比乌斯函数的记号。具体定义如下：
如果一个数包含平方因子，那么miu(n) = 0。例如：miu(4), miu(12), miu(18) = 0。
如果一个数不包含平方因子，并且有k个不同的质因子，那么miu(n) = (-1)^k。例如：miu(2), miu(3), miu(30) = -1,miu(1), miu(6), miu(10) = 1。
给出一个区间[a,b]，S(a,b) = miu(a) + miu(a + 1) + …… miu(b)。
例如：S(3, 10) = miu(3) + miu(4) + miu(5) + miu(6) + miu(7) + miu(8) + miu(9) + miu(10)
= -1 + 0 + -1 + 1 + -1 + 0 + 0 + 1 = -1。

Input
输入包括两个数a, b，中间用空格分隔(2 <= a <= b <= 10^10)
Output
输出S(a, b)。

Input示例
3 10

Output示例
-1

李陶冶 (题目提供者)

Sol:
杜教筛，筛积性函数的前缀和
哈希表来存一下，真是一道神题.
自己瞎写的，不好见谅

#include<cstdio>
#include<cstring>
#define N 20000000
#define P 1300000
#define Mod 2980001
int cnt,sum;int pr[P+10];short mu[N+10];bool bo[N+10];
long long a,b,n;
struct map{long long k;int v,nxt;}hash[Mod];int head[Mod];
void Pre(){
    mu[1]=1;
    for(int i=2;i<=N;i++){
        if(!bo[i]) pr[++cnt]=i,mu[i]=-1;
        for(int j=1;j<=cnt&&i*pr[j]<=N;j++){
            bo[i*pr[j]]=1;if(!(i%pr[j])) break;
            mu[i*pr[j]]=-mu[i];
        }
    }
    for(int i=2;i<=N;i++) mu[i]+=mu[i-1];
}
inline void add(long long k,int v){
    int ha=k%Mod;hash[++sum].k=k,hash[sum].v=v;
    hash[sum].nxt=head[ha],head[ha]=sum;
}
int calc(long long x,int res=0){
    if(x<=N) return mu[x];int ha=x%Mod;
    for(int i=head[ha];i;i=hash[i].nxt) if(hash[i].k==x) return hash[i].v;
    for(long long i=2,j;i<=x;i=j+1) j=x/(x/i),res+=calc(x/i)*(j-i+1);
    res=1-res;add(x,res);return res;
}
int main(){
    Pre();scanf("%lld%lld",&a,&b);
    printf("%d\n",calc(b)-calc(a-1));
    return 0;
}