OFFSET
1,2
LINKS
Robert Price, Table of n, a(n) for n = 1..1459
MAPLE
with(numtheory): A259309:=n->`if`(isprime(1 + sigma(n)^4), n, NULL): seq(A259309(n), n=1..500); # Wesley Ivan Hurt, Jul 09 2015
MATHEMATICA
Select[ Range[10000], PrimeQ[ 1 + DivisorSigma[1, #]^4] & ]
Select[ Range[10000], PrimeQ[ Cyclotomic[8, DivisorSigma[1, #]]] &]
PROG
(Magma) [n: n in [1..250] | IsPrime(1+SumOfDivisors(n)^4)]; // Vincenzo Librandi, Jun 24 2015
(PARI) is(n)=my(s=sigma(n)); isprime(s^4+1) \\ Charles R Greathouse IV, May 22 2017
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Robert Price, Jun 24 2015
STATUS
approved