51nod_1040:最大公约数之和

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

LL Eular(LL n)
{
    LL ret=n;
    for(LL i=2; i*i<= n; i++)
        if(n%i==0)
        {
            ret-=ret/i;
            while(n%i==0) n/= i;
        }
    if(n>1) ret-=ret/n;
    return ret;
}
LL cal(LL n)
{
    LL ret=0,i;
    for(i=1;i*i<n;i++)
           if(n%i==0)
        {
            ret=ret+(n/i)*Eular(i);
            ret=ret+i*Eular(n/i);
            }
    if(i*i==n) ret+=i*Eular(i);
    return ret;
}

int main()
{
    int n;
    while(cin>>n)
        cout<<cal(n)<<endl;
}