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

 

Logo
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.

License Agreements, Terms of Use, Privacy Policy. .

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