We'll prove that for every $\varepsilon>0$, there exists $N$ such that for $n>N$, $d(n)<n^\varepsilon$. For each $k\in\mathbb Z_0^+$, let $p_k$ be the $k$-th prime, and let $f_k:\mathbb R_0^+\to\mathbb R$ be the function such that for $x\in\mathbb R_0^+$, $$f_k(x)=\frac{x+1}{p_k^{\varepsilon x}}.$$ We'll make five key observations, which will together imply our main theorem.
- For every $k\in\mathbb Z^+$, $f_k(x)$ is increasing for $x<\frac{1}{\varepsilon\log\left(p_k\right)}-1$, and decreasing otherwise.
We can verify that $$\frac{\mathrm d}{\mathrm dx}\,f_k(x)=\frac{1-\varepsilon\log(p_k)(x+1)}{p_k^{\varepsilon x}}.$$ The result follows easily.
- For every $k\in\mathbb Z^+$, there exists $m\in\mathbb Z_0^+$ such that for every $n\in\mathbb Z_0^+$, $f_k(m)\ge f_k(n)$.
Because of 1, $f_k(x)$ attains a maximum value at $x=\frac{1}{\varepsilon\log(p_k)}-1$, is increasing to the left of it, and decreasing to the right of it. Clearly, either the floor or the ceiling of this number will satisfy the required condition. Let the maximum value of $f_k$ over the non-negative integers be $M_k$.
- For sufficiently large $k$, $M_k=1$.
If $p_k\ge e^{1/\varepsilon}$, because of $1$, $f_k(x)$ will be decreasing, and its maximum value will be $f_k(0)=1$. We can therefore define $$M=\prod_{i=1}^\infty M_i,$$ as the product will only have finitely many terms different from $1$.
- For all $k\in\mathbb Z^+$, there exists $\beta_k'\in\mathbb Z_0^+$ such that for $x\geq\beta_k'$, $f_k(x)<\frac{1}{M}$.
By L'Hopital, $$\lim_{x\to\infty}\frac{x+1}{p_k^{\varepsilon x}}=\lim_{x\to\infty}\frac{1}{\varepsilon\log(p_k)p_k^{\varepsilon x}}=0.$$ The result follows easily. Let the least possible value of $\beta_k'$ be $\beta_k$.
- For sufficiently large $k$, $\beta_k\le 1$.
If $p_k>\max(e^{1/\varepsilon},(2M)^{1/\varepsilon})$, for all $x\ge 1$, $$f_k(x)\le f_k(1)<\frac{1}{M},$$ so that $\beta_k\le 1$. We can therefore define $$B=\prod_{i=1}^\infty\beta_i,$$ as the product will have only finitely many terms different from either $0$ or $1$.
Now, consider $n$ such that $d(n)\ge n^\varepsilon$, and let $n=\prod_{i=1}^\infty p_i^{\alpha_i}$ be its prime factorization. Since $d(n)=\prod_{i=1}^\infty(\alpha_i+1)$, this implies $$\prod_{i=1}^\infty f_i(\alpha_i)\ge1.$$ By 5, it's clear that for any $i$, $\alpha_i$ can only take integer values from $0$ to $\beta_i-1$, since otherwise, the product would be less than $M\cdot\frac{1}{M}=1$. There are therefore at most $B$ possible choices for the exponents $\alpha_i$ – taking $N$ greater than all the values of $n$ they generate is enough to finish the problem. $\blacksquare$