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