P2398 GCD SUM

最新推荐文章于 2021-06-12 17:12:03 发布

guapi2333

最新推荐文章于 2021-06-12 17:12:03 发布

阅读量108

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/guapi2333/article/details/115706578

版权

数论专栏收录该内容

25 篇文章 2 订阅

订阅专栏

Label

欧拉函数

Description

给定整数 $n(1\leq n\leq 10^5)$ ，求：

$\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)$

Solution

看上去此题所求又是一个反演的形式（其实此题可用反演解），但我们可以利用欧拉函数来分析。

有了poj2480的基础，我们已经知道： $\sum_{i=1}^{n}gcd(i,n)=\sum_{d|n}d\varphi(\frac{n}{d})$ ，则易得：

$\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)=2\sum_{i=1}^{j}\sum_{j=1}^{n}gcd(i,j)-\sum_{i=1}^{n}i=2\sum_{i=1}^{n}\sum_{d|i}d\varphi(\frac{i}{d})-\sum_{i=1}^{n}i$

而 $\sum_{i=1}^{n}\sum_{d|i}d\varphi(\frac{i}{d})=\sum_{d=1}^{n}\sum_{i=1}^{n}\varphi(\frac{i}{d})=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\varphi(i)$

（此步化简：考虑所有可能的 $d$ 的取值，进而考虑每种 $d$ 对答案的贡献）

所以， $\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)=2(\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\varphi(i))-\frac{n(n+1)}{2}$

由于 $\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\varphi(i)$ 可前缀和预处理，故算法整体复杂度为 $O (n)$ 。

Code

#include<cstdio>
#include<iostream>
#define ri register int
#define ll long long
using namespace std;

const int MAXN=1e5+20;
int N,prime[MAXN>>2],phi[MAXN],cnt;
ll ans,sum[MAXN];
bool isprime[MAXN];

void Phi()
{
	isprime[1]=true; phi[1]=1;
	for(ri i=2;i<=N;++i)
	{
		if(!isprime[i])
		{
			prime[++cnt]=i;
			phi[i]=i-1;	
		}
		for(ri j=1;j<=cnt&&i*prime[j]<=N;++j)
		{
			isprime[i*prime[j]]=true;
			if(i%prime[j]==0)	phi[i*prime[j]]=phi[i]*prime[j];
			else	phi[i*prime[j]]=phi[i]*phi[prime[j]];
			if(i%prime[j]==0)	break;
		}
	}
}

int main()
{
	scanf("%d",&N);
	Phi();
	for(ri i=1;i<=N;++i)	sum[i]=sum[i-1]+(ll)phi[i];
	for(ri d=1;d<=N;++d)	ans+=(ll)d*sum[N/d];
	ans=ans*2LL-(ll)N*(ll)(N+1)/2LL;
	cout<<ans;
	return 0;
}