OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
3^4-3^3-3^2-3-1 = 41 is prime. Thus, 3 is a member of this sequence.
MATHEMATICA
Select[Prime[Range[600]], PrimeQ[#^4-#^3-#^2-#-1]&] (* Harvey P. Dale, Nov 23 2024 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(p) for p in range(10**4) if isprime(p**4-p**3-p**2-p-1) and isprime(p)}
(PARI) s=[]; forprime(p=2, 4000, if(isprime(p^4-p^3-p^2-p-1), s=concat(s, p))); s \\ Colin Barker, Apr 01 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Mar 31 2014
STATUS
approved