OFFSET
1,2
COMMENTS
No other terms below 10^16. - Max Alekseyev, Nov 02 2018
MATHEMATICA
Select[Range[1, 500000], Mod[13^#, #] == 2 &] (* G. C. Greubel, Nov 19 2017 *)
Join[{1}, Select[Range[5000000], PowerMod[13, #, #] == 2 &]] (* Robert Price, Apr 10 2020 *)
PROG
(PARI) isok(n) = Mod(13, n)^n == 2; \\ Michel Marcus, Nov 19 2017
CROSSREFS
Cf. A116609.
Solutions to b^n == 2 (mod n): A015919 (b=2), A276671 (b=3), A130421 (b=4), A124246 (b=5), A277401 (b=7), this sequence (b=13), A333269 (b=17).
KEYWORD
more,nonn
AUTHOR
Zak Seidov, Feb 19 2006
EXTENSIONS
One more term from Ryan Propper, Jun 11 2006
Term a(1)=1 is prepended and a(7)-a(8) are added by Max Alekseyev, Jun 29 2011
STATUS
approved