OFFSET
0,3
COMMENTS
Equals the leftmost column in the limit of a triangular matrix B to the n-th power as n increases, where B(m,k) = (k+1)*binomial(m,k) for m > k with zeros along the main diagonal except for B(0,0) = 1.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..250
FORMULA
a(n) ~ c * n^(2*n + 3/2) / exp(2*n), where c = 4.57091094257826787312511057805555932113507180436658388555642557963... . - Vaclav Kotesovec, Dec 09 2016
EXAMPLE
E.g.f.: A(x) = 1 + x + 5*x^2/2! + 52*x^3/3! + 931*x^4/4! + 25516*x^5/5! + 992799*x^6/6! + 52032702*x^7/7! + 3533592843*x^8/8! + 301810098928*x^9/9! + 31663565386063*x^10/10! + ...
PROG
(PARI) {a(n) = if(n==0, 1, sum(k=0, n-1, (k+1)*binomial(n, k) * a(k) ))}
for(n=0, 20, print1(a(n), ", "))
(PARI) {a(n) = my(B=matrix(n+1, n+1, m, k, if(m==1&&k==1, 1, if(m>k, k*binomial(m-1, k-1))))); (B^(n+1))[n+1, 1]}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 08 2016
STATUS
approved