OFFSET
0,2
PROG
(PARI) {a(n)=local(M=matrix(n+3, n+3, m, j, if(m==j, 1, if(m==j+1, -m+1, -(m-j-1)*polcoeff(log(sum(i=0, m, (i+1)!/1!*x^i)), m-j-1))))); if(n<0, 0, sum(i=1, #M, (M^0-M)^i/i)[n+3, 3])}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 06 2005
STATUS
approved