|
PROG
|
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0, n, log(1+x*subst(A, x, m*x+x*O(x^n)))^m/m!)); n!*polcoeff(A, n)}
for(n=0, 20, print1(a(n), ", "))
(PARI) {Stirling1(n, k)=n!*polcoeff(binomial(x, n), k)}
{a(n)=local(A=1+x); for(i=1, n, A=sum(m=0, n, sum(k=0, m, Stirling1(m, k)*subst(A, x, k*x+x*O(x^n))^m)*x^m/m!)); n!*polcoeff(A, n)}
for(n=0, 20, print1(a(n), ", "))
|