OFFSET
1,1
LINKS
W. Keller and J. Richstein Fermat quotients that are divisible by p.
MATHEMATICA
Do[p = Prime[n]; q = 2; While[PowerMod[q, p-1, p^8] != 1, q = NextPrime[q]]; Print[q], {n, 100}] (* Ryan Propper, Apr 01 2007 *)
PROG
(PARI) { a(n) = local(p, x, y); if(n==1, return(257)); p=prime(n); x=znprimroot(p^8)^(p^7); vecsort( vector(p-1, i, y=lift(x^i); while(!isprime(y), y+=p^8); y ) )[1] } \\ Max Alekseyev, May 30 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Nov 29 2006
EXTENSIONS
More terms from Ryan Propper, Apr 01 2007
More terms from Max Alekseyev, May 30 2007
STATUS
approved