OFFSET
0,4
LINKS
N. J. A. Sloane, Table of n, a(n) for n=0..200
N. J. A. Sloane, Transforms
FORMULA
Take Euler transform of A034825 and shift right. (Christian G. Bower).
MAPLE
For Maple program see link in A000235.
with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: shr:= proc(p) n->`if`(n=0, 1, p(n-1)) end: b[0]:= etr(n->1): for j from 1 to 7 do b[j]:= etr(shr(b[j-1])) od: a:= shr(b[7]): seq(a(n), n=0..40); # Alois P. Heinz, Sep 08 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]; shr[p_] = If[# == 0, 1, p[#-1]]&; b[0] = etr[1&]; For[j = 1, j <= 7, j++, b[j] = etr[shr[b[j-1]]]]; a = shr[b[7]]; Table[a[n], {n, 0, 31}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved