poj 2480 Longge's problem(積分関数&Euler関数)

http://poj.org/problem?id=2480
Σgcd(i,N)1<=i<=Nを解く.
最大公約数については、gcd(n,m 1*m 2)=gcd(nm 1)*gcd(n,m 2)例えばgcd(6,12)=gcd(6,3)*gcd(6,4)の性質がある.
f(N)=Σgcd(i,N)とすると,それも積関数であるべきである.(積性関数の積は積性関数、積性関数の和は積性関数)
ここでpはnと互いに素を結ぶ数字とする.
また、積関数の性質を設定します.
      (*)
同時に、次のことを知っています.
オラ関数には、次のような計算方法があります.
結果:
持ち込み(*)中:
書き込みコード:

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;

int main()
{
    LL n,ans;
    while(cin>>n){
        int len=sqrt(n*1.0);
        ans=n;
        LL a=0,p=1;
        for(int i=2;i<=len;i++){
            a=0;
            if(n%i==0){
               a++;
               n/=i;
               p=i;
               while(n%i==0) {
                   a++;
                   n/=i;
               }
               ans=ans/p*(p+a*p-a);
            }
        }
        if(n>1){
            ans=ans/n*(n+n-1);
        }
        printf("%lld
",ans);
    }
    return 0;
}