完数与盈数（时间复杂度极低的算法）

题目描述
一个数如果恰好等于它的各因子（该数本身除外）子和，如：6=3+2+1，则称其为“完数”；
若因子之和大于该数，则称其为“盈数”。求出2 到60 之间所有“完数”和“盈数”，并以如
下形式输出： E: e1 e2 e3 …(ei 为完数) G: g1 g2 g3 …(gi 为盈数)
输入
无

输出
按描述要求输出（注意EG后面的冒号之后有一个空格）。

如果题目范围为2~n，当n很大时，普通的方法会需要很长时间对因子进行求和。以下算法更好。

```cpp
#include <stdio.h>
#include <math.h>

using namespace std;

const int maxn = 61;
int prime[maxn], pnum = 0;
bool p[maxn] = {false};

void findPrime(){
    for(int i = 2; i < maxn; i++){
        if(p[i] == false){
            prime[pnum++] = i;
            for(int j = i + i; j < maxn; j += i){
                p[j] = true;
            }
        }
    }
}

struct factor{
    int x, cnt;
}fac[10];

int main()
{
    int E[maxn], G[maxn], E_cnt = 0, G_cnt = 0;

    findPrime();

    for(int i = 2; i < maxn; i++){
        int n = i;
        int sqr = (int)sqrt(1.0 * n);
        int num = 0;//num of fac member

        for(int j = 0; j < pnum && prime[j] <= sqr; j++){
            if(n % prime[j] == 0){
                fac[num].x = prime[j];
                fac[num].cnt = 0;

                while(n % prime[j] == 0){
                    fac[num].cnt++;
                    n /= prime[j];
                }

                num++;
            }
        }

        if(n != 1){
            fac[num].x = n;
            fac[num].cnt = 1;
            num++;
        }

        int sum = 1;
        for(int k = 0; k < num; k++){
            sum *= ( (1 - (int)pow(1.0 * fac[k].x, 1.0 * fac[k].cnt + 1)) / (1 - fac[k].x) );
        }
        sum -= i;

        if(sum == i){
            E[E_cnt++] = i;
        }
        else if(sum > i){
            G[G_cnt++] = i;
        }

    }

    printf("E:");
    for(int i = 0; i < E_cnt; i++) printf(" %d", E[i]);
    printf(" G:");
    for(int i = 0; i < G_cnt; i++) printf(" %d", G[i]);
    printf("\n");


    return 0;
}