%I #9 Nov 17 2014 14:41:02

%S 2,3,11,23,29,31,41,59,71,113,151,163,191,239,241,269,359,431,433,499,

%T 503,521,541,563,661,683,701,751,773,829,883,983,1039,1259,1483,1499,

%U 1511,1549,1571,1609,1693,1721,1759,1913

%N Primes p such that p^3 - p - 1 is prime.

%C Primes in A126421.

%H Harvey P. Dale, <a href="/A236168/b236168.txt">Table of n, a(n) for n = 1..1000</a>

%e 269 is prime and 269^3 - 269 - 1 is also prime. So, 269 is a member of this sequence.

%t Select[Prime[Range[300]],PrimeQ[#^3-#-1]&] (* _Harvey P. Dale_, Nov 17 2014 *)

%o (Python)

%o import sympy

%o from sympy import isprime

%o {print(p) for p in range(10**4) if isprime(p) and isprime(p**3-p-1)}

%o (PARI)

%o s=[]; forprime(p=2, 2000, if(isprime(p^3-p-1), s=concat(s, p))); s \\ _Colin Barker_, Jan 20 2014

%Y Cf. A126421.

%K nonn

%O 1,1

%A _Derek Orr_, Jan 19 2014