Divisor summatory function



In number theory, the divisor summatory function is a function that is a sum over the divisor function. It frequently occurs in the study of the asymptotic behaviour of the Riemann zeta function. The various studies of the behaviour of the divisor function are sometimes called divisor problems.

Definition
The divisor summatory function is defined as


 * $$D(x)=\sum_{n\le x} d(n) = \sum_{j,k \atop jk\le x} 1$$

where


 * $$d(n)=\sigma_0(n) = \sum_{j,k \atop jk=n} 1$$

is the divisor function. The divisor function counts the number of ways that the integer n can be written as a product of two integers. More generally, one defines


 * $$D_k(x)=\sum_{n\le x} d_k(n)= \sum_{m\le x}\sum_{mn\le x} d_{k-1}(n)$$

where dk(n) counts the number of ways that n can be written as a product of k numbers. This quantity can be visualized as the count of the number of lattice points fenced off by a hyperbolic surface in k dimensions. Thus, for k=2, D(x) = D2(x) counts the number of points on a square lattice bounded on the left by the vertical-axis, on the bottom by the horizontal-axis, and to the upper-right by the hyperbola jk = x. Roughly, this shape may be envisioned as a hyperbolic simplex. This allows us to provide an alternative expression for D(x), and a simple way to compute it in $$O(\sqrt{x})$$ time:


 * $$D(x)=\sum_{k=1}^x \left\lfloor\frac{x}{k}\right\rfloor = 2 \sum_{k=1}^u \left\lfloor\frac{x}{k}\right\rfloor - u^2$$, where $$u = \left\lfloor \sqrt{x}\right\rfloor$$

If the hyperbola in this context is replaced by a circle then determining the value of the resulting function is known as the Gauss circle problem.

Sequence of D(n): 0, 1, 3, 5, 8, 10, 14, 16, 20, 23, 27, 29, 35, 37, 41, 45, 50, 52, 58, 60, 66, 70, 74, 76, 84, 87, 91, 95, 101, 103, 111, ...

Dirichlet's divisor problem
Finding a closed form for this summed expression seems to be beyond the techniques available, but it is possible to give approximations. The leading behavior of the series is given by


 * $$D(x) = x\log x + x(2\gamma-1) + \Delta(x)\ $$

where $$\gamma$$ is the Euler–Mascheroni constant, and the error term is


 * $$\Delta(x) = O\left(\sqrt{x}\right).$$

Here, $$O$$ denotes Big-O notation. This estimate can be proven using the Dirichlet hyperbola method, and was first established by Dirichlet in 1849. The Dirichlet divisor problem, precisely stated, is to improve this error bound by finding the smallest value of $$\theta$$ for which


 * $$\Delta(x) = O\left(x^{\theta+\epsilon}\right)$$

holds true for all $$\epsilon > 0$$. As of today, this problem remains unsolved. Progress has been slow. Many of the same methods work for this problem and for Gauss's circle problem, another lattice-point counting problem. Section F1 of Unsolved Problems in Number Theory surveys what is known and not known about these problems.


 * In 1904, G. Voronoi proved that the error term can be improved to $$O(x^{1/3}\log x).$$
 * In 1916, G. H. Hardy showed that $$\inf \theta \ge 1/4$$. In particular, he demonstrated that for some constant $$K$$, there exist values of x for which $$\Delta(x) > Kx^{1/4}$$ and values of x for which $$\Delta(x) < -Kx^{1/4}$$.
 * In 1922, J. van der Corput improved Dirichlet's bound to $$\inf \theta \le 33/100 = 0.33$$.
 * In 1928, van der Corput proved that $$\inf \theta \le 27/82 = 0.3\overline{29268}$$.
 * In 1950, Chih Tsung-tao and independently in 1953 H. E. Richert proved that $$\inf \theta \le 15/46 = 0.32608695652...$$.
 * In 1969, Grigori Kolesnik demonstrated that $$\inf \theta \le 12/37 = 0.\overline{324}$$.
 * In 1973, Kolesnik demonstrated that $$\inf \theta \le 346/1067 = 0.32427366448...$$.
 * In 1982, Kolesnik demonstrated that $$\inf \theta \le 35/108 = 0.32\overline{407}$$.
 * In 1988, H. Iwaniec and C. J. Mozzochi proved that $$\inf \theta \leq 7/22 = 0.3\overline{18}$$.
 * In 2003, M.N. Huxley improved this to show that $$\inf \theta \leq 131/416 = 0.31490384615...$$.

