项目Euler＃23-无法找出代码中的错误

Question

链接： http ： //projecteuler.net/problem=23

这是一个警告，如果您不想知道答案，请不要再看了。

我被困了很长时间，我知道网上有解决方案。 但是我只是无法在我的代码中找出问题所在。 我已经附加了我创建的代码和factors模块的代码。 任何帮助表示赞赏。

实际答案=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):
    return sum(factors(num))

Answer 1

factors(4)应该给出[1,2] ，但是产量[1] 。 修正您的factors功能。