OFFSET
1,1
LINKS
Daniel Starodubtsev, Table of n, a(n) for n = 1..10000
EXAMPLE
110685^3+110685+1 (1356020665779811), 110685^3+110685-1 (1356020665779809), 110685^3-110685+1 (1356020665558441) and 110685^3-110685-1 (1356020665558439) are all prime. Thus 110685 is a member of this sequence.
PROG
(Python)
import sympy
from sympy import isprime
{print(n) for n in range(10**6) if isprime(n**3+n+1) and isprime(n**3-n+1) and isprime(n**3+n-1) and isprime(n**3-n-1)}
(PARI) for(n=1, 800000, if(isprime(n^3+n+1)&&isprime(n^3-n+1)&&isprime(n^3+n-1)&&isprime(n^3-n-1), print1(n, ", "))) \\ Colin Barker, Jan 31 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jan 30 2014
STATUS
approved