51nod 1238 最小公倍数之和 V3

背景：

$\text{luogu P1829 }$[国家集训队]Crash的数字表格 / JZPTAB的加强版。

题目传送门：

http://www.51nod.com/Challenge/Problem.html#!#problemId=1238

题意：

你可以认为以下的$n=m$。
求${\sum }_{i=1}^n{\sum }_{j=1}^m\text{lcm}(i,j)$。

思路：

$$\sum _{i=1}^n\sum _{j=1}^m\text{lcm}(i,j)$$

$$=\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{\gcd (i,j)}$$

你设一个$k$枚举$gcd(i,j)$。
$$=\sum _{k=1}^{\min (n,m)}\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{\gcd (i,j)}[\gcd (i,j)=k]$$

搞掉那一个$k$，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }{k}^2\frac{ij}{k}[\gcd (i,j)=1]$$

化简一下，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }ijk[\gcd (i,j)=1]$$

再套一个${\sum }_{d|n}{\mu }_d=[n=1]$，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }ijk\sum _{d|\gcd (i,j)}{\mu }_d$$

搞一下$d$的枚举方式，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }ijk\sum _{d=1}^{\min (\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{\mu }_d[d|\gcd (i,j)]$$

前置一下$d$，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{\mu }_d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }ijk[d|\gcd (i,j)]$$

更换一下枚举项，从而消除判断条件，得：
$$=\sum _{k=1}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{\mu }_d\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }ijk{d}^2$$

$k,d^2$置：
$$=\sum _{k=1}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{\mu }_d{d}^{2}k\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }ij$$

设$S={\sum }_{i=1}^n{\sum }_{i=1}^{n}ij$
因为${\sum }_{i=1}^n{\sum }_{i=1}^{n}ij=\frac{n(n+1)}{2}$
所以有：$S=\frac{n(n+1)}{2}$

则上式变为：
$$=\sum _{k=1}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{\mu }_d{d}^{2}k{S}_{\lfloor \frac{n}{kd}\rfloor }{S}_{\lfloor \frac{m}{kd}\rfloor }$$
现在就可以$\Theta (n\sqrt{n})$出解了。


你发现式子太长了，不和谐，于是用$T=kd$代替一下，同时$i=k,j=d$，得：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\sum _{ij=T}i{j}^2{\mu }_j$$

因为$ij=T$，所以有：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\sum _{ij=T}Tj{\mu }_j$$

$woc$后面的式子$i$不见了，那就等价于$j|T$，得：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\sum _{j|T}Tj{\mu }_j$$

等价于：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\cdot i{d}_T\cdot \sum _{j|T}i{d}_j\cdot {\mu }_j$$

由于一个数$\times 1$得原数，因此等价于：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\cdot i{d}_T\cdot \sum _{j|T}i{d}_j\cdot {\mu }_j\cdot I_{\frac{T}{j}}$$

等价于：
$$=\sum _{T=1}^{\min (n,m)}{S}_{\lfloor \frac{n}{T}\rfloor }{S}_{\lfloor \frac{m}{T}\rfloor }\cdot (id\cdot (id\cdot \mu \ast I))(T)$$


考虑$id\cdot (id\cdot \mu \ast I)$怎么杜教筛。
$$id\cdot ((id\cdot \mu )\ast I)$$

若$A$为完全积性函数，有$A\cdot (B\ast C)=A\cdot B\ast A\cdot C$，所以将$(id\cdot \mu )$看为一个整体，有：
$$=(id\cdot \mu \cdot id)\ast (I\cdot id)$$

$$=(i{d}^2\cdot \mu )\ast (I\cdot id)$$

我们选择$i{d}^2$进行狄利克雷卷积。

设$f=(i{d}^2\cdot \mu )\ast (I\cdot id)$，$g=i{d}^2$，$Sum=\sum f$，则有：
考虑$f\ast g=(i{d}^2\cdot \mu )\ast (I\cdot id)\ast i{d}^2=f\ast g=(i{d}^2\cdot \mu )\ast id\ast i{d}^2$。
因为狄利克雷卷积满足交换律，因此有：$f\ast g=(i{d}^2\cdot \mu )\ast i{d}^2\ast id$

$((i{d}^2\cdot \mu )\ast i{d}^2)(i)$
${\sum }_{j|i}{j}^2\cdot {\mu }_j\cdot (\frac{n}{j}{)}^2$
${\sum }_{j|i}{n}^2\cdot {\mu }_j$
$=n^2{\sum }_{j|i}\cdot {\mu }_j$
$=n^2[n=1]$
$=1$

