OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
2^3 +- 3 -> ( 5,11) (primes);
3^3 +- 4 -> (23,31) (primes).
MAPLE
a:=proc(n): if isprime(ithprime(n)^3+ithprime(n)+1)=true and isprime(ithprime(n)^3-ithprime(n)-1)=true then ithprime(n) else end if end proc: seq(a(n), n= 1..4000); # Emeric Deutsch, Apr 29 2008
MATHEMATICA
Select[Prime[Range[900]], PrimeQ[ #^3-(#+1)]&&PrimeQ[ #^3+(#+1)]&]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vladimir Joseph Stephan Orlovsky, Apr 21 2008
EXTENSIONS
More terms from Stefan Steinerberger and Emeric Deutsch, Apr 26 2008
STATUS
approved