OFFSET
1,2
COMMENTS
This sequence includes all powers of 3. a(2) = 2, a(3) = 3, a(6) = 11 and a(45) = 1006003 are the only known primes in this sequence.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..400 (terms 1..100 from Harvey P. Dale)
MATHEMATICA
Do[f=PowerMod[3, n, n]; g=PowerMod[3, n, n^2]; If[f==g, Print[n]], {n, 1, 1100000}]
Select[Range[600000], PowerMod[3, #, #]==PowerMod[3, #, #^2]&] (* Harvey P. Dale, Feb 21 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Dec 07 2006
STATUS
approved