51nod 最大公约数之和

给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6
1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起 = 15
收起
输入
1个数N(N <= 10^9)
输出
公约数之和
输入样例
6
输出样例
15

基本思路:与N的最大公约数,肯定是N的因子,求最大公约式之和,及求因子x对和的贡献(与N最大公约数为x的个数)及GCD(k,n)=x,k的个数,即GCD(k/x,n/x)=1,则求与n/x互质的数的个数,用欧拉函数求解即可。

#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll enlar(ll n)//欧拉函数求解与n互质的个数
{
	ll ans=n;
  for(ll i=2;i*i<=n;i++)
  {
  	if(n%i==0)//第一次找到的必定为素因子
  	{
  		ans/=i;//ans=ans*(1-1/i)可画简为ans=ans/i*(i-1),先除后乘避免数据溢出
  		ans*=(i-1);
  		while(n%i==0)//把该素因子全部约掉
  		  n/=i;
	  }
  }
 if(n>1)
 {
 	ans/=n;
 	ans*=(n-1);
 }
	
  return ans; 
}
int main()
{
	ll x;
	cin>>x;
    ll ans=0;
for(ll i=1;i*i<=x;i++)	
 {
 	 if(x%i==0)//同时进行两次运算,增加效率
 	  {
 	  	
 	  	 ans+=i*enlar(x/i);
 	  	  if(x!=i*i)
 	  	   ans+=x/i*enlar(i);
 	  	 
	   } 

 }
cout<<ans<<endl;

 } 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值