login
A088242
Values of x, where x^2 + xy + y^2=p (x<y) is a prime of the form 6n + 1 (=A002476).
6
1, 1, 2, 1, 3, 1, 4, 2, 1, 3, 3, 2, 5, 6, 3, 5, 1, 3, 4, 7, 2, 1, 6, 5, 1, 9, 7, 6, 1, 3, 10, 8, 3, 9, 4, 7, 11, 8, 1, 11, 5, 7, 1, 2, 7, 9, 4, 13, 5, 8, 1, 3, 9, 5, 14, 11, 9, 8, 11, 3, 13, 12, 7, 10, 1, 15, 2, 6, 14, 13, 4, 10, 3, 13, 7, 17, 3, 7, 9, 13, 8, 11, 16, 15, 6, 3, 12, 17, 7, 9, 1, 3, 16
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 *)
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Nov 03 2003
EXTENSIONS
More terms from Ray Chandler, Nov 04 2003
STATUS
approved