OFFSET
1,1
COMMENTS
Odd primes that divide 7^p-1 for some prime p [after Robert Israel].
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= n -> isprime(n) and isprime(numtheory:-order(7, n)):
select(filter, [2, 3, 5, seq(i, i=11..10000, 2)]); # Robert Israel, Jan 03 2025
MATHEMATICA
Select[Prime@Range@4000, PrimeQ@MultiplicativeOrder[7, #]&]
PROG
(Magma) [p: p in PrimesInInterval(2, 4000) | IsPrime(Modorder(7, p))];
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Dec 17 2024
STATUS
approved