洛谷 P2398 GCD SUM

【题目链接】

洛谷 P2398 GCD SUM

【题目考点】

1. 数论：欧拉函数

欧拉函数$\varphi (x)$表示小于等于$x$的与$x$互质的数的个数。

欧拉函数具有如下性质：
$n={\sum }_{d|n}\varphi (d)$，即正整数n等于n的所有约数的欧拉函数值的加和。

2. 线性筛：求1到n的欧拉函数值

已知：对于质数$p$
如果$p\mid x$，则$\varphi (px)=p\varphi (x)$
如果$p\nmid x$，则$\varphi (px)=(p-1)\varphi (x)$

在线性筛求质数时，i是当前访问的数值，prime[j]是当前遍历到的质数，phi[i]表示欧拉函数$\varphi (i)$。
当确定i*prime[j]的最小质因数为prime[j]时，判断
如果i%prime[j] == 0，则phi[i*prime[j]] = prime[j]*phi[i]
否则如果i%prime[j] != 0，则phi[i*prime[j]] = (prime[j]-1)*phi[i]
递推可以求出1到n的欧拉函数。

【解题思路】

本题求：
$$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)$$
根据欧拉函数的性质$n=\sum _{d|n}\varphi (d)$，有：
$gcd(i,j)=\sum _{d|gcd(i,j)}\varphi (d)$
已知：$d|a,d|b$等价于$d|gcd(a,b)$，即公约数整除最大公约数。
那么：
$gcd(i,j)=\sum _{d|i,d|j}\varphi (d)$
$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)=\sum _{i=1}^n\sum _{j=1}^n\sum _{d|i,d|j}\varphi (d)$

即：$i、j$分别枚举区间$[1,n]$的值，如果$d$是$i$与$j$的公约数，那么对结果贡献为$\varphi (d)$
由于$i、j$能取到的最大值为$n$，所以$i、j$的公约数小于等于$n$，因此$i、j$的公约数$d$的取值范围为$[1,n]$。
因此可以在$[1,n]$范围内枚举$d$的值，对于每个$d$的取值考察其对结果的贡献。
只有当$d|i,d|j$，即$i$是$d$的倍数且$j$是$d$的倍数时，才有贡献$\varphi (d)$。
$i$的取值范围为$[1,n]$，其中$d$的倍数有$d,2d,...,\lfloor \frac{n}{d}\rfloor d$，共有$\lfloor \frac{n}{d}\rfloor$ 个$d$的倍数。
$j$的取值范围为$[1,n]$，同样其中共有$\lfloor \frac{n}{d}\rfloor$ 个$d$的倍数。
$i$取$d$的倍数的情况有$\lfloor \frac{n}{d}\rfloor$种，$j$取$d$的倍数的情况有$\lfloor \frac{n}{d}\rfloor$种，根据乘法原理，$i、j$同时取$d$的倍数的情况有$\lfloor \frac{n}{d}{\rfloor }^2$种。
因此，当$i、j$的公约数为$d$时，对结果的贡献为$\lfloor \frac{n}{d}{\rfloor }^2\varphi (d)$
因此：$$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)=\sum _{d=1}^n\lfloor \frac{n}{d}{\rfloor }^2\varphi (d)$$

在$[1,n]$范围内枚举$d$的值，将$\lfloor \frac{n}{d}{\rfloor }^2\varphi (d)$加和即为结果。

该计算过程需要用到1到n的欧拉函数值，可以使用线性筛预处理求出$\varphi (1)\sim \varphi (n)$的值。
该算法时间复杂度为$O(n)$。

【题解代码】

#include <bits/stdc++.h>
using namespace std;
#define N 100005
bool isPrime[N];
int prime[N], pn, phi[N];
long long n, ans;
void eulerSeive(int n)
{
	memset(isPrime, 1, sizeof(isPrime));
	phi[1] = 1;
	for(int i = 2; i <= n; ++i)
	{
		if(isPrime[i])
		{
			prime[++pn] = i;
			phi[i] = i-1;
		}
		for(int j = 1; j <= pn && i*prime[j] <= n; ++j)
		{
			isPrime[i*prime[j]] = false;
			if(i%prime[j] == 0)
			{
				phi[i*prime[j]] = prime[j]*phi[i];
				break;
			}
			else
				phi[i*prime[j]] = (prime[j]-1)*phi[i];
		}
	}
}
int main()
{
	cin >> n;
	eulerSeive(n);
	for(int d = 1; d <= n; ++d) 
		ans += (n/d)*(n/d)*phi[d];	
	cout << ans;
	return 0;
}