spoj5971 lcmsum

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

考虑将$lcm$转化为$gcd$，原式变成$$\sum _{i=1}^n\frac{i\ast n}{gcd(i,n)}$$

蜜汁变形（$i$和$n-i$其实是等价的）：
$$\frac{{\sum }_{i=1}^{n-1}\frac{i\ast n}{gcd(i,n)}+{\sum }_{i=1}^{n-1}\frac{(n-i)\ast n}{gcd(n-i,n)}}{2}+n$$

然后发现：$gcd(i,n)=gcd(n-i,n)$

那么只需求：
$$\sum _{i=1}^{n-1}\frac{n^2}{gcd(i,n)}$$

枚举$gcd$（设枚举的$gcd$为$d$）：
$$\sum _{d|n}\sum _{i=1}^{n-1}\frac{n^2}{d}\ast (gcd(i,n)==d)$$

把$\frac{n^2}{d}$提出来，把d除下去，得到：
$$\sum _{d|n}\frac{n^2}{d}\ast \sum _{i=1}^{\lfloor \frac{n-1}{d}\rfloor }(gcd(i,\frac{n}{d})==1)$$

把$\frac{n}{d}$看做D。（其实$\frac{n}{d}$和d是等价的）变成：
$$n\ast \sum _{D|n}D\ast \sum _{i=1}^{D-1}(gcd(i,D)==1)$$

这里可以把$\lfloor \frac{n-1}{d}\rfloor$看做是$D-1$。那么${\sum }_{i=1}^{D-1}(gcd(i,D)==1)$就可以理解为小于D的与D互质的数的个数。这就是欧拉函数嘛。于是原式变为：
$$n\ast \sum _{D|n,D!=1}D\ast \varphi (D)$$

我们发现，只要得到每个数的约数就能算出答案了。但是对于每个数去求一次约数，由于有3e5个询问，时限显然不够。那么可以反过来考虑——考虑每一个数对它倍数的贡献。$D!=1$是因为如果$D=1$，那么$d=n$，那么$i=n$，由于把$i=n$排除了，所以要排除$D=1$的情况。

for(int i=1;i<maxn;++i)
		for(int j=2;((j<maxn)&&((i*j)<maxn));++j)
		ans[i*j]+=(ll)j*phi[j];
	for(int i=1;i<maxn;++i) ans[i]=(ll)(ans[i]/2)*i+i;

这里$j$代表的就是刚才的$D$，枚举$i$的所有倍数，把答案记在它的倍数上。然后代进之前推的完整的式子，就可以$O(1)$查询了。
下面是完整代码。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e6+10;
int phi[maxn],mark[maxn],P[maxn],cnt=0,T,x;
ll ans[maxn];
inline void pre(){
	phi[1]=1;
	for(int i=2;i<=maxn;++i){
		if(!mark[i]) P[++cnt]=i,phi[i]=i-1;
		for(int j=1;((j<=cnt)&&(P[j]*i<maxn));++j){
			mark[P[j]*i]=1;
			if(i%P[j]==0){
				phi[i*P[j]]=phi[i]*P[j];
				break;
			}phi[i*P[j]]=phi[i]*phi[P[j]];
		}
	}
	for(int i=1;i<maxn;++i)
		for(int j=2;((j<maxn)&&((i*j)<maxn));++j)
		ans[i*j]+=(ll)j*phi[j];
	for(int i=1;i<maxn;++i) ans[i]=(ll)(ans[i]/2)*i+i;
}
inline void print(ll x){
	if(x>9) print(x/10);
	putchar(x%10+'0');
}
int main(){
	pre(),scanf("%d",&T);
	while(T--)scanf("%d",&x),print(ans[x]),putchar(10);
}