P2398 GCD SUM （三种思路）

题目链接：点击这里

题目大意：
给定正整数 $n$ 求：
$$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)$$

题目分析：

• 方法一：
  用莫比乌斯反演化简式子(具体推导：点击这里)可以得到：
  $$\sum _{i=1}^n\sum _{j=1}^n\gcd (i,j)$$
  $$=\sum _{x=1}^n\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{x}\rfloor \phi (x)$$
  然后通过整除分块，线筛欧拉函数前缀和解决，时间复杂度为 $O(n)$

int cnt,pri[maxn],mu[maxn],phi[maxn],sum[maxn];
bool vis[maxn];
void get_phi()
{
	phi[1] = sum[1] = vis[1] = 1;
	for(int i = 2;i < maxn;i++)
	{
		if(!vis[i])
		{
			pri[++cnt] = i;
			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 phi[i*pri[j]] = phi[i]*(pri[j]-1);
		}
		sum[i] = sum[i-1]+phi[i];
	}
}
signed main()
{
	get_phi();
	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);
	return 0;
}

• 方法二：
  考虑枚举所有 $\gcd$ 值计算，令 $\gcd (i,j)=1$ ，则有 $\gcd (ik,jk)=k$ ，所以 $\gcd (i,j)=k$ 的个数为：
  $$2\ast \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\phi (i)-1$$
  减1是因为 $(x,y)$ 与 $(y,x)$ 均可对答案产生贡献，但 $(1,1)$ 只能对答案产生一次贡献
  线筛求一下欧拉函数前缀和即可，时间复杂度为 $O(n)$

int cnt,pri[maxn],mu[maxn],phi[maxn],sum[maxn];
bool vis[maxn];
void get_phi()
{
	phi[1] = sum[1] = vis[1] = 1;
	for(int i = 2;i < maxn;i++)
	{
		if(!vis[i])
		{
			pri[++cnt] = i;
			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 phi[i*pri[j]] = phi[i]*(pri[j]-1);
		}
		sum[i] = sum[i-1]+phi[i];
	}
}
signed main()
{
	get_phi();
	int n = read();
	ll ans = 0;
	for(int i = 1;i <= n;i++)
		ans += (sum[n/i]*2-1)*i;
	printf("%lld\n",ans);
	return 0;
}

• 方法三：
  设 $$f(n)=\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)=n]$$
  $$F(x)=\sum _{x|d}f(d)=\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{x}\rfloor$$
  $$=f(x)+f(2x)+...+f(kx)=\sum _{i=1}^{\lfloor \frac{n}{x}\rfloor }f(ix)$$
  所以有：
  $$f(x)=\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{x}\rfloor -\sum _{i=2}^{\lfloor \frac{n}{x}\rfloor }f(ix)$$
  时间复杂度为 $O(nlogn)$ （调和级数：$n(\frac{1}{1}+\frac{1}{2}+\frac{1}{3}+...+\frac{1}{n})\approx nlogn$）

ll f[maxn];
signed main()
{
	int n = read();
	ll ans = 0;
	for(int i = n;i;i--)
	{
		f[i] = n/i*(n/i);
		for(int j = 2*i;j <= n;j += i)
			f[i] -= f[j];
		ans += f[i]*i;
	}
	printf("%lld\n",ans);
	return 0;
}