OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..600
N. J. A. Sloane, Transforms
FORMULA
a(n) ~ c * d^n / n^(3/2), where d = 1 + A051491 = 3.955765285651994974714..., c = 2.0132071688691451819... . - Vaclav Kotesovec, Aug 25 2014
MAPLE
with(numtheory):
b:= proc(n) option remember; local d, j; `if`(n<=1, n,
(add(add(d*b(d), d=divisors(j))*b(n-j), j=1..n-1))/(n-1))
end:
a:= n-> add(b(k+1)*binomial(n, k), k=0..n):
seq(a(n), n=0..30); # Alois P. Heinz, May 17 2013
MATHEMATICA
b[n_] := b[n] = Module[{d, j}, If[n <= 1, n, Sum[Sum[d*b[d], {d, Divisors[j]}]*b[n-j], {j, 1, n-1}]/(n-1)]]; a[n_] := Sum[b[k+1]*Binomial[n, k], {k, 0, n}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Oct 07 2014, translated from Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved