[杜教筛] 51Nod 1244: 莫比乌斯函数之和

最新推荐文章于 2020-08-30 19:31:23 发布

Lynstery

最新推荐文章于 2020-08-30 19:31:23 发布

阅读量683

点赞数

分类专栏： --------------数学相关------------ 杜教筛

本文链接：https://blog.csdn.net/CHHNZ/article/details/74276079

版权

102 篇文章 0 订阅

订阅专栏

杜教筛

6 篇文章 0 订阅

订阅专栏

题意

求 $\sum_{i=L}^R\mu(i)$ ， $L,R\le10^{11}$

题解

杜教筛模板题。

g (1) S (n) = \sum i = 1 n (f * g) (i) - \sum i = 2 n g (i) S (⌊ n i ⌋)

$g(1)S(n)=\sum_{i=1}^n(f*g)(i)-\sum_{i=2}^ng(i)S(\lfloor\frac{n}{i}\rfloor)$

μ∗1=ϵ $\mu∗1=\epsilon$ ，

g $g$ 取常数函数

1 $1$ 就好了:

S (n) = \sum i = 1 n - \sum i = 2 n g (i) S (⌊ n i ⌋)

$S(n)=\sum_{i=1}^n -\sum_{i=2}^ng(i)S(\lfloor\frac{n}{i}\rfloor)$

#include<cstdio>
#include<map>
#include<cstring>
#include<algorithm>
#include <tr1/unordered_map> 
using namespace std;
using namespace std::tr1; 
const int maxn=10000005, N=10000000;
typedef long long LL;
int Q,p[maxn],mu[maxn];
bool vis[maxn];
void get_mu(){
    memset(vis,1,sizeof(vis));
    mu[1]=1;
    for(int i=2;i<=N;i++){
        if(vis[i]) p[++p[0]]=i, mu[i]=-1;
        for(int j=1;j<=p[0]&&p[j]*i<=N;j++){
            vis[p[j]*i]=false;
            if(i%p[j]==0){ mu[p[j]*i]=0; break; }
            mu[p[j]*i]=-mu[i];
        }
    } 
    for(int i=1;i<=N;i++) mu[i]+=mu[i-1];
}
unordered_map< LL, int > lst; 
int Sum(LL n){
    if(n<=N) return mu[n];
    if(lst.find(n)!=lst.end()) return lst[n];
    int res=1;
    for(LL i=2,nxt;i<=n;i=nxt+1){
        nxt=n/(n/i);
        res-=Sum(n/i)*(nxt-i+1);
    }
    return lst[n]=res;
}
int main(){
    freopen("51nod1244.in","r",stdin);
    freopen("51nod1244.out","w",stdout);
    get_mu();
    LL L,R; scanf("%lld%lld",&L,&R); 
    printf("%d\n",Sum(R)-Sum(L-1));
    return 0;
}