OFFSET
1,1
COMMENTS
Subset of A187023.
If p is a prime, then p^p-1 has at least a prime factor that is congruent to 1 modulo p.
Also smallest prime factor of (p^p - 1)/(p - 1). - Jianing Song, Nov 03 2019
EXAMPLE
a(4) = 29 because prime(4) = 7 and 7^7 -1 = 823542 = 2 * 3 * 29 * 4733 => 29 == 1 (mod 7).
MAPLE
with(numtheory): for n from 1 to 34 do:i:=0:p:=ithprime(n):x:=p^p -1:y:=factorset(x):n1:=nops(y):for k from 1 to n1 while(i=0) do:z:=y[k]:if irem(z, p)=1 then i:=1: printf ( "%d %d \n", n, z):else fi:od:od:
MATHEMATICA
Table[p=First/@FactorInteger[Prime[n]^Prime[n]-1]; Select[p, Mod[#1, Prime[n]] == 1 &, 1][[1]], {n, 1, 10}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, May 20 2012
STATUS
approved