A061201 Partial sums of A007425: (tau<=)_3(n).
1, 4, 7, 13, 16, 25, 28, 38, 44, 53, 56, 74, 77, 86, 95, 110, 113, 131, 134, 152, 161, 170, 173, 203, 209, 218, 228, 246, 249, 276, 279, 300, 309, 318, 327, 363, 366, 375, 384, 414, 417, 444, 447, 465, 483, 492, 495, 540, 546, 564, 573, 591, 594, 624, 633, 663
Offset: 1
References
- M. N. Huxley, Area, Lattice Points and Exponential Sums, Oxford, 1996; p. 239.
Links
- Vaclav Kotesovec, Table of n, a(n) for n = 1..100000 (terms 1..1000 from Harry J. Smith)
- Benoit Cloitre, On the divisor and circle problems
- Vaclav Kotesovec, Graph - The asymptotic ratio (1000000 terms)
Crossrefs
Programs
-
Magma
[&+[NumberOfDivisors(k)*Floor(n/k): k in [1..n]]: n in [1..56]]; // Bruno Berselli, Apr 13 2011
-
Maple
b:= proc(k, n) option remember; uses numtheory; `if`(k=1, 1, add(b(k-1, d), d=divisors(n))) end: a:= proc(n) option remember; `if`(n=0, 0, b(3, n)+a(n-1)) end: seq(a(n), n=1..76); # Alois P. Heinz, Oct 23 2023 -
Mathematica
a[n_] := Sum[ DivisorSigma[0, k]*Floor[n/k], {k, 1, n}]; Table[a[n], {n, 1, 56}] (* Jean-François Alcover, Sep 20 2011, after Benoit Cloitre *) (* Asymptotics: *) n*(Log[n]^2/2 + (3*EulerGamma - 1)*Log[n] + 3*EulerGamma^2 - 3*EulerGamma - 3*StieltjesGamma[1] + 1) (* Vaclav Kotesovec, Sep 09 2018 *) Accumulate[a[n_]:=DivisorSum[n, DivisorSigma[0, #]&]; Array[a, 60]] (* Vincenzo Librandi, Jan 12 2020 *) -
PARI
a(n)=sum(k=1,n,numdiv(k)*floor(n/k)) \\ Benoit Cloitre, Apr 19 2007
-
PARI
{ for (n=1, 1000, write("b061201.txt", n, " ", sum(k=1, n, numdiv(k)*(n\k))) ) } \\ Harry J. Smith, Jul 18 2009 -
PARI
my(N=60, x='x+O('x^N)); Vec(sum(k=1, N, numdiv(k)*x^k/(1-x^k))/(1-x)) \\ Seiichi Manyama, Jul 24 2022 -
Python
from math import isqrt from sympy import integer_nthroot def A061201(n): return (m:=integer_nthroot(n,3)[0])**3+3*sum(-(s:=isqrt(r:=n//i))**2+(sum(r//k for k in range(1,s+1))<<1)-sum(n//(i*j) for j in range(1,m+1)) for i in range(1,m+1)) # Chai Wah Wu, Oct 23 2023
Formula
(tau<=)k(n) = Sum{i=1..n} tau_k(i).
a(n) = n * ( log(n)^2/2 + (3*g-1)*log(n) + 3*g^2-3*g-3*g1+1 ) + O(sqrt(n)), where g is the Euler-Mascheroni number ~ 0.57721... (see A001620), and g1 is the first Stieltjes constant ~ -0.072816 (see A082633). The determination of the precise size of the error term is an unsolved problem - see references. - Andrew Lelechenko, Apr 15 2011 [corrected by Vaclav Kotesovec, Sep 09 2018]
a(n) = Sum_{k=1..n} A000005(k)*floor(n/k). - Benoit Cloitre, Apr 19 2007
To compute a(n) for huge n (see A180365) in sublinear use a(n) = 3*Sum_{i=1..n3} A006218(n/i) - Sum_{j=1..n3} floor(n/(i*j)) + n3^3, where n3 = floor(n^(1/3)). - Andrew Lelechenko, Apr 15 2011
a(n) = Sum_{k=1..n} Sum_{i=1..n} floor(n/(i*k)). - Wesley Ivan Hurt, Sep 14 2017
G.f.: (1/(1-x)) * Sum_{k>=1} A000005(k) * x^k/(1 - x^k). - Seiichi Manyama, Jul 24 2022
Comments