OFFSET
1,1
EXAMPLE
33^3-33^2+1=34849 prime
34^3-34^2+1=38149 prime
35^3-35^2+1=41651 prime
36^3-36^2+1=45361 prime so n(1)=33
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