OFFSET
0,16
FORMULA
MATHEMATICA
nmax = 75; CoefficientList[Series[Product[1/(1 + x^(3 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = (1/n) Sum[DivisorSum[k, (-1)^(k/#) # &, Mod[#, 3] == 2 &] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 75}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jun 27 2024
STATUS
approved