|
|
A260377
|
|
Primes p such that A001221(p-1)^(p-1) == 1 (mod p^2).
|
|
3
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
No further terms up to 10^9.
|
|
LINKS
|
|
|
EXAMPLE
|
A001221(1092) == 4 and 4^1092 == 1 (mod 1093^2), so 1093 is a term of the sequence.
|
|
PROG
|
(PARI) forprime(p=1, 1e9, if(Mod(omega(p-1), p^2)^(p-1)==1, print1(p, ", ")))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|