OFFSET
1,2
LINKS
Robert Price, Table of n, a(n) for n = 1..1786
MAPLE
with(numtheory): A258976:=n->`if`(isprime(1+sigma(n)^2), n, NULL): seq(A258976(n), n=1..500); # Wesley Ivan Hurt, Jul 09 2015
MATHEMATICA
Select[ Range[10000], PrimeQ[ 1 + DivisorSigma[1, #]^2] & ]
Select[ Range[10000], PrimeQ[ Cyclotomic[4, DivisorSigma[1, #]]] &]
PROG
(Magma) [n: n in [1..250] | IsPrime(1 + DivisorSigma(1, n)^2)]; // Vincenzo Librandi, Jun 16 2015
(PARI) lista(nn) = for (n=1, nn, if (isprime(1+sigma(n)^2), print1(n, ", "))); \\ Michel Marcus, Jun 17 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Price, Jun 15 2015
STATUS
approved