A236071 - OEIS

%I #12 Jan 20 2019 15:29:18

%S 2,5,7,11,13,23,53,61,71,79,137,139,193,229,239,251,293,317,373,433,

%T 523,599,601,683,727,859,877,887,911,991,1009,1163,1229,1297,1303,

%U 1429,1481,1483,1789,1801,1871,1999,2011

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

%C Primes in A126424.

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

%e 139 is prime and 139^4 - 139 - 1 is prime, so 139 is a member of this sequence.

%t Select[Prime[Range[400]],PrimeQ[#^4-#-1]&] (* _Harvey P. Dale_, Jan 20 2019 *)

%o (Python)

%o import sympy

%o from sympy import isprime

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

%o (PARI)

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

%Y Cf. A049408.

%K nonn

%O 1,1

%A _Derek Orr_, Jan 19 2014