tsy‘s number（莫比乌斯反演+狄利克雷卷积+欧拉筛+数论分块）

再化简公式前，先介绍一下强化版的欧拉函数积性公式
$\varphi (nm)=\varphi (n)\varphi (m)\frac{GCD(n,m)}{\varphi (GCD(n,m))}$

如下是详细的化简过程
原式 $=\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j{k}^2\varphi (gcd(i,j,k))$

$=\sum _{d=1}^n\varphi (d)\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j{k}^2[gcd(i,j,k)=d]$

$=\sum _{d=1}^n\varphi (d)d^3\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{d}\rfloor }j{k}^2[gcd(i,j,k)=1]$

设$f(x)=\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{d}\rfloor }j{k}^2[gcd(i,j,k)=x]$

又设$g(x)=\sum _{x|d}f(d)=\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{d}\rfloor }j{k}^2[x|gcd(i,j,k)]$

则$g(x)=x^3\sum _{i=1}^{\lfloor \frac{n}{xd}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{xd}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{xd}\rfloor }j{k}^2[1|gcd(i,j,k)]$

$=x^3\sum _{i=1}^{\lfloor \frac{n}{xd}\rfloor }1\sum _{j=1}^{\lfloor \frac{n}{xd}\rfloor }j\sum _{k=1}^{\lfloor \frac{n}{xd}\rfloor }{k}^2$

$=x^3\lfloor \frac{n}{xd}\rfloor \frac{(1+\lfloor \frac{n}{xd}\rfloor )\lfloor \frac{n}{xd}\rfloor }{2}\frac{(1+\lfloor \frac{n}{xd}\rfloor )\lfloor \frac{n}{xd}\rfloor (2\lfloor \frac{n}{xd}\rfloor +1)}{6}$

莫比乌斯反演：
$f(x)=\sum _{x|d}\mu (\frac{d}{x})\ast g(d)$

$f(1)=\sum _{1|d}\mu (\frac{d}{1})g(d)=\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\mu (d)g(d)$

原式$=\sum _{d=1}^n\varphi (d)d^3\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\mu (i)i^3\lfloor \frac{n}{id}\rfloor \frac{(1+\lfloor \frac{n}{id}\rfloor )\lfloor \frac{n}{id}\rfloor }{2}\frac{(1+\lfloor \frac{n}{id}\rfloor )\lfloor \frac{n}{id}\rfloor (2\lfloor \frac{n}{id}\rfloor +1)}{6}$

设T = id ，将 i 的求和变为 T 的求和：
$\sum _{T=1}^n[\lfloor \frac{n}{T}\rfloor \frac{(1+\lfloor \frac{n}{T}\rfloor )\lfloor \frac{n}{T}\rfloor }{2}\frac{(1+\lfloor \frac{n}{T}\rfloor )\lfloor \frac{n}{T}\rfloor (2\lfloor \frac{n}{T}\rfloor +1)}{6}\sum _{d|T}(\frac{T}{d}{)}^3\mu (d)d^3\varphi (\frac{T}{d})]$

先利用欧拉筛法求出右边的狄利克雷卷积，再用数论分块的方式求出整体的值。筛法的时间复杂度是o(n)，求和的时间复杂度是o(n^1/2)，总的时间复杂度是o(n)。

先讲一下线性筛的思路：
令$f(x)=\sum _{d|x}(\frac{x}{d}{)}^3\varphi (\frac{x}{d})\mu (d)d^3$

1. $x$为素数时，$f(x)=x^3(x-2)$
2. $x$为合数时
   1. $x$只含一个素因子$p$，$f(x)=f(\frac{x}{p}p)=f(\frac{x}{p})f(p)$
   2. $x$含两个素因子$p$
      1.$\frac{x}{p^2}=1$时， $f(x)=f(\frac{x}{p^2}{p}^2)=f(\frac{x}{p^2})f(p^2)=f(p^2)=p^6(p-1{)}^2$
      2. $\frac{x}{p^2}\ne 1$时，$f(x)=f(\frac{x}{p^2}{p}^2)=f(\frac{x}{p^2})f(p^2)$
   3. $x$含三个及以上的素因子$p$时，$f(x)=\sum _{d|x}(\frac{x}{d}{)}^3\varphi (\frac{x}{d})\mu (d)d^3=\sum _{d|x}(\frac{\frac{x}{p}p}{d}{)}^3\varphi (\frac{\frac{x}{p}p}{d})\mu (d)d^3=p^4\sum _{d|\frac{x}{p}}(\frac{\frac{x}{p}}{d}{)}^3\varphi (\frac{\frac{x}{p}}{d})\mu (d)d^3=p^{4}f(\frac{x}{p})$

注意：上述表达中，用函数表达的都是已经算过的，直接拿来用就行了

接下来再讲一下数论分块的思路
数论分块有一道经典的例题，计算$\sum _{i=1}^{n}k(\mathrm{m}\mathrm{o}\mathrm{d}i)$
以$n=k=10$为例，$k(\mathrm{m}\mathrm{o}\mathrm{d}i)$分布如下
$10532211111$
这个序列是非递增序列，且有很多等值区间，可以证明最多有$2\sqrt{k}$个不同的数。那么可以化加法为乘法，一次性求出等值区间的和，将时间复杂度从$o(n)$，降为$o(n^{\frac{1}{2}})$。

有了这个基础之后，就可以类似的思想对和式进行分块，可以降低时间复杂度。
注意：在模下作除法要谨慎，因为2、6与模不互素，求不了逆元，所以要先作除法再取模，因为这个wa了很多次555…

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const ll MAXN=1e7+10;
const ll MOD=1ll<<30;
ll FG[MAXN]={0,1},pri[MAXN],jud[MAXN],cnt,loop,n,sum;
ll QM(ll a,ll b){
    ll res=1,rem=a;
    while(b){
        if(b&1) res=(res*rem)%MOD;
        rem=(rem*rem)%MOD;
        b>>=1;
    }
    return res;
}
void euler(){
    memset(jud,0,sizeof(jud));
    cnt=0;
    for(ll i=2;i<MAXN;i++){
        if(!jud[i]){
            pri[cnt++]=i;
            FG[i]=QM(i,3)*(i-2)%MOD;
        }
        for(ll j=0;j<cnt&&i*pri[j]<MAXN;j++){
            jud[i*pri[j]]=1;
            if(i%pri[j]==0){
                if((i/pri[j])%pri[j]==0) FG[i*pri[j]]=FG[i]*QM(pri[j],4)%MOD;
                else if(i/pri[j]==1) FG[i*pri[j]]=QM(pri[j],6)*QM(pri[j]-1,2)%MOD;
                else FG[i*pri[j]]=FG[i/pri[j]]*FG[pri[j]*pri[j]]%MOD;
                break;
            }
            FG[i*pri[j]]=FG[i]*FG[pri[j]]%MOD;
        }
    }
    for(ll i=2;i<MAXN;i++) FG[i]=(FG[i]+FG[i-1])%MOD;
}
int main(){
    euler();
    scanf("%lld",&loop);
    while(loop--){
        scanf("%lld",&n);
        sum=0;
        for(ll i=1,r=0;i<=n;i=r+1){
            ll t=n/i;
            r=n/t;
            ll A=(1+t)*t/2;
            ll B=(2*t+1)%3?A/3%MOD*(2*t+1)%MOD:(2*t+1)/3%MOD*A%MOD;
            sum=(sum+(FG[r]-FG[i-1]+MOD)%MOD*t%MOD*(A%MOD)%MOD*B%MOD)%MOD;
        }
        printf("%lld\n",sum);
    }
    return 0;
}