OFFSET
1,3
COMMENTS
n^Phi(n)-1 is never prime. This sequence is an outgrowth of Euler's generalization to Fermat's little theorem.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
FORMULA
MATHEMATICA
a = {}; Do[ a = Append[a, FactorInteger[ n^EulerPhi[n] - 1, FactorComplete -> False][[1, 1]]], {n, 1, 100}]; a
PROG
(PARI)
A020639(n) = if(1==n, n, forprime(p=2, , if(!(n%p), return(p))));
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jan 15 2002
STATUS
approved