OFFSET
0,1
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, 2558, 0, 0, 0, 0, 0, 0, 0, -1).
FORMULA
G.f.: (26 +27*x +53*x^2 +80*x^3 +373*x^4 +453*x^5 +826*x^6 +1279*x^7 +826*x^8 -453*x^9 +373*x^10 -80*x^11 +53*x^12 -27*x^13 +26*x^14 -x^15)/(1 -2558*x^8 +x^16). - Vincenzo Librandi, Nov 22 2013
a(n) = 2558*a(n-8) - a(n-16). - Vincenzo Librandi, Nov 22 2013
MATHEMATICA
Numerator[Convergents[Sqrt[710], 30]] (* or *) CoefficientList[Series[(26 + 27 x + 53 x^2 + 80 x^3 + 373 x^4 + 453 x^5 + 826 x^6 + 1279 x^7 + 826 x^8 - 453 x^9 + 373 x^10 - 80 x^11 + 53 x^12 - 27 x^13 + 26 x^14 - x^15)/(1 - 2558 x^8 + x^16), {x, 0, 30}], x] (* Vincenzo Librandi, Nov 22 2013 *)
PROG
(Magma) I:=[26, 27, 53, 80, 373, 453, 826, 1279, 67334, 68613, 135947, 204560, 954187, 1158747, 2112934, 3271681]; [n le 16 select I[n] else 2558*Self(n-8)-Self(n-16): n in [1..30]]; // Vincenzo Librandi, Nov 22 2013
CROSSREFS
KEYWORD
nonn,cofr,frac,easy
AUTHOR
STATUS
approved