A236044 - OEIS

%I #6 Sep 27 2014 17:50:14

%S 2,5,11,47,89,107,131,191,197,239,347,641,701,839,941,977,1049,1061,

%T 1289,1709,1847,1871,1949,1979,2129,2789,2897,3371,3557,3719,3761,

%U 3917,4001,4091,4211,4289,4337,4397,4547,4751

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

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

%e 941 is prime and 941^4 + 941 + 1 is prime, thus 941 is a member of this sequence.

%t Select[Prime[Range[700]],PrimeQ[#^4+#+1]&] (* _Harvey P. Dale_, Sep 27 2014 *)

%o (Python)

%o import sympy

%o from sympy import isprime

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

%Y Cf. A049408.

%K nonn

%O 1,1

%A _Derek Orr_, Jan 18 2014