OFFSET
2,1
MATHEMATICA
cr = {}; ck = {}; 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[ck, GCD[d4, n]]; AppendTo[cr, n]]], {n, 2, 200000}]; aa = {}; Do[If[ck[[n]] == 3, AppendTo[aa, cr[[n]]]], {n, 1, Length[cr]}]; aa
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 11 2011
STATUS
approved