OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from G. C. Greubel)
MATHEMATICA
Join[{1}, Select[Range[1, 750], Mod[13^#, #] == 1 &]] (* G. C. Greubel, Nov 19 2017 *)
Join[{1}, Select[Range[320], PowerMod[13, #, #] == 1 &]] (* Robert Price, Apr 10 2020 *)
PROG
(PARI) isok(n) = Mod(13, n)^n == 1; \\ Michel Marcus, Nov 19 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Feb 19 2006
EXTENSIONS
1 prepended by Max Alekseyev, Jun 28 2011
STATUS
approved