OFFSET
1,1
LINKS
Abhiram R Devesh, Table of n, a(n) for n = 1..10000
EXAMPLE
Prime p = 3 is in this sequence as p^5 + p^3 + p - 4 = 269 (prime).
Prime p = 5 is in this sequence as p^5 + p^3 + p - 4 = 3251 (prime).
MATHEMATICA
Select[Prime[Range[500]], PrimeQ[#^5+#^3+#-4]&] (* Harvey P. Dale, Jul 03 2015 *)
PROG
(Python)
import sympy.ntheory as snt
p=1
while p>0:
....p=snt.nextprime(p)
....pp=p+(p**3)+(p**5)-4
....if snt.isprime(pp) == True:
........print(p, pp)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Abhiram R Devesh, Jun 14 2014
STATUS
approved