OFFSET
1,2
COMMENTS
From Alexander Adamchuk, Jan 11 2007: (Start)
2 divides a(n) for n>1. 2^2 divides a(n) for n>2. 5 divides a(n) for n>3.
11 divides a(n) for n = {5,7,8,9,10,12,13,14,15,16,17,18,19,20,22,23,24,26,27, 28,29,30,31,31,33,34,35,...}.
11^2 divides a(n) for n = {7,12,14,15,26,27,29,30,31,33,34,...}.
Prime factors of a(n) in order of their appearance in a(n) are {2,5,11,61,23,67,1181,661,47,1321,367,3851,5501,727,461,269,...}. (End)
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..300
MATHEMATICA
Select[Range[30000], IntegerQ[(PowerMod[3, #, #^2 ]-1)/#^2 ]&]
Join[{1}, Select[Range[335*10^4], PowerMod[3, #, #^2]==1&]] (* Harvey P. Dale, Oct 02 2019 *)
PROG
(PARI) is(k) = Mod(3, k^2)^k == 1; \\ Amiram Eldar, May 21 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Jan 05 2007
EXTENSIONS
More terms from Ryan Propper and Alexander Adamchuk, Jan 05 2007
STATUS
approved