OFFSET
1,1
COMMENTS
LINKS
Danko Alorvor, Table of n, a(n) for n = 1..10000
EXAMPLE
2 is a term since the proper divisors of 20 (1, 2, 4, 5, and 10) are all deficient numbers and their sum equals 22, hence the abundance of 20 is 22 - 20 = 2.
PROG
(Python) from sympy import is_abundant, proper_divisors, is_deficient
def primitives_abundances(stop):
abundances = []
for n in range(1, stop + 1):
if is_abundant(n):
divisors = proper_divisors(n)
if all(is_deficient(divisor) for divisor in divisors):
abundances.append(sum(divisors) - n)
return abundances
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Danko Alorvor, Sep 11 2024
STATUS
approved