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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 78960998, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1).
MATHEMATICA
Numerator[Convergents[Sqrt[1000], 30]] (* Harvey P. Dale, Oct 29 2013 *)
PROG
(PARI) A42936=contfracpnqn(c=contfrac(sqrt(1000)), #c)[1, ][^-1] \\ Discards possibly incorrect last term. NB: a(n)=A42936[n+1]. Could be extended using: {A42936=concat(A42936, 78960998*A42936[-18..-1]-A42936[-36..-19])} \\ M. F. Hasler, Nov 01 2019
(PARI) A042936(n)={[A42936[n%18+i]|i<-[1, 19]]*([0, -1; 1, 78960998]^(n\18))[, 1]} \\ Faster but longer with n=divrem(n, 18), M. F. Hasler, Nov 01 2019
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
STATUS
approved
