51Nod 1040 最大公约数之和

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;
typedef long long LL;

LL Eular(LL m)
{
    LL res = m;
    for(LL i=2; i*i<=m; i++)
    {
        if(m%i==0)
        {
            res -= res/i;
            while(m%i==0)
                m /= i;
        }
    }
    if(m > 1)
        res -= res/m;
    return res;
}

int main()
{
    LL m;
    while(~scanf("%lld",&m))
    {
        LL ans = 0;
        for(LL i=1; i*i<=m; i++)
        {
            if(m%i == 0){
                ans += Eular(m/i)*i;
                if(m != i*i)
                    ans += m/i*Eular(i);
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}