完全数

如果一个数恰好等于它的因子之和，则称该数为“完全数” 。各个小于它的约数（真约数,列出某数的约数，去掉该数本身，剩下的就是它的真约数）的和等于它本身的自然数叫做完全数（Perfect number），又称完美数或完备数。（第一个完全数是6，第二个完全数是28，第三个完全数是496，后面的完全数还有8128、33550336等等）

梅森素数：
古希腊数学家欧几里得在名著《几何原本》中证明了素数有无穷多个，并论述完全数时提出：如果2 ^ P-1是素数(其中指数P也是素数)，则2 ^ (P-1) * (2^P-1)是完全数。

#include<iostream>
#include<vector>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
typedef long long ll;
int INF = 0x3f3f3f3f;
int prime[25] = {2,3,5,7,11,13,17,19,23};
bool is_prime(ll n) {
    for (int i = 2; (ll)i * i <= n; i++)
        if (n % i == 0)
            return false;
    return true;
}

ll mod_pow(ll x, ll n, ll mod){
    if(n == 0) return 1;
    ll res = mod_pow(x * x % mod, n/2, mod);
    if(n & 1) res = res * x % mod;
    return res;
}

void solve() {
    for (int i = 0; i < 9; i++) {
        ll res = mod_pow(2, prime[i], INF);
        if (is_prime(res - 1)) {
            int x = prime[i] - 1;
            ll ans = mod_pow(2, x, INF);
            cout << ans * (res - 1)<< endl;
        }
    }
}
int main() {
    solve();
}

暴力做法很难求出第5个完全数。

题目：求出前5个完全数。
官方题解：

#include<stdio.h>
int main()
{
    int x[5]= {1,2,4,6,12};
    int i;
    for(i =0; i <5; i++)
    {
        printf("%d\n",(1<< x[i])*((1<< x[i]+1)-1));
    }
    return 0;
}