示例#1
0
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
示例#2
0
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)