A042528 Numerators of continued fraction convergents to sqrt(793).

28, 169, 704, 4393, 246712, 1484665, 6185372, 38596897, 2167611604, 13044266521, 54344677688, 339112332649, 19044635306032, 114606924168841, 477472331981396, 2979440916057217, 167326163631185548, 1006936422703170505

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,8786,0,0,0,-1).

Formula

G.f.: (28 +169*x +704*x^2 +4393*x^3 +704*x^4 -169*x^5 +28*x^6 -x^7)/(1 -8786*x^4 +x^8). - Vincenzo Librandi, Nov 26 2013

a(n) = 8786*a(n-4) - a(n-8). - Vincenzo Librandi, Nov 26 2013

Mathematica

Numerator[Convergents[Sqrt[793], 30]] (* or *) CoefficientList[Series[(28 + 169 x + 704 x^2 + 4393 x^3 + 704 x^4 - 169 x^5 + 28 x^6 - x^7)/(1 - 8786 x^4 + x^8), {x, 0, 30}], x] (* Vincenzo Librandi, Nov 25 2013 *)

Prog

(Magma) I:=[28, 169, 704, 4393, 246712, 1484665, 6185372, 38596897]; [n le 8 select I[n] else 8786*Self(n-4)-Self(n-8): n in [1..30]]; // Vincenzo Librandi, Nov 26 2013

Crossrefs

Cf. A042529.

Sequence in context: A305270 A366954 A129136 * A345215 A219821 A341160

Adjacent sequences: A042525 A042526 A042527 * A042529 A042530 A042531

Keyword

nonn,cofr,frac,easy

Author

N. J. A. Sloane.

Status

approved