OFFSET
1,1
MATHEMATICA
op40Q[n_]:=Module[{s=Sqrt[n], cf, len}, cf=If[IntegerQ[s], {0, 0}, ContinuedFraction[ s][[2]]]; len=Length[cf]; OddQ[len]&&cf[[(len+1)/2]] ==40]; Select[Range[500000], op40Q] (* Harvey P. Dale, Dec 06 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Corrected and extended by Harvey P. Dale, Dec 06 2015
STATUS
approved