OFFSET
1,1
COMMENTS
x^2+n*y^2=(+/-)2^s where s is 0 or 1.
Definition: Unity is singular when GCD[n,y]<>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, n]]], {n, 2, 330}]; cr
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 10 2011
STATUS
approved