A143638 E.g.f. satisfies: A(x) = exp(x*A(((x+1)^7-1)/7)).

1, 1, 3, 34, 605, 16416, 644647, 33690574, 2252245353, 187575203080, 19000833293771, 2295318297423834, 325536649109809117, 53508774130762119508, 10080999100649218887615, 2156137639664134179951166, 519200838601168582073365073, 139740129055162031424178122096

Offset: 0

Author

Alois P. Heinz, Aug 27 2008

Keywords

• nonn

Links

• Alois P. Heinz, Table of n, a(n) for n = 0..100

Crossrefs

Cf. 7th column of A143632.

Programs

Maple

A:= proc(n,k::nonnegint) option remember; if n<=0 or k=0 then 1 else A(n-1,k)(((x+1)^k-1)/k) fi; unapply(convert(series(exp(x*%), x,n+1), polynom), x) end: a:= n-> coeff(A(n,7)(x), x,n)*n!: seq(a(n), n=0..20);

Mathematica

A[n_, k_] := Module[{f}, f[x_] = If[n <= 0 || k == 0, 1, A[n-1, k][((x+1)^k-1)/k]]; Normal[Series[Exp[x*f[x]], { x, 0, n+1}]] /. x -> #]&; a[n_] := Coefficient[A[n, 7][x], x, n]*n!; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Feb 14 2014, after Maple *)