OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..900
EXAMPLE
2 is prime and 2^4+2+1 (19) and 2^4+2-1 (17) are twin primes. Thus, 2 is a member of this sequence.
MATHEMATICA
Select[Prime[Range[100000]], PrimeQ[#^4 + # - 1] && PrimeQ[#^4 + # + 1]&] (* Vincenzo Librandi, Feb 14 2014 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(p) for p in range(10**6) if isprime(p) and isprime(p**4+p+1) and isprime(p**4+p-1)}
(Magma) [p: p in PrimesUpTo(90000) | IsPrime(p^4+p-1) and IsPrime(p^4+p+1)]; // Vincenzo Librandi, Feb 14 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Feb 01 2014
STATUS
approved