OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,646,0,0,0,-1).
FORMULA
G.f.: -(x^2-x-1)*(x^4+18*x^2+1) / (x^8-646*x^4+1). - Colin Barker, Nov 19 2013
a(n) = 646*a(n-4) - a(n-8) for n>7. - Vincenzo Librandi, Dec 21 2013
MATHEMATICA
Denominator[Convergents[Sqrt[322], 30]] (* Harvey P. Dale, Apr 25 2012 *)
CoefficientList[Series[-(x^2 - x - 1) (x^4 + 18 x^2 + 1)/(x^8 - 646 x^4 + 1), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 21 2013 *)
PROG
(Magma) I:=[1, 1, 17, 18, 629, 647, 10981, 11628]; [n le 8 select I[n] else 646*Self(n-4)-Self(n-8): n in [1..70]]; // Vincenzo Librandi, Dec 21 2013
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
EXTENSIONS
More terms from Colin Barker, Nov 19 2013
STATUS
approved