链接:http ://projecteuler.net/problem=23
这是一个警告,如果您不想知道答案,请不要再看下去。
我已经卡了很长一段时间了。我知道网上有解决方案。但我只是无法找出我的代码中的问题。我附上了我创建的代码和因素模块的代码。任何帮助表示赞赏。
实际答案 = 4179871。我的答案 = 4190404
from math import *
from time import *
from prime import *
from factors import *
abundant = list(n for n in xrange(12,28124) if n < sum_of_factors(n))
sums = {}
for i in abundant:
for j in abundant:
if (i+j) > 28123:
break
else:
sums[i+j] = 1
non_abun = [i for i in range(1,28124)]
print sum(non_abun) - sum(sums)
因子代码 -
from time import * from math import *
def factors(num): factors_array = [1,]
n = num for k in range(2,int(ceil(sqrt(num)))): if(n%k == 0): if(k not in factors_array): factors_array.append(k) factors_array.append(n/k) return factors_array
def sum_of_factors(num): 返回总和(因子(数量))