OFFSET
0,2
COMMENTS
a(n) = A008955(2*n,n).
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 0..120
FORMULA
a(n) ~ c * d^n * n!^2 / n^(3/2), where d = 30.472614208903344955058219429421999299236787591960717... and c = 0.27436634881777520262458169083560234658... - Vaclav Kotesovec, Aug 28 2017
MAPLE
b:= proc(n, k) option remember; `if`(k=0, 1,
add(b(j-1, k-1)*j^2, j=1..n))
end:
a:= n-> b(2*n, n):
seq(a(n), n=0..14); # Alois P. Heinz, Feb 19 2022
MATHEMATICA
Flatten[{1, Table[Coefficient[Expand[Product[1 + k^2*x, {k, 0, 2*n}]], x^n], {n, 1, 15}]}] (* Vaclav Kotesovec, Aug 28 2017 *)
Table[Sum[(-1)^(n-j) * StirlingS1[2*n+1, 2*n+1-j] * StirlingS1[2*n+1, j+1], {j, 0, 2*n}], {n, 0, 15}] (* Vaclav Kotesovec, Aug 28 2017 *)
PROG
(Haskell)
a234324 n = a008955 (2 * n) n
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Dec 24 2013
STATUS
approved