%I #7 Apr 01 2014 17:00:25
%S 3,11,17,41,59,71,101,113,179,233,293,347,389,449,461,503,521,617,641,
%T 683,797,953,1319,1439,1487,1493,1823,1877,1973,2087,2339,2351,2633,
%U 2663,2789,2801,2909,2927,2957,2963,2999,3011,3167,3467,3527,3677,3851,3881,3923
%N Primes p such that p^4-p^3-p^2-p-1 is prime.
%e 3^4-3^3-3^2-3-1 = 41 is prime. Thus, 3 is a member of this sequence.
%o (Python)
%o import sympy
%o from sympy import isprime
%o {print(p) for p in range(10**4) if isprime(p**4-p**3-p**2-p-1) and isprime(p)}
%o (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
%Y Cf. A173179.
%K nonn
%O 1,1
%A _Derek Orr_, Mar 31 2014
|