OFFSET
1,1
COMMENTS
Terms must be divisible by the square of a Mirimanoff prime p (or base-3 Wieferich prime, A014127) such that the multiplicative order of 3 modulo p is not divisible by 3. So far, the only known Mirimanoff primes are 11 and 1006003. The multiplicative order of 3 modulo 11 is 5, not a multiple of 3, while the multiplicative order of 3 modulo 1006003 is 1006002, which is a multiple of 3. As a result, all known terms are divisible by 3*11^2 = 363.
EXAMPLE
726 is a term because 726 divides 3^726 - 3 and 726 = 2 * 3 * 11^2.
PROG
(PARI) forstep(n=3, 10^9, 3, if(Mod(3, n)^n==3 && !issquarefree(n), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jianing Song, Feb 16 2019
EXTENSIONS
More terms from Jinyuan Wang, Feb 18 2019
STATUS
approved
