OFFSET
1,1
COMMENTS
(97*m)^2+2*m for m >= 1 are terms of the sequence (see comment in A031749). The term 89453959 is not of this form. - Chai Wah Wu, Jun 19 2016
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
MATHEMATICA
cf97Q[n_]:=Module[{s=Sqrt[n]}, If[IntegerQ[s], 1, Min[ContinuedFraction[s][[2]]]]==97]; Select[Range[738*10^4], cf97Q] (* Harvey P. Dale, Nov 20 2018 *)
PROG
(Python)
from sympy import continued_fraction_periodic
A031775_list = [n for n, d in ((n, continued_fraction_periodic(0, 1, n)[-1]) for n in range(1, 10**5)) if isinstance(d, list) and min(d) == 97] # Chai Wah Wu, Jun 10 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved