A139505 Primes of the form x^2 + 25x*y + y^2 for x and y nonnegative.
151, 163, 307, 397, 409, 541, 547, 601, 673, 811, 823, 859, 967, 997, 1153, 1231, 1237, 1327, 1567, 1669, 1741, 1879, 2083, 2143, 2281, 2293, 2557, 2677, 2707, 2833, 2971, 3037, 3259, 3313, 3433, 3877, 4003, 4129, 4153, 4603, 4639, 4861, 4957, 5101, 5227
Offset: 1
Keywords
Links
- Harvey P. Dale, Table of n, a(n) for n = 1..5000
- N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
Crossrefs
Programs
-
Mathematica
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*) With[{nn=80},Select[Union[#[[1]]^2+25#[[1]]#[[2]]+#[[2]]^2&/@Tuples[ Range[ 0,nn],2]],PrimeQ[#]&
Harvey P. Dale, Feb 10 2020 *)