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

0, 1, 1, -1, 0, -1, 1, -1, -1, -1, 2, 0, 0, -3, 0, 0, 3, -3, -1, -1, 4, -4, -1, -5, 6, 2, 2, -17, 4, 4, 13, -16, -7, -11, 30, -14, -7, -34, 42, 7, 16, -80, 27, 6, 95, -117, -21, -60, 223, -97, -48, -265, 326, 53, 129, -800, 242, 93, 752, -948, -105, -499, 1806, -853, -240, -2189, 2750, 124

Offset

1,11

Links

Robert Israel, Table of n, a(n) for n = 1..2500

Maple

N:= 20: # for a(1)..a(N)
P:= 1: a:= Vector(N):
for n from 1 to N do
  c:= coeff(P, x, n);
  if isprime(n) then a[n]:= 1-c  else a[n]:= -c fi;
  P:= series(P/(1-a[n]*x^n), x, N+1);
od:
convert(a, list); # Robert Israel, Mar 01 2022

Crossrefs

Cf. A000040, A002121, A010051, A030010, A147557, A305871, A305882, A308298, A328777, A348128.

Sequence in context: A280164 A049597 A210951 * A233440 A280728 A175676

Adjacent sequences: A348124 A348125 A348126 * A348128 A348129 A348130

Keyword

sign

Author

Ilya Gutkovskiy, Oct 01 2021

Status

approved