login
A020418
Numbers k such that the continued fraction for sqrt(k) has period 79.
1
1621, 3001, 3637, 4153, 5113, 5653, 7114, 8026, 9161, 12538, 12569, 13261, 14177, 18058, 19373, 21625, 25733, 25981, 26297, 28277, 28394, 29674, 30122, 32201, 32393, 32861, 33773, 34801, 35746, 36137, 36185, 36257, 36994, 39049, 39701, 41429, 43721
OFFSET
1,1
MATHEMATICA
cf79Q[n_]:=Module[{c=Sqrt[n]}, !IntegerQ[c]&&Length[ContinuedFraction[c][[2]]]==79]; Select[Range[50000], cf79Q] (* Harvey P. Dale, Dec 23 2012 *)
PROG
(Python)
from sympy import continued_fraction_periodic
A020418_list = [n for n, s in ((i, continued_fraction_periodic(0, 1, i)[-1]) for i in range(1, 10**5)) if isinstance(s, list) and len(s) == 79] # Chai Wah Wu, Jun 08 2017
CROSSREFS
Sequence in context: A269212 A269095 A269274 * A347039 A224722 A045288
KEYWORD
nonn
STATUS
approved