So, $$\inf \theta$$ lies somewhere between 1/4 and 131/416 (approx. 0.3149); it is widely conjectured to be 1/4. Theoretical evidence lends credence to this conjecture, since $$\Delta(x)/x^{1/4}$$ has a (non-Gaussian) limiting distribution. The value of 1/4 would also follow from a conjecture on exponent pairs.

Piltz divisor problem
In the generalized case, one has


 * $$D_k(x) = xP_k(\log x)+\Delta_k(x) \,$$

where $$P_k$$ is a polynomial of degree $$k-1$$. Using simple estimates, it is readily shown that


 * $$\Delta_k(x)=O\left(x^{1-1/k} \log^{k-2} x\right)$$

for integer $$k\ge 2$$. As in the $$k=2$$ case, the infimum of the bound is not known for any value of $$ k$$. Computing these infima is known as the Piltz divisor problem, after the name of the German mathematician Adolf Piltz (also see his German page). Defining the order $$\alpha_k$$ as the smallest value for which $$\Delta_k(x)=O\left(x^{\alpha_k+\varepsilon}\right)$$ holds, for any $$\varepsilon>0$$, one has the following results (note that $$\alpha_2$$ is the $$\theta$$ of the previous section):


 * $$\alpha_2\le\frac{131}{416}\ ,$$


 * $$\alpha_3 \le\frac{43}{96}\ ,$$ and



\begin{align} \alpha_k & \le \frac{3k-4}{4k}\quad(4\le k\le 8) \\[6pt] \alpha_9 & \le\frac{35}{54}\ ,\quad \alpha_{10}\le\frac{41}{60}\ ,\quad \alpha_{11}\le\frac{7}{10} \\[6pt] \alpha_k & \le \frac{k-2}{k+2}\quad(12\le k\le 25) \\[6pt] \alpha_k & \le \frac{k-1}{k+4}\quad(26\le k\le 50) \\[6pt] \alpha_k & \le \frac{31k-98}{32k}\quad(51\le k\le 57) \\[6pt] \alpha_k & \le \frac{7k-34}{7k}\quad(k\ge 58) \end{align} $$
 * E. C. Titchmarsh conjectures that $$\alpha_k =\frac{k-1}{2k}\ .$$

Mellin transform
Both portions may be expressed as Mellin transforms:


 * $$D(x)=\frac{1}{2\pi i} \int_{c-i\infty}^{c+i\infty}

\zeta^2(w) \frac {x^w}{w}\, dw$$

for $$c>1$$. Here, $$\zeta(s)$$ is the Riemann zeta function. Similarly, one has


 * $$\Delta(x)=\frac{1}{2\pi i} \int_{c^\prime-i\infty}^{c^\prime+i\infty}

\zeta^2(w) \frac {x^w}{w} \,dw$$

with $$0<c^\prime<1$$. The leading term of $$D(x)$$ is obtained by shifting the contour past the double pole at $$w=1$$: the leading term is just the residue, by Cauchy's integral formula. In general, one has


 * $$D_k(x)=\frac{1}{2\pi i} \int_{c-i\infty}^{c+i\infty}

\zeta^k(w) \frac {x^w}{w} \,dw$$

and likewise for $$\Delta_k(x)$$, for $$k\ge 2$$.