OFFSET
0,3
COMMENTS
For n>0, a(n) is divisible by n: a(n)/n = A111535(n).
FORMULA
a(n) = [x^n] Log( Sum_{m=0..n} (n-1+m)!/(n-1)!*x^m ).
MATHEMATICA
T[n_, k_] := T[n, k] = Which[n<0 || k<0, 0, k==0 || k==1, 1, n==0, k!, True, (T[n-1, k+1]-T[n-1, k])/n - Sum[T[n, j] T[n-1, k-j], {j, 1, k-1}]];
a[n_] := T[n, n];
Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Aug 09 2018 *)
PROG
(PARI) {a(n)=if(n<0, 0, if(n==0, 1, polcoeff(log(sum(m=0, n, (n-1+m)!/(n-1)!*x^m)), n)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 06 2005
STATUS
approved