1. 程式人生 > >【杜教篩】51Nod1244[莫比烏斯函式之和]題解

【杜教篩】51Nod1244[莫比烏斯函式之和]題解

題目概述

ni=1μ(i)

解題報告

杜教篩可以用來求積性函式的字首和,具體想法是用另外一個函式卷待求函式,如下:

i=1n(fg)(i)=i=1nd|if(id)g(d)=d=1ng(d)i=1ndf(i)=i=1ng(i)S(ni)g(1)S(n)=i=1n(fg)(i)i=2ng(i)S(ni)S(n)=ni=1(fg)(i)ni=2g(i)S(ni)g(1)
我們發現這變成了一個遞迴的過程,由於 i2 開始,所以如果我們能求出 g 的字首和,就變成了除法分塊!

對於莫比烏斯函式,我們可以用常數函式 1μ

,這樣就變成了:

S(n)=i=1ne(i)i=2nS(ni)=1i=2nS(ni)
所以除法分塊+記憶化就可以啦,效率是 O(n34)證明?我不會啊!如果預處理前 O(n23) 個,複雜度就變成 O(n23) 了。證明?我不會啊!

解題報告

#include<cstdio>
#include<map>
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;map
<LL,int>
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; }