A249580 -id:A249580 - cp's OEIS frontend

A249579 List of quadruples (r,s,t,u): the matrix M = [[4,12,9][2,5,3][1,2,1]] is raised to successive powers, then (r,s,t,u) are the square roots of M[3,1], M[3,3], M[1,1], M[1,3] respectively.

0, 1, 1, 0, 1, 1, 2, 3, 3, 4, 7, 9, 10, 13, 23, 30, 33, 43, 76, 99, 109, 142, 251, 327, 360, 469, 829, 1080, 1189, 1549, 2738, 3567, 3927, 5116, 9043, 11781, 12970, 16897, 29867, 38910, 42837, 55807, 98644, 128511, 141481, 184318, 325799, 424443, 467280

Offset: 0

Russell Walsmith, Nov 02 2014

			M^0 = [[1,0,0][0,1,0][0,0,1]]: r = sqrt(M[3,1]) = a(0) = 0, s = sqrt(M[3,3]) = a(1) = 1, t = sqrt(M[1,1]) = a(2) = 1, u = sqrt(M[1,3])u = a(3) = 0.
M^2 = [[49, 126, 81][21, 55, 36][9, 24, 16]]: r = sqrt(M[3, 1]) = a(8) = 3, s = sqrt(M[3, 3]) = a(9) = 4, t = sqrt(M[1, 1]) = a(10) = 7, u = sqrt(M[1, 3]) = a(11) = 9.

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,0,3,0,0,0,1).

Cf. A249580
a(4n) = A006190
a(4n+2) = A052924.

I:=[0,1,1,0,1,1,2,3]; [n le 8 select I[n] else 3*Self(n-4)+Self(n-8): n in [1..50]]; // Vincenzo Librandi, Nov 14 2014

CoefficientList[Series[- x (3 x^6 - x^5 - 2 x^4 + x^3 + x + 1) / (x^8 + 3 x^4 - 1), {x, 0, 50}], x] (* Vincenzo Librandi, Nov 14 2014 *)

concat(0, Vec(-x*(3*x^6-x^5-2*x^4+x^3+x+1)/(x^8+3*x^4-1) + O(x^100))) \\ Colin Barker, Nov 13 2014

Some identities:
a4(n-1) + a(4n) = a(4n+1),
a(4n) + a(4n+1) = a(4n+2),
3a(4n) = a(4n+3).
a(n) = 3*a(n-4)+a(n-8). - Colin Barker, Nov 13 2014
G.f.: -x*(3*x^6-x^5-2*x^4+x^3+x+1) / (x^8+3*x^4-1). - Colin Barker, Nov 13 2014

More terms from Colin Barker, Nov 13 2014

A249581 List of quadruples (r,s,t,u): the matrix M = [[9,24,16][3,10,8][1,4,4]] is raised to successive powers, then (r,s,t,u) are the square roots of M[3,1], M[3,3], M[1,1], M[1,3] respectively.

Original entry on oeis.org

0, 1, 1, 0, 1, 2, 3, 4, 5, 8, 13, 20, 23, 36, 59, 92, 105, 164, 269, 420, 479, 748, 1227, 1916, 2185, 3412, 5597, 8740, 9967, 15564, 25531, 39868, 45465, 70996, 116461, 181860, 207391, 323852, 531243, 829564, 946025, 1477268, 2423293, 3784100

Offset: 0

Russell Walsmith, Nov 03 2014

The general form of these matrices is [[t^2,2tu,u^2][rt,st+ru,su][r^2,2rs,s^2]]. Different symmetries have different properties.

Iff |r * u - s * t| = 1 then terms to the left of a(0) are all integers.

			M^0 = [[1,0,0][0,1,0][0,0,1]]. r = sqrt(M[3,1]) = a(0) = 0; s = sqrt(M[3,3]) = a(1) = 1; t = sqrt(M[1,1]) = a(2) = 1; u = sqrt(M[1,3]) = a(3) = 0.
M^1 = [[9,24,16][3,10,8][1,4,4]]. r = sqrt(M[3,1]) = a(4) = 1; s = sqrt(M[3,3]) = a(5) = 2; t = sqrt(M[1,1]) = a(6) = 3; u = sqrt(M[1,3]) = a(7) = 4.

Colin Barker, Table of n, a(n) for n = 0..1000
Russell Walsmith, DCL-Chemy III: Hyper-Quadratics
Index entries for linear recurrences with constant coefficients, signature (0,0,0,5,0,0,0,-2).

Cf. A249579, A249580, A147722, A052984.

LinearRecurrence[{0,0,0,5,0,0,0,-2},{0,1,1,0,1,2,3,4},50] (* Harvey P. Dale, Aug 01 2016 *)

concat(0, Vec(x*(4*x^6-2*x^5-3*x^4+x^3+x+1)/(2*x^8-5*x^4+1) + O(x^100))) \\ Colin Barker, Nov 04 2014

a(4n) + a(4n + 1) = a(4n + 2).
a(4n + 1) + a(4n + 2) + a(4n + 3) - a(4n) = a(4n + 5)
4a(4n) = a(4n+3).
a(4n+1) = A147722(n), a(4n+2) = A052984(n).
a(n) = 5*a(n-4)-2*a(n-8). - Colin Barker, Nov 04 2014
G.f.: x*(4*x^6-2*x^5-3*x^4+x^3+x+1) / (2*x^8-5*x^4+1). - Colin Barker, Nov 04 2014

cp's OEIS Frontend

A249579 List of quadruples (r,s,t,u): the matrix M = [[4,12,9][2,5,3][1,2,1]] is raised to successive powers, then (r,s,t,u) are the square roots of M[3,1], M[3,3], M[1,1], M[1,3] respectively.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions