【51nod 1238】最小公倍数之和 V3

传送门

---

problem

给出一个数 $n$，计算：

$$\sum _{i=1}^n\sum _{j=1}^n\mathrm{l}\mathrm{c}\mathrm{m}(i,j)$$

由于结果很大，输出答案模 $10^9+7$ 的结果。

数据范围：$n\le 10^{10}$。

---

solution

接下来的式子可能有点多，也有点复杂，建议在纸上好好推一下。

$\mathrm{l}\mathrm{c}\mathrm{m}$ 并不好处理，我们先把它化成 $\gcd$：

$$\begin{array}{rl}\sum _{i=1}^n\sum _{j=1}^n\mathrm{l}\mathrm{c}\mathrm{m}(i,j)& =\sum _{i=1}^n\sum _{j=1}^n\frac{ij}{\gcd (i,j)}\\ & =\sum _{k=1}^n\frac{1}{k}\sum _{i=1}^n\sum _{j=1}^{n}ij[\gcd (i,j)=k]\\ & =\sum _{k=1}^{n}k\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }ij[\gcd (i,j)=1]\end{array}$$

注意到枚举 $i,j$ 的上界都是 $\lfloor \frac{n}{k}\rfloor$，我们可以不用反演，有一种更简单的方法，即讨论一下 $i<j$，$i=j$，$i>j$ 三种情况（其中 $i>j$ 和 $i<j$ 可以合并）：

$$\begin{array}{rl}ans& =\sum _{k=1}^{n}k(2\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{i}ij[\gcd (i,j)=1]-\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }{i}^2[\gcd (i,i)=1])\\ & =\sum _{k=1}^{n}k(2\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }i(\sum _{j=1}^{i}j[\gcd (i,j)=1])-1)\end{array}$$

在进行下一步化简之前，我们先来证明一个式子： $\sum _{i=1}^{n}i[\gcd (i,n)=1]=\frac{n\phi (n)}{2}$。
根据辗转相减法，我们有 $\gcd (i,n)=\gcd (n-i,n)$，这就意味着，当 $\gcd (i,n)=1$ 时，必然有 $\gcd (n-i,n)=1$。因此，当答案要加上 $i$ 的贡献时，必然也会有 $n-i$ 的贡献，而 $i+(n-i)=n$。
又发现 $\sum _{i=1}^n[\gcd (i,n)=1]=\phi (n)$，所以原式 $=\frac{n\phi (n)}{2}$。
注意一下，当 $n=1$ 时，原式 $=1$。
证毕。

所以上式可以继续化简（那个 $[i=1]$ 就是特判 $i=1$ 的情况）：

$$\begin{array}{rl}ans& =\sum _{k=1}^{n}k(2\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }i(\frac{i\phi (i)+[i=1]}{2})-1)\\ & =\sum _{k=1}^{n}k\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }{i}^2\phi (i)\end{array}$$

于是我们令 $f(n)=n^2\phi (n)$（容易发现它是个积性函数），设 $F(n)$ 为它的前缀和。

那么答案就是：

$$ans=\sum _{k=1}^{n}k\times F(\lfloor \frac{n}{k}\rfloor )$$

这个显然可以整除分块，接下来就是考虑如何计算 $F(n)$。

考虑给 $f$ 卷上个 $I{d}^2$（$I{d}^2(n)=n^2$），得到：

$$\begin{array}{rl}(f\ast I{d}^2)(n)& =\sum _{d|n}f(d)I{d}^2(\frac{n}{d})\\ & =\sum _{d|n}{d}^2\phi (d)\frac{n^2}{d^2}\\ & =n^2\sum _{d|n}\phi (d)\\ & =n^3\end{array}$$

那么肯定也有 ${\sum }_{d|n}I{d}^2(d)f(\frac{n}{d})=n^3$，所以：

$$\begin{array}{rl}\sum _{i=1}{i}^3& =\sum _{i=1}^i\sum _{d|i}I{d}^2(d)f(\frac{i}{d})\\ & =\sum _{d=1}^n{d}^2\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)\\ & =\sum _{d=1}^n{d}^2\times F(\lfloor \frac{n}{d}\rfloor )\end{array}$$

又因为 $\sum _{i=1}^n{i}^3=\frac{n^2(n+1{)}^2}{4}$，所以：

$$F(n)=\frac{n^2(n+1{)}^2}{4}-\sum _{d=2}^n{d}^2\times F(\lfloor \frac{n}{d}\rfloor )$$

然后又一个整除分块就可以了。

时间复杂度 $O(n^{\frac{2}{3}})$。

---

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<unordered_map>
#define N 5000005
#define P 1000000007
#define ll long long
using namespace std;
ll n,inv2,inv4,inv6,f[N];
int sum,prime[N],mark[N],phi[N];
ll power(ll a,ll b,ll ans=1){
	for(;b;b>>=1,a=a*a%P)
		if(b&1)  ans=ans*a%P;
	return ans;
}
void linear_sieves(){
	f[1]=1;
	for(int i=2;i<N;++i){
		if(!mark[i])  prime[++sum]=i,phi[i]=i-1;
		for(int j=1;j<=sum&&i*prime[j]<N;++j){
			mark[i*prime[j]]=1;
			if(i%prime[j]==0)  {phi[i*prime[j]]=phi[i]*prime[j];break;}
			phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
		f[i]=(ll)i*i%P*phi[i]%P,f[i]=(f[i]+f[i-1])%P;
	}
	inv2=power(2,P-2),inv4=power(4,P-2),inv6=power(6,P-2);
}
ll Squ(ll x)  {return x%P*x%P;}
ll sum2(ll n)  {return n%P*((n+1)%P)%P*((2*n%P+1)%P)%P*inv6%P;}
ll sum3(ll n)  {return Squ(n%P)*Squ((n+1)%P)%P*inv4%P;}
unordered_map<ll,ll>Sum;
ll F(ll n){
	if(n<N)  return f[n];
	if(Sum[n])  return Sum[n];
	ll ans=sum3(n);
	for(ll i=2,j;i<=n;i=j+1){
		j=n/(n/i);
		ans=(ans-((sum2(j)-sum2(i-1)+P)%P*F(n/i)%P)+P)%P;
	}
	return Sum[n]=ans;
}
int main(){
	linear_sieves();
	scanf("%lld",&n);
	int ans=0;
	for(ll i=1,j;i<=n;i=j+1){
		j=n/(n/i);
		ans=(ans+((i+j)%P*((j-i+1)%P)%P*inv2%P)*F(n/i)%P)%P;
	}
	printf("%d",ans);
	return 0;
}