A353606 - OEIS

A353606

Product_{n>=1} (1 + a(n)*x^n) = 1 + x + Sum_{n>=2} prime(n-1)*x^n.

1, 2, 1, 4, 1, 0, -1, 11, -3, -7, 4, -10, 5, -15, -4, 151, -9, -50, -3, -63, 19, -176, 6, 591, -27, -637, 13, -999, 50, -1957, -49, 27250, -162, -7887, 83, -12821, 468, -27226, -40, 127341, -1215, -99166, -526, -174140, 2640, -362870, 1673, 1419061, -4516, -1344620

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

LINKS

Table of n, a(n) for n=1..50.

MATHEMATICA

nn = 50; f[x_] := Product[(1 + a[n] x^n), {n, 1, nn}]; sol = SolveAlways[0 == Series[f[x] - 1 - x - Sum[Prime[k - 1] x^k, {k, 2, nn}], {x, 0, nn}], x]; Table[a[n], {n, 1, nn}] /. sol // Flatten

CROSSREFS

Cf. A008578, A147557, A353605.

Sequence in context: A295881 A371890 A117971 * A220420 A190616 A238018

Adjacent sequences: A353603 A353604 A353605 * A353607 A353608 A353609

KEYWORD

sign

AUTHOR

Ilya Gutkovskiy, May 07 2022

STATUS

approved