OFFSET
0,2
COMMENTS
Convolution inverse of A353065.
MATHEMATICA
nmax = 42; CoefficientList[Series[Product[(1 - x^k)^Prime[k + 1], {k, 1, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = -(1/n) Sum[Sum[d Prime[d + 1], {d, Divisors[k]}] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 42}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Apr 28 2022
STATUS
approved