题意:给出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;
}