OFFSET
1,1
COMMENTS
Primes of the form 2*k^2 + 2*k - 39.
2*p + h is not verified if h is an odd prime that belongs to A055025 because (2*h-1)/2 is a multiple of 2.
EXAMPLE
a(1) = 73 because 2*73 + 79 = 225, which is a square.
MATHEMATICA
Select[Prime[Range[50000]], IntegerQ[Sqrt[2 # + 79]] &]
PROG
(Magma) [p: p in PrimesUpTo(600000) | IsSquare(2*p+79)];
(PARI) lista(nn) = {forprime(p=2, nn, if(issquare(2*p + 79), print1(p, ", "))); } \\ Altug Alkan, Mar 24 2016
(Python)
from sympy import isprime
from gmpy2 import is_square
for p in range(0, 1000000):
if(is_square(2*p+79) and isprime(p)):print(p)
# Soumil Mandal, Apr 03 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Mar 24 2016
STATUS
approved