login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A123943 The (1,5)-entry in the 5 X 5 matrix M^n, where M={{5, 3, 2, 1, 1}, {3, 2, 1, 1, 0}, {2, 1, 1, 0, 0}, {1, 1, 0, 0, 0}, {1, 0, 0, 0, 0}} (n>=0). 0
0, 1, 5, 40, 315, 2490, 19681, 155563, 1229604, 9719061, 76821600, 607214857, 4799560053, 37936780428, 299860673343, 2370164848026, 18734305316497, 148080078051971, 1170457572108040, 9251554605638681, 73126326541645648 (list; graph; refs; listen; history; internal format)
OFFSET

0,3

REFERENCES

See A123942 for references.

FORMULA

a(n)=8a(n-1)-6a(n-3)+a(n-5) for n>=5 (follows from the minimal polynomial of the matrix M).

MAPLE

with(linalg): M[1]:=matrix(5, 5, [5, 3, 2, 1, 1, 3, 2, 1, 1, 0, 2, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0]): for n from 2 to 22 do M[n]:=multiply(M[1], M[n-1]) od: 0, seq(M[n][1, 5], n=1..22);

a[0]:=0: a[1]:=1: a[2]:=5: a[3]:=40: a[4]:=315: for n from 5 to 22 do a[n]:=8*a[n-1]-6*a[n-3]+a[n-5] od: seq(a[n], n=0..22);

MATHEMATICA

M = {{5, 3, 2, 1, 1}, {3, 2, 1, 1, 0}, {2, 1, 1, 0, 0}, {1, 1, 0, 0, 0}, {1, 0, 0, 0, 0}}; v[1] = {0, 0, 0, 0, 1}; v[n_] := v[n] = M.v[n - 1]; a1 = Table[v[n][[1]], {n, 1, 50}]

CROSSREFS

Cf. A122099, A122100.

Sequence in context: A144069 A073505 A145841 * A067412 A078846 A027259

Adjacent sequences:  A123940 A123941 A123942 * A123944 A123945 A123946

KEYWORD

nonn

AUTHOR

Roger Bagula and Gary Adamson (rlbagulatftn(AT)yahoo.com), Oct 25 2006

EXTENSIONS

Edited by N. J. A. Sloane (njas(AT)research.att.com), Dec 04 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 17 14:50 EST 2012. Contains 206050 sequences.