OFFSET
3,1
COMMENTS
No quadratic number with a fully periodical continued fraction of period 1 can be written as (a+b*sqrt(n))/c with n allowed to have square factors.
Subsequence of the squarefree terms is A087642.
MATHEMATICA
cr = {}; Do[If[IntegerQ[Sqrt[n]], , kk = NumberFieldFundamentalUnits[Sqrt[n]]; d1 = kk[[1]][[2]][[1]]; d2 = kk[[1]][[1]] kk[[1]][[2]][[2]]; d3 = Expand[(d1 + d2) (d1 - d2)]; If[d3 == 1, AppendTo[cr, n]]], {n, 2, 1000}]; cr (* Artur Jasinski, Oct 10 2011 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Thomas Baruchel, Sep 16 2003
EXTENSIONS
Edited by Max Alekseyev, May 03 2010
STATUS
approved