A258977 Primes in A258974.
2, 17, 37, 197, 577, 577, 401, 1297, 577, 1601, 3137, 2917, 3137, 8101, 7057, 15377, 2917, 14401, 14401, 8101, 7057, 15877, 5477, 15377, 15877, 7057, 50177, 14401, 32401, 8101, 14401, 24337, 44101, 78401, 12101, 57601, 44101, 32401, 50177, 24337, 30977
Offset: 1
Keywords
Links
- Robert Price, Table of n, a(n) for n = 1..1786
- OEIS Wiki, Cyclotomic Polynomials at x=n, n! and sigma(n)
Programs
-
Magma
[m: n in [1..200] | IsPrime(m) where m is 1 + DivisorSigma(1, n)^2]; // Vincenzo Librandi, Jun 16 2015
-
Maple
with(numtheory): A258977:=n->`if`(isprime(1+sigma(n)^2), 1+sigma(n)^2, NULL): seq(A258977(n), n=1..300); # Wesley Ivan Hurt, Jul 09 2015
-
Mathematica
Select[Table[1 + DivisorSigma[1, n]^2, {n, 10000}], PrimeQ] Select[Table[Cyclotomic[4, DivisorSigma[1, n]], {n, 10000}], PrimeQ]
-
PARI
for(n=1,100,p=1+sigma(n)^2;if(isprime(p),print1(p,", "))) \\ Derek Orr, Jun 18 2015
Comments