OFFSET
0,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 0..10000
FORMULA
EXAMPLE
a(7) = 93 because 7^2+1 = 50 and the sum of the 6 divisors { 1, 2, 5, 10, 25, 50} is 93.
MAPLE
with(numtheory):for n from 0 to 110 do:x:=divisors(n^2+1):n1:=nops(x):s:=0:for m from 1 to n1 do: s:=s+x[m]:od: printf(`%d, `, s):od:
MATHEMATICA
Table[Total[Divisors[n^2 + 1]], {n, 0, 100}] (* T. D. Noe, Jul 28 2011 *)
DivisorSigma[1, Range[0, 50]^2+1] (* Harvey P. Dale, Aug 03 2020 *)
PROG
(PARI) a(n) = sigma(n^2+1); \\ Michel Marcus, Mar 17 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jul 28 2011
STATUS
approved