OFFSET
1,1
COMMENTS
These primes are neither sorted nor uniqued. They are listed in the order found in A258978.
LINKS
Robert Price, Table of n, a(n) for n = 1..932
MAPLE
with(numtheory): A258980:=n->`if`(isprime(1 + sigma(n) + sigma(n)^2 + sigma(n)^3 + sigma(n)^4), 1 + sigma(n) + sigma(n)^2 + sigma(n)^3 + sigma(n)^4, NULL): seq(A258980(n), n=1..300); # Wesley Ivan Hurt, Jul 09 2015
MATHEMATICA
Select[Table[1 + DivisorSigma[1, n] + DivisorSigma[1, n]^2 + DivisorSigma[1, n]^3 + DivisorSigma[1, n]^4, {n, 10000}], PrimeQ]
Select[Table[Cyclotomic[5, DivisorSigma[1, n]], {n, 10000}], PrimeQ]
PROG
(Magma) [m: n in [1..200] | IsPrime(m) where m is (1 + DivisorSigma(1, n) + DivisorSigma(1, n)^2 + DivisorSigma(1, n)^3 + DivisorSigma(1, n)^4)]; // Vincenzo Librandi, Jun 16 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Price, Jun 15 2015
STATUS
approved