Explicit bounds for Dickman's function

技术摘要：Dickman 函数的显式界限

问题陈述
Dickman 函数 $\rho(u)$ 是延迟微分方程 $u\rho'(u) + \rho(u-1) = 0$ 的连续解，其初始条件为在 $0 \le u \le 1$ 时 $\rho(u) = 1$。它代表了最大素因子不超过 $n^{1/u}$ 的整数 $n$ 的自然密度。虽然 $\rho(u)$ 在 $u \to \infty$ 时的渐近行为已通过 de Bruijn 和 Alladi 的工作得到充分理解，但这些结果通常涉及隐式误差项（例如 $O(1/u)$）。本文解决的主要挑战是为 $\rho(u)$ 建立数值显式的上界和下界，从而实现无需数值求解延迟微分方程的近似计算。

方法论
作者将鞍点法应用于 $\rho(u)$ 的逆拉普拉斯积分表示。该方法改进了 Tenenbaum 的证明技术，后者确立了具有 $O(u^{-1})$ 相对误差的 Alladi 结果。为了实现此处所需的更紧凑的显式界限，该方法论包括：

1. 精确误差追踪： 追踪所有的常数因子，而不是将其吸收进大 $O$ 符号中。
2. 精细积分： 通过将积分范围划分为额外的区间，以更精确地隔离各项贡献。
3. 高阶展开： 在鞍点附近利用七阶泰勒多项式（而非三阶）来捕捉展开式中的高阶项。
4. 显式引理： 推导关于积分函数 $I(s) = \int_0^s \frac{e^z-1}{z} dz$ 导数的特定不等式，并使用与不完全伽马函数相关的引理来限定积分的尾部贡献。

作为主要近似函数的 $\tilde{\rho}(u)$ 定义为：
$$\tilde{\rho}(u) := \frac{\exp\{\gamma - u\xi + I(\xi)\}}{\sqrt{2\pi I''(\xi)}}$$
其中 $\xi = \xi(u)$ 是唯一的正解，满足 $e^\xi = 1 + u\xi$，$\gamma$ 为欧拉常数。

主要贡献与结果
本文建立了两个关于 $\rho(u)$ 的显式界限和估计的主定理：

• 定理 1（显式界限）： 对于 $u > 1$，$\rho(u)$ 被限制在：
  $$\tilde{\rho}(u) \left(1 - \frac{1}{12u}\right) < \rho(u) < \tilde{\rho}(u) \left(1 - \frac{1}{12u(1 + \frac{1}{\log u})}\right)$$
  这提供了一个仅使用初等函数和 Lambert $W$ 函数（用于求解 $\xi$）即可计算 $\rho(u)$ 的便捷范围。

• 定理 2（带有误差项的精细估计）： 对于 $u > 1$，该函数可以表示为：
  $$\rho(u) = \tilde{\rho}(u) \left(1 + \alpha(u) + \frac{\theta(u)}{u^2}\right)$$
  其中 $\alpha(u)$ 是一个涉及 $I(\xi)$ 导数的显式表达式：
  $$\alpha(u) := \frac{1}{8}\frac{I^{(4)}(\xi)}{(I''(\xi))^2} - \frac{5}{24}\frac{(I^{(3)}(\xi))^2}{(I''(\xi))^3}$$
  至关重要的是，误差项被显式限定：对于所有 $u \ge 5$，有 $|\theta(u)| < 0.005$。这导致相对误差小于 $0.005/u^2$。

论文还提供了一种用于验证 $\rho(u)$ 数值计算的方法：

• 对于 $1 \le u \le 50$，使用 Marsaglia, Zaman, and Marsaglia 的泰勒多项式算法。
• 对于较大的 $u$，数值计算逆拉普拉斯积分，并注意到其被积函数收敛于标准正态密度函数。

意义与主张
本文声称，这些结果使得“无需数值求解延迟微分方程即可对 $\rho(u)$ 进行近似计算”。通过提供显式常数和误差界限，这项工作使 de Bruijn 和 Alladi 的渐近估计变得“完全显式”。

作者指出，尽管其他著作（如 Smida, Xuan）提供了误差量级更小的估计，但缺乏数值显式性。此处推导的显式界限在涉及 $\Psi(x, y)$（即最大素因子 $\le y$ 的不超过 $x$ 的整数计数）的应用中特别有用。具体而言，结合 Gorodetsky 的条件结果（假设黎曼假设），这些界限可以产生“易于数值评估”的 $\Psi(x, y)$ 条件下界。

论文通过将结果与 van de Lune 和 Wattel 的值进行对比来验证其正确性，后者使用梯形法则将 $\rho(u)$ 近似到至少五个有效数字（对于 $u \le 1000$）。定理 1 和定理 2 导出的值与这些已建立的数值基准是一致的。