BZOJ2705Longge的问题

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

phi可以在根号时间内求出。。
1A。。
TM还写了快速读入。。
我是不是闲的。。
附上本蒟蒻的代码：

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
long long n,ans;
int m;

long long read()
{
    long long w=0,c=1; char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

long long phi(long long x)
{
    long long t=x,i;
    for (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()
{
    int i;
    n=read(),m=sqrt(n);
    for (i=1;i<=m;i++)
      if (n%i==0)
        {
            ans+=(long long)i*phi(n/i);
            if (i*i<n) ans+=(long long)(n/i)*phi(i);
        }
    printf("%lld",ans);
    return 0;
}