login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A083100 a(n) = 2*a(n-1) + 7*a(n-2). 21
1, 9, 25, 113, 401, 1593, 5993, 23137, 88225, 338409, 1294393, 4957649, 18976049, 72655641, 278143625, 1064876737, 4076758849, 15607654857, 59752621657, 228758827313, 875786006225, 3352883803641, 12836269650857, 49142725927201 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

a(n) = a(n-1) + 8*A015519(n). a(n)/A015519(n+1) converges to sqrt(8).

a(n-1) is the number of compositions of n when there is 1 type of 1 and 8 types of other natural numbers. - Milan Janjic, Aug 13 2010

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (2, 7).

FORMULA

G.f.: (1+7*x)/(1-2*x-7*x^2).

a(n) = (1/2)*(1-2*sqrt(2))^n + sqrt(2)*(1+2*sqrt(2))^n - (1-2*sqrt(2))^n*sqrt(2) + (1/2)*(1+2*sqrt(2))^n, with n >= 0. - Paolo P. Lava, Jun 10 2008

a(n) = binomial transform of 1,8,8,64,64,512. - Al Hakanson (hawkuu(AT)gmail.com), Aug 17 2009

If p[1]=1, and p[i]=8,(i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=det A. - Milan Janjic, Apr 29 2010

G.f.: G(0)/(2*x) - 1/x, where G(k)= 1 + 1/(1 - x*(8*k-1)/(x*(8*k+7) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 28 2013

MATHEMATICA

CoefficientList[Series[(1 + 7 x)/(1 - 2 x - 7 x^2), {x, 0, 25}], x] (* Or *) a[n_] := Simplify[((1 + Sqrt[8])^n + (1 - Sqrt[8])^n)/2]; Array[a, 25, 0] (* Or *) LinearRecurrence[{2, 7}, {1, 1}, 28] (* Or *) Table[ MatrixPower[{{1, 2}, {4, 1}}, n][[1, 1]], {n, 0, 25}] (* Robert G. Wilson v, Sep 18 2013 *)

PROG

(PARI) a(n)=([0, 1; 7, 2]^n*[1; 9])[1, 1] \\ Charles R Greathouse IV, Apr 06 2016

(PARI) x='x+O('x^30); Vec((1+7*x)/(1-2*x-7*x^2)) \\ G. C. Greubel, Jan 08 2018

(MAGMA) I:=[1, 1]; [n le 2 select I[n] else 2*Self(n-1) + 7*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 08 2018

CROSSREFS

Essentially a duplicate of A084058.

The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.

Sequence in context: A146674 A083672 A193644 * A084058 A108570 A092769

Adjacent sequences:  A083097 A083098 A083099 * A083101 A083102 A083103

KEYWORD

easy,nonn

AUTHOR

Mario Catalani (mario.catalani(AT)unito.it), Apr 23 2003

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 11 10:25 EST 2021. Contains 348842 sequences. (Running on oeis4.)