OFFSET
1,3
MATHEMATICA
Reap[For[n = 1, n <= 200, n++, If[PrimeQ[p = 6 n + 1], s = Solve[x^2 + x y + y^2 == p && 0 < x < y, {x, y}, Integers]; Sow[x /. s[[1]]]]]][[2, 1]] (* Jean-François Alcover, Mar 07 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Nov 03 2003
EXTENSIONS
More terms from Ray Chandler, Nov 04 2003
STATUS
approved