[BZOJ4804][积性函数][线性筛][莫比乌斯反演]欧拉心算

400题留念

推一推答案就是

$$\sum_{T=1}^n \lfloor\frac{n}{T}\rfloor^2f(T)$$
$$f(n)=\sum_{d|n}\phi(d)\mu(\dfrac{n}{d})$$

$f(n)$是个两个积性函数的卷积，同样也是积性函数，那么可以考虑用线性筛求

• 若n为素数 $f(n)=n-2$
• 若n是素数p的幂，$f(n)=\phi(n)-\phi(n/p)$
• 否则记录一下n的最小质因数p及次数c，则$f(n)=f(\dfrac{n}{p^c})*f(p_{now}*p^c)$

线性筛一下 再分块

#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 10000010

using namespace std;

typedef long long ll;

ll f[N];
int p[N],low[N];
int n,t;

inline void rea(int &x){
    char c=getchar(); x=0;
    for(;c>'9'||c<'0';c=getchar());
    for(;c>='0'&&c<='9';x=x*10+c-'0',c=getchar());
}

inline void First(){
    f[1]=1;
    for(int i=2;i<=N-10;i++){
        if(!p[i]) p[++*p]=i,f[i]=i-2,low[i]=i;
        for(int j=1;j<=*p&&(ll)i*p[j]<=N-10;j++)
            if(p[i*p[j]]=1,i%p[j]) f[i*p[j]]=f[i]*f[p[j]],low[i*p[j]]=p[j];
            else{
                low[i*p[j]]=low[i]*p[j];
                if(low[i]==i) f[i*p[j]]=(i*p[j]-i)-(i*p[j]-i)/p[j];
                else f[i*p[j]]=f[i/low[i]]*f[low[i]*p[j]];
                break;
            }
    }
    for(int i=2;i<=N-10;i++) f[i]+=f[i-1];
}

int main(){
    rea(t); First();
    while(t--){
        rea(n); ll Ans=0;
        for(int i=1,j=1;i<=n;i=j+1){
            j=n/(n/i);
            Ans+=(f[j]-f[i-1])*(n/i)*(n/i);
        }
        printf("%lld\n",Ans);
    }
}