2

链接: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): 返回总和(因子(数量))

4

1 回答 1

4

factors(4)应该给予[1,2],但屈服[1]。修复你的factors功能。

于 2013-07-04T04:21:16.410 回答