OFFSET
1,1
LINKS
MAPLE
q:=n->(irem(2^(2^padic:-ordp(ithprime(n)-1, 2))-1, ithprime(n)) = 0):
select(q, [$1..10^5])[]; # Lorenzo Sauras Altuzarra, Feb 20 2023
PROG
(PARI) is_a023394(p)=p>2 && Mod(2, p)^lift(Mod(2, znorder(Mod(2, p)))^p)==1 && isprime(p) \\ after Charles R Greathouse IV in A023394
my(i=1); forprime(p=1, , if(is_a023394(p), print1(i, ", ")); i++) \\ Felix Fröhlich, Apr 30 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Lorenzo Sauras Altuzarra, Apr 28 2021
EXTENSIONS
More terms from Michel Marcus, Apr 29 2021
More terms from Amiram Eldar, Apr 29 2021
STATUS
approved