OFFSET
1,1
MATHEMATICA
nn = 29; f[x_] := Product[(1 + a[n] x^n), {n, 1, nn}]; sol = SolveAlways[0 == Series[f[x] - 1 - Sum[Prime[k + 1] x^k, {k, 1, nn}], {x, 0, nn}], x]; Table[a[n], {n, 1, nn}] /. sol // Flatten
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, May 07 2022
STATUS
approved