如果一个数恰好等于它的因子之和,则称该数为“完全数” 。各个小于它的约数(真约数,列出某数的约数,去掉该数本身,剩下的就是它的真约数)的和等于它本身的自然数叫做完全数(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;
}