cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A356533 a(n) = sigma_2(n)^2.

Original entry on oeis.org

1, 25, 100, 441, 676, 2500, 2500, 7225, 8281, 16900, 14884, 44100, 28900, 62500, 67600, 116281, 84100, 207025, 131044, 298116, 250000, 372100, 280900, 722500, 423801, 722500, 672400, 1102500, 708964, 1690000, 925444, 1863225, 1488400, 2102500, 1690000, 3651921
Offset: 1

Views

Author

Vaclav Kotesovec, Aug 11 2022

Keywords

Crossrefs

Cf. A001157, A127473, A035116, A072861, A356535 (partial sums).

Programs

  • Mathematica
    Table[DivisorSigma[2, n]^2, {n, 1, 40}]
  • PARI
    a(n) = sigma(n, 2)^2; \\ Michel Marcus, Aug 11 2022

Formula

Dirichlet g.f.: zeta(s) * zeta(s-2)^2 * zeta(s-4) / zeta(2*s-4).
Multiplicative with a(p^e) = ((p^(2*e+2)-1)/(p^2-1))^2. - Amiram Eldar, Aug 11 2022
a(n) = A001157(n)^2. - R. J. Mathar, Aug 18 2022

A356536 a(n) = Sum_{k=1..n} sigma_3(k)^2.

Original entry on oeis.org

1, 82, 866, 6195, 22071, 85575, 203911, 546136, 1119185, 2405141, 4179365, 8357301, 13188505, 22773721, 35220505, 57132266, 81279662, 127696631, 174756231, 259359435, 352134859, 495847003, 643907227, 912211627, 1160305628, 1551633152, 1969426752, 2600039296
Offset: 1

Views

Author

Vaclav Kotesovec, Aug 11 2022

Keywords

Comments

Partial sums of A356534.
In general, for m>0, Sum_{k=1..n} sigma_m(k)^2 ~ zeta(2*m+1) * zeta(m+1)^2 * n^(2*m+1) / ((2*m+1) * zeta(2*m+2)).

Crossrefs

Programs

  • Mathematica
    Table[Sum[DivisorSigma[3, k]^2, {k, 1, n}], {n, 1, 40}]
    Accumulate[DivisorSigma[3,Range[40]]^2] (* This program is much more efficient than the first program above. *) (* Harvey P. Dale, Feb 27 2023 *)
  • PARI
    a(n) = sum(k=1, n, sigma(k, 3)^2); \\ Michel Marcus, Aug 11 2022

Formula

a(n) ~ zeta(7) * n^7 / 6.
Showing 1-2 of 2 results.