南昌网络赛tsy's number(莫比乌斯反演&线性递推)

南昌网络赛tsy’s number(莫比乌斯反演&线性递推)

题目大意

给出n，求
$$\sum _{i=1}^n\sum _{j=1}^n\sum _{j=1}^n\frac{\phi (i)\phi (j^2)\phi (k^3)}{\phi (i)\phi (j)\phi (k)}\phi (gcd(i,j,k))$$

解题思路

据欧拉函数的性质，原式可以化作
$$\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j{k}^2\phi (gcd(i,j,k))$$
也就可以化作
$$\begin{gathered} \sum _{d=1}^n\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j{k}^2\phi (d)[gcd(i,j,k)==d] \\ \to \sum _{d=1}^n{d}^3\phi (d)\sum _{i=1}^{[\frac{n}{d}]}\sum _{j=1}^{[\frac{n}{d}]}\sum _{k=1}^{[\frac{n}{d}]}j{k}^2[gcd(i,j,k)==1] \end{gathered}$$
令一新的函数为
$$f(s)=\sum _{i=1}^{[\frac{n}{d}]}\sum _{j=1}^{[\frac{n}{d}]}\sum _{k=1}^{[\frac{n}{d}]}j{k}^2[gcd(i,j,k)==s]$$
再令一$F(s)={\sum }_{s|d}f(d)$

则$F(s)$可以写作
$$\begin{gathered} F(s)=\sum _{i=1}^{[\frac{n}{d}]}\sum _{j=1}^{[\frac{n}{d}]}\sum _{k=1}^{[\frac{n}{d}]}j{k}^2[s|gcd(i,j,k)] \\ \to =s^3\sum _{i=1}^{[\frac{n}{sd}]}\sum _{i=1}^{[\frac{n}{sd}]}\sum _{i=1}^{[\frac{n}{sd}]}j{k}^2[1|gcd(i,j,k)] \\ \to =s^3\sum _{i=1}^{[\frac{n}{sd}]}\sum _{i=1}^{[\frac{n}{sd}]}\sum _{i=1}^{[\frac{n}{sd}]}j{k}^2 \end{gathered}$$
据莫比乌斯反演$f(s)={\sum }_{s|d}\mu (\frac{d}{s})F(d)$

原式就可以化成
$$=\sum _{d=1}^n\phi (d)d^3\sum _{t=1}^{[\frac{n}{d}]}\mu (t)t^3\sum _{i=1}^{[\frac{n}{td}]}\sum _{i=1}^{[\frac{n}{td}]}\sum _{i=1}^{[\frac{n}{td}]}j{k}^2$$
再令$T=td$将式子转化就有
$$ans[n]=\sum _{T=1}^n{T}^3\sum _{d|T}\phi (d)\mu (\frac{T}{d})\sum _{i=1}^{[\frac{n}{td}]}\sum _{i=1}^{[\frac{n}{td}]}\sum _{i=1}^{[\frac{n}{td}]}j{k}^2$$
其中${\sum }_{i=1}^{[\frac{n}{td}]}{\sum }_{i=1}^{[\frac{n}{td}]}{\sum }_{i=1}^{[\frac{n}{td}]}j{k}^2$可以整除分块得到答案

而${\sum }_{d|T}\phi (d)\mu (\frac{T}{d})$可以通过线性递推得到

令$g(d)={\sum }_{i|d}\phi (i)\mu (\frac{d}{i})$

根据欧拉函数与莫比乌斯函数的性质，我们可以进行如下分类

当d为素数时易证$g(d)=d-2$

当d为合数时：

设d由$q\ast p$(p为素数,q为合数)递推来

$g(d)$可以写作$g(d)={\sum }_{i|q}\phi (ip)\mu (\frac{q}{i})+{\sum }_{i|q}\phi (i)\mu (\frac{qp}{i})$

设$A(d)={\sum }_{i|q}\phi (ip)\mu (\frac{q}{i})$,$B(d)={\sum }_{i|q}\phi (i)\mu (\frac{qp}{i})$

当q与p互素时：

由于欧拉函数为积性函数因此$A(d)=\phi (p)\ast A(q)$

由于增加了一个互素的素因子，因此所有的莫比乌斯函数取相反数$B(d)=-A(q)$

(或根据g(d)为欧拉函数与莫比乌斯函数的迪利克雷卷积，可以得其也为积性函数)

当p为q的因数时：

由于qp种有重复因数，根据莫比乌斯函数$B(d)=0$

设$q=o\ast P^m$

当$m==1$时
$$\begin{gathered} A(d)=\sum _{i|o}\phi (ip)\mu (\frac{op}{i})+\sum _{i|o}\phi (i{p}^2)\mu (\frac{o}{i}) \\ \to =-\sum _{i|o}\phi (ip)\mu (\frac{o}{i})+p\sum _{i|o}\phi (ip)\mu (\frac{o}{i}) \\ \to =(p-1)\ast \sum _{i|o}\phi (ip)\mu (\frac{o}{i}) \\ \to =(p-1)\ast phi(p)\ast \sum _{i|o}\phi (i)\mu (\frac{o}{i}) \\ \to =(p-1)\ast phi(p)\ast g[o] \end{gathered}$$

当$m\ge 2$时

由于欧拉函数的性质$\phi (qp)=\phi (q)\ast p$

因此$A(d)=p\ast g[q]$

由此线性递推出所有的答案即可

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int maxn=1e7+5;
const LL mod=1<<30;
bool check[maxn];
int prime[maxn];
int g[maxn];
int f[maxn];
int phi[maxn];
int tot=0;
int sum[maxn];
void init()
{
	for(int i=1;i<maxn;i++)
	{
		f[i]=i;
	}
	for(int i=1;i<maxn;i++)
	{
		g[i]=(1LL*g[i-1]+1LL*i)%mod;
		f[i]=(1LL*f[i]*g[i])%mod;
	}
	for(int i=1;i<maxn;i++)
	{
		g[i]=(1LL*g[i-1]+1LL*i*i)%mod;
		f[i]=(1LL*f[i]*g[i])%mod;
	}
	g[1]=phi[1]=1;
	for(int i=2;i<maxn;i++) check[i]=true;
	for(int i=2;i<maxn;i++)
	{
		if(check[i])
		{
			phi[i]=i-1;
			prime[tot++]=i;
			g[i]=i-2;//for mu(1)=1
		}
		for(int j=0;j<tot&&prime[j]*i<maxn;j++)
		{
			check[i*prime[j]]=false;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j];
				if((i/prime[j])%prime[j]==0) g[i*prime[j]]=g[i]*prime[j]%mod;
				else g[i*prime[j]]=(prime[j]-1)*phi[prime[j]]*g[i/prime[j]];
				break;
			}
			else
			{
				phi[i*prime[j]]=phi[i]*phi[prime[j]];
				g[i*prime[j]]=((-g[i]+1LL*phi[prime[j]]*g[i])%mod+mod)%mod;
			}
		}
	}
	for(int i=1;i<maxn;i++)
	{
		sum[i]=(sum[i-1]+((((1LL*g[i]*i%mod)*i)%mod)*i)%mod)%mod;
	}
}
int solve(int n)
{
	int ans=0;
	for(int l=1,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		ans=(ans+((((sum[r]-sum[l-1])%mod+mod)%mod)*f[n/l])%mod)%mod;
	}
	return ans;
}
int main()
{
	int t;
	init();
	scanf("%d",&t);
	while(t--)
	{
		int n;
		scanf("%d",&n);
		printf("%d\n",solve(n));
	}
}