A380915 E.g.f. A(x) satisfies A(x) = exp(x / (1 - x*A(x)^3)) / (1 - x*A(x)^3).
1, 2, 19, 421, 14453, 676741, 40225525, 2901397997, 246222420841, 24038780973913, 2654362957336481, 327087730518759937, 44498835149618922253, 6624743172003104909957, 1071295799491745519081629, 186999332904147675923216341, 35044146207707289182759039825
Offset: 0
Keywords
Programs
-
PARI
a(n, q=1, r=3, s=0, t=1, u=1) = q*n!*sum(k=0, n, (r*n+(s-r)*k+q)^(k-1)*binomial((r*u+1)*n+((s-r)*u+t-1)*k+q*u-1, n-k)/k!);
Formula
a(n) = n! * Sum_{k=0..n} (3*n-3*k+1)^(k-1) * binomial(4*n-3*k,n-k)/k!.