OFFSET
1,3
EXAMPLE
MAPLE
a:=series(log(1+add(x^ithprime(k)/ithprime(k), k=1..100)), x=0, 24): seq(n!*coeff(a, x, n), n=0..23); # Paolo P. Lava, Mar 26 2019
MATHEMATICA
nmax = 23; Rest[CoefficientList[Series[Log[1 + Sum[x^Prime[k]/Prime[k], {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!]
a[n_] := a[n] = Boole[PrimeQ[n]] (n - 1)! - Sum[k Binomial[n, k] Boole[PrimeQ[n - k]] (n - k - 1)! a[k], {k, 1, n - 1}]/n; a[0] = 0; Table[a[n], {n, 23}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jun 06 2018
STATUS
approved