OFFSET
1,2
FORMULA
a(n)=n!*T(n,1), T(n,m)=m/n*sum(k=1..n-m, sum(i=k..n-m, T(n-m,i)*stirling1(i,k)*k!/i!)*binomial(n+k-1,n-1)), n>m, T(n,n)=1.
PROG
(Maxima) array(TL, fixnum, 30, 30); T(n, m):=if n=m then 1 else if TL[n, m]=0 then TL[n, m]:m/n*sum(sum(T(n-m, i)*stirling1(i, k)*k!/i!, i, k, n-m)*binomial(n+k-1, n-1), k, 1, n-m) else TL[n, m]; makelist(n!*T(n, 1), n, 1, 25);
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, May 07 2012
STATUS
approved