给出一个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
枚举因子,gcd(n,k)=x , 那么gcd( n/x , k/x)=1,即互质
所以 利用欧拉函数求解Euler(n/x) ,再利用公式 即可求解
欧拉函数
#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll Eular(ll n)
{
ll res=n;
for(ll i=2;i*i<=n;i++)
if(!(n%i))
{
res-=res/i;
while(!(n%i))
n/=i;
}
if(n>1)
res-=res/n;
return res;
}
int main()
{
ios::sync_with_stdio(false);
ll n,ans=0;
cin>>n;
for(ll i=1;i*i<=n;i++)
if(n%i==0)
{
ans+=i*Eular(n/i);
if(i*i!=n)
ans+=n/i*Eular(i);
}
cout<<ans<<endl;
return 0;
}