|
%I #23 Apr 10 2020 21:43:01
%S 1,2,3,4,6,8,9,12,14,16,18,20,24,27,28,32,36,40,42,48,54,56,60,64,68,
%T 72,80,81,84,96,98,100,108,112,120,126,128,136,140,144,160,162,168,
%U 180,183,192,196,200,204,216,220,224,240,243,252,256,272,280,288,294,300,320
%N Positive integers n such that 13^n == 1 (mod n).
%H Seiichi Manyama, <a href="/A116621/b116621.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from G. C. Greubel)
%t Join[{1}, Select[Range[1, 750], Mod[13^#, #] == 1 &]] (* _G. C. Greubel_, Nov 19 2017 *)
%t Join[{1}, Select[Range[320], PowerMod[13, #, #] == 1 &]] (* _Robert Price_, Apr 10 2020 *)
%o (PARI) isok(n) = Mod(13, n)^n == 1; \\ _Michel Marcus_, Nov 19 2017
%Y Column k=13 of A333432.
%Y Cf. A116609, A116611, A116620.
%K nonn
%O 1,2
%A _Zak Seidov_, Feb 19 2006
%E 1 prepended by _Max Alekseyev_, Jun 28 2011
|
