OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1823
N. J. A. Sloane, Transforms
FORMULA
G.f.: Sum_{j>=1} j*(Sum_{k>=1} prime(k)*x^k)^j. - Ilya Gutkovskiy, Apr 21 2019
a(n) = Sum_{k=0..n} k * A340991(n,k). - Alois P. Heinz, Feb 01 2021
MAPLE
b:= proc(n) option remember; `if`(n=0, [1, 0], (p->
p+[0, p[1]])(add(ithprime(j)*b(n-j), j=1..n)))
end:
a:= n-> b(n)[2]:
seq(a(n), n=1..27); # Alois P. Heinz, Sep 11 2019
MATHEMATICA
b[n_] := b[n] = If[n==0, {1, 0}, #+{0, #[[1]]}&[Sum[Prime[j] b[n-j], {j, 1, n}]]];
a[n_] := b[n][[2]];
Array[a, 27] (* Jean-François Alcover, Nov 09 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
STATUS
approved