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,0,0,0,0,15486,0,0,0,0,0,0,0,-1).
FORMULA
G.f.: -(x^14 -x^13 +2*x^12 -7*x^11 +86*x^10 -265*x^9 +351*x^8 -616*x^7 -351*x^6 -265*x^5 -86*x^4 -7*x^3 -2*x^2 -x -1) / (x^16 -15486*x^8 +1). - Colin Barker, Nov 14 2013
a(n) = 15486*a(n-8) - a(n-16). - Vincenzo Librandi, Dec 14 2013
MATHEMATICA
Denominator[Convergents[Sqrt[158], 30]] (* Harvey P. Dale, Aug 17 2013 *)
CoefficientList[Series[-(x^14 - x^13 + 2 x^12 - 7 x^11 + 86 x^10 - 265 x^9 + 351 x^8 - 616 x^7 - 351 x^6 - 265 x^5 - 86 x^4 - 7 x^3 - 2 x^2 - x - 1)/(x^16 - 15486 x^8 + 1), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 14 2013*)
PROG
(Magma) I:=[1, 1, 2, 7, 86, 265, 351, 616, 15135, 15751, 30886, 108409, 1331794, 4103791, 5435585, 9539376]; [n le 16 select I[n] else 15486*Self(n-8)-Self(n-16): n in [1..40]]; // Vincenzo Librandi, Dec 14 2013
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
EXTENSIONS
More terms from Colin Barker, Nov 14 2013
STATUS
approved