%I #11 Feb 07 2019 02:12:16
%S 2,3,8,9,15,18,21,29,41,48,60,69,71,104,108,126,155,168,186,191,221,
%T 239,261,276,360,366,381,414,428,429,480,518,570,630,701,749,755,825,
%U 836,879,1140,1170,1173,1220,1290,1395,1449,1455,1466,1485,1605,1611,1655
%N Positive numbers n such that n^3 +- (n+1) are primes.
%H Vincenzo Librandi, <a href="/A137471/b137471.txt">Table of n, a(n) for n = 1..1000</a>
%e 2^3 +- 3 -> (5, 11) (primes);
%e 3^3 +- 4 -> (23, 31) (primes).
%p a:=proc(n) if isprime(n^3+n+1)=true and isprime(n^3-n-1)=true then n else end if end proc: seq(a(n),n=1..1200); # _Emeric Deutsch_, Apr 29 2008
%t Select[Range[900],PrimeQ[ #^3-(#+1)]&&PrimeQ[ #^3+(#+1)]&]
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Apr 21 2008
%E More terms from _Stefan Steinerberger_ and _Emeric Deutsch_, Apr 26 2008
|