対数積分

同義語:logarithmic integral

概要

対数積分(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)$ だが、差の符号は無限回変わる。

$$\newcommand{AA}[0]{\mathscr{A}} \newcommand{abs}[1]{\left\lvert#1\right\rvert} \newcommand{Arg}[0]{\operatorname{Arg}} \newcommand{BB}[0]{\mathscr{B}} \newcommand{C}[0]{\mathbb{C}} \newcommand{CC}[0]{\mathscr{C}} \newcommand{floor}[1]{\left\lfloor#1\right\rfloor} \newcommand{ind}[0]{\operatorname{ind}} \newcommand{Ker}[0]{\operatorname{Ker}} \newcommand{Li}[0]{\operatorname{Li}} \newcommand{li}[0]{\operatorname{li}} \newcommand{mmod}[1]{\ \left(\mathrm{mod}\ #1\right)} \newcommand{Mod}[1]{\ \left(\mathrm{mod}\ #1\right)} \newcommand{N}[0]{\mathbb{N}} \newcommand{ord}[0]{\operatorname{ord}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{rank}[0]{\mathrm{rank}} \newcommand{SS}[0]{\mathscr{S}} \newcommand{TT}[0]{\mathscr{T}} \newcommand{UU}[0]{\mathscr{U}} \newcommand{wenvert}[1]{\left\lvert\left\lvert#1\right\rvert\right\rvert} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 広義積分, 対数関数, Landau記号, 素数計数関数

定義

対数積分

正の実数 $x\ne1$ に対し、対数積分(logarithmic integral)$\operatorname{li}(x)$ を次で定義する。

  • $0< x<1$ のとき、$\operatorname{li}(x):=\displaystyle\int_0^x\frac{dt}{\log t}$(広義積分。被積分関数は $(0,x]$ で連続かつ負であり、$t\to0+$ で $0$ に収束するので収束する)。
  • $x>1$ のとき、被積分関数 $1/\log t$ は $t=1$ で特異点をもつので、Cauchyの主値
    $$\operatorname{li}(x):=\lim_{\varepsilon\to0+}\left(\int_0^{1-\varepsilon}\frac{dt}{\log t}+\int_{1+\varepsilon}^{x}\frac{dt}{\log t}\right)$$
    で定義する(極限の存在は prop-logarithmic-integral-well-defined)。
    また、特異点を避けて $2$ から積分した
    $$\operatorname{Li}(x):=\int_2^x\frac{dt}{\log t}=\operatorname{li}(x)-\operatorname{li}(2)\qquad(x>1)$$
    を Euler の対数積分、あるいはオフセット対数積分(offset logarithmic integral)という。$\operatorname{li}(2)=1.04516\ldots$ である。
主値の存在

$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/log 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. $\pi(x)\sim\operatorname{li}(x)$($x\to\infty$)である(素数定理)。
  2. ある正の定数 $c,\alpha$ が存在して $\pi(x)=\operatorname{li}(x)+O\!\left(x\exp(-c\log^{\alpha}x)\right)$ が成り立つ。$\alpha=1/2$ にとれる。
  3. $\pi(x)=\operatorname{li}(x)+O\!\left(x\exp\!\left(-0.2098\,\dfrac{(\log x)^{3/5}}{(\log\log x)^{1/5}}\right)\right)$ が成り立つ。
  4. $x\ge229$ のとき $\bigl|\pi(x)-\operatorname{li}(x)\bigr|\le0.2795\,\dfrac{x}{(\log x)^{3/4}}\exp\!\left(-\sqrt{\dfrac{\log x}{6.455}}\right)$ が成り立つ。
  5. Riemann予想を仮定すると $\pi(x)=\operatorname{li}(x)+O(\sqrt{x}\log x)$ が成り立つ。
素数定理の出典

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)。いずれも証明は各文献に譲る。

Chebyshev 関数の評価からの導出

$\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$

差 li(x)−π(x) の符号

Littlewood の定理と Skewes 数
  1. 差 $\operatorname{li}(x)-\pi(x)$ は $x\to\infty$ で無限回符号を変える。すなわち $\pi(x)>\operatorname{li}(x)$ となる $x$ が無限に存在する。
  2. $\pi(x)>\operatorname{li}(x)$ となる $x$ が $1.3982\times10^{316}$ より小さい範囲に存在する。
Littlewood の定理と Skewes 数の出典

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)よりはるかに大きい。

関連項目

参考文献

[3]
Bruce C. Berndt, Ramanujan's Notebooks, Part IV, Springer-Verlag, 1994, pp. 126–131(対数積分の Ramanujan の級数)
[6]
J. E. Littlewood, Sur la distribution des nombres premiers, C. R. Acad. Sci. Paris, 1914, 1869--1872

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する