Python识别完美数

01 完美数

完美数(perfect number，又称完全数)指，它所有的真因子(即除了自身以外的因子)和，恰好等于它自身。

第一个完美数：6，

第二个完美数：28，

第三个完美数：496，

第四个完美数：8128，

第五个完美数：33550336,

.......

02 探索

在茫茫数海中，第五个完美数(33550336)要大得多，居然藏在千万位数的深处！它在十五世纪被人们发现，计算机问世后，借助这一有力工具，数论爱好者们继续探索。

笛卡尔曾公开预言：“能找出的完美数是不会多的，好比人类一样，要找一个完美人亦非易事。”

时至今日，人们一直没有发现有奇完美数的存在。于是是否存在奇完美数成为数论中的一大难题。只知道即便有，这个数也是非常之大，并且需要满足一系列苛刻的条件。

经过不少数学家研究，到2013年为止，一共找到了48个完美数。

0 3 有趣性质

1. 目前发现的完美数都是以6或8结尾，会不会有奇完全数存在？如果存在，它必须大于10^300，至今无人能回答这些问题。

2. 所有的完美数都是三角形数。例如：

 

3. 所有完美数的倒数都是调和数。例如：

 

4. 可以表示成连续奇立方数之和。除6以外的完全数，都可以表示成连续奇立方数之和，并规律式增加。例如：

 

04 判断

如何判断是为否完美数呢？在计算机数值型可以表达的的范围内，我们可以尝试找一找。

这是一道leetcode题(No.507)，我前段时间写过一个解，在leetcode平台上已通过：

 

已知完美数都以6或8结尾，所以才有了上面的方法，注意这不是寻找一个数所有因子的方法。使用6或8结尾这个小trick，实现更高效（ > 95.26%），但不严谨。

很遗憾，今天我发现这是一个错误的解法，虽然在Leetcode上已经通过。原因如下，我们试图打印尽可能多的完美数：

 

 

第1个完美数: 6

第2个完美数: 28

第3个完美数: 120

第4个完美数: 496

第5个完美数: 2016

第6个完美数: 8128

第7个完美数: 32640

第8个完美数: 130816

第9个完美数: 523776

第10个完美数: 2096128

很明显，120不是一个完美数，因此，可以确定Leetcode平台遗漏了这些cases，已经将此问题提交到Leetocode，如下所示：

05 正解

如果遍历所有的小于num的数，check是否为其因子，时间复杂度为o(n),在平台上提交会超时。

一种更好的解法，时间复杂度为O(sqrt(n)), 因为num的两个因子：num_i和num_j，假设num_i < num_j ，则 num_i 的最大值为 sqrt(num), 所以我们只需要遍历到sqrt(num)即可。

代码如下：

 

06 更多完美数

6. 8,589,869,056

7. 137,438,691,328

8. 2,305,843,008,139,952,128

9. 2,658,455,991,569,831,744,654,692,615,953,842,176

10.191,561,942,608,236,107,294,793,378,084,303,638,130,997,321,548,169,216

11.13,164,036,458,569,648,337,239,753,460,458,722,910,223,472,318,386,943,117,783,728,128

12.14,474,011,154,664,524,427,946,373,126,085,988,481,573,677,491,474,835,889,066,354,349,131,199,152,128

……

……

47 ……2^42643800 X (2^42643801-1)

48 ……2^57885160 X (2^57885161-1)

由于后面数字位数较多，例子只列到12个，第13个有314位。

到第39个完全数有25674127位数，据估计它以四号字打出时需要一本字典大小的书。