bzoj 2705: [SDOI2012]Longge的问题【欧拉函数+gcd】

bzoj 2705: [SDOI2012]Longge的问题
题意：给你一个整数n，求 $\sum_{i=1}^{n}gcd(i,n)$ .
分析：
$k = gcd(i, n) => 1 = gcd(i / k, n / k) => euler(N)_{N = n / k} * k$，所以$ans = \sum_{d|n}\phi(n/d) * d$.
折半枚举最大公约数k，$ans += i * euler(n / i) ;$

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

LL euler(LL x) {
    LL ans = x;
    for(LL i = 2; i * i <= x; i++) {
        if(x % i) continue;
        ans = ans / i * (i - 1);
        while(x % i == 0) x /= i;
    }
    if(x > 1) ans = ans / x * (x - 1);
    return ans;
}

int main() {
    LL n, ans = 0;
    scanf("%lld", &n);
    for(LL i = 1; i * i <= n; i++) {
        if(n % i) continue;
        ans += i * euler(n / i);
        if(i * i != n) ans += n / i * euler(i);
    }
    printf("%lld\n", ans);
    return 0;
}