【杜教筛】51Nod1244[莫比乌斯函数之和]题解

题目概述

求 ∑ni=1μ(i) 。

解题报告

杜教筛可以用来求积性函数的前缀和，具体想法是用另外一个函数卷待求函数，如下：

∑i=1n(f∗g)(i)=∑i=1n∑d|if(id)g(d)=∑d=1ng(d)∑i=1⌊nd⌋f(i)=∑i=1ng(i)S(⌊ni⌋)⇔g(1)S(n)=∑i=1n(f∗g)(i)−∑i=2ng(i)S(⌊ni⌋)⇔S(n)=∑ni=1(f∗g)(i)−∑ni=2g(i)S(⌊ni⌋)g(1)

我们发现这变成了一个递归的过程，由于 i 从 2 开始，所以如果我们能求出 g 的前缀和，就变成了除法分块！

对于莫比乌斯函数，我们可以用常数函数 1 卷 μ ，这样就变成了：

S(n)=∑i=1ne(i)−∑i=2nS(⌊ni⌋)=1−∑i=2nS(⌊ni⌋)

所以除法分块+记忆化就可以啦，效率是 O(n34) 。 证明？我不会啊！如果预处理前 O(n23) 个，复杂度就变成 O(n23) 了。 证明？我不会啊！

解题报告

#include
#include
using namespace std;
typedef long long LL;
const int maxn=4650000;

int p[maxn+5],mu[maxn+5];bool pri[maxn+5];
LL L,R;mapint> f;

void Make()
{
    pri[1]=true;mu[1]=1;
    for (int i=2;i<=maxn;i++)
    {
        if (!pri[i]) p[++p[0]]=i,mu[i]=-1;
        for (int j=1,t;j<=p[0]&&(t=i*p[j])<=maxn;j++)
        {
            pri[t]=true;mu[t]=-mu[i];
            if (!(i%p[j])) {mu[t]=0;continue;}
        }
    }
    for (int i=2;i<=maxn;i++) mu[i]+=mu[i-1];
}
LL Sum(LL n)
{
    if (n<=maxn) return mu[n];
    if (f.count(n)) return f[n];LL ans=1;
    for (LL l=2,r;l<=n;l=r+1)
        r=n/(n/l),ans-=Sum(n/l)*(r-l+1);
    return f[n]=ans;
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    scanf("%lld%lld",&L,&R);f.clear();Make();
    return printf("%lld\n",Sum(R)-Sum(L-1)),0;
}