OFFSET
1,1
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 + 2y /. s[[1]]]]]][[2, 1]] (* Jean-François Alcover, Mar 09 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Nov 03 2003
EXTENSIONS
More terms from Ray Chandler, Nov 04 2003
STATUS
approved