OFFSET
0,3
COMMENTS
(a(n)/n!)^(1/n) tends to 15.1198... - Vaclav Kotesovec, Aug 06 2019
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..310
FORMULA
a(n) = A215861(2n,n).
MAPLE
b:= proc(n, k) option remember; `if`(k<0 or k>n, 0,
`if`(n=0, 1, add(binomial(n-1, i)*b(n-1-i, k-1)*
`if`(i<2, 1, i!/2 +(i+1)^(i-1)), i=0..n-k)))
end:
a:= n-> b(2*n, n):
seq(a(n), n=0..20);
MATHEMATICA
b[n_, k_] := b[n, k] = If[k < 0 || k > n, 0,
If[n == 0, 1, Sum[Binomial[n - 1, i]*b[n - 1 - i, k - 1]*
If[i<2, 1, i!/2 + (i+1)^(i-1)], {i, 0, n-k}]]];
a[n_] := b[2n, n];
a /@ Range[0, 20] (* Jean-François Alcover, Dec 29 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jul 22 2019
STATUS
approved