并非盈数之和
完全数是指真因数之和等于自身的那些数。例如,28的真因数之和为1 + 2 + 4 + 7 + 14 = 28,因此28是一个完全数。
一个数n被称为亏数,如果它的真因数之和小于n;反之则被称为盈数。
由于12是最小的盈数,它的真因数之和为1 + 2 + 3 + 4 + 6 = 16,所以最小的能够表示成两个盈数之和的数是24。通过数学分析可以得出,所有大于28123的数都可以被写成两个盈数的和;尽管我们知道最大的不能被写成两个盈数的和的数要小于这个值,但这是通过分析所能得到的最好上界。
找出所有不能被写成两个盈数之和的正整数,并求它们的和。
代码演示
#include <iostream>
using namespace std;
#define MAX_N 28123
int is_prime[MAX_N + 5] = {0};
int prime[MAX_N + 5] = {0};
int f[MAX_N + 5] = {0};
int num[MAX_N + 5] = {0};
void init() {
f[1] = 0;
for (int i = 2; i <= MAX_N; i++) {
if (!is_prime[i]) {
is_prime[i] = i;
f[i] = i + 1;
prime[++prime[0]] = i;
}
for (int j = 1; j <= prime[0] && prime[j] * i <= MAX_N; j++) {
if (i % prime[j] == 0) {
is_prime[i * prime[j]] = is_prime[i] * prime[j];
f[i * prime[j]] = f[i] *
(is_prime[i] * prime[j] * prime[j] - 1) /
(is_prime[i] * prime[j] - 1);
break;
} else {
is_prime[i * prime[j]] = prime[j];
f[i * prime[j]] = f[prime[j]] * f[i];
}
}
}
return ;
}
int main() {
f[0] = 0;
init();
for (int i = 0; i <= MAX_N; i++) {
f[i] -= i;
if (f[i] > i) f[++f[0]] = i;
}
for (int i = 1; i <= f[0]; i++) {
for (int j = 1; j <= f[0]; j++) {
if (f[i] + f[j] > MAX_N) break;
num[f[i] + f[j]] = 1;
}
}
int sum = 0;
for (int i = 0; i <= MAX_N; i++) {
if (!num[i]) sum += i;
}
cout << sum << endl;
return 0;
}