洛谷P1390 公约数的和

最新推荐文章于 2022-10-14 11:58:11 发布

Cr.rech

最新推荐文章于 2022-10-14 11:58:11 发布

阅读量244

点赞数

CC 4.0 BY-SA版权

分类专栏：做题笔记

本文链接：https://blog.csdn.net/bcr_233/article/details/88879911

做题笔记专栏收录该内容

18 篇文章

订阅专栏

本文介绍了一种使用数论技巧求解特定数学问题的方法。通过枚举和递推公式，我们能够有效地计算∑∑gcd(i,j)从i=1到n，j=1到n的值。文章提供了一个O(nlogn)复杂度的C++代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

设 $f[d]=∑∑gcd⁡(i,j)=d,F[d]=∑∑d∣gcd⁡(i,j)f[d]=\sum\sum \gcd(i,j)=d,F[d]=\sum\sum d|\gcd(i,j)$
不难看出 $F[d]=(n/d)^2$
那么 $f [d] = F [d] - \sum f [k d]$

通过枚举 $d$ 的方式，我们可以求出 $∑i=1n∑j=1ngcd⁡(i,j)\sum\limits_{i=1}^n\sum\limits_{j=1}^n\gcd(i,j)$
然后稍作处理就是题目要求的答案了。

好像很简单的样子啊qwq
所以倒着 $O(nlog⁡n)\mathcal{O(n\log n)}$ 扫一遍就珂以了QAQ

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#define ll long long //不开ll见祖宗
using namespace std;

ll f[2000010];

int main() {
     ll n,ans=0;
     scanf("%lld",&n);
     for (int i=n; i>=1; --i) { //枚举d
     	f[i]=(n/i)*(n/i);
     	for (int j=(i<<1); j<=n; j+=i) { //枚举kd
     		f[i]-=f[j];
     	}
     	ans+=(ll)(f[i]*i);
     }
     ans=ans-n*(n+1)/2; //除去gcd(i,i)形式的结果
     ans/=2;  //除去重复部分
     printf("%lld\n",ans);
     return 0;
}