51Nod2026 Gcd and Lcm-题解

【题目地址】

题意简述

我们令$f(n)={\sum }_{d|n}\mu (d)d$，求下面式子在模$998244353$后的值：

$$\sum _{i=1}^n\sum _{j=1}^{n}f(gcd(i,j))\times f(lcm(i,j))$$

其中$gcd$为最大公约数，$lcm$为最小公倍数。

---

我们先来看一个例子：

求：
$$\sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)\times lcm(i,j)$$

很容易知道，这个可以这样变换：

$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)\frac{ij}{gcd(i,j)} \\ =\sum _{i=1}^n\sum _{j=1}^{n}ij \\ =\left(\sum _{i=1}^{n}i\right)\left(\sum _{j=1}^{n}j\right) \\ ={\left(\sum _{i=1}^{n}i\right)}^2 \end{gathered}$$

所以我们看，原来的式子是否也能这样变换。

首先我们可以知道$f=(\mu \cdot id)\ast 1$（其中$\ast$为狄利克雷卷积）。

由于积性函数的乘积以及狄利克雷卷积仍旧为积性函数，所以$f$也为积性函数。

那么我们考虑，对于一对$i,j$，我们根据唯一分解定理，假设：

$$\begin{gathered} i=\prod p_x^{a_x} \\ j=\prod p_y^{b_y} \end{gathered}$$

那么$f(gcd(i,j))\times f(lcm(i,j))$可以写成：

$$\prod f(p_w^{min(a_w,b_w)})\prod f(p_w^{max(a_w,b_w)})$$
（其中是因为积性函数$f(xy)=f(x)f(y)$，在$x\perp y$时（互质），所以将$i,j$拆开）

那么我们将属于$i$的放在一起，属于$j$的放在一起，可以得到：

$$\begin{gathered} \prod f(p_x^{a_x})\prod f(p_y^{a_y}) \\ =f(i)f(j) \end{gathered}$$

所以就可以按照前面的例子的做法，将其分开变换后可以得到，原式等于：

$${\left(\sum _{i=1}^{n}f(i)\right)}^2$$

那么现在我们考虑如何求$f$的前缀和，由于$n\le 10^9$，所以要用杜教筛。

我们看，$f=(\mu \cdot id)\ast 1$，那么我们令$g=id$，则$f\ast g$就等于：

$$\begin{array}{rl}f\ast g=& (\mu \cdot id)\ast 1\ast id\\ =& (\mu \cdot id)\ast id\ast 1\\ =& \left(\sum _{d|n}\mu (d)d\times \frac{n}{d}\right)\ast 1\\ =& \left(n\sum _{d|n}\mu (d)\right)\ast 1\\ =& (id\cdot ϵ)\ast 1\\ =& ϵ\ast 1\\ =& 1\end{array}$$

那么我们将其带入杜教筛的公式，可以得到前缀和等于：

$$S(n)=\sum _{i=1}^{n}1-\sum _{i=2}^{n}iS(\lfloor \frac{n}{i}\rfloor )$$

我们预处理$n^{\frac{2}{3}}$的前缀和，然后杜教筛就好啦。

关于$f$的线性筛，方法如下：

因为当$\mu (p^c)$中$c>1$时$\mu (p^c)=0$，所以原式为：

$$f(p^c)=\sum _{i=0}^c\mu (p^i)p^i=1-p$$
根据积性函数定义：
$$\begin{gathered} f(p^c)=1-p \\ f(x\times p^c)=f(p^c)f(x)[x\perp p^c] \end{gathered}$$

所以这个题就可以在$O(n^{\frac{2}{3}})$时间内解决了。

#include<map>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int M=1e6+10;
const ll Mod=1e9+7;
ll prime[M],mu[M],cnt;
bool vis[M];
ll F[M],c[M],f[M],inv_2;
ll fpow(ll a,ll b){
	ll res=1;
	for(;b;b>>=1,a=(a*a)%Mod){
		if(b&1)res=(res*a)%Mod;
	}
	return res;
}
void init(){
	inv_2=fpow(2,Mod-2);
	F[1]=1;
	for(ll i=2;i<M;i++){
		if(!vis[i]){
			prime[++cnt]=i;
			F[i]=((1-i)%Mod+Mod)%Mod;
			c[i]=i;f[i]=i;
		}
		for(ll j=1,v;j<=cnt&&i*prime[j]<M;j++){
			v=i*prime[j];
			vis[v]=1;
			if(!(i%prime[j])){
				c[v]=c[i]*f[i];f[v]=f[i];
				F[v]=F[v/c[v]]*F[prime[j]]%Mod;
				break;
			}
			F[v]=F[i]*F[prime[j]]%Mod;
			c[v]=prime[j];f[v]=prime[j];
		}
	}
	for(ll i=2;i<M;i++){
		F[i]=(F[i-1]+F[i])%Mod;
	}
}
ll S(ll x){return ((x*(x+1))%Mod)*inv_2%Mod;}
ll Sarea(ll L,ll R){return ((S(R)-S(L-1))%Mod+Mod)%Mod;}
map <ll,ll> mp;
ll calc(ll x){
	if(x<M) return F[x];
	if(mp.count(x)) return mp[x];
	ll ans=x;
	for(ll i=2,j;i<=x;i=j+1){
		j=(x/(x/i));
		ans=(ans-(Sarea(i,j)*calc(x/i))%Mod)%Mod;
	}
	if(ans<0)ans+=Mod;
	return mp[x]=ans;
}
ll n;
ll Sqr(ll a){return a*a%Mod;}
int main(){
	scanf("%lld",&n);
	init();
	printf("%lld\n",Sqr(calc(n)));
	return 0;
}

---

End