

A260377


Primes p such that A001221(p1)^(p1) == 1 (mod p^2).


3




OFFSET

1,1


COMMENTS

No further terms up to 10^9.
All currently (as of 20150723) known terms of A001220 and A019434 are in the sequence, as is A014127(2), but A014127(1) is not.


LINKS

Table of n, a(n) for n=1..8.


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(p1), p^2)^(p1)==1, print1(p, ", ")))


CROSSREFS

Cf. A001220, A014127, A019434, A244546.
Sequence in context: A023394 A176689 A256510 * A056130 A273871 A078726
Adjacent sequences: A260374 A260375 A260376 * A260378 A260379 A260380


KEYWORD

nonn,hard,more


AUTHOR

Felix FrÃ¶hlich, Jul 23 2015


STATUS

approved



