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.
EXAMPLE
a(1)=6 because unity of quadratic field Q(6) is 5+2*Sqrt[6] and GCD[2,6]=2 <>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 (*Artur Jasinski*)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 10 2011
STATUS
approved