OFFSET
0,3
MATHEMATICA
nmax = 70; CoefficientList[Series[Product[1/(1 + MoebiusMu[k] x^k), {k, 1, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = (1/n) Sum[Sum[d (-MoebiusMu[d])^(k/d), {d, Divisors[k]}] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 70}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Oct 02 2022
STATUS
approved