from math import sqrt from mathe import sum_proper_divisor amicable = [] for x in range(1, 10001): if x == sum_proper_divisor(sum_proper_divisor(x)) and x != sum_proper_divisor(x): amicable.append(x) print sum(amicable) print "\n ----------------- \n" for x in range(1,10001): if x == sum_proper_divisor(x): print x
from mathe import sum_proper_divisor abundant = [] vals = [] for num in xrange(2, 28123): if num < sum_proper_divisor(num): abundant.append(num) for num in xrange(1, 28123): val = True print num for abun in abundant: if abun > num: break if (num - abun) in abundant: val = False break if val: vals.append(num) print sum(vals)