OFFSET
1,1
COMMENTS
Numbers k such that gcd(k, 6) = 1 and if m is the least positive integer such that k divides 3^m - 2^m, then (3^m - 2^m)/k is a prime number.
The corresponding primes are 5, 71, 19, 2002867877, 29927, 29, 7, 5, ...
LINKS
FORMULA
13 is a term since ord(3/2, 13) = 4 and (3^4 - 2^4)/13 = 5 is a prime number.
MATHEMATICA
ord[n_] := Module[{k = 1}, While[! Divisible[PowerMod[3, k, n] - PowerMod[2, k, n], n], k++]; k]; f[k_] := 3^k - 2^k; Select[Range[1000], CoprimeQ[6, #] && PrimeQ[f[ord[#]]/#] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 24 2021
STATUS
approved