This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A010845 #52 Jul 02 2025 16:01:55
%S A010845 1,4,25,226,2713,40696,732529,15383110,369194641,9968255308,
%T A010845 299047659241,9868572754954,355268619178345,13855476147955456,
%U A010845 581929998214129153,26186849919635811886,1256968796142518970529
%N A010845 a(n) = 3*n*a(n-1) + 1, a(0) = 1.
%C A010845 a(n)/(A000142*A000244) is an increasingly good approximation to cube root of e.
%C A010845 Related to Incomplete Gamma Function at 1/3. - _Michael Somos_, Mar 26 1999
%C A010845 For positive n, a(n) equals 3^n times the permanent of the n X n matrix with (4/3)'s along the main diagonal, and 1's everywhere else. - _John M. Campbell_, Jul 10 2011
%D A010845 M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 262.
%H A010845 Vincenzo Librandi, <a href="/A010845/b010845.txt">Table of n, a(n) for n = 0..200</a>
%H A010845 M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 262.
%H A010845 Roland Bacher, <a href="https://doi.org/10.37236/2522">Counting Packings of Generic Subsets in Finite Groups</a>, Electr. J. Combinatorics, 19 (2012), #P7. - From _N. J. A. Sloane_, Feb 06 2013
%H A010845 M. Z. Spivey and L. L. Steil, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL9/Spivey/spivey7.html">The k-Binomial Transforms and the Hankel Transform</a>, J. Integ. Seqs. Vol. 9 (2006), #06.1.1.
%F A010845 E.g.f.: exp(x)/(1-3*x).
%F A010845 a(n) = floor( n!*e^(1/3)*3^n ) = n! * (Sum_{k=0..n} 3^(n-k) / k!) = n! * (e^(1/3) * 3^n - Sum_{k>n} 3^(n-k) / k!). - _Michael Somos_, Mar 26 1999
%F A010845 a(n) = Sum_{k=0..n} P(n, k)*3^k. - _Ross La Haye_, Aug 29 2005
%F A010845 Binomial transform of A032031. - _Carl Najafi_, Sep 11 2011
%F A010845 Conjecture: a(n) +(-3*n-1)*a(n-1) +3*(n-1)*a(n-2)=0. - _R. J. Mathar_, Feb 16 2014
%F A010845 a(n) = hypergeometric_U(1,n+2,1/3)/3. - _Peter Luschny_, Nov 26 2014
%F A010845 From _Peter Bala_, Mar 01 2017: (Start)
%F A010845 a(n) = Integral_{x >= 0} (3*x + 1)^n*exp(-x) dx.
%F A010845 The e.g.f. y = exp(x)/(1 - 3*x) satisfies the differential equation (1 - 3*x)*y' = (4 - 3*x)*y. Mathar's recurrence above follows easily from this.
%F A010845 The sequence b(n) := (3^n)*n! also satisfies Mathar's recurrence with b(0) = 1, b(1) = 3. This leads to the continued fraction representation a(n) = (3^n)*n!*( 1 + 1/(3 - 3/(7 - 6/(10 - ... - (3*n - 3)/(3*n + 1) )))) for n >= 2. Taking the limit as n -> oo gives the continued fraction representation exp(1/3) = 1 + 1/(3 - 3/(7 - 6/(10 - ... - (3*n - 3)/((3*n + 1) - ... )))). Cf. A010844. (End)
%e A010845 1 + 4*x + 25*x^2 + 226*x^3 + 2713*x^4 + 40696*x^5 + 732529*x^6 + ...
%t A010845 Table[ Gamma[ n, 1/3 ]*Exp[ 1/3 ]*3^(n-1), {n, 1, 24} ]
%t A010845 a[ n_] := If[ n<0, 0, Floor[ n! E^(1/3) 3^n ]] (* _Michael Somos_, Sep 04 2013 *)
%t A010845 Range[0, 20]! CoefficientList[Series[Exp[x]/(1 - 3 x), {x, 0, 20}], x] (* _Vincenzo Librandi_, Feb 17 2014 *)
%o A010845 (PARI) {a(n) = if( n<0, 0, n! * sum(k=0, n, 3^(n-k) / k!))} /* _Michael Somos_, Sep 04 2013 */
%Y A010845 Cf. A000522, A010844, A056545, A056546, A056547 for analogs.
%K A010845 nonn,easy
%O A010845 0,2
%A A010845 _Simon Plouffe_
%E A010845 Better description and formulas from _Michael Somos_
%E A010845 More terms from _James Sellers_, Jul 04 2000