OFFSET
0,2
COMMENTS
Weigh transform of odd primes.
MATHEMATICA
nmax = 33; CoefficientList[Series[Product[(1 + x^k)^Prime[k + 1], {k, 1, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = (1/n) Sum[Sum[(-1)^(k/d + 1) d Prime[d + 1], {d, Divisors[k]}] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 33}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 28 2022
STATUS
approved