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, 23326, 0, 0, 0, 0, 0, 0, 0, -1).
FORMULA
G.f.: (30 +31*x +185*x^2 +216*x^3 +1481*x^4 +1697*x^5 +9966*x^6 +11663*x^7 +9966*x^8 -1697*x^9 +1481*x^10 -216*x^11 +185*x^12 -31*x^13 +30*x^14 -x^15)/(1 -23326*x^8 +x^16). - Vincenzo Librandi, Dec 07 2013
a(n) = 23326*a(n-8) - a(n-16). - Vincenzo Librandi, Dec 07 2013
MATHEMATICA
Numerator[Convergents[Sqrt[952], 30]] (* or *)CoefficientList[Series[(30 + 31 x + 185 x^2 + 216 x^3 + 1481 x^4 + 1697 x^5 + 9966 x^6 + 11663 x^7 + 9966 x^8 - 1697 x^9 + 1481 x^10 - 216 x^11 + 185 x^12 - 31 x^13 + 30 x^14 - x^15)/(1 - 23326 x^8 + x^16), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 07 2013 *)
PROG
(Magma) I:=[30, 31, 185, 216, 1481, 1697, 9966, 11663, 709746, 721409, 4316791, 5038200, 34545991, 39584191, 232466946, 272051137]; [n le 16 select I[n] else 23326*Self(n-8)-Self(n-16): n in [1..30]]; // Vincenzo Librandi, Dec 07 2013
CROSSREFS
KEYWORD
nonn,cofr,frac,easy
AUTHOR
EXTENSIONS
More terms from Vincenzo Librandi, Dec 07 2013
STATUS
approved