Project Euler problem 23

又是一个跟亲和数有关的东西。

不过这次稍微变形了下。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <stack>
#include <cmath>
#include <map>
#include <ctime>
#define MAXN 111111
#define INF 100000007
using namespace std;
int up = 28123;
int a[33333];
int ok[33333];
int num[33333];
int cnt;
int main()
{
    for(int i = 1; i <= up; i++)
        for(int j = i + i; j <= up; j += i)
            a[j] += i;
    for(int i = 1; i <= up; i++)
        if(a[i] > i)
            num[cnt++] = i;
    for(int i = 0; i < cnt; i++)
        for(int j = i; j < cnt; j++)
            if(num[i] + num[j] <= up)
                ok[num[i] + num[j]] = 1;
    int sum = 0;
    for(int i = 1; i <= up; i++)
        if(!ok[i]) sum += i;
    cout << sum << endl;
    return 0;
}