OFFSET
1,1
COMMENTS
2nd prime p such that p^2 divides prime(n)^(p-1) - 1.
2nd prime p such that p divides the Fermat quotient q_p(p_n) = ((p_n)^(p-1) - 1)/p, where p_n = prime(n).
a(5) is unknown: 71 is the only known prime p that divides q_p(11).
If a(5) is found, the sequence continues a(6) = 863, a(7) = 3, a(8) = 7, a(9) = 2481757.
LINKS
EXAMPLE
PROG
(PARI) {default(primelimit, 10^7); for(n=1, 9, a=prime(n); c=0; forprime(p=2, 10^7, if(Mod(a, p^2)^(p-1)==1, c++; if(c==2, print1(p, ", "); next(2)))); print1(">10^7, "))} \\ Jens Kruse Andersen, Jun 18 2014
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Jonathan Sondow, Jun 20 2010, Jun 24 2010
STATUS
approved