BZOJ2705: [SDOI2012]Longge的问题

Description
Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i
<=N)。
Input
一个整数，为N。
Output
一个整数，为所求的答案。
Sample Input
6
Sample Output
15
HINT
【数据范围】
对于60%的数据，$0<N<=2^{16}$。
对于100%的数据，$0<N<=2^{32}$。
Source
round1 day1

枚举$n$的约数$k$，设$s(k)$为满足$gcd(m,n)=k(1≤m≤n)$的个数，则$ans=\sum(k*s(k))$。因为$gcd(m,n)=k$,所以$gcd(m/k,n/k)=1$,于是$s(k)=φ(n/k)$

#include<cstdio>
#include<iostream>
#include<math.h>
using namespace std;
typedef long long ll;
ll n,ans;
void read(ll &x)
{
    char t=getchar();int f=1;x=0;
    while((t<48)or(t>57)){if(t=='-')f=-1;t=getchar();}
    while((t>=48)and(t<=57)){x=(x*10)+t-48;t=getchar();}
    x*=f;
}
ll phi(ll x)
{
    ll ans=1;
    for (int i=2;i*i<=x;++i)
    if (x%i==0)
    {
        ans*=(i-1);
        x/=i;
        while (x%i==0)
        {
            ans*=i;
            x/=i;
        }
    }
    if (x>1) ans*=(x-1);
    return ans;
}
int main()
{
    read(n);
    for (int i=1;i*i<=n;++i)
    if (n%i==0)
    {
        ans+=i*phi(n/i);
        if (i*i!=n) ans+=(n/i)*phi(i);//对应的另一约数
    }
    cout<<ans;
    return 0;
}