给出一个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;
}