OFFSET
1,6
FORMULA
Numerators of coefficients in expansion of Sum_{k>=1} x^(k*(k+1)) / (k * (1 - x^k)).
EXAMPLE
0, 1, 1, 1, 1, 3/2, 1, 3/2, 1, 3/2, 1, 11/6, 1, 3/2, 4/3, 3/2, 1, 11/6, ...
MATHEMATICA
nmax = 80; CoefficientList[Series[Sum[x^(k (k + 1))/(k (1 - x^k)), {k, 1, nmax}], {x, 0, nmax}], x] // Numerator // Rest
PROG
(PARI) a(n) = numerator(sumdiv(n, d, if (d^2 < n, 1/d))); \\ Michel Marcus, May 14 2024
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Ilya Gutkovskiy, May 14 2024
STATUS
approved