BZOJ 2705 [SDOI2012]Longge的问题

Description
Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。

---

【题目分析】
因为数据范围过大。不能存下来，只能使用根号n的算法每次都要算。

---

【代码】

#include <cstdio>
#include <cmath>
using namespace std;
long long m,t,n,ans;
long long phi(long long x)
{
    t=x;
    for (long long i=2;i<=m;++i)
    if (x%i==0){t=t/i*(i-1);while (x%i==0) x/=i;}
    if (x>1) t=t/x*(x-1);
    return t;
}
int main()
{
    scanf("%lld",&n);
    m=sqrt(n);
    for (long long i=1;i<=m;++i)
    if (n%i==0)
    {
        ans+=i*phi(n/i);
        if (i*i<n) ans+=(n/i)*phi(i);
    }
    printf("%lld",ans);
}