1040 最大公约数之和

1040 最大公约数之和

基准时间限制：1 秒 空间限制：131072 KB 分值: 80 难度：5级算法题

给出一个n，求1-n这n个数，同n的最大公约数的和。比如：n = 6
1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6，加在一起 = 15

Input
1个数N(N <= 10^9)

Output
公约数之和

Input示例
6

Output示例
15

与n的公约数，肯定是n的因子
那么我们枚举因子x就好了，显然这个因子的对答案的贡献就是gcd(n,i)=x的个数
gcd(n,i)=x的个数，就是gcd(n/x,i/x)=1的个数，那么就是求phi(n/x)啦

#include<stdio.h>
#include<iostream>
using namespace std;
#define ll long long
ll eular(ll n)
{
    int res=n,a=n;
    for(int i=2;i*i<=a;i++)
    {
        if(a%i==0)
        {
            res=res/i*(i-1);
            while(a%i==0)
            a/=i;
        }
    }if(a>1)res=res/a*(a-1);
    return res;
}
int main()
{
    ll n;
    ll sum=0;
    scanf("%d",&n);
    for(int i=1;i*i<=n;i++)
    {
        if(n%i==0)
        {
            sum+=1ll*i*eular(n/i);
            if(i*i!=n)
            sum+=1ll*n/i*eular(i);
        }
    }
    printf("%lld\n",sum);

    return 0;
}