OFFSET
1,2
LINKS
Robert Price, Table of n, a(n) for n = 1..249
MATHEMATICA
k=1; Do[ p=Prime[k]; If[ IntegerQ[ (PowerMod[ p+1, n^2, n^3 ] - 1 )/n^3 ], Print[ {k, p, n} ]], {n, 1, 100000} ]
Join[{1}, Select[Range[3000000], PowerMod[3, #^2, #^3] == 1 &]] (* Robert Price, Mar 31 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Apr 03 2007
STATUS
approved