python/021-amicable_numbers.py


1
2
3
4
5
6
7
8
9
10
11
12
13

from helper.numbers import proper_divisors


amicable_nums = []
for a in range(1, 10_000):
    if a in amicable_nums:
        continue

    b = sum(proper_divisors(a))
    if a != b and sum(proper_divisors(b)) == a:
        amicable_nums.extend([a, b])

print(sum(amicable_nums))