cp's OEIS Frontend

A118589 Expansion of e.g.f. exp(x + x^2 + x^3).

%I A118589 #39 Dec 17 2024 18:17:08
%S A118589 1,1,3,13,49,261,1531,9073,63393,465769,3566611,29998101,262167313,
%T A118589 2394499693,23249961099,233439305401,2439472944961,26649502709073,
%U A118589 300078056044963,3498896317045789,42244252226263281,524289088799352661
%N A118589 Expansion of e.g.f. exp(x + x^2 + x^3).
%C A118589 Equals row sums of triangle A118588.
%H A118589 Seiichi Manyama, <a href="/A118589/b118589.txt">Table of n, a(n) for n = 0..575</a> (terms 0..200 from Vincenzo Librandi)
%F A118589 a(n) = n! * Sum_{k=1..n}(Sum_{j=0..k}(binomial(j,n-3*k+2*j)*binomial(k,j))/k!), n>0. - _Vladimir Kruchinin_, Sep 01 2010
%F A118589 Recurrence equation: a(n) = a(n-1) + 2*(n-1)*a(n-2) + 3*(n-1)*(n-2)*a(n-3) with initial conditions a(0) = a(1) = 1 and a(2) = 3. - _Peter Bala_, May 14 2012
%F A118589 E.g.f.: 1 + x*(E(0)-1)/(x+1) where E(k) =  1 + (1+x+x^2)/(k+1)/(1-x/(x+1/E(k+1) )); (recursively defined continued fraction). - _Sergei N. Gladkovskii_, Jan 26 2013
%F A118589 a(n) ~ 3^(n/3-1/2) * n^(2*n/3) * exp(7/9*(n/3)^(1/3) + (n/3)^(2/3) - 2*n/3 - 14/81) * (1 + 419/(4374*(n/3)^(1/3)) + 16229573/(191318760*(n/3)^(2/3))). - _Vaclav Kotesovec_, Oct 09 2013
%t A118589 Range[0, 21]!*CoefficientList[ Series[ Exp[x*(1-x^3)/(1 - x)], {x, 0, 21}], x] (* _Zerinvary Lajos_, Mar 23 2007 *)
%t A118589 With[{nn=30},CoefficientList[Series[Exp[x+x^2+x^3],{x,0,nn}],x] Range[0,nn]!] (* _Harvey P. Dale_, Dec 17 2024 *)
%o A118589 (PARI) {a(n)=n!*polcoeff(exp(x+x^2+x^3 +x*O(x^n)),n,x)}
%o A118589 (Maxima) a(n):=n!*sum(sum(binomial(j,n-3*k+2*j)*binomial(k,j),j,0,k)/k!,k,1,n); /* _Vladimir Kruchinin_, Sep 01 2010 */
%Y A118589 Cf. A118588 (triangle), A119013 (eigenvector).
%K A118589 nonn
%O A118589 0,3
%A A118589 _Paul D. Hanna_, May 08 2006