%I #18 Jan 07 2025 16:32:33
%S 2,5,9,11,12,26,44,47,62,69,71,89,125,140,147,179,219,254,264,285,294,
%T 312,317,326,341,344,384,407,461,495,659,680,714,740,837,845,861,866,
%U 867,957,989,1071,1079,1152,1215,1310,1437,1481,1499,1511,1530,1577
%N Integers k such that (k^3 + k^2) -+ 1 are primes.
%H Robert Israel, <a href="/A154732/b154732.txt">Table of n, a(n) for n = 1..10000</a>
%e 2^3 + 2^2 = 12 -+ 1 = 11 and 13 (both prime).
%p select(k -> andmap(isprime,[k^3+k^2-1,k^3+k^2+1]), [$1..10000]); # _Robert Israel_, Jan 07 2025
%t lst={};Do[k=n^3+n^2;If[PrimeQ[k-1]&&PrimeQ[k+1],AppendTo[lst,n]],{n,8!}];lst
%t Select[Range[3000], PrimeQ[#^3 + #^2 - 1] && PrimeQ[#^3 + #^2 + 1] &] (* _Vincenzo Librandi_, Dec 26 2015 *)
%o (Magma) [n: n in [1..5*10^3] |IsPrime(n^3+n^2-1) and IsPrime(n^3+n^2+1)]; // _Vincenzo Librandi_, Dec 26 2015
%o (PARI) isok(n) = isprime(n^3+n^2+1) && isprime(n^3+n^2-1); \\ _Michel Marcus_, Dec 27 2015
%Y Cf. A111503, A154731.
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Jan 14 2009