OFFSET
1,2
COMMENTS
No other terms below 10^15. - Max Alekseyev, Nov 24 2017
Some larger terms: 1440926367749746685, 76025040962646716305439353859479569558065. - Max Alekseyev, Jun 29 2011
MATHEMATICA
Join[{1, 2}, Select[Range[1, 5000], Mod[13^#, #] == 3 &]] (* G. C. Greubel, Nov 19 2017 *)
Join[{1, 2}, Select[Range[10000000], PowerMod[13, #, #] == 3 &]] (* Robert Price, Apr 10 2020 *)
PROG
(PARI) isok(n) = Mod(13, n)^n == 3; \\ Michel Marcus, Nov 19 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Feb 19 2006
EXTENSIONS
Two more terms from Ryan Propper, Jan 09 2008
Terms 1,2 are prepended and a(9)-a(14) are added by Max Alekseyev, Jun 29 2011; Nov 24, 2017
STATUS
approved