1. 程式人生 > >51nod 1040 求1-n這n個數,同n的最大公約數的和(尤拉函式)

51nod 1040 求1-n這n個數,同n的最大公約數的和(尤拉函式)

題目:給出一個n,求1-n這n個數,同n的最大公約數的和。比如:n = 6

          1,2,3,4,5,6 同6的最大公約數分別為1,2,3,2,1,6,加在一起 = 15

思路:一個數與n的最大公約數肯定是n的因子中的一個,所以只需要列舉n的每一個因子x,然後看有多少個滿足條件gcd(k,n)==x即gcd(k/x,n/x)==1的k就可以了,求這個可以用尤拉函式

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
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;
}
ll n;
int main()
{
    scanf("%lld",&n);
    ll ans = 0;
    for(ll i=1;i*i<=n;i++)
    if(n%i==0)
    {
        ans+=i*eular(n/i);
        if(i!=n/i)
            ans+=(n/i)*eular(i);
    }
    printf("%lld\n",ans);
}