1,1

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

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

Select[Prime[Range[700]], PrimeQ[#^4+#+1]&] (* Harvey P. Dale, Sep 27 2014 *)

(Python)

import sympy

from sympy import isprime

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

Cf. A049408.

Sequence in context: A153989 A160859 A106887 * A285809 A089609 A277041

Adjacent sequences: A236041 A236042 A236043 * A236045 A236046 A236047

nonn

Derek Orr, Jan 18 2014

approved