A318755 -id:A318755 - cp's OEIS frontend

A061502 a(n) = Sum_{k<=n} tau(k)^2, where tau = number of divisors function A000005.

1, 5, 9, 18, 22, 38, 42, 58, 67, 83, 87, 123, 127, 143, 159, 184, 188, 224, 228, 264, 280, 296, 300, 364, 373, 389, 405, 441, 445, 509, 513, 549, 565, 581, 597, 678, 682, 698, 714, 778, 782, 846, 850, 886, 922, 938, 942, 1042, 1051, 1087

Offset: 1

N. J. A. Sloane, Jun 14 2001

nonn

R. Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; Chapter II, Problem 56.

N. J. A. Sloane, Table of n, a(n) for n = 1..1024
Adrian Dudek, On the Success of Mishandling Euclid's Lemma, arXiv:1602.03555 [math.HO], 2016. See B(n) p. 2.
Chaohua Jia and Ayyadurai Sankaranarayanan, The mean square of the divisor function, Acta Arithmetica 164 (2014), 181-208.
Michaela Cully-Hugill and Timothy Trudgian, Two explicit divisor sums, arXiv:1911.07369 [math.NT], Nov 19 2019
Vaclav Kotesovec, Graph - The asymptotic ratio (1000000 terms)
Florian Luca and László Tóth, The r-th moment of the divisor function: an elementary approach, Journal of Integer Sequences 20 (2017), Article 17.7.4, 8 pp.
Adolf Piltz, Über das Gesetz, nach welchem die mittlere Darstellbarkeit der natürlichen Zahlen als Produkte einer gegebenen Anzahl Faktoren mit der Grösse der Zahlen wächst, 1881.
Ramanujan's Papers, Some formulas in the analytic theory of numbers, Messenger of Mathematics, XLV, 1916, 81-84, Formula (3).
D. Suryanarayana and R. Rama Chandra Rao, On an Asymptotic Formula of Ramanujan, Mathematica Scandinavica, 32, 258-264, 1973.
B. M. Wilson, Proofs of some formulas enunciated by Ramanujan, Proc. London Math. Soc. (2) 21 (1922) 235-255.

Cf. A000005, A035116, A061503, A318755.
Cf. A092742 (A), A245074 (B), A319090 (C), A319091 (D).
Cf. A057434, A072379, A074789.

[&+[NumberOfDivisors(k^2)*Floor(n/k): k in [1..n]]: n in [1..60]]; // Vincenzo Librandi, Sep 10 2016

Table[Sum[DivisorSigma[0, k^2]*Floor[n/k], {k, 1, n}], {n, 1, 50}] (* Vaclav Kotesovec, Aug 30 2018 *)
Accumulate[Table[DivisorSigma[0, n]^2, {n, 1, 50}]] (* Vaclav Kotesovec, Sep 10 2018 *)

for (n=1, 1024, write("b061502.txt", n, " ", sum(k=1, n, numdiv(k)^2)) ) \\ Harry J. Smith, Jul 23 2009

vector(60, n, sum(k=1, n, numdiv(k)^2)) \\ Michel Marcus, Jul 23 2015

first(n)=my(v=vector(n),s); forfactored(k=1,n, v[k[1]] = s += numdiv(k)^2); v; \\ Charles R Greathouse IV, Nov 28 2018

a(n) = Sum_{k=1..n} tau(k^2)*floor(n/k).
Asymptotic to A*n*log(n)^3 + B*n*log(n)^2 + C*n*log(n) + D*n + O(n^(1/2+eps)) where A = 1/Pi^2 and B = (12*gamma-3)/Pi^2 - 36*zeta'(2)/Pi^4. [corrected by Vaclav Kotesovec, Aug 30 2018]
C = 36*gamma^2/Pi^2 - (288*z1/Pi^4 + 24/Pi^2)*gamma + (864*z1^2/Pi^6 + 72*z1/Pi^4 - 72/Pi^4*z2 + 6/Pi^2) - 24*g1/Pi^2 and D = 24*gamma^3/Pi^2 - (432*z1 /Pi^4+ 36/Pi^2)*gamma^2 + (3456*z1^2/Pi^6 + 288*(z1-z2)/Pi^4 + 24/Pi^2 - 72*g1/Pi^2)*gamma + g1*(288*z1/Pi^4 + 24/Pi^2)-10368*z1^3/Pi^8 - 864*z1^2/Pi^6 + 1728*z2*z1/Pi^6 + 72*(z2-z1)/Pi^4- 48*z3/Pi^4 + (12*g2-6)/Pi^2, where gamma is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994, z3 = Zeta'''(2) = A201995 and g1, g2 are the Stieltjes constants, see A082633 and A086279. - Vaclav Kotesovec, Sep 10 2018
See Cully-Hugill & Trudgian, Theorem 2, for an explicit version of the asymptotic given above. - Charles R Greathouse IV, Nov 19 2019

Definition corrected by N. J. A. Sloane, May 25 2008

A320896 a(n) = Sum_{k=1..n} k * tau(k)^2, where tau is A000005.

Original entry on oeis.org

1, 9, 21, 57, 77, 173, 201, 329, 410, 570, 614, 1046, 1098, 1322, 1562, 1962, 2030, 2678, 2754, 3474, 3810, 4162, 4254, 5790, 6015, 6431, 6863, 7871, 7987, 9907, 10031, 11183, 11711, 12255, 12815, 15731, 15879, 16487, 17111, 19671, 19835, 22523, 22695, 24279

Offset: 1

Vaclav Kotesovec, Oct 23 2018

nonn

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ramanujan's Papers, Some formulas in the analytic theory of numbers, Messenger of Mathematics, XLV, 1916, 81-84.

Cf. A061502, A318755, A320897.

Cf. A000005, A006218, A143127, A319085, A320895.

Accumulate[Table[k*DivisorSigma[0, k]^2, {k, 1, 50}]]

a(n) = sum(k=1, n, k*numdiv(k)^2); \\ Michel Marcus, Oct 23 2018

a(n) ~ n^2 * (3*(Pi^6*(-1 - 24*g^2 + 32*g^3 + g*(8 - 96*s1) + 16*s1 + 16*s2) - 13824*z1^3 + 576*Pi^2*z1*((-1 + 8*g)*z1 + 4*z2) - 8*Pi^4*(3*(1 - 8*g + 24*g^2 - 16*s1)*z1 - 6*z2 + 48*g*z2 + 8*z3)) + 6*(Pi^6*(1 - 8*g + 24*g^2 - 16*s1) + 576*Pi^2*z1^2 - 24*Pi^4*(-z1 + 8*g*z1 + 2*z2))*log(n) + 6*((-1 + 8*g)*Pi^6 - 24*Pi^4*z1)*log(n)^2 + 4*Pi^6*log(n)^3) / (8*Pi^8), where g is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994, z3 = Zeta'''(2) = A201995 and s1, s2 are the Stieltjes constants, see A082633 and A086279.

