A357364 Primes p such that either p^(q-1) == 1 (mod q^2) or q^(p-1) == 1 (mod p^2), where q = A151800(A151800(A151800(p))).

11, 23, 41, 107, 389, 1987673, 35603983

Offset: 1

Author

Felix Fröhlich, Sep 25 2022

Keywords

• nonn
• hard
• more

Crossrefs

Cf. A151800, A357362, A357363, A357365.

Programs

PARI

is(n) = my(b=precprime(precprime(precprime(n-1)-1)-1)); Mod(b, n^2)^(n-1)==1 || Mod(n, b^2)^(b-1)==1
forprime(p=7, , if(is(p), print1(p, ", ")))