A259309 Numbers k such that 1 + sigma(k)^4 is prime.
1, 3, 5, 12, 14, 15, 19, 23, 28, 33, 34, 35, 39, 40, 47, 53, 57, 58, 73, 76, 79, 86, 88, 89, 104, 112, 116, 118, 126, 131, 133, 134, 138, 139, 145, 147, 148, 154, 163, 165, 173, 175, 179, 183, 185, 189, 193, 197, 201, 204, 206, 207, 213, 216, 219, 220, 224
Offset: 1
Links
- Robert Price, Table of n, a(n) for n = 1..1459
- OEIS Wiki, Cyclotomic Polynomials at x=n, n! and sigma(n)
Programs
-
Magma
[n: n in [1..250] | IsPrime(1+SumOfDivisors(n)^4)]; // Vincenzo Librandi, Jun 24 2015
-
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, #]]] &]
-
PARI
is(n)=my(s=sigma(n)); isprime(s^4+1) \\ Charles R Greathouse IV, May 22 2017
Comments