A320897 a(n) = Sum_{k=1..n} k^2 * tau(k)^2, where tau is A000005.

Original entry on oeis.org

1, 17, 53, 197, 297, 873, 1069, 2093, 2822, 4422, 4906, 10090, 10766, 13902, 17502, 23902, 25058, 36722, 38166, 52566, 59622, 67366, 69482, 106346, 111971, 122787, 134451, 162675, 166039, 223639, 227483, 264347, 281771, 300267, 319867, 424843, 430319, 453423

Offset: 1

Vaclav Kotesovec, Oct 23 2018

nonn

In general, for m>=0, Sum_{k=1..n} k^m * tau(k)^2 ~ n^(m+1) * (log(n))^3 / ((m+1) * Pi^2).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ramanujan's Papers, Some formulas in the analytic theory of numbers, Messenger of Mathematics, XLV, 1916, 81-84.

Cf. A061502, A318755, A320896.

Cf. A000005, A006218, A143127, A319085, A320895.

Accumulate[Table[k^2*DivisorSigma[0, k]^2, {k, 1, 50}]]

a(n) = sum(k=1, n, k^2*numdiv(k)^2); \\ Michel Marcus, Oct 23 2018

a(n) ~ n^3 * (2*Pi^6*(-1 + 12*g - 54*g^2 + 108*g^3 + 36*s1 - 324*g*s1 + 54*s2) - 93312*z1^3 + 2592*Pi^2*z1*(-z1 + 12*g*z1 + 6*z2) - 72*Pi^4*(z1 - 12*g*z1 + 54*g^2*z1 - 36*s1*z1 - 3*z2 + 36*g*z2 + 6*z3) + 6*(Pi^6*(1 - 12*g + 54*g^2 - 36*s1) + 1296*Pi^2*z1^2 - 36*Pi^4*(-z1 + 12*g*z1 + 3*z2))*log(n) + 9*((-1 + 12*g)*Pi^6 - 36*Pi^4*z1)*log(n)^2 + 9*Pi^6*log(n)^3) / (27*Pi^8), where g is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994, z3 = Zeta'''(2) = A201995 and s1, s2 are the Stieltjes constants, see A082633 and A086279.

A319089 a(n) = tau(n)^3, where tau is A000005.

Original entry on oeis.org

1, 8, 8, 27, 8, 64, 8, 64, 27, 64, 8, 216, 8, 64, 64, 125, 8, 216, 8, 216, 64, 64, 8, 512, 27, 64, 64, 216, 8, 512, 8, 216, 64, 64, 64, 729, 8, 64, 64, 512, 8, 512, 8, 216, 216, 64, 8, 1000, 27, 216, 64, 216, 8, 512, 64, 512, 64, 64, 8, 1728, 8, 64, 216, 343

Offset: 1

Vaclav Kotesovec, Sep 10 2018

Eric Weisstein's World of Mathematics, Divisor Function.

Cf. A000005, A006218, A035116, A061502, A318755 (partial sums).

with(numtheory): seq(tau(n)^3, n=1..100); # Ridouane Oudra, Mar 07 2023

Mathematica
```
DivisorSigma[0, Range[100]]^3
```

a(n) = numdiv(n)^3; \\ Altug Alkan, Sep 10 2018

for(n=1, 100, print1(direuler(p=2, n, (1 + 4*X + X^2)/(1 - X)^4)[n], ", ")) \\ Vaclav Kotesovec, Mar 09 2023

Multiplicative with a(p^e) = (e+1)^3. - Amiram Eldar, Dec 31 2022
a(n) = Sum_{d1|n} Sum_{d2|n} tau(d1*d2). - Ridouane Oudra, Mar 07 2023
From Vaclav Kotesovec, Mar 09 2023: (Start)
Dirichlet g.f.: Product_{p prime} p^(2*s) * (1 + 4*p^s + p^(2*s)) / (p^s - 1)^4.
Dirichlet g.f.: zeta(s)^8 * Product_{p prime} (1 - 9/p^(2*s) + 16/p^(3*s) - 9/p^(4*s) + 1/p^(6*s)), (with a product that converges for s=1). (End)

cp's OEIS Frontend

A061502 a(n) = Sum_{k<=n} tau(k)^2, where tau = number of divisors function A000005.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

PARI

Formula

Extensions

A320896 a(n) = Sum_{k=1..n} k * tau(k)^2, where tau is A000005.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A320897 a(n) = Sum_{k=1..n} k^2 * tau(k)^2, where tau is A000005.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A319089 a(n) = tau(n)^3, where tau is A000005.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula