GCD - Extreme (II) UVA - 11426 （欧拉函数与gcd的关系）

思路：若gcd(n,x)=i，则gcd(n/i,x/i)=1;

假设b[n]表示1到n-1与n的gcd的和，那么G[n]=G[n-1]+b[n];

a[i]表示与gcd（n, x）= i 的x的个数；b[n]=sum( a[i] * i ) , 所以我们只需求a[i]即可；根据gcd(n, x)=i ----->gcd(n/i, x/i) = 1,

因此仅仅要求出欧拉函数phi(n / i)，就能够得到与n / i互质的个数，从而求出gcd(x , n) = i的个数，这样总体就能够求解了

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int32_t N = 4e6 + 10;
int t, n, phi[N], b[N], g[N];
void euler() {
    for(int i = 2; i < N; i++) {
        if(!phi[i])
            for(int j = i; j < N; j += i) {
                if(!phi[j])
                    phi[j] = j;
                phi[j] = phi[j] / i * (i - 1);
            }
    }
}
int32_t main() {
    euler();
    for(int i = 1; i <= N; i++) {
        for(int j = i * 2; j <= N; j += i) {
            b[j] += phi[j / i] * i;
        }
    }
    for(int i = 2; i <= N; i++)
        g[i] = g[i - 1] + b[i];
    while(cin >> n && n) {
        cout << g[n] << endl;
    }
    return 0;
}