即$((i{d}^2\cdot \mu )\ast i{d}^2)(i)=I$，带回去，得到：
$f\ast g=(i{d}^2\cdot \mu )\ast i{d}^2\ast id=I\ast id=id$

杜教筛是这么处理积性函数的：
$$g(1)\ast Sum(n)=\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)Sum(\lfloor \frac{n}{i}\rfloor )$$

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

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

现在可以整除分块了，自然就可以杜教筛了。


考虑如何线性筛。
$(f)(1)=(i{d}^2\cdot \mu )\ast id$
$=1$

设$p\in prime$。
$(f)(p)=(i{d}^2\cdot \mu )\ast id$
$={\sum }_{i|p}i{d}_i^2\cdot {\mu }_i\cdot i{d}_{\frac{p}{i}}$
$={\sum }_{i|p}{i}^2\cdot {\mu }_i\cdot \frac{p}{i}$
$={\sum }_{i|p}ip{\mu }_i$
$=p-p^2$
$=p(1-p)$

$(f)(p^k)=(i{d}^2\cdot \mu )\ast id$
$={\sum }_{i|p^k}i{d}_i^2\cdot {\mu }_i\cdot i{d}_{\frac{p}{i}}$
$={\sum }_{i|p^k}{i}^2\cdot {\mu }_i\cdot \frac{p^k}{i}$
$={\sum }_{i|p^k}i{p}^k{\mu }_i$
$=p^k-p^{k+1}$
$=p^k(1-p)$

可以$\Theta (1)$求，就可以线性筛了。

代码：

#include<cstdio>
#include<cstring>
#include<map>
#include<algorithm>
#define LL long long
#define MAXN 10000000
#define mod 1000000007
#define inv2 500000004
#define inv6 166666668
using namespace std;
	LL n;
	int t=0;
	int prime[MAXN],low[MAXN];
	LL f[MAXN];
	bool bz[MAXN];
void init(int ma)
{
	f[0]=0,f[1]=1;
	bz[0]=bz[1]=true;
	for(int i=2;i<=ma;i++)
	{
		if(!bz[i]) prime[++t]=low[i]=i,f[i]=((LL)i*(1-i)%mod+mod)%mod;
		for(int j=1;j<=t&&i*prime[j]<=ma;j++)
		{
			bz[i*prime[j]]=true;
			if(!(i%prime[j]))
			{
				low[i*prime[j]]=low[i]*prime[j];
				if(low[i]==i)
					f[i*prime[j]]=((((1ll-prime[j])*i%mod+mod)%mod)*prime[j]%mod+mod)%mod;
				else
					f[i*prime[j]]=(f[i/low[i]]*f[low[i]*prime[j]]%mod+mod)%mod;
				break;
			}
			f[i*prime[j]]=f[i]*f[prime[j]]%mod;
			low[i*prime[j]]=prime[j];
		}
	}
	for(int i=1;i<=ma;i++)
		f[i]=(f[i]+f[i-1]+mod)%mod;
}
LL calc_2(LL x)
{
	x%=mod;
	return x*(x+1)%mod*inv2%mod;
}
LL calc_6(LL x)
{
	x%=mod;
	return x*(x+1)%mod*(2*x+1)%mod*inv6%mod;
}
map<LL,LL> Sum;
LL calc(LL n)
{
	if(n<=MAXN-10) return f[n];
	{
		map<LL,LL>::iterator t=Sum.find(n);
		if(t!=Sum.end()) return t->second;
	}
	LL ans=calc_2(n);
	for(LL l=2,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		ans=(ans-(calc_6(r)-calc_6(l-1)+mod)%mod*calc(n/l)%mod+mod)%mod;
	}
	return Sum[n]=ans;
}
LL work()
{
	LL ans=0;
	for(LL l=1,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		ans=(ans+calc_2(n/l)*calc_2(n/l)%mod*(calc(r)-calc(l-1)+mod)%mod)%mod;
	}
	return ans;
}
LL dg(LL x,LL k)
{
	if(!k) return 1;
	LL op=dg(x,k>>1);
	if(k&1) return op*op%mod*x%mod; else return op*op%mod;
}
int main()
{
	scanf("%lld",&n);
	init(MAXN-10);
	printf("%lld",work());
}