A154597 a(n) = 15*a(n-1) + a(n-2) with a(0) = 0, a(1) = 1.

0, 1, 15, 226, 3405, 51301, 772920, 11645101, 175449435, 2643386626, 39826248825, 600037119001, 9040383033840, 136205782626601, 2052127122432855, 30918112619119426, 465823816409224245, 7018275358757483101

Offset

0,3

Comments

Limit_{n -> infinity} a(n)/a(n-1) = (15 + sqrt(229))/2. - Klaus Brockhaus, Oct 07 2009

For more information about this type of recurrence follow the Khovanova link and see A054413, A086902 and A178765. - Johannes W. Meijer, Jun 12 2010

For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 15's along the main diagonal, and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011

a(n) equals the number of words of length n - 1 on alphabet {0,1,...,15} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015

From Michael A. Allen, Apr 30 2023: (Start)

Also called the 15-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.

a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 15 kinds of squares available. (End)

Links

G. C. Greubel, Table of n, a(n) for n = 0..840 (a(0) = 0 added by Jianing Song)

Michael A. Allen and Kenneth Edwards, Fence tiling derived identities involving the metallonacci numbers squared or cubed, Fib. Q. 60:5 (2022) 5-17.

Tanya Khovanova, Recursive sequences. [From Johannes W. Meijer, Jun 12 2010]

Index entries for linear recurrences with constant coefficients, signature (15,1).

Formula

G.f.: x/(1 - 15*x - x^2). - Klaus Brockhaus, Jan 12 2009, corrected Oct 07 2009

a(n) = ((15 + sqrt(229))^n - (15 - sqrt(229))^n)/(2^n*sqrt(229)).

From Johannes W. Meijer, Jun 12 2010: (Start)

Limit_{k -> infinity} a(n+k)/a(k) = (A090301(n) + a(n)*sqrt(229))/2.

Limit_{n -> infinity} A090301(n)/a(n) = sqrt(229).

a(2n+1) = 15*A098245(n-1).

a(3n+1) = A041427(5n), a(3n+2) = A041427(5n+3), a(3n+3) = 2*A041427(5n+4). (End)

Mathematica

a=0; lst={}; s=0; Do[a=s-(a-1); AppendTo[lst, a]; s+=a*15, {n, 3*4!}]; lst (* Vladimir Joseph Stephan Orlovsky, Oct 27 2009 *)
CoefficientList[Series[x/(1 - 15*x - x^2), {x, 0, 50}], x] (* G. C. Greubel, Apr 16 2017 *)

Prog

(Magma) Z<x>:=PolynomialRing(Integers()); N<r>:=NumberField(x^2-229); S:=[ ((15+r)^n-(15-r)^n)/(2^n*r): n in [1..17] ]; [ Integers()!S[j]: j in [1..#S] ]; // Klaus Brockhaus, Jan 12 2009
(PARI) x='x+O('x^50); Vec(x/(1 - 15*x - x^2)) \\ G. C. Greubel, Apr 16 2017 *)

Crossrefs

Row n=15 of A073133, A172236 and A352361 and column k=15 of A157103.

First bisection is A098247.

Cf. A166125 (decimal expansion of sqrt(229)), A166126 (decimal expansion of (15 + sqrt(229))/2).

Cf. also A041427, A090301, A098245.

Sequences with g.f. 1/(1-k*x-x^2) or x/(1-k*x-x^2): A000045 (k=1), A000129 (k=2), A006190 (k=3), A001076 (k=4), A052918 (k=5), A005668 (k=6), A054413 (k=7), A041025 (k=8), A099371 (k=9), A041041 (k=10), A049666 (k=11), A041061 (k=12), A140455 (k=13), A041085 (k=14), this sequence (k=15), A041113 (k=16), A178765 (k=17), A041145 (k=18), A243399 (k=19), A041181 (k=20).

Sequence in context: A012643 A364857 A067222 * A041422 A129836 A075262

Adjacent sequences: A154594 A154595 A154596 * A154598 A154599 A154600

Keyword

nonn,easy

Author

Al Hakanson (hawkuu(AT)gmail.com), Jan 12 2009

Extensions

Extended beyond a(7) by Klaus Brockhaus and Philippe Deléham, Jan 12 2009

Name from Philippe Deléham, Jan 12 2009

Edited by Klaus Brockhaus, Oct 07 2009

Missing a(0) added by Jianing Song, Jan 29 2019

Status

approved