A283515 Numbers k such that sigma(k^(k-1)) is a prime.
2, 3, 4, 16, 19, 31, 7547
Offset: 1
Examples
sigma(4^3) = sigma(64) = 127 (prime).
Programs
-
Magma
[n: n in [1..500] | IsPrime(SumOfDivisors(n^(n-1)))];
-
Mathematica
fQ[n_] := PrimeQ[DivisorSigma[1, n^(n - 1)]]; Select[Range@1000, fQ] (* Robert G. Wilson v, Mar 10 2017 *)
-
PARI
isok(n) = isprime(sigma(n^(n-1))); \\ Michel Marcus, Mar 10 2017
Extensions
a(7) from Giovanni Resta, Mar 10 2017
Comments