OFFSET
0,4
FORMULA
MATHEMATICA
nmax = 22; CoefficientList[Series[-Log[1 - Sum[x^(k (k + 1)/2)/(k (k + 1)/2)!, {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = Boole[IntegerQ[(8 n + 1)^(1/2)]] + Sum[Binomial[n, k] Boole[IntegerQ[(8 (n - k) + 1)^(1/2)]] k a[k], {k, 1, n - 1}]/n; a[0] = 0; Table[a[n], {n, 0, 22}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 09 2019
STATUS
approved