UVA 11426 GCD - Extreme (II) (欧拉函数)

题意：给出n，求下列值。

题解：欧拉函数
gcd(n, m) = i 的个数为p[m / i] (n≤m)，令f(m)=i * p[m / i]，将1~n的f(i)累加即可，这里预处理一下，查询直接O(1)。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<map>
#define ll long long
using namespace std;
const int maxn = 5e6;
int n;
ll ans;
int p[maxn + 10];
ll s[maxn + 10], f[maxn + 10];
bool isPrime[maxn + 10];
void Eular() {
    for (int i = 1; i <= maxn; i++) p[i] = i;
    memset(isPrime, true, sizeof(isPrime));
    isPrime[0] = isPrime[1] = false;
    p[1] = 0;
    for (int i = 2; i <= maxn; i++) {
        if (isPrime[i]) {
            for (int j = i; j <= maxn; j += i) {
                isPrime[j] = false;
                p[j] -= p[j] / i;
            }
        }
    }
    for (int i = 1; i <= maxn; i++) {
        for (int j = 2 * i; j <= maxn; j += i) {
            f[j] += 1ll * p[j / i] * i;
        }
    }
    for (int i = 1; i <= maxn; i++) {
        s[i] = s[i - 1] + f[i];
    }
}
int main() {
    Eular();
	while (~scanf("%d", &n) && n) {
        printf("%lld\n", s[n]);
	}
	return 0;
}