A380641 Expansion of e.g.f. exp(x/(1 - 3*x)^3).

1, 1, 19, 379, 8857, 244801, 7904251, 292980619, 12257946289, 570627408097, 29212843607011, 1629314013114811, 98250285167099209, 6365331315043185889, 440712959779710869707, 32460639303987670526731, 2533396174719346231613281, 208776665140069914314618689

Offset

0,3

Links

Table of n, a(n) for n=0..17.

Formula

a(n) = n! * Sum_{k=0..n} 3^k * binomial(3*n-2*k-1,k)/(n-k)!.

E.g.f.: exp( Sum_{k>=1} k*(k+1)/2 * 3^(k-1) * x^k ).

From Vaclav Kotesovec, Jan 29 2025: (Start)

a(n) = (12*n-11)*a(n-1) - 6*(n-1)*(9*n-19)*a(n-2) + 108*(n-3)*(n-2)*(n-1)*a(n-3) - 81*(n-4)*(n-3)*(n-2)*(n-1)*a(n-4).

a(n) ~ 3^n * n^(n - 1/8) / (2 * exp(n - 4*n^(3/4)/3 - sqrt(n)/6 + n^(1/4)/72 + 1/81)) * (1 + 16957/(207360*n^(1/4))). (End)

Mathematica

Table[n! * Sum[3^k * Binomial[3*n-2*k-1, k]/(n-k)!, {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jan 29 2025 *)

Prog

(PARI) a(n) = n!*sum(k=0, n, 3^k*binomial(3*n-2*k-1, k)/(n-k)!);

Crossrefs

Cf. A091695, A380637.

Sequence in context: A023283 A387896 A075839 * A158592 A072359 A222835

Adjacent sequences: A380638 A380639 A380640 * A380642 A380643 A380644

Keyword

nonn,easy

Author

Seiichi Manyama, Jan 28 2025

Status

approved