OFFSET
1,1
COMMENTS
Primes in A126421.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
269 is prime and 269^3 - 269 - 1 is also prime. So, 269 is a member of this sequence.
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[#^3-#-1]&] (* Harvey P. Dale, Nov 17 2014 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(p) for p in range(10**4) if isprime(p) and isprime(p**3-p-1)}
(PARI)
s=[]; forprime(p=2, 2000, if(isprime(p^3-p-1), s=concat(s, p))); s \\ Colin Barker, Jan 20 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jan 19 2014
STATUS
approved