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
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Oct 01 2021
STATUS
approved