cp's OEIS Frontend

A236071 Primes p such that p^4 - p - 1 is prime.

%I A236071 #12 Jan 20 2019 15:29:18
%S A236071 2,5,7,11,13,23,53,61,71,79,137,139,193,229,239,251,293,317,373,433,
%T A236071 523,599,601,683,727,859,877,887,911,991,1009,1163,1229,1297,1303,
%U A236071 1429,1481,1483,1789,1801,1871,1999,2011
%N A236071 Primes p such that p^4 - p - 1 is prime.
%C A236071 Primes in A126424.
%H A236071 Harvey P. Dale, <a href="/A236071/b236071.txt">Table of n, a(n) for n = 1..1000</a>
%e A236071 139 is prime and 139^4 - 139 - 1 is prime, so 139 is a member of this sequence.
%t A236071 Select[Prime[Range[400]],PrimeQ[#^4-#-1]&] (* _Harvey P. Dale_, Jan 20 2019 *)
%o A236071 (Python)
%o A236071 import sympy
%o A236071 from sympy import isprime
%o A236071 {print(p) for p in range(10**4) if isprime(p**4-p-1) and isprime(p)}
%o A236071 (PARI)
%o A236071 s=[]; forprime(p=2, 3000, if(isprime(p^4-p-1), s=concat(s, p))); s \\ _Colin Barker_, Jan 19 2014
%Y A236071 Cf. A049408.
%K A236071 nonn
%O A236071 1,1
%A A236071 _Derek Orr_, Jan 19 2014