OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
975^4-975+1 (903687889651) and 975^4-975-1 (903687889649) are twin primes. Thus, 975 is a member of this sequence.
MATHEMATICA
Select[Range[3000], AllTrue[#^4-#+{1, -1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 27 2014 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(n) for n in range(10**4) if isprime(n**4-n+1) and isprime(n**4-n-1)}
(PARI)
s=[]; for(n=1, 3000, if(isprime(n^4-n+1)&&isprime(n^4-n-+1), s=concat(s, n))); s \\ Colin Barker, Jan 31 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jan 30 2014
STATUS
approved