OFFSET
1,4
COMMENTS
Essentially the same as A001065, but with a(1)=1.
Note that if n is a noncomposite number then a(n)=1.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..4096
MATHEMATICA
Array[Boole[# == 1] + DivisorSigma[1, #] - # &, 80] (* Michael De Vlieger, Oct 30 2017 *)
PROG
(PARI) a(n) = if (n==1, 1, sigma(n) - n); \\ Michel Marcus, Oct 30 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Omar E. Pol, Nov 22 2010
STATUS
approved