OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
N. J. A. Sloane, Transforms
MAPLE
with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: c:= etr(k-> `if`(k<=1, k, ithprime(k-1))): a:= n-> c(n+1): seq(a(n), n=0..40); # Alois P. Heinz, Sep 07 2008
MATHEMATICA
etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n-j], {j, 1, n}]/n]; b]; c = etr[If[# <= 1, #, Prime[#-1]]&]; a = c[#+1]&; Table[a[n], {n, 0, 35}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved