OFFSET
1,1
EXAMPLE
53 is in the sequence because it is prime, nextprime(53) = 59 and (53+59)/2 = 56 = 8*7, an oblong number.
PROG
(PARI){for(i=3, 10^5, if(isprime(i), k=(i+nextprime(i+1))/4; if(issquare(8*k+1), print1(i, ", "))))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, May 12 2014
STATUS
editing