欧拉计划第23题

并非盈数之和

完全数是指真因数之和等于自身的那些数。例如，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;
}