[51Nod1244]莫比乌斯函数之和

最新推荐文章于 2021-01-12 08:04:42 发布

清疚

最新推荐文章于 2021-01-12 08:04:42 发布

阅读量167

点赞数

分类专栏：数论文章标签：杜教筛 51nod

本文链接：https://blog.csdn.net/qq_38064038/article/details/79317906

版权

数论专栏收录该内容

13 篇文章 0 订阅

订阅专栏

题意

给定 $l,r$ ，求 $\sum_{i=l}^r \mu(i)$

解法

和上一题差不多，只是推的式子不一样而已
$[n=1]=\sum_{d|n}\mu(d) · · · ①$
由①得： $\mu(n)=[n=1]-\sum_{d|n,d<n}\mu(d) · · · ②$
令 $M(n)=\sum_{i=1}^n\mu(i)$ ，将②带入得到：

$M (n) = \sum i = 1 n [[i = 1] - \sum d | i, d < i μ (d)] = 1 - \sum i = 2 n \sum d | i, d < i μ (d)$ $M(n)=\sum_{i=1}^n[[i=1]-\sum_{d|i,d<i}\mu(d)]=1-\sum_{i=2}^n\sum_{d|i,d<i}\mu(d)$
采用和上一题一样的套路，变化得到：
$M (n) = 1 - \sum i = 2 n \sum d = 1 ⌊ n i ⌋ μ (d) = 1 - \sum i = 2 n M (⌊ n i ⌋)$ $M(n)=1-\sum_{i=2}^n\sum_{d=1}^{\lfloor \frac{n}{i} \rfloor}\mu(d)=1-\sum_{i=2}^n M(\lfloor \dfrac{n}{i} \rfloor)$
然后就可以直接杜教筛了

复杂度

O( $k*n^\frac{2}{3}$ )， k <script type="math/tex" id="MathJax-Element-1249">k</script>为常数

代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<map>
#define Rint register int
#define Lint long long int
using namespace std;
const int N=5000010;
bool vis[N];
int pri[N/10],miu[N],sum[N];
int cnt;
map<Lint,int> f;
void Prepare()
{
    miu[1]=1;
    for(int i=2;i<N;i++)
    {
        if( !vis[i] )   pri[++cnt]=i,miu[i]=-1;
        for(int j=1;j<=cnt;j++)
        {
            int x=pri[j]*i;
            if( x>=N )   break ;
            vis[x]=1;
            if( i%pri[j] )   miu[x]=-miu[i];
            else   break ;
        }
    }
    for(int i=1;i<N;i++)   sum[i]=sum[i-1]+miu[i];
}
int cal(Lint n)
{
    if( n<N )   return sum[n];
    if( f.count(n) )   return f[n];
    Lint x=2;int ret=1;
    while( x<=n )
    {
        Lint y=n/(n/x);
        ret-=(y-x+1)*cal(n/x),x=y+1;
    }
    return f[n]=ret;
}
int main()
{
    Lint l,r;
    Prepare();
    scanf("%lld%lld",&l,&r);
    printf("%d\n",cal(r)-cal(l-1));
    return 0;
}