P1390 公约数的和（莫比乌斯反演）

题目传送门

题意：

给你一个$n(1\le n\le 2\ast 10^6)$，求$\sum _{i=1}^n\sum _{j=i+1}^{n}gcd(i,j)$。

思路：

稍微对这个式子改写一下，变成了$\frac{\sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)-\frac{(1+n)\ast n}{2}}{2}$。

我们只研究式子$\sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)$就行了。

$\begin{array}{rl}& \\ & \sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)\\ =& \sum _{g=1}^n\sum _{i=1}^n\sum _{j=1}^n[gcd(i,j)=g]\ast g\\ =& \sum _{g=1}^n\sum _{i=1}^{\lfloor \frac{n}{g}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{g}\rfloor }[gcd(i,j)=1]\ast g\\ =& \sum _{g=1}^n\sum _{i=1}^{\lfloor \frac{n}{g}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{g}\rfloor }\sum _{d|(i,j)}\mu (d)\ast g\\ =& \sum _{g=1}^n\sum _{d|(i,j)}\mu (d)\lfloor \frac{n}{gd}{\rfloor }^2\ast g\end{array}$

然后就可以整除分块求答案了。

$Code$

// Author : ACfunhsl
// Time : 2021/5/17 14:13:11
#define int long long
const int N = 2e6+50;
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;
bool ok[N];
int p[N],cnt=0,mu[N];
void euler()
{
	mu[1] = 1;
	for(int i=2; i<N; i++)
	{
		if(!ok[i])
		{
			p[++cnt] = i;
			mu[i] = -1;
		}
		for(int j=1; j<=cnt&&i*p[j]<N; j++)
		{
			ok[i*p[j]] = 1;
			if(i%p[j]==0) break;
			mu[i*p[j]] = -mu[i];
		}
	}
	for(int i=1; i<N; i++)
		mu[i] += mu[i-1];
}
int cal(int n,int k)
{
	int res = 0;
	n/=k;
	for(int l=1,r; l<=n; l=r+1)
	{
		r = n/(n/l);
		res += (n/l)*(n/l)*(mu[r] - mu[l-1]);
	}
	return res;
}
signed main()
{
	euler();
	int n;
	cin>>n;
	int res = 0;
	for(int i=1;i<=n;i++)
		res += cal(n,i)*i;
	res -= (1+n)*n/2;
	res /= 2;
	cout<<res<<endl;
	

	return 0;
}