The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A056545 a(n) = 4*n*a(n-1) + 1 with a(0)=1. 11
 1, 5, 41, 493, 7889, 157781, 3786745, 106028861, 3392923553, 122145247909, 4885809916361, 214975636319885, 10318830543354481, 536579188254433013, 30048434542248248729, 1802906072534894923741, 115385988642233275119425 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,2 COMMENTS For positive n, a(n) equals 4^n times the permanent of the n X n matrix with (5/4)'s along the main diagonal and 1's everywhere else. - John M. Campbell, Jul 10 2011 LINKS Harvey P. Dale, Table of n, a(n) for n = 0..365 Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1. FORMULA a(n) = floor(e^(1/4)*4^n*n!). From Philippe Deléham, Mar 14 2004: (Start) a(n) = n!*Sum_{k=0..n} (4^(n-k)/k!. E.g.f.: exp(x)/(1 - 4*x). (End) a(n) = Sum_{k=0..n} P(n, k)*4^k. - Ross La Haye, Aug 29 2005 a(n) = hypergeometric_U(1, n+2 , 1/4)/4. - Peter Luschny, Nov 26 2014 a(n) = exp(1/4)*4^n*Gamma(n+1, 1/4). a(n) ~ sqrt(2*Pi)*4^n*n^(n+1/2)*exp(1/4-n). - Vladimir Reshetnikov, Oct 14 2016 From Peter Bala, Mar 01 2017: (Start) a(n) = Integral_{x = 0..inf} (4*x + 1)^n*exp(-x) dx. The e.g.f. y = exp(x)/(1 - 4*x) satisfies the differential equation (1 - 4*x)*y' = (5 - 4*x)*y. a(n) = (4*n + 1)*a(n-1) - 4*(n - 1)*a(n-2). The sequence b(n) := 4^n*n! also satisfies the same recurrence with b(0) = 1, b(1) = 4. This leads to the continued fraction representation a(n) = 4^n*n!*( 1 + 1/(4 - 4/(9 - 8/(13 - ... - (4*n - 4)/(4*n + 1) )))) for n >= 2. Taking the limit gives the continued fraction representation exp(1/4) = 1 + 1/(4 - 4/(9 - 8/(13 - ... - (4*n - 4)/((4*n + 1) - ... )))). Cf. A010844. (End) EXAMPLE a(2) = 4*2*a(1) + 1 = 8*5 + 1 = 41. MATHEMATICA Round@Table[Exp[1/4] 4^n Gamma[n + 1, 1/4], {n, 0, 20}] (* Round is equivalent to FullSimplify here, but is much faster; Vladimir Reshetnikov, Oct 14 2016 *) nxt[{n_, a_}]:={n+1, 4a(n+1)+1}; NestList[nxt, {0, 1}, 20][[All, 2]] (* Harvey P. Dale, Mar 19 2019 *) CROSSREFS Cf. A000522, A010844, A010845, A056546, A056547, A001907 for analogs. A056545/(A000142*A000302) is an increasingly good approximation to 4th root of e. Sequence in context: A240996 A346982 A143415 * A325888 A347951 A275787 Adjacent sequences: A056542 A056543 A056544 * A056546 A056547 A056548 KEYWORD nonn,easy AUTHOR Henry Bottomley, Jun 20 2000 EXTENSIONS More terms from James A. Sellers, Jul 04 2000 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 27 17:36 EST 2023. Contains 359845 sequences. (Running on oeis4.)