51NOD--1040 最大公约数之和--数论欧拉函数

给出一个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与i的最大公约数一定是n的约数x。枚举约数x，若gcd(n,i)为x，即gcd(n/x,i/x)为1。这里用欧拉定理处理。处理出与n/x互质的个数。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=50000+66;
const ll mod=1e9+7;

ll eular(ll n)
{
    ll ans=n;
    for(ll i=2;i*i<=n;i++)
    if(n%i==0)
    {
        ans-=ans/i;
        while(n%i==0)
            n/=i;
    }
    if(n>1) ans-=ans/n;
    return ans;
}
int main()
{
    int n;
    scanf("%d",&n);
    ll ans=0;
    for(int i=1;i*i<=n;i++)
    {
        if(n%i==0)
        {
            int o1=i;//好比2、3是6的倍数，i只到根号6,那么就要求出来6/2
            ans+=o1*eular(n/o1);
            if(i!=n/i)
            {
                int o2=n/i;
                ans+=o2*eular(n/o2);
            }
        }
    }
    printf("%lld\n",ans);
}