OFFSET

1,3

COMMENTS

MATHEMATICA

nn = 100; t = Table[0, {nn}]; Do[If[Mod[Prime[n], 4] == 3, t[[n]] = -1], {n, nn}]; n = 0; While[Times @@ t == 0, n++; s = FactorInteger[n^2 + 1][[-1, 1]]; p = PrimePi[s]; If[p <= nn && t[[p]] == 0, t[[p]] = n]]; Do[If[Mod[Prime[n], 4] == 3, t[[n]] = 0], {n, nn}]; t

CROSSREFS

KEYWORD

nonn

AUTHOR

T. D. Noe, Apr 03 2013

STATUS

approved