login
A137472
Primes p such that p^3 +- (p+1) are primes.
1
2, 3, 29, 41, 71, 191, 239, 701, 2309, 2531, 4229, 5261, 6491, 6971, 8009, 10559, 11279, 11393, 11789, 12203, 12329, 14519, 17021, 18041, 18719, 20369, 23561, 23633, 28499, 33329, 34499, 37409, 38453, 41051, 41843, 42359, 44159
OFFSET
1,1
LINKS
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
Sequence in context: A059453 A235481 A214889 * A065932 A143882 A108657
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Stefan Steinerberger and Emeric Deutsch, Apr 26 2008
STATUS
approved