P1390 公约数的和 （莫比乌斯反演+欧拉函数）

题目链接：点击这里

写在前面：

先在这里记录一下莫比乌斯函数的三条性质和一个式子：
1.对于给定的正整数 $n$ 有：
$$\sum _{d|n}\mu (d)=[n=1]$$
2.
$$\sum _{d|n}\mu (\frac{n}{d})d=\phi (n)$$
3.对于给定的正整数 $n$ 有：
$$\sum _{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$$
4.
$$\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)=1]$$
$$=\sum _{i=1}^n\sum _{j=1}^n\sum _{d|\gcd (i,j)}\mu (d)$$
$$=\sum _{d=1}^n\mu (d)\sum _{d|i}\sum _{d|j}$$
$$=\sum _{d=1}^n\mu (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{n}{d}\rfloor$$
题目大意：

给定 $n$ ，求
$$\sum _{i=1}^n\sum _{j=i+1}^n\gcd (i,j)$$
题目分析：
$$\sum _{i=1}^n\sum _{j=i+1}^n\gcd (i,j)$$
$$=\frac{{\sum }_{i=1}^n{\sum }_{j=1}^n\gcd (i,j)-{\sum }_{i=1}^n\gcd (i,i)}{2}$$
其中 ${\sum }_{i=1}^n\gcd (i,i)=\frac{n(n+1)}{2}$
$$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)$$
$$=\sum _{k=1}^{n}k\sum _{i=1}^n\sum _{j=1}^n[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 }\sum _{d|\gcd (i,j)}\mu (d)$$
$$=\sum _{k=1}^{n}k\sum _{d=1}^n\mu (d)\lfloor \frac{n}{kd}\rfloor \lfloor \frac{n}{kd}\rfloor$$
令 $x=kd$
$$=\sum _{x=1}^n\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{x}\rfloor \sum _{k|x}k\mu (\frac{x}{k})$$
利用性质2：
$$=\sum _{x=1}^n\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{x}\rfloor \phi (x)$$
此时就可以通过整除分块求 $\lfloor \frac{n}{x}\rfloor$ ，线筛求 $\phi (x)$ 及其前缀和来求出这个式子的值了
时间复杂度：$O(n)$

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define ll long long
#define int long long
#define inf 0x3f3f3f3f
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 2e6+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int cnt,pri[maxn],mu[maxn],phi[maxn],sum[maxn];
bool vis[maxn];
void get_mu()
{
	mu[1] = phi[1] = sum[1] = vis[1] = 1;
	for(int i = 2;i < maxn;i++)
	{
		if(!vis[i])
		{
			pri[++cnt] = i;
			mu[i] = -1;
			phi[i] = i-1;
		}
		for(int j = 1;j <= cnt && i*pri[j] < maxn;j++)
		{
			vis[i*pri[j]] = true;
			if(i%pri[j] == 0) 
			{
				phi[i*pri[j]] = phi[i]*pri[j];
				break;
			}
			else {
				mu[i*pri[j]] -= mu[i];
				phi[i*pri[j]] = phi[i]*(pri[j]-1);
			}
		}
		sum[i] = sum[i-1]+phi[i];
	}
}
signed main()
{
	get_mu();
	int n = read();
	ll ans = 0;
	for(int l = 1,r;l <= n;l = r+1)
	{
		r = n/(n/l);
		ans += 1ll*(sum[r]-sum[l-1])*(n/l)*(n/l);
	}
	printf("%lld\n",(ans-1ll*n*(n+1)/2)/2);
	return 0;
}