login
Values of y, where x^2 + xy + y^2=p (x<y) is a prime of the form 6n + 1 (=A002476).
6

%I #10 Mar 07 2020 07:58:12

%S 2,3,3,5,4,6,5,7,8,7,8,9,7,7,10,9,12,11,11,9,13,14,11,12,15,10,12,13,

%T 17,16,11,13,17,13,17,15,12,15,20,13,18,17,21,21,18,17,21,14,21,19,24,

%U 23,19,22,15,18,20,21,19,25,18,19,23,21,27,17,27,25,19,20,27,23,28,21,26

%N Values of y, where x^2 + xy + y^2=p (x<y) is a prime of the form 6n + 1 (=A002476).

%t 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[y /. s[[1]]]]]][[2, 1]] (* _Jean-François Alcover_, Mar 07 2020 *)

%Y Cf. A002476, A088242, A088243, A088296, A088298, A088299, A088977.

%K nonn

%O 1,1

%A _Lekraj Beedassy_, Nov 03 2003

%E More terms from _Ray Chandler_, Nov 04 2003