OFFSET
1,1
EXAMPLE
33^3 - 33^2 + 1 = 34849 is prime,
34^3 - 34^2 + 1 = 38149 is prime,
35^3 - 35^2 + 1 = 41651 is prime,
36^3 - 36^2 + 1 = 45361 is prime, so 33 is a term.
MATHEMATICA
fQ[n_] := Block[{j = {0, 1, 2, 3}}, Union@PrimeQ[(n + j)^3 - (n + j)^2 + 1] == {True}]; t = {}; Do[ If[ fQ[n], AppendTo[t, n]], {n, 1005807}] (* Robert G. Wilson v *)
slv[n_]:=x/.Solve[x^3-x^2+1==n, x][[1]]; slv/@(Transpose[Select[ Partition[ Table[ n^3-n^2+1, {n, 980000}], 4, 1], And@@PrimeQ[#]&]][[1]]) (* Harvey P. Dale, May 03 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Nov 16 2005
EXTENSIONS
a(6)-a(32) from Robert G. Wilson v, Nov 18 2005
STATUS
approved