対数積分(logarithmic integral)とは、$1/\log t$ を積分して定まる特殊関数 $\operatorname{li}(x)$ のことであり、$x>1$ では $t=1$ の特異点を Cauchy の主値で処理する。$2$ から積分した $\operatorname{Li}(x)$ もよく用いられる。部分積分により漸近展開 $x/\log x+x/\log^2x+\cdots$ をもち、特に $x/\log x$ と漸近的に等しい。素数定理は $x$ 以下の素数の個数 $\pi(x)$ が $\operatorname{li}(x)$ で近似されることを述べ、これは $x/\log x$ よりはるかによい近似である。計算できる範囲では常に $\operatorname{li}(x)>\pi(x)$ だが、差の符号は無限回変わる。
前提知識: 広義積分, 対数関数, Landau記号, 素数計数関数
正の実数 $x\ne1$ に対し、対数積分(logarithmic integral)$\operatorname{li}(x)$ を次で定義する。
$x>1$ に対し、def-logarithmic-integral の極限は存在し、有限である。さらに $x\mapsto\operatorname{li}(x)$ は $(0,1)$ と $(1,\infty)$ のそれぞれで微分可能(微分)で、$\operatorname{li}'(x)=1/\log x$ である。
$t>0$、$t\ne1$ に対し
$$g(t):=\frac{1}{\log t}-\frac{1}{t-1}=\frac{(t-1)-\log t}{(t-1)\log t}$$
とおく。$\log t=(t-1)-\tfrac12(t-1)^2+O((t-1)^3)$($t\to1$)により分子は $\tfrac12(t-1)^2+O((t-1)^3)$、分母は $(t-1)^2+O((t-1)^3)$ なので、$g(t)\to\tfrac12$($t\to1$)である。よって $g$ は $g(1):=\tfrac12$ と定めれば $(0,\infty)$ で連続であり、$t\to0+$ で $g(t)\to0-(-1)=1$ なので $(0,x]$ で有界である。したがって $\int_0^x g(t)\,dt$ は通常の広義積分として収束する。一方、
$$\int_0^{1-\varepsilon}\frac{dt}{t-1}+\int_{1+\varepsilon}^{x}\frac{dt}{t-1}=\Bigl[\log|t-1|\Bigr]_0^{1-\varepsilon}+\Bigl[\log(t-1)\Bigr]_{1+\varepsilon}^{x}=\log\varepsilon-0+\log(x-1)-\log\varepsilon=\log(x-1)$$
は $\varepsilon$ によらない。$1/\log t=1/(t-1)+g(t)$ を代入すれば、極限は存在して
$$\operatorname{li}(x)=\log(x-1)+\int_0^x g(t)\,dt$$
である。最後の主張について、$x>1$ のときは上の表示を微分して $\operatorname{li}'(x)=1/(x-1)+g(x)=1/\log x$ であり、$0< x<1$ のときは定義から直接、微積分学の基本定理により $\operatorname{li}'(x)=1/\log x$ である。$\square$
指数積分 $\operatorname{Ei}(u)$($u\ne0$ に対する $\int_{-\infty}^{u}e^{s}\,ds/s$ の主値)を用いると、$t=e^{s}$ の置換積分により $\operatorname{li}(x)=\operatorname{Ei}(\log x)$ である(DLMF §6.2(i))。$\operatorname{li}(x)$ は $x\to1$ で $-\infty$ に発散し、$x=1.45136\ldots$(Ramanujan–Soldner 定数)でただ一度 $0$ になり、$x>1$ で単調増加(単調関数)である。文献によっては $\operatorname{Li}(x)$ を $\operatorname{li}(x)$ と書くこともあるが、両者の差は定数 $\operatorname{li}(2)$ なので、漸近的な議論では区別する必要がない。
$x>e$ のとき
$$\frac{x}{\log x}-e+\frac{x-e}{\log^2x}<\operatorname{li}(x)-\operatorname{li}(e)<\frac{x}{\log x}+\frac{4x}{\log^2x}+\sqrt{x}-e$$
が成り立つ。したがって
$$\operatorname{li}(x)=\frac{x}{\log x}+O\!\left(\frac{x}{\log^2x}\right),\qquad \operatorname{li}(x)\sim\frac{x}{\log x}\quad(x\to\infty)$$
である。ここで $O$ は Landau記号 である。
$x>e$ とする。部分積分により
$$\operatorname{li}(x)-\operatorname{li}(e)=\int_e^x\frac{dt}{\log t}=\left[\frac{t}{\log t}\right]_e^x+\int_e^x\frac{dt}{\log^2t}=\frac{x}{\log x}-e+\int_e^x\frac{dt}{\log^2t}$$
である。$e\le t< x$ で $1/\log^2t>1/\log^2x$ なので $\int_e^x dt/\log^2t>(x-e)/\log^2x$ であり、下からの評価を得る。
上からの評価について、$x\ge e^2$ のときは積分を $\sqrt{x}$ で分ける。$t\ge e$ で $1/\log^2t\le1$ なので $\int_e^{\sqrt{x}}dt/\log^2t\le\sqrt{x}-e$ であり、$t\ge\sqrt{x}$ で $1/\log^2t\le1/\log^2\sqrt{x}=4/\log^2x$ なので $\int_{\sqrt{x}}^xdt/\log^2t\le4(x-\sqrt{x})/\log^2x<4x/\log^2x$ である。合わせて主張の右辺を得る。$e< x< e^2$ のときは $\int_e^xdt/\log^2t\le x-e$ であり、$\log^2x<4$ より $x<4x/\log^2x$、よって $x-e<4x/\log^2x+\sqrt{x}-e$ となり、やはり右辺を得る。
$\operatorname{li}(e)$ は定数であり、$\sqrt{x}=O(x/\log^2x)$ なので、両側の評価から $\operatorname{li}(x)=x/\log x+O(x/\log^2x)$ が従い、これを $x/\log x$ で割れば $\operatorname{li}(x)/(x/\log x)\to1$ である。$\square$
非負整数 $r$ に対し、$x>1$ で
$$\operatorname{li}(x)=\frac{x}{\log x}\sum_{k=0}^{r}\frac{k!}{\log^kx}+(r+1)!\int_2^x\frac{dt}{\log^{r+2}t}+c_r,\qquad c_r:=\operatorname{li}(2)-2\sum_{k=0}^{r}\frac{k!}{\log^{k+1}2}$$
が成り立ち、剰余項について
$$(r+1)!\int_2^x\frac{dt}{\log^{r+2}t}+c_r=O\!\left(\frac{x}{\log^{r+2}x}\right)\qquad(x\to\infty)$$
である($O$ の定数は $r$ に依存する)。したがって
$$\operatorname{li}(x)\sim\frac{x}{\log x}\left(1+\frac{1!}{\log x}+\frac{2!}{\log^2x}+\frac{3!}{\log^3x}+\cdots\right)$$
は $\operatorname{li}(x)$ の漸近展開である。
$m\ge1$ に対し $I_m(x):=\int_2^x dt/\log^mt$ とおく。$\dfrac{d}{dt}\dfrac{t}{\log^mt}=\dfrac{1}{\log^mt}-\dfrac{m}{\log^{m+1}t}$ なので、部分積分により
$$I_m(x)=\frac{x}{\log^mx}-\frac{2}{\log^m2}+m\,I_{m+1}(x)$$
である。これを $m=1,2,\dots,r+1$ について順に代入すると、$r$ に関する数学的帰納法により
$$I_1(x)=\sum_{k=0}^{r}k!\left(\frac{x}{\log^{k+1}x}-\frac{2}{\log^{k+1}2}\right)+(r+1)!\,I_{r+2}(x)$$
を得る($r=0$ は上の等式そのもの、$r\to r+1$ は $I_{r+2}$ に上の等式を適用する)。$\operatorname{li}(x)=\operatorname{li}(2)+I_1(x)$ なので第 1 の等式を得る。
剰余項について、$x\ge4$ とし、積分を $\sqrt{x}$ で分ける。$t\ge2$ で $1/\log^{r+2}t\le1/\log^{r+2}2$ なので $\int_2^{\sqrt{x}}dt/\log^{r+2}t\le\sqrt{x}/\log^{r+2}2=O(\sqrt{x})$ であり、$\sqrt{x}=O(x/\log^{r+2}x)$ である。$t\ge\sqrt{x}$ では $1/\log^{r+2}t\le2^{r+2}/\log^{r+2}x$ なので $\int_{\sqrt{x}}^xdt/\log^{r+2}t\le2^{r+2}x/\log^{r+2}x$ である。定数 $c_r$ も $O(x/\log^{r+2}x)$ なので、剰余項の評価を得る。$\square$
thm-logarithmic-integral-asymptotic-expansion の級数 $\sum_{k\ge0}k!/\log^kx$ は、どのような $x>1$ に対しても収束しない。$k!/\log^kx$ は $k>\log x$ で増大するからである。漸近展開とは、各 $r$ について有限和で打ち切った誤差が最後の項より小さい位数であることを意味し、級数の収束を意味しない。実際、$x$ を固定して $r$ を増やすと剰余項 $(r+1)!\,I_{r+2}(x)$ はいずれ増大する。
$x>1$ に対し、Eulerの定数 $\gamma$ を用いて
$$\operatorname{li}(x)=\gamma+\log\log x+\sum_{n=1}^{\infty}\frac{(\log x)^n}{n\cdot n!}$$
が成り立つ。また Ramanujan による次の級数表示が成り立つ。
$$\operatorname{li}(x)=\gamma+\log\log x+\sqrt{x}\sum_{n=1}^{\infty}\frac{(-1)^{n-1}(\log x)^n}{n!\,2^{n-1}}\sum_{k=0}^{\lfloor(n-1)/2\rfloor}\frac{1}{2k+1}.$$
第 1 の級数は $\operatorname{li}(x)=\operatorname{Ei}(\log x)$ と指数積分の級数 $\operatorname{Ei}(u)=\gamma+\log u+\sum_{n\ge1}u^n/(n\cdot n!)$($u>0$)から従う(DLMF §6.6、式 6.6.2)。第 2 の級数は Ramanujan によるもので、証明は Ber94 pp. 126–131 に譲る(式の形は MathWorldLi にもある)。Ramanujan の級数は第 1 の級数より収束が速く、数値計算に適する。
$x=10^n$ に対する $\operatorname{li}(x)$、$x/\log x$、および $x$ 以下の素数の個数 $\pi(x)$(素数計数関数)を比べると次の通りである。$\operatorname{li}(x)$ は thm-logarithmic-integral-series の級数から小数第 1 位まで求め、$\pi(x)$ の値は OEIS6880 による。
| $x$ | $\operatorname{li}(x)$ | $x/\log x$ | $\pi(x)$ |
|---|---|---|---|
| $10$ | $6.2$ | $4.3$ | $4$ |
| $10^2$ | $30.1$ | $21.7$ | $25$ |
| $10^3$ | $177.6$ | $144.8$ | $168$ |
| $10^4$ | $1246.1$ | $1085.7$ | $1229$ |
| $10^5$ | $9629.8$ | $8685.9$ | $9592$ |
| $10^6$ | $78627.5$ | $72382.4$ | $78498$ |
| $10^7$ | $664918.4$ | $620420.7$ | $664579$ |
| $10^8$ | $5762209.4$ | $5428681.0$ | $5761455$ |
$\operatorname{li}(x)$ は $x/\log x$ よりはるかによく $\pi(x)$ を近似し、この範囲では常に $\operatorname{li}(x)>\pi(x)$ である。
$x\to\infty$ のとき
$$\operatorname{li}(x)-\frac{x}{\log x}\sim\frac{x}{\log^2x}$$
であり、特に $\operatorname{li}(x)-x/\log x\to\infty$ である。
thm-logarithmic-integral-asymptotic-expansion で $r=1$ とすると $\operatorname{li}(x)=x/\log x+x/\log^2x+O(x/\log^3x)$ であり、これを $x/\log^2x$ で割ればよい。$\square$
1 は 2 から従い、また prop-logarithmic-integral-basic-bound により $\pi(x)\sim x/\log x$ とも同値である。2 で $\alpha=1/2$ がとれることは de la Vallée Poussin による(Dav80 Chapter 18)。3 は Ford(For02)、4 は Trudgian(Tru16)による明示的評価である。5 は von Koch による(MV07 §13.1)。いずれも証明は各文献に譲る。
$\theta(x):=\sum_{p\le x}\log p$ を Chebyshev関数 とする。$x\ge2$ で $1\le E(x)\le x$ を満たす増加関数 $E$ について
$$\theta(x)=x+O(E(x))$$
が成り立つならば、
$$\pi(x)=\operatorname{li}(x)+O\!\left(\frac{E(x)}{\log x}+\int_2^x\frac{E(t)}{t\log^2t}\,dt\right)$$
が成り立つ。
$x\ge2$ とする。Abel の総和公式(Abelの級数変形法とその応用)を $a_n=\log n$($n$ が素数のとき)、$a_n=0$(それ以外)、$f(t)=1/\log t$ に適用すると、$A(t)=\sum_{n\le t}a_n=\theta(t)$ であり
$$\pi(x)=\sum_{p\le x}\log p\cdot\frac{1}{\log p}=\frac{\theta(x)}{\log x}+\int_2^x\frac{\theta(t)}{t\log^2t}\,dt$$
となる。ここに $\theta(t)=t+O(E(t))$ を代入すると
$$\pi(x)=\frac{x}{\log x}+\int_2^x\frac{dt}{\log^2t}+O\!\left(\frac{E(x)}{\log x}+\int_2^x\frac{E(t)}{t\log^2t}\,dt\right)$$
である。一方、thm-logarithmic-integral-asymptotic-expansion の証明の等式($m=1$)により
$$\operatorname{Li}(x)=\int_2^x\frac{dt}{\log t}=\frac{x}{\log x}-\frac{2}{\log2}+\int_2^x\frac{dt}{\log^2t}$$
なので、$\operatorname{li}(x)=\operatorname{Li}(x)+\operatorname{li}(2)$ と合わせて $\pi(x)=\operatorname{li}(x)+O(1)+O\!\left(E(x)/\log x+\int_2^xE(t)\,dt/(t\log^2t)\right)$ を得る。$E\ge1$ より $\int_2^x E(t)\,dt/(t\log^2t)\ge1/\log2-1/\log x$ であり、これは $x\ge4$ で正の定数以上、$2\le x<4$ では $E(x)/\log x\ge1/\log4$ なので、誤差項全体は下から正の定数で抑えられ、$O(1)$ はこれに吸収される。$\square$
1 は Littlewood による(Lit14)。ex-logarithmic-integral-values のように計算できる範囲では常に $\operatorname{li}(x)>\pi(x)$ であるが、それが永久に続くわけではない。最初に符号が変わる $x$ の上界は Skewes が Riemann 予想のもとで $10^{10^{10^{34}}}$ を与えたのが最初であり(Ske33)、この種の上界を Skewes 数という。2 の上界は Bays–Hudson による(BH00)。いずれも証明は各論文に譲る。
$\operatorname{li}(x)$ の数値計算には thm-logarithmic-integral-series の級数が使える。$\operatorname{Li}(x)=\operatorname{li}(x)-\operatorname{li}(2)$ は $\pi(x)$ の近似としては $\operatorname{li}(x)$ より約 $1.045$ 小さいだけであり、thm-logarithmic-integral-prime-number-theorem の誤差項に比べて無視できる。一方、$x/\log x$ による近似は prop-logarithmic-integral-difference により $\operatorname{li}(x)$ より $x/\log^2x$ 程度悪く、これは $\pi(x)-\operatorname{li}(x)$ の大きさ(thm-logarithmic-integral-prime-number-theorem の 2)よりはるかに大きい。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する