OFFSET
1,1
COMMENTS
2^3-2^2-1=3;2^3+2^2+1=13, 3^3-3^2-1=17;3^3+3^2+1=37, ...
MATHEMATICA
lst={}; Do[c=n^3; s=n^2; p1=c-s-1; p2=c+s+1; If[PrimeQ[p1]&&PrimeQ[p2], AppendTo[lst, s]], {n, 7!}]; lst
#^2&/@Select[Range[1600], AllTrue[#^3+{-#^2-1, #^2+1}, PrimeQ]&] (* Harvey P. Dale, Feb 24 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 30 2009
STATUS
approved