2020-10-10

完全数（续）

参考完全数
前文得到6，28，496、8128、33550336五个完全数，效率很低。
古希腊数学家欧几里得给出了完全数的通项公式：2${}^{n-1}$(2${}^n$-1)。通项公式中的n和2${}^n$-1都是质数（素数）。利用这个公式，可以求得更多完全数。

from math import ceil, sqrt
from numba import jit
import time


@jit(nopython=True)
def prime(p):
    x = ceil(sqrt(p))
    for i in range(3, x+1, 2):
        if p % i == 0:
            return False
    return True


def calc(m):
    for n in range(2, m):
        if prime(n) and prime(2 ** n - 1):
            print(2 ** (n - 1) * (2 ** n - 1))


t0 = time.time()
nn = int(input("Please input N:"))
calc(nn)
print("耗时{}s。".format(time.time() - t0))

D:\Python\Python38\python.exe D:/Python/study/test26.py
Please input N:63
6
28
496
8128
33550336
8589869056
137438691328
2305843008139952128
2658455991569831744654692615953842176
耗时11.4685800075531s。