完全数(Perfect number),又称完美数或完备数,是一些特殊的自然数。它所有的真因子(即除了自身以外的约数)的和(即因子函数),恰好等于它本身。如果一个数恰好等于它的因子之和,则称该数为“完全数”。第一个完全数是 6,它有约数 1、2、3、6,除去它本身 6 外,其余 3 个数相加,1+2+3=6。第二个完全数是 28,它有约数 1、2、4、7、14、28,除去它本身 28 外,其余 5 个数相加,1+2+4+7+14=28。第三个完全数是 496,有约数 1、2、4、8、16、31、62、124、248、496,除去其本身 496 外,其余 9 个数相加,1+2+4+8+16+31+62+124+248=496。
#coding: utf-8
#问题:1000以内的完美数
#答案:[6,28,496]
import time
#方法一:常规做法
def solution1():
start_time = time.time()
for n in range(1000):
num_list = []
for num in range(1,1000):
sum_num = 0
for i in range(1,num):
if num % i == 0:
sum_num += i
else:
continue
if sum_num == num:
num_list.append(num)
print(num_list)
print(time.time()-start_time)
#运行1000次的时间:50.213872
#方法二:大于这个数1/2的公因子是不存在的
def solution2():
start_time = time.time()
for n in range(1000):
num_list = []
for num in range(1,1000):
sum_num = 0
for i in range(1,int(num/2)+1):
if num % i == 0:
sum_num += i
else:
continue
if sum_num == num:
num_list.append(num)
print(num_list)
print(time.time()-start_time)
#运行1000次的时间:25.8714797
#方法三:一对因子总会发布在这个数的根号两边,或者就是这平方根值
def solution3():
start_time = time.time()
for n in range(1000):
num_list = []
for num in range(2,1000):#从1开始,会把1也写入数组
sum_num = 1#每个数都会有个因子是1
for i in range(2,int(pow(num,1/2)+1)):#从2开始,防止把num写入数组
if num % i == 0:
res = int(num/i)
if res == i:#判断是不是恰好是平方根
sum_num += i#正好是平方根,只写入一个数据
else:
sum_num += (i+res)#把一对因子写入
else:
continue
if sum_num == num:
num_list.append(num)
print(num_list)
print(time.time()-start_time)
#运行1000次的时间:4.66926717758
if __name__ == '__main__':
solution3()
PS:好像有公式和素数、质数相关?欧拉公式?我数学不好,不是很懂。如果能直接生成完美数就好了。