OFFSET
0,3
COMMENTS
The limiting ratio a(n)/a(n-1) alternates between 2.2433981132057377... and 4.554247641507091...
LINKS
Eric Weisstein's World of Mathematics, Symplectic Group.
Index entries for linear recurrences with constant coefficients, signature (0,11,0,-8).
FORMULA
Given the four matrices m(1) = {{1, 0}, {0, 1}}, m(2) = {{1, 1}, {0, 1}}, m(3) = {{1, 0}, {1, 1}}, and m(4) = m(1)+m(2)+m(3), repeatedly multiply the vector {0,1} with m(1), m(2) etc in a round-robin periodic fashion. a(n) is the top value of the vector after each second multiplication.
G.f.: x*(1 + 5*x + x^2) / (1 - 11*x^2 + 8*x^4). - R. J. Mathar, Feb 10 2011
a(n) = 5*a(n-1) - a(n-2) if n is even. a(n) = 4*a(n-1) - 8*a(n-2) if n is odd. - R. J. Mathar, Jun 18 2014
EXAMPLE
Multiplication with m1 updates the vector to {1,1}, then with m2 to {5,7}, then with m3 to {5,7}, then with m4 to {12,7}, then with m1 {12,19}, then with m2 {55,69}, ... The leading components are 1, 1, 5, 5, 12, 12, 55, 55, ... and each second defines the sequence.
MAPLE
with(LinearAlgebra) ; Am := [ [1, 0, 0, 1], [1, 1, 0, 1], [1, 0, 1, 1], [3, 1, 1, 3] ]: v := Vector[column](2, [0, 1]) :
for l from 0 to 10 do
A := Matrix(2, 2, op(1+(l mod 4), Am) ) ; v := A.v;
if type(l, 'even') then print(v[1]) ; end if;
end do: # R. J. Mathar, Feb 10 2011
MATHEMATICA
a[1] = {{1, 0}, {0, 1}};
a[2] = {{1, 1}, {0, 1}};
a[3] = {{1, 0}, {1, 1}};
a[4] = a[1] + a[2] + a[3]
v[0] = {0, 1};
v[n_] := v[n] = a[1 + Mod[n, 4]].v[n - 1]
Table[v[n][[1]], {n, 0, 60, 2}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Roger L. Bagula, Aug 05 2010
STATUS
approved