OFFSET
1,12
COMMENTS
Sum of divisors d of n with sigma(d) > 2*d.
a(n) = n when n is a primitive abundant number (A091191). - _Alonso del Arte_, Jan 19 2013
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
a(12) = 12 because the divisors of 12 are 1, 2, 3, 4, 6, 12, but of those only 12 is abundant.
a(13) = 0 because the divisors of 13 are 1 and 13, neither of which is abundant.
MAPLE
MATHEMATICA
Table[Total@ Select[Divisors@ n, DivisorSigma[1, #] > 2 # &], {n, 96}] (* _Michael De Vlieger_, Jul 16 2016 *)
PROG
(PARI) a(n)=sumdiv(n, d, (sigma(d, -1)>2)*d) \\ _Charles R Greathouse IV_, Jan 15 2013
(Python)
from sympy import divisors, divisor_sigma
def A187795(n): return sum(d for d in divisors(n, generator=True) if divisor_sigma(d) > 2*d) # _Chai Wah Wu_, Sep 22 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
_Timothy L. Tiffin_, Jan 06 2013
STATUS
approved