OFFSET
1,4
COMMENTS
If a(n) = 0, then n is a multiply-perfect number (A007691). - Alonso del Arte, Mar 30 2014
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537 (first 1000 terms from T. D. Noe, terms 1001..20000 from Alois P. Heinz).
Walter Nissen, Abundancy : Some Resources.
FORMULA
a(n) = sigma(n) mod n.
a(p) = 1 for p prime.
EXAMPLE
a(12) = 4 because sigma(12) = 28 and 28 == 4 (mod 12).
a(13) = 1 because 13 is prime.
a(14) = 10 because sigma(14) = 24 and 24 == 10 (mod 14).
MAPLE
with(numtheory): seq(sigma(i) mod i, i=1..100);
MATHEMATICA
Table[Mod[DivisorSigma[1, n], n], {n, 80}] (* Alonso del Arte, Mar 30 2014 *)
PROG
(Haskell)
a054024 n = mod (a000203 n) n -- Reinhard Zumkeller, Mar 23 2013
(PARI) a(n)=sigma(n)%n \\ Charles R Greathouse IV, Nov 04 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Asher Auel, Jan 19 2000
STATUS
approved