OFFSET
0,2
LINKS
Harvey P. Dale, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: 1 + Sum_{k>=1} k*(k^2 + 1)*x^k/(1 - x^k). - Ilya Gutkovskiy, Mar 17 2017
MATHEMATICA
nxt[{n_, a_}]:={n+1, DivisorSigma[1, n+1]+DivisorSigma[3, n+1]}; NestList[nxt, {0, 1}, 50][[All, 2]] (* Harvey P. Dale, Apr 14 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Mar 20 2004
STATUS
approved