OFFSET
1,1
COMMENTS
sigma(k) is the sum of the divisors of k (A000203).
a(8) > 10^4.
Corresponding values of k^(k-1): 2, 9, 64, 1152921504606846976, ...
Corresponding values of sigma(k^(k-1)): 3, 13, 127, 2305843009213693951, ...
Subsequence of A280257 (numbers k such that tau(k^(k-1)) is prime).
Prime terms are in A088790.
For k < 1000, sigma(k^(k+1)) is prime only for k = 5: sigma(5^6) = sigma(15625) = 19531 (prime).
EXAMPLE
sigma(4^3) = sigma(64) = 127 (prime).
MATHEMATICA
fQ[n_] := PrimeQ[DivisorSigma[1, n^(n - 1)]]; Select[Range@1000, fQ] (* Robert G. Wilson v, Mar 10 2017 *)
PROG
(Magma) [n: n in [1..500] | IsPrime(SumOfDivisors(n^(n-1)))]
(PARI) isok(n) = isprime(sigma(n^(n-1))); \\ Michel Marcus, Mar 10 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jaroslav Krizek, Mar 10 2017
EXTENSIONS
a(7) from Giovanni Resta, Mar 10 2017
STATUS
approved