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

2, 5, 7, 11, 13, 23, 53, 61, 71, 79, 137, 139, 193, 229, 239, 251, 293, 317, 373, 433, 523, 599, 601, 683, 727, 859, 877, 887, 911, 991, 1009, 1163, 1229, 1297, 1303, 1429, 1481, 1483, 1789, 1801, 1871, 1999, 2011

Offset: 1

Author

Derek Orr, Jan 19 2014

Keywords

• nonn

Comments

Primes in A126424.

Examples

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

Links

• Harvey P. Dale, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A049408.

Programs

Mathematica

Select[Prime[Range[400]],PrimeQ[#^4-#-1]&] (* Harvey P. Dale, Jan 20 2019 *)

PARI

s=[]; forprime(p=2, 3000, if(isprime(p^4-p-1), s=concat(s, p))); s \\ Colin Barker, Jan 19 2014

Python

import sympy
from sympy import isprime
{print(p) for p in range(10**4) if isprime(p**4-p-1) and isprime(p)}