bzoj2705[SDOI2012]Longge的问题 欧拉函数

题意：求$\sum^{i=1}_{n} gcd(i,n)$
这种题目居然也是省选题qwq..
考虑到n<=1e9，于是不可能直接暴力搞。
直接求出每一种约数的贡献然后直接乘起来就可以了
问题就转化成求每一个约数出现的次数，即$\sum gcd(i,j)=d$
继续转化就是gcd(i/d,j/d)=1,那么贡献为phi(max(i,j)/d).
但是时间复杂度很鬼畜啊，感觉根号套根号不久退化成O(n)了，为啥能过了。。。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
typedef long long ll;
ll n,m,ans;
ll phi(ll x)
{
    ll y=x;
    fo(i,2,m)
    if (x%i==0)
    {
        y=y/i*(i-1);
        while(x%i==0)x/=i;
    }
    if (x>1)y=y/x*(x-1);
    return y;
}
int main()
{
    scanf("%lld",&n);
    m=1ll*sqrt(n);
    fo(i,1,m)
    {
        if (n%i==0)
        {
            ans+=i*phi(n/i);
            if (i*i<n)ans+=(n/i)*phi(i);
        }
    }
    printf("%lld\n",ans);
}