OFFSET
0,7
LINKS
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).
FORMULA
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
EXAMPLE
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.
MATHEMATICA
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 204 *)
PROG
(PARI) 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
(Magma) 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
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
Russell Walsmith, Nov 02 2014
EXTENSIONS
More terms from Colin Barker, Nov 13 2014
STATUS
approved