OFFSET
1,1
MATHEMATICA
cr = {}; Do[If[IntegerQ[Sqrt[n]], , kk = NumberFieldFundamentalUnits[Sqrt[n]]; d1 = kk[[1]][[2]][[1]]; d2 = kk[[1]][[1]] kk[[1]][[2]][[2]]; d4 = Numerator[d2/Sqrt[n]]; If[GCD[d4, n] == 1, , AppendTo[cr, GCD[d4, n]]]], {n, 2, 20000}]; cr
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 11 2011
STATUS
approved