login
A357363
Primes p such that either p^(q-1) == 1 (mod q^2) or q^(p-1) == 1 (mod p^2), where q = A151800(A151800(p)).
3
5, 19, 263, 1667
OFFSET
1,1
PROG
(PARI) is(n) = my(b=precprime(precprime(n-1)-1)); Mod(b, n^2)^(n-1)==1 || Mod(n, b^2)^(b-1)==1
forprime(p=5, , if(is(p), print1(p, ", ")))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Felix Fröhlich, Sep 25 2022
STATUS
approved