bzoj 2818 Gcd

#include<bits/stdc++.h>

using namespace std;

#define ll long long
const int maxn = 1e7 + 5;
int n, pri[maxn];
ll f[maxn];
bool Notpri[maxn]; int tot = 0;

void pre(int x) {
	f[1] = 1;
	for (int i = 2; i <= n; ++i) {
		if(!Notpri[i]) pri[++tot] = i, f[i] = i - 1;
		for (int j = 1; j <= tot && i * pri[j] <= n; ++j) {
			Notpri[i * pri[j]] = true;
			if(i % pri[j] == 0) {
				f[i * pri[j]] = f[i] * pri[j];
				break;
			} f[i * pri[j]] = f[i] * (pri[j] - 1);
		}
	}
	for (int i = 1; i <= n; ++i) f[i] += f[i - 1];
}

int main() {
	scanf("%d", &n);
	pre(n);
	ll ans = 0;
	for (int i = 1; i <= tot; ++i) {
		ans += f[n / pri[i]] * 2; ans--;
	} printf("%lld\n", ans);
	return 0;
}