洛谷P1390 公约数的和

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]=gcd(i,j)=d,F[d]=dgcd(i,j)
不难看出 F [ d ] = ( n / d ) 2 F[d]=(n/d)^2 F[d]=(n/d)2
那么 f [ d ] = F [ d ] − ∑ f [ k d ] f[d]=F[d]-∑f[kd] f[d]=F[d]f[kd]

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

好像很简单的样子啊qwq
所以倒着 O ( n log ⁡ n ) \mathcal{O(n\log n)} O(nlogn) 扫一遍就珂以了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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值