A258974 a(n) = 1 + sigma(n)^2.
2, 10, 17, 50, 37, 145, 65, 226, 170, 325, 145, 785, 197, 577, 577, 962, 325, 1522, 401, 1765, 1025, 1297, 577, 3601, 962, 1765, 1601, 3137, 901, 5185, 1025, 3970, 2305, 2917, 2305, 8282, 1445, 3601, 3137, 8101, 1765, 9217, 1937, 7057, 6085, 5185, 2305
Offset: 1
Links
- Robert Price, Table of n, a(n) for n = 1..10000
- OEIS Wiki, Cyclotomic Polynomials at x=n, n! and sigma(n)
Crossrefs
Programs
-
Magma
[(1 + DivisorSigma(1, n)^2): n in [1..50]]; // Vincenzo Librandi, Jun 16 2015
-
Maple
with(numtheory): A258974:=n->1+sigma(n)^2: seq(A258974(n), n=1..100); # Wesley Ivan Hurt, Jul 09 2015
-
Mathematica
Table[1 + DivisorSigma[1, n]^2, {n, 10000}] Table[Cyclotomic[4, DivisorSigma[1, n]], {n, 10000}]
-
PARI
a(n)=sigma(n)^2+1 \\ Charles R Greathouse IV, Jun 18 2015
Formula
a(n) = 1 + A000203(n)^2.
a(n) = 1 + A072861(n). - Omar E. Pol, Jun 19 2015