OFFSET
1,4
COMMENTS
Logarithmic transform of A010051.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..472
EXAMPLE
MAPLE
a:= proc(n) option remember; (t-> `if`(n=0, 0, t(n)-add(a(j)*t(n-j)*
j*binomial(n, j), j=1..n-1)/n))(i-> `if`(isprime(i), 1, 0))
end:
seq(a(n), n=1..25); # Alois P. Heinz, Dec 04 2018
MATHEMATICA
nmax = 24; 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]] - Sum[k Binomial[n, k] Boole[PrimeQ[n - k]] a[k], {k, 1, n - 1}]/n; a[0] = 0; Table[a[n], {n, 24}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jun 06 2018
STATUS
approved