Question

以下是问题的链接：link。我的代码似乎以指数方式减慢，但我无法确定原因或更有效的算法。我的方法是识别所有数量达到极限的数量，并通过从原始数字中减去每个数字并查看差异是否在丰富数字列表中来确定不是数字总和的数字。关于发生了什么以及更好地处理这个问题的任何想法？

以下是我使用的代码：

import numpy as np 
import math
import itertools

def divisors(n): return sorted(np.unique(np.array([[x,n/x] for x in range(1,int(round(math.sqrt(n))+1)) if n%x == 0]).flatten()).tolist())[0:-1]



ubound = 28123
abundant_numbers = [x for x in range(1,ubound) if x < sum(divisors(x))] 


def is_sum_of_abundant(n):
    isob = False
    for i in abundant_numbers:
        if (n - i) <=0: 
            continue
        else:
            if (n - i) in abundant_numbers:
                isob = True 
    return isob

s = 0

for x in range(1,ubound):
    print "%.2f percent\n" % ((float(x)/ubound)*100)
    if is_sum_of_abundant(x):
        print "{} is abundant".format(x)
    else:
        s+=x
        print "{} is not abundant".format(x)
print s

Answer 1

您尝试的一件事是计算除数之和的更好方法 - 请参阅sigma函数here的定义。基本上你找到了主要因素并使用

这一事实

sigma(ab) = sigma(a) * sigma(b)
sigma(p^n) = (p^(n+1)-1)/(p-1)

其中sigma（n）是正除数的总和。

项目Euler 23：为什么我的代码这么慢，我该怎么做才能加快它？

1 个答案: