OFFSET
1,1
COMMENTS
Primes in the sequence A236072.
EXAMPLE
6971 is prime, 6971^4 - 6971 - 1 is prime, and 6971^4 + 6971 + 1 is prime. So 6971 is a member of this sequence.
PROG
(Python)
import sympy
from sympy import isprime
{print(p) for p in range(10**5) if isprime(p**4+p+1) and isprime(p**4-p-1) and isprime(p)}
(PARI) s=[]; forprime(p=2, 55871, if(isprime(p^4+p+1)&&isprime(p^4-p-1), s=concat(s, p))); s \\ Colin Barker, Jan 19 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jan 19 2014
STATUS
approved