%I #8 Feb 10 2020 13:56:37
%S 151,163,307,397,409,541,547,601,673,811,823,859,967,997,1153,1231,
%T 1237,1327,1567,1669,1741,1879,2083,2143,2281,2293,2557,2677,2707,
%U 2833,2971,3037,3259,3313,3433,3877,4003,4129,4153,4603,4639,4861,4957,5101,5227
%N Primes of the form x^2 + 25x*y + y^2 for x and y nonnegative.
%H Harvey P. Dale, <a href="/A139505/b139505.txt">Table of n, a(n) for n = 1..5000</a>
%H N. J. A. Sloane et al., <a href="https://oeis.org/wiki/Binary_Quadratic_Forms_and_OEIS">Binary Quadratic Forms and OEIS</a> (Index to related sequences, programs, references)
%t a = {}; w = 25; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)
%t With[{nn=80},Select[Union[#[[1]]^2+25#[[1]]#[[2]]+#[[2]]^2&/@Tuples[ Range[ 0,nn],2]],PrimeQ[#]&&#<nn^2&]] (* _Harvey P. Dale_, Feb 10 2020 *)
%Y Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.
%K nonn
%O 1,1
%A _Artur Jasinski_, Apr 24 2008
|