OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
m=131; 131=2*65+1; 3^65 mod 131 = 1 and 65^65 mod 131 = 1.
MAPLE
filter:= proc(n) local k;
k:= (n-1)/2;
3 &^ k mod n = 1 and k &^ k mod n = 1
end proc:
select(filter, [seq(i, i=3..3000, 4)]); # Robert Israel, Nov 28 2019
MATHEMATICA
2#+1&/@Select[Range[1, 1200, 2], PowerMod[3, #, 2#+1]==PowerMod[ #, #, 2#+1] == 1&] (* Harvey P. Dale, May 05 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alzhekeyev Ascar M, May 26 2016
STATUS
approved