Lcm

IV.Lcm

既然上一道题中的DZY都能自定义函数，那我们为什么不能呢？

定义$f(x)$为$x$中是否含有平方项。没有则为$1$，有则为$0$。显然，它是积性函数。而我们要求的，就是

$\sum_{i=1}^n\sum_{j=1}^m\dfrac{ij}{\gcd(i,j)}f(\gcd(i,j))$

都是老一套了。做多了就发现都是这个德行。

$\begin{aligned}&\sum_{i=1}^n\sum_{j=1}^m\dfrac{ij}{\gcd(i,j)}f(\gcd(i,j))\\=&\sum_{i=1}^n\sum_{j=1}^m\sum_{d=1}^{\min(n,m)}[\gcd(i,j)=d]\dfrac{ij}{d}f(d)\\=&\sum_{d=1}^{\min(n,m)}\sum_{i=1}^{n/d}\sum_{j=1}^{m/d}dij[\gcd(i,j)=1]f(d)\\=&\sum_{d=1}^{\min(n,m)}f(d)d\sum_{i=1}^{n/d}\sum_{j=1}^{m/d}ij[\gcd(i,j)=1]\\=&\sum_{d=1}^{\min(n,m)}f(d)d\sum_{i=1}^{n/d}\sum_{j=1}^{m/d}ij\sum_{x|i,x|j}\mu(x)\\=&\sum_{d=1}^{\min(n,m)}f(d)d\sum_{x=1}^{\min(n/d,m/d)}x^2\mu(x)\sum_{i=1}^{n/(dx)}\sum_{j=1}^{m/(dx)}ij\end{aligned}$

我们设$sum(x)=\dfrac{x(x+1)}{2}$。

$\begin{aligned}&\sum_{d=1}^{\min(n,m)}f(d)d\sum_{x=1}^{\min(n/d,m/d)}x^2\mu(x)\sum_{i=1}^{n/(dx)}\sum_{j=1}^{m/(dx)}ij\\=&\sum_{d=1}^{\min(n,m)}f(d)d\sum_{x=1}^{\min(n/d,m/d)}x^2\mu(x)sum(\dfrac{n}{dx})sum(\dfrac{m}{dx})\\=&\sum_{T=1}^{\min(n,m)}sum(\dfrac{n}{T})sum(\dfrac{m}{T})\sum_{d|T}f(d)d(\dfrac{T}{d})^2\mu(\dfrac{T}{d})\\=&\sum_{T=1}^{\min(n,m)}sum(\dfrac{n}{T})sum(\dfrac{m}{T})T\sum_{d|T}f(d)(\dfrac{T}{d})\mu(\dfrac{T}{d})\end{aligned}$

我们设$g(T)=T\sum_{d|T}f(d)(\dfrac{T}{d})\mu(\dfrac{T}{d})$。因为$f$是积性函数，$\mu(x)*id(x)$（注意这里是数乘不是卷积）也是积性函数，因此$g$是积性函数。

我们考虑线性筛$g$。

设$p$为一质数,则$g(p)=p*(f(1)*p*\mu(p)+f(p)*1*\mu(1))=p*(1-p)$

考虑如何求出$g(xp)$，其中$p$为一质数，$x$为任意数（保证大于$p$）。

当$\gcd(x,p)=1$，即$x,p$互质时，我们按照积性函数性质，有$g(xp)=g(x)*g(p)$。

否则，如果$xp$中含有三个及以上的因数$p$，则$f(xp)=0$，可以直接得出$g(xp)=0$。

否则，我们将$xp$分割成$p^2$和$\dfrac{x}{p}$两个必然互质的部分，运用积性函数性质，则有$g(xp)=g(\dfrac{x}{p})*g(p^2)$。

考虑如何求出$g(p^2)$。按照性质递推，会发现它等于$-p^3$。

当然咯，为了减少计算量，因为$g(T)=T\sum_{d|T}f(d)(\dfrac{T}{d})\mu(\dfrac{T}{d})$，因此这个因子$T$我们统一留到最后加上。这时，有$g(p)=(1-p),g(p^2)=-p$（复辟）。

然后就直接整除分块即可。

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=4e6;
int n,m,T,pri[N+5],g[N+5],res;
void init(){
	g[1]=1;
	for(int i=2;i<=N;i++){
		if(!pri[i])pri[++pri[0]]=i,g[i]=1-i;
		for(int j=1;j<=pri[0]&&i*pri[j]<=N;j++){
			pri[i*pri[j]]=true;
			if(!(i%pri[j])){
				int r=i/pri[j];
				if(r%pri[j])g[i*pri[j]]=-pri[j]*g[r];
				break;
			}
			g[i*pri[j]]=g[i]*g[pri[j]];
		}
	}
	for(int i=1;i<=N;i++)g[i]=g[i-1]+g[i]*i;
}
int main(){
	scanf("%d",&T),init();
	while(T--){
		scanf("%d%d",&n,&m),res=0;
		for(int l=1,r;l<=min(n,m);l=r+1)r=min(n/(n/l),m/(m/l)),res+=((1ll*(n/l)*(n/l+1)/2)*(1ll*(m/l)*(m/l+1)/2)*(g[r]-g[l-1]));
		printf("%d\n",res&(~((1<<30)|(1<<31))));
	}
	return 0;
}