OFFSET
1,2
COMMENTS
We allowed k to vary up to 10^7.
MATHEMATICA
nn = 60; t = Table[0, {nn}]; ps = Prime[Range[nn]]; Do[num = n^3 + 1; j = 0; lastP = 0; While[num > 0 && j < nn, j++; p = ps[[j]]; While[Mod[num, p] == 0, lastP = j; num = num/p]]; If[num == 1 && t[[lastP]] == 0, t[[lastP]] = n; Print[{lastP, n}]], {n, 10^7}]; t
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Apr 03 2013
STATUS
approved