本記事は「 素数計数関数の恒等式 」シリーズの第7回です。この記事だけでも読めるよう、これまでのおさらいから始めます。詳細は各記事を参照してください。本記事の主要な結果として、$N$以下の素数の個数を数える素数計数関数$\pi(N)$の近似式を得たことです。
$\pi(N)\approx\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$
\begin{align}
\tilde{\pi}(N)&= \sum_{\ell=1}^{\lfloor \log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{3}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr), \\
\pm\delta\tilde{\pi}(N)&=\pm\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\biggl(-\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{1}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr).
\end{align}
ここで$\mathrm{Li}(x)$は補正対数積分関数であり、
\begin{align}
\mathrm{Li}(x)&\coloneqq\mathrm{li}(x)-\mathrm{li}(2), \\
\mathrm{li}(x)&\coloneqq\int_{0}^{x}\frac{dt}{\ln t}
\end{align}
で定義されます。おさらいの後で、この近似式の数値的な諸結果を示します。最後に4つの近似の下で、素数計数関数の恒等式から目的の近似式を導出します。
まず
第1回
で素数計数関数$\pi(N)$の恒等式を考えました。
第2回
では、それを複素数の$s$乗に一般化しました。
結果は、与えられた2以上の実数$N$について、
\begin{equation}
H_s(N)-1=\displaystyle \sum_{\ell=1}^{\lfloor \log_2{N} \rfloor} \sum_{m=1}^{\lfloor N/2^\ell \rfloor} \frac{\eta_{\ell m}}{m^s}q_{\ell s} \Bigl(\sqrt[\ell]{N/m} \Bigr),\quad \eta_{\ell m}=\frac{1}{\Omega(m)+\ell}
\end{equation}
が成り立つことでした。ここで$\eta_{\ell m}$を展開係数と呼び、
\begin{align}
H_s(N) &\coloneqq \displaystyle \sum_{k=1}^{\lfloor N \rfloor} \frac{1}{k^s}, \\
q_s(N) &\coloneqq \displaystyle \sum_{p \leq N} \frac{1}{p^s}
\end{align}
と定義し、$p$は素数を表します。
第3回
では、素数計数関数の恒等式を反転させ
\begin{equation}
\pi(N)=\displaystyle \sum_{\ell=1}^{\lfloor \log_2{N} \rfloor} \sum_{m=1}^{\lfloor N/2^\ell \rfloor} \mu_{\ell m} \Bigl( \Bigl\lfloor \sqrt[\ell]{N/m} \Bigr\rfloor -1\Bigr)
\end{equation}
を得て、反転係数$\mu_{\ell m}$の計算規則も求めました。
第4回
は第3回の後編であり、本記事と密接に関係する回です。
まず反転式の床関数$\lfloor x \rfloor$を$x,x-1,x-1/2$で近似しました(近似1)。新たに得た近似式をそれぞれ$\tilde{\pi}_0(N),\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$と定義し、アーベルの総和公式を用いて
\begin{align} \tilde{\pi}_0(N)&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\biggl(S_\ell\biggl( \frac{N}{2^\ell}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{1}{\ell}\int_{1}^{ N/2^\ell }S_\ell(t)t^{1/\ell-1}dt\biggr), \\ \tilde{\pi}_1(N)&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\frac{2}{\ell}\int_{1}^{ N/2^\ell }S_\ell(t)t^{1/\ell-1}dt, \\ \tilde{\pi}_{1/2}(N)&=\frac{\tilde{\pi}_0(N)+\tilde{\pi}_1(N)}{2} \end{align}
と表しました。ここで、$S_\ell(x)$は
\begin{equation}
S_\ell(x)\coloneqq\sum_{m=1}^{\lfloor x \rfloor} \frac{\mu_{\ell m}}{\sqrt[\ell]{m}}
\end{equation}
と定義しています。
ここで定義した$\tilde{\pi}_0(N),\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$は、 第4回 では$\pi_0(N),\pi_1(N),\pi_{1/2}(N)$と定義していました。しかし、後に出てくる$\pi_0(N)\coloneqq (\pi(x+0)+\pi(x-0))/2$との混乱を避けるため、以降は近似1の近似式をチルダ付きで表すことにします。
第4回 の図2として示した$[2,128)$における近似式の計算結果を図1に再掲します。
素数計数関数と3種の近似式の比較
$\tilde{\pi}_0(N)$が素数計数関数$\pi(N)$をよく近似していることが分かります。また
第4回
の命題5では素数定理から
\begin{align}
S_1(N)\approx\frac{1}{\ln 2N}
\end{align}
と予想しました。この近似の妥当性を確認するため、$S_1(N)$と$1/\ln 2N$を比較した
第4回
の図3を図2として再掲します。
$S_1(N)$の近似の検証
第4回
で素数定理から予想したこの近似を本記事では導出します。
さらに命題7では、素数計数関数の反転式を複素数の$s$乗へと一般化し、
\begin{equation}
q_s(N)=\displaystyle \sum_{\ell=1}^{\lfloor \log_2{N} \rfloor} \sum_{m=1}^{\lfloor N/2^\ell \rfloor} \frac{\mu_{\ell m}}{m^s}\Bigl(H_{\ell s} \Bigl(\sqrt[\ell]{N/m} \Bigr)-1\Bigr)
\end{equation}
も得ました。
素数計数関数の恒等式は離散的で有限の理論であったため、
第5回
で測度を用いて連続性を議論し、
第6回
で無限極限を考えました。特に、本記事でも無限極限の式から一意に定まる反転係数$\mu_{\ell m}$を使用します。
$n \ge 1$に対し、$m_n=p_1\cdots p_n$と定義する。
ここで、$p_1,\dots,p_n$は互いに相異なる素数である。また$m_0=1$と定める。
以下の展開係数と対応する反転係数は、素因数分解の順序交換を行っても変わらない。
\begin{align} &\eta_{1,m_{n}}=\frac{1}{n+1},\quad \eta_{\ell,1}=\frac{1}{\ell}, \\ &\mu_{1,m_{n}}=B_n,\quad \mu_{\ell,1}=\frac{\mu(\ell)}{\ell}. \end{align}
ここで、$B_n$は第一種ベルヌーイ数、$\mu(\cdot)$はメビウス関数である。
本記事の近似式のうち、特別な場合は
\begin{align}
\tilde{\pi}_{1/2}(N)&\approx \tilde{\pi}(N)=\sum_{\ell=1}^{\lfloor \log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{3}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr), \\
\tilde{\pi}_{0}(N)&\approx\tilde{\pi}(N)-\delta\tilde{\pi}(N)=\sum_{\ell=1}^{\lfloor \log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{2}\frac{\sqrt[\ell]{N}}{\ln\sqrt[\ell]{N}}+\frac{1}{2}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr), \\
\tilde{\pi}_{1}(N)&\approx\tilde{\pi}(N)+\delta\tilde{\pi}(N)=\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr),
\end{align}
と表されます。これらの近似式と素数計数関数$\pi(N)$を比較していきます。$N=128$までを図3に、$N=1000$までを図4に示します。
素数計数関数$\pi(N)$との比較($N=128$まで)
図3では、青色の素数計数関数$\pi(N)$が緑色の$\tilde{\pi}_0(N)$より大きい値を取っています。また、橙色の$\tilde{\pi}_1(N)$が$\pi(N)$の概ね中間を通っており、$\pi(N)$をよく近似しています。近似1のみを使った図1では$\tilde{\pi}_{1/2}(N)$が$\pi(N)$をよく近似していたため、相違が見られます。
素数計数関数$\pi(N)$との比較($N=1000$まで)
図4でも$\tilde{\pi}_1(N)$が$\pi(N)$に最も近く、$\pi(N)$は$\tilde{\pi}_1(N)$を中心として正負に揺らいでいます。
次に、$N=10^1\sim10^8$まで$\pi(N)$と近似式$\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$の数値を比較します。相対誤差を
\begin{align}
\text{相対誤差}\coloneqq\biggl(\frac{(\text{推定値}\pm\text{不確かさ})-\text{真の値}}{\text{真の値}}\biggr)\times100 \%
\end{align}
と定義し、近似式が真の値$\pi(N)$を含むかどうかを判定します。これらの結果を表1にまとめました。
表1 近似の様子(+相対誤差、判定)
| $N$ | $\pi(N)$ | $\tilde{\pi}(N)$ | $\delta\tilde{\pi}(N)$ | 相対誤差 / % | 判定 |
|---|---|---|---|---|---|
| $10^1$ | 4 | 3.822 | 0.596 | -4.439$\pm$14.904 | ○ |
| $10^2$ | 25 | 23.846 | 1.796 | -4.617$\pm$7.184 | ○ |
| $10^3$ | 168 | 161 | 8 | -4.276$\pm$4.481 | ○ |
| $10^4$ | 1,229 | 1,188 | 39 | -3.325$\pm$3.157 | × |
| $10^5$ | 9,592 | 9,354 | 233 | -2.481$\pm$2.433 | × |
| $10^6$ | 78,498 | 76,971 | 1,556 | -1.945$\pm$1.983 | ○ |
| $10^7$ | 664,579 | 653,553 | 11,114 | -1.659$\pm$1.672 | ○ |
| $10^8$ | 5,761,455 | 5,678,192 | 83,360 | -1.445$\pm$1.447 | ○ |
この表から、近似式の値の桁数と最高位の数字は$\pi(N)$と一致していることが読み取れます。そして相対誤差は$N$の増加とともに減少し、幅も同程度になっています。$10^8$までは、ほとんどの近似式の値が真の値を含みますが、$10^4,10^5$では含まないことが分かります。
近似の様子について調べたので、次は既存の近似式と比較します。$N=10^8$における$N/\ln N,\mathrm{li}(N),R(N),\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$の相対誤差を表2にまとめました。ここで、$R(N)$はリーマンの素数計数関数であり、
\begin{align}
R(N)\coloneqq \sum_{\ell=1}^\infty \frac{\mu(\ell)}{\ell}\mathrm{li}\bigl(N^{1/\ell}\bigr)
\end{align}
と定義されます。
表2 $N=10^8$における既存の近似式との比較
| 近似式 | 相対誤差 / % |
|---|---|
| $N/\ln N $ | -5.776 |
| $\mathrm{li}(N)$ | 0.013 |
| $R(N)$ | 0.002 |
| $\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$ | -1.445$\pm$1.447 |
表2より、$\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$の相対誤差は、素数定理に現れる$N/\ln N$より小さくなっています。$\tilde{\pi}_1(N)(=\tilde{\pi}(N)+\delta\tilde{\pi}(N))$は$R(N)$と定数項を除いて等しく、相対誤差はどちらも0.002%です。そのため$\tilde{\pi}_1(N)$であれば$\mathrm{li}(N)$よりも相対誤差が小さいと分かります。したがって近似式$\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$は$N/\ln N$よりも精度がよく、特別な場合に$R(N)$と等しいため、それらの中間に位置づけられます。
ここからは$\pi(N)$の近似式を理論的に導出します。 第4回 では$\tilde{\pi}_0(N),\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$が$S_\ell(x)$で表されていたので、その近似式から求めていきます。
$\ell\ge2$の場合はリーマンの近似公式を参考にして、$S_\ell(x)$を$S_1(x)$で近似します。
リーマンの近似公式から導かれる仮定の下、2以上の自然数$\ell$、1以上の実数$x$について、$S_\ell(x)$を
\begin{equation}
S_\ell(x)\approx\frac{\mu(\ell)}{\ell}S_1\bigl(\sqrt[\ell]{x}\bigr)
\end{equation}
と近似できる。
まず$S_\ell(x)$は次のように定義される。
\begin{align}
S_\ell(x)\coloneqq \sum_{m=1}^{\lfloor x \rfloor}\frac{\mu_{\ell m}}{\sqrt[\ell]{m}}=\mu_{\ell 1}+\frac{\mu_{\ell2}}{\sqrt[\ell]{2}}+\cdots.
\end{align}
prinfの命題8より、反転係数$\mu_{\ell1}$は展開の仕方に依らず
\begin{align}
\mu_{\ell 1}=\frac{\mu(\ell)}{\ell}
\end{align}
と定まる。そこで新たな関数$S_\ell^*(x)$を導入する。
\begin{align}
S_\ell(x)=:\frac{\mu(\ell)}{\ell}S_\ell^*(x).
\end{align}
これをprin2の定理3で示した$\tilde{\pi}_{1/2}(N)$の表式に代入する。
\begin{align}
\tilde{\pi}_{1/2}(N)&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\biggl(\frac{1}{2}S_\ell\biggl( \frac{N}{2^\ell}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{3}{2\ell}\int_{1}^{ N/2^\ell }S_\ell(t)t^{1/\ell-1}dt\biggr)\\
&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{2}S_\ell^*\biggl( \frac{N}{2^\ell}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{3}{2\ell}\int_{1}^{ N/2^\ell }S_\ell^*(t)t^{1/\ell-1}dt\biggr) \\
&=:\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\frac{\mu(\ell)}{\ell}f_\ell(N).
\end{align}
最後の行で$\mu(\ell)/\ell$以外の部分を$f_\ell(N)$とした。
ところで、リーマンは$\pi_0(x)$の有名な近似式
\begin{align}
\pi_0(x)\approx R(x)=\sum_{\ell=1}^\infty \frac{\mu(\ell)}{\ell}\mathrm{li}\bigl(x^{1/\ell}\bigr)
\end{align}
を得たRiesel。ここで$\pi_0(x)$は
\begin{align}
\pi_0(x)\coloneqq \frac{\pi(x+0)+\pi(x-0)}{2}
\end{align}
と定義され、$x$が素数のとき$\pi(x)-1/2$、それ以外では$\pi(x)$と等しい関数である。
このリーマンの近似式と$\tilde{\pi}_{1/2}(N)$を比較して
\begin{align}
f_\ell(N) \approx f_1\bigl(\sqrt[\ell]{N}\bigr)
\end{align}
と仮定することができる。
ここからは$f_\ell(N)$と$f_1\bigl(\sqrt[\ell]{N}\bigr)$を具体的に計算する。
\begin{align}
f_\ell(N)&=\frac{1}{2}S_\ell^*\biggl( \frac{N}{2^\ell}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{3}{2\ell}\int_{1}^{ N/2^\ell }S_\ell^*(t)t^{1/\ell-1}dt, \\
f_1\bigl(\sqrt[\ell]{N}\bigr)&=\frac{1}{2}S_1^*\biggl( \frac{\sqrt[\ell]{N}}{2}\biggr)\frac{\sqrt[\ell]{N}}{2}+\frac{3}{2}\int_{1}^{ \sqrt[\ell]{N}/2 }S_1^*(t)dt.
\end{align}
1項目を比較すると、1以上の実数$x$について
\begin{align}
S_\ell^*\biggl(\frac{N}{2^\ell}\biggr)&\approx S_1^*\biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\biggr)\\
\to S_\ell^*(x)&\approx S_1^*\bigl(\sqrt[\ell]{x}\bigr)
\end{align}
が得られる。
2項目を比較するため、$f_\ell(N)$の積分変数を$u=t^{1/\ell}$と変数変換する。$du=t^{1/\ell-1}dt/\ell$であり、$u$は$1$から$\sqrt[\ell]{N/2^\ell}$まで動くので、
\begin{align}
\frac{1}{\ell}\int_{1}^{ N/2^\ell }S_\ell^*(t)t^{1/\ell-1}dt=\int_1^\sqrt[\ell]{N/2^\ell}S_\ell^*\bigl(u^\ell\bigr)du
\end{align}
となる。被積分関数を比較して
\begin{align}
S_\ell^*\bigl(t^\ell\bigr)&\approx S_1^*(t) \\
\to S_\ell^*(t)&\approx S_1^*\bigl(\sqrt[\ell]{t}\bigr)
\end{align}
が得られる。すなわち$f_\ell(N) \approx f_1\bigl(\sqrt[\ell]{N}\bigr)$の仮定から
\begin{align}
S_\ell^*(x)\approx S_1^*\bigl(\sqrt[\ell]{x}\bigr)
\end{align}
が導かれる。これは$\tilde{\pi}_{1/2}(N)$だけでなく$\tilde{\pi}_0(N)$や$\tilde{\pi}_1(N)$として計算しても同様である。
$S_\ell(x)$の表式に戻すと
\begin{align}
S_\ell(x)=\frac{\mu(\ell)}{\ell}S_\ell^*(x)\approx\frac{\mu(\ell)}{\ell}S_1\bigl(\sqrt[\ell]{x}\bigr)
\end{align}
となり、目的の近似式が得られる。ここで、$S_1^*\bigl(\sqrt[\ell]{x}\bigr)=S_1\bigl(\sqrt[\ell]{x}\bigr)$を用いた。
$S_\ell(x)$を$S_1(x)$で近似できました。以降の議論では、$S_1(x)$のみを考えればよいことが分かります。
$S_1(x)$について調べるため、$N$以下の素数の逆数和$q(N)$と自然数の$\ell$乗の逆数和$H_\ell(N)$の恒等式の反転式を考えます。その準備として$H_\ell(N)$の近似式を求めます。
十分大きい$N$に対し、自然数の$\ell$乗の逆数和$H_\ell(N)$は
\begin{align}
H(N)&\approx \ln N + \gamma, \\
H_{\ell\ge2}(N)&\approx\zeta(\ell)-\frac{1}{(\ell-1)N^{\ell-1}}
\end{align}
と近似できる。ここで$\gamma$はオイラー・マスケローニ定数で、$\zeta(\ell)$はリーマンゼータ関数である。
$\ell=1$のときは、Harmonicより$N$番目の調和数は漸近的に
\begin{align}
H(N)\approx\ln N+\gamma+\frac{1}{2N}+\cdots
\end{align}
で与えられる。よって十分大きい$N$に対し、
\begin{align}
H(N)\approx\ln N+\gamma
\end{align}
と近似できる。
$\ell \ge2$のときは、定義から$H_\ell(N)$を次のように書き直せる。
\begin{align}
H_\ell(N)&=\sum_{n=1}^N\frac{1}{n^\ell}\\
&=\sum_{n=1}^\infty\frac{1}{n^\ell}-\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell} \\
&=\zeta(\ell)-\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}.
\end{align}
2項目の評価を行う。各$n$について
\begin{align}
\int_{n}^{n+1}\frac{dx}{x^\ell}\le\frac{1}{n^\ell}\le \int_{n-1}^n\frac{dx}{x^\ell}
\end{align}
の不等式が成り立つ。$\lfloor N \rfloor+1$から$\infty$まで足し上げて
\begin{align}
\int_{\lfloor N \rfloor+1}^\infty\frac{dx}{x^\ell}\le\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}\le\int_{\lfloor N \rfloor}^\infty \frac{dx}{x^\ell}
\end{align}
となり、計算して
\begin{align}
\frac{1}{(\ell-1)(\lfloor N \rfloor+1)^{\ell-1}}\le\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}\le\frac{1}{(\ell-1)\lfloor N \rfloor^{\ell-1}}
\end{align}
を得る。両辺に$(\ell-1)N^{\ell-1}$を掛けると
\begin{align}
\biggl(\frac{N}{\lfloor N \rfloor+1}\biggr)^{\ell-1}\le(\ell-1)N^{\ell-1}\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}\le\biggl(\frac{N}{\lfloor N \rfloor}\biggr)^{\ell-1}
\end{align}
となる。$N\to \infty$で両辺とも1に収束するため、はさみうちの原理から
\begin{align}
\lim_{N \to \infty}(\ell-1)N^{\ell-1}\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}=1
\end{align}
が成り立ち、十分大きい$N$に対して
\begin{align}
\sum_{n=\lfloor N \rfloor+1}^\infty \frac{1}{n^\ell}\approx\frac{1}{(\ell-1)N^{\ell-1}}
\end{align}
と近似できる。
これを$H_\ell(N)$の2項目に代入し、十分大きい$N$での近似式
\begin{align}
H_\ell(N)\approx\zeta(\ell)-\frac{1}{(\ell-1)N^{\ell-1}}
\end{align}
を得る。
この結果と近似3を用いて、素数の逆数和$q(N)$についての近似式を得ることができます。 第4回 の一般化した恒等式の反転式から始めます。
素数の逆数和$q(N)$と自然数の$\ell$乗の逆数和$H_\ell(N)$の恒等式の反転式が、$H_\ell(N)$を近似式で置き換えても近似的に成り立つと仮定する(近似3)。
この仮定の下で、$q(N)$は次のように近似できる。
\begin{align}
q(N)&\approx (\gamma-1+\ln2)S_1\biggl(\frac{N}{2}\biggr)+\int_1^{N/2}\frac{S_1(t)}{t}dt \\
&+\sum_{\ell=2}^{\lfloor \log_2{N}\rfloor}\Biggl[(\zeta(\ell)-1)\Biggl(S_\ell\biggl(\frac{N}{2^\ell}\biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}+\frac{\ell-1}{\ell}\int_1^{N/2^\ell}\frac{S_\ell(t)}{t\cdot t^{\frac{\ell-1}{\ell}}}dt\Biggr)-\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}S_\ell\biggl(\frac{N}{2^\ell}\biggr)\Biggr].
\end{align}
prin2の命題7で示した$s$乗に一般化された恒等式の反転式の$s=1$の場合を考える。2以上の$N$に対し、
\begin{align}
q(N)&=\sum_{\ell=1}^{\lfloor \log_2N\rfloor}\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m}\Bigl(H_\ell\Bigl(\sqrt[\ell]{N/m}\Bigr)-1\Bigr) \\
&=\sum_{m=1}^{\lfloor N/2 \rfloor}\frac{\mu_{1m}}{m}(H(N/m)-1)+\sum_{\ell =2}^{\lfloor \log_2 N\rfloor}\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m}\Bigl(H_\ell\Bigl(\sqrt[\ell]{N/m}\Bigr)-1\Bigr)
\end{align}
が成り立ち、後の都合のため$\ell=1$と$\ell\ge2$を分けた。
補題2の$H(N)\approx\ln N+\gamma$,$H_{\ell\ge2}(N)\approx\zeta(\ell)-\frac{1}{(\ell-1)N^{\ell-1}}$で置き換えても、この反転式が近似的に成り立つと仮定する(近似3)。
$\ell=1,\ell\ge2$の場合に分け、具体的に計算する。
・$\ell=1$のとき
\begin{align}
&\sum_{m=1}^{\lfloor N/2 \rfloor}\frac{\mu_{1m}}{m}(\ln N-\ln m+\gamma-1) \\
&=(\ln N+\gamma-1)\sum_{m=1}^{\lfloor N/2\rfloor}\frac{\mu_{1m}}{m}-\sum_{m=1}^{\lfloor N/2 \rfloor}\frac{\mu_{1m}}{m}\ln m.
\end{align}
1項目の$\sum_{m=1}^{\lfloor N/2 \rfloor}\mu_{1m}/m$は定義から$S_1(N/2)$である。
2項目には
第4回
と
第5回
で登場したアーベルの総和公式を用いる。
Apostolのp.77,p.78より任意の数論的関数$a(n)$に対し、
\begin{equation} A(x)=\sum_{n\le x}a(n) \end{equation}
とする。ここで$x<1$では$A(x)=0$である。
$0< y< x$のとき、$f$が$[y,x]$の範囲で連続な導関数を持つと仮定して、
\begin{equation} \sum_{y< n\le x}a(n)f(n)=A(x)f(x)-A(y)f(y)-\int_y^xA(t)f'(t)dt \end{equation}
が成り立つ。
初項を分離してから適用することで、
\begin{align}
\sum_{m=1}^{\lfloor N/2 \rfloor}\frac{\mu_{1m}}{m}\ln m &=\frac{\mu_{11}}{1}\ln 1+\sum_{1< m\le N/2}\frac{\mu_{1m}}{m}\ln m\\
&=S_1\biggl(\frac{N}{2}\biggr)\ln\biggl(\frac{N}{2}\biggr)-\int_1^{N/2}\frac{S_1(t)}{t}dt
\end{align}
と書ける。よって$\ell=1$の項は
\begin{align}
&(\ln N+\gamma-1)S_1\biggl(\frac{N}{2}\biggr)-\biggl(S_1\biggl(\frac{N}{2}\biggr)\ln\biggl(\frac{N}{2}\biggr)-\int_1^{N/2}\frac{S_1(t)}{t}dt\biggr) \\
&=(\gamma-1+\ln2)S_1\biggl(\frac{N}{2}\biggr)+\int_1^{N/2}\frac{S_1(t)}{t}dt
\end{align}
となる。
・$\ell \ge2$のとき
\begin{align}
&\sum_{m=1}^{\lfloor N/2^\ell\rfloor} \frac{\mu_{\ell m}}{m}\Bigl(H_\ell\Bigl(\sqrt[\ell]{N/m}\Bigr)-1\Bigr) \\
&=\sum_{m=1}^{\lfloor N/2^\ell\rfloor} \frac{\mu_{\ell m}}{m}\Biggl(\zeta(\ell)-\frac{1}{(\ell-1)(N/m)^{\frac{\ell-1}{\ell}}}-1\Biggr) \\
&=(\zeta(\ell)-1)\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m}-\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m}m^{\frac{\ell-1}{\ell}}.
\end{align}
1項目の総和は$S_\ell(x)\coloneqq \sum_{m=1}^{\lfloor x \rfloor}\mu_{\ell m}/\sqrt[\ell]{m}$とアーベルの総和公式を用いて
\begin{align}
\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m^{1/\ell}}\frac{1}{m^{\frac{\ell-1}{\ell}}}&=\frac{\mu_{11}}{1}\frac{1}{1}+\sum_{1< m\le N/2^\ell}\frac{\mu_{\ell m}}{m^{1/\ell}}m^{-\frac{\ell-1}{\ell}}\\
&=S_\ell\biggl(\frac{N}{2^\ell}\biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}+\frac{\ell-1}{\ell}\int_1^{N/2^\ell}S_\ell(t)t^{-\frac{2\ell-1}{\ell}}dt
\end{align}
と書ける。ここで、$f(m)=m^{-(\ell-1)/\ell}、f'(m)=-\frac{\ell-1}{\ell}m^{-(2\ell-1)/\ell}$を用いた。
2項目の総和は
\begin{align}
\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m^{1/\ell}}\frac{m^{\frac{\ell-1}{\ell}}}{m^{\frac{\ell-1}{\ell}}}=\sum_{m=1}^{\lfloor N/2^\ell \rfloor}\frac{\mu_{\ell m}}{m^{1/\ell}}=S_\ell\biggl(\frac{N}{2^\ell}\biggr)
\end{align}
と簡単になる。よって$\ell\ge2$の項は
\begin{align}
(\zeta(\ell)-1)\Biggl(S_\ell\biggl(\frac{N}{2^\ell}\biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}+\frac{\ell-1}{\ell}\int_1^{N/2^\ell}\frac{S_\ell(t)}{t\cdot t^\frac{\ell-1}{\ell}}dt\Biggr)-\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}S_\ell\biggl(\frac{N}{2^\ell}\biggr)
\end{align}
となる。$\ell=1$の式と合わせて目的の式が得られる。
$q(N)$と$H_\ell(N)$の恒等式の反転式は、$H_\ell(N)$を近似式で置き換えても、近似的に成り立つと仮定した(近似3)。ここで$m=\lfloor N/2 ^\ell \rfloor$の項では
\begin{align}
H_\ell\Bigl(\sqrt[\ell]{N/\lfloor N/2^\ell \rfloor}\Bigr)=H_\ell(2)
\end{align}
となる。しかし、補題2の$H_\ell(N)$の近似式は、十分大きい$N$に対する近似である。よって$m=\lfloor N/2^\ell \rfloor$の項では、近似式で置き換えることは適切でない。$N$を大きくしても、そのような項は必ず存在し、解決しない。
素数の逆数和$q(N)$についての近似式が得られました。
次に$S_1(x)$の近似式を求めます。そのための準備として、対数積分$\mathrm{li}(x)$の近似式を導きます。
$x\to 0^+$で対数積分$\mathrm{li}(x)$は
\begin{align}
\mathrm{li}(x)\approx\frac{x}{\ln x}, \quad x\to0^+,
\end{align}
と近似できる。
DLMFの$\S$6.2の式6.2.1より指数積分$E_1(z)$を導入する。
\begin{align}
E_1(z)\coloneqq \int_z^\infty \frac{e^{-t}}{t}dt, \quad z\neq0.
\end{align}
この関数について漸近展開が知られており、DLMFの$\S$6.12の式6.12.1
\begin{align}
E_1(z)\sim\frac{e^{-z}}{z}\biggl(1-\frac{1!}{z}+\frac{2!}{z^2}-\frac{3!}{z^3}+\cdots\biggr), \quad z\to \infty,\:|\mathrm{ph}\:z|\le \frac{3}{2}\pi-\delta\biggl(<\frac{3}{2}\pi\biggr)
\end{align}
である。これを利用するため、対数積分$\mathrm{li}(x)$を指数積分$E_1(x)$で表す。
\begin{align}
\mathrm{li}(x)&=\int_0^x\frac{dt}{\ln t}\\
&=-\int_{-\ln x}^\infty\frac{e^{-u}}{u}du \\
&=-E_1(-\ln x).
\end{align}
1行目から2行目の式変形で$-\ln t=u$と変数変換した。$dt=-e^{-u}du$の関係があり、その負符号で積分区間を入れ替えることにより、$u$は$-\ln x$から$+\infty$まで動く。
$x\to0^+$で$-\ln x\to+\infty$となるので、さきほどの漸近展開が使えて
\begin{align}
\mathrm{li}(x)&=-E_1(-\ln x) \\
&\sim-\frac{e^{\ln x}}{-\ln x}\biggl(1-\frac{1!}{-\ln x}+\frac{2!}{(-\ln x)^2}-\frac{3!}{(-\ln x)^3}+\cdots\biggr)\\
&=\frac{x}{\ln x}\biggl(1+\frac{1!}{\ln x}+\frac{2!}{(\ln x)^2}+\frac{3!}{(\ln x)^3}+\cdots\biggr)\\
&\approx\frac{x}{\ln x}
\end{align}
を得る。
メルテンスの定理と命題3の$q(N)$の近似式から$S_1(x)$の近似式を導出します。
1以上の実数$x$に対して、$S_1(x)$は
\begin{align}
S_1(x)\approx\frac{1}{\ln 2x}
\end{align}
と近似できる。
まず発見論的に$S_1(x)$の近似式を推測する。
Mertensのメルテンスの第二定理より、素数の逆数和$q(N)$は
\begin{align}
q(N)\approx\ln\ln N+b
\end{align}
と書くことができ、$b$はMeissel-Mertens定数である。
十分大きい$N$について、$q(N)$は積分を用いて
\begin{align}
q(N)\approx \int_2^N \frac{dx}{x\ln x} +\ln \ln2 +b
\end{align}
と表せる。一方で、命題3で示した$q(N)$の近似式の$S_1(t)$を含む積分に注目すると
\begin{align}
\int_{1}^{N/2}\frac{S_1(t)}{t}dt&=\int_{2}^N \frac{S_1(u/2)}{u/2}\frac{du}{2}\\
&=\int_{2}^{N}\frac{S_1(u/2)}{u}du
\end{align}
となる。ここで$u=2t$と変数変換した。このとき$du=2dt$の関係があり、$u$は2から$N$まで動く。比較することで
\begin{align}
S_1(x) \approx \frac{1}{\ln 2x}
\end{align}
と推測できる。
この近似式が、命題3の$q(N)$の近似式の解となるかを確かめる。比較した項は
\begin{align}
\int_2^N \frac{S_1(u/2)}{u}du \approx\int_2^N \frac{1}{u\ln u}du=\ln\ln N-\ln\ln2
\end{align}
となる。それ以外の$\ell=1$の残りの項と、$\ell\ge2$の場合に分けて計算を進める。
・$\ell=1$の残りの項
\begin{align}
(\gamma-1+\ln2)S_1\biggl(\frac{N}{2}\biggr)\approx(\gamma-1+\ln2)\frac{1}{\ln N}=O\biggl(\frac{1}{\ln N}\biggr).
\end{align}
ここでHardyのp.5より$O(\phi)$はランダウの記号であり、ある関数$f(x)$に対して
\begin{align}
|f(x)|< K\phi(x)
\end{align}
を満たすものである。$K$は正の定数で、$\phi$は$x$の正の関数とする。
・$\ell\ge2$の場合
$q(N)$の近似式に命題1の$S_\ell(x)\approx\mu(\ell)/\ell ・S_1\bigl(\sqrt[\ell]{x}\bigr)$の近似を用いる。
\begin{align}
\sum_{\ell=2}^{\lfloor \log_2{N}\rfloor}\frac{\mu(\ell)}{\ell}\Biggl[(\zeta(\ell)-1)\Biggl(S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}+\frac{\ell-1}{\ell}\int_1^{N/2^\ell}\frac{S_1\bigl(\sqrt[\ell]{t}\bigr)}{t\cdot t^{\frac{\ell-1}{\ell}}}dt\Biggr)-\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)\Biggr].
\end{align}
2項目は$u=2\sqrt[\ell]{t}$と変数変換する。$t^{\frac{1-\ell}{\ell}}dt/\ell=du/2$の関係があり、$u$は2から$\sqrt[\ell]{N}$まで動く。すなわち
\begin{align}
\sum_{\ell=2}^{\lfloor \log_2{N}\rfloor}\frac{\mu(\ell)}{\ell}\Biggl[(\zeta(\ell)-1)\Biggl(S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}+(\ell-1)\int_2^{\sqrt[\ell]{N}}\frac{S_1(u/2)}{(u/2)^\ell}\frac{du}{2}\Biggr)-\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)\Biggr]
\end{align}
となる。次に、各項に対し$S_1(x)\approx 1/\ln2x$の近似をし、式を整理する。
・1項目
\begin{align}
\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}&\approx
\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)\frac{1}{\ln \sqrt[\ell]{N}}\biggl(\frac{N}{2^\ell}\biggr)^{-\frac{\ell-1}{\ell}}\\
&=\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)\frac{\ell\cdot2^{\ell-1}}{N^{\frac{\ell-1}{\ell}}\ln N}
\end{align}
・2項目
\begin{align}
&\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\int_{2}^{\sqrt[\ell]{N}}\frac{S_1(u/2)}{u^\ell}du\\
&\approx \frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\int_{2}^{\sqrt[\ell]{N}} \frac{1}{u^\ell \ln u}du.
\end{align}
ここで積分変数を$s=u^{1-\ell}$と変数変換する。$ds=(1-\ell)u^{-\ell}du$の関係があり、$s$は$2^{1-\ell}$から$N^\frac{1-\ell}{\ell}$まで動く。すると
\begin{align}
&\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\int_{2^{1-\ell}}^{N^{\frac{1-\ell}{\ell}}} \frac{1}{\ln s^{1/(1-\ell)}}\frac{ds}{1-\ell} \\
&=\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\int_{2^{1-\ell}}^{N^{\frac{1-\ell}{\ell}}} \frac{1}{\ln s}ds \\
&=\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\bigl(\mathrm{li}\bigl(N^{\frac{1-\ell}{\ell}}\bigr)-\mathrm{li}\bigl(2^{1-\ell}\bigr)\bigr)
\end{align}
のように対数積分$\mathrm{li}(x)$を用いて表せる。
補題4より$x\to0^+$で$\mathrm{li}(x)\approx x/\ln x$と近似できるので
\begin{align}
&\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\biggl(\frac{1}{N^{\frac{\ell-1}{\ell}} \ln N^{\frac{1-\ell}{\ell}}} -\mathrm{li}\bigl(2^{1-\ell}\bigr)\biggr) \\
&=-\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)\frac{\ell \cdot 2^{\ell-1}}{N^{\frac{\ell-1}{\ell}} \ln N}
-\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\mathrm{li}\bigl(2^{1-\ell}\bigr)
\end{align}
となる。
・3項目
\begin{align}
-\frac{\mu(\ell)}{\ell}\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}S_1\Biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\Biggr)&\approx-\frac{\mu(\ell)}{\ell}\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}}\frac{1}{\ln\sqrt[\ell]{N}} \\
&=-\frac{\mu(\ell)}{\ell-1}\frac{1}{N^{\frac{\ell-1}{\ell}}\ln N}.
\end{align}
以上の結果をまとめて、$q(N)$の近似式は
\begin{align}
&\ln\ln N-\ln\ln2+O\biggl(\frac{1}{\ln N}\biggr)+\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\Biggl[\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)\Biggl(\frac{\ell\cdot2^{\ell-1}}{N^{\frac{\ell-1}{\ell}}\ln N}-\frac{\ell\cdot2^{\ell-1}}{N^{\frac{\ell-1}{\ell}}\ln N}-(\ell-1)2^{\ell-1}\mathrm{li}\bigl(2^{1-\ell}\bigr)\Biggr)-\frac{\mu(\ell)}{\ell-1}\frac{1}{N^{\frac{\ell-1}{\ell}}\ln N} \Biggr]\\
&=\ln\ln N-\ln\ln2+O\biggl(\frac{1}{\ln N}\biggr)+\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\Biggl[-\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\mathrm{li}\bigl(2^{1-\ell}\bigr)-\frac{\mu(\ell)}{\ell-1}\frac{1}{N^{\frac{\ell-1}{\ell}}\ln N} \Biggr]
\end{align}
と書ける。
最後に$\ell$の総和について評価を行う。メビウス関数の大きさ$|\mu(\ell)|\le1$と$\zeta(\ell)-1\approx 2^{-\ell}$、補題4より$\mathrm{li}\bigl(2^{1-\ell}\bigr)\approx 2^{1-\ell}/\ln2^{1-\ell}$が成り立つ。これらを用いて、1項目の絶対値は
\begin{align}
&\Biggl|\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}-\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\mathrm{li}\bigl(2^{1-\ell}\bigr) \Bigg|\\
&\approx \sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{1}{\ell}2^{-\ell}(\ell-1)2^{\ell-1}\frac{2^{1-\ell}}{|1-\ell|\ln2}\\
&=\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{1}{\ell\cdot 2^\ell \ln2}
\end{align}
と整理できる。これは
\begin{align}
\sum_{\ell=2}^{\infty}\frac{1}{2^\ell}=\frac{1}{2}
\end{align}
で上から抑えられるため、1項目は$O(1)$と評価できる。
2項目の絶対値は
\begin{align}
\Biggl| \sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell-1}\frac{1}{N^{\frac{\ell-1}{\ell}}\ln N}\Biggr|<\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}\ln N}
\end{align}
で抑えられる。さらに
\begin{align}
\frac{1}{N^{\frac{\ell-1}{\ell}} \ln N} \le \frac{1}{\sqrt{N} \ln N}
\end{align}
を用いて
\begin{align}
\sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{1}{(\ell-1)N^{\frac{\ell-1}{\ell}}\ln N}
&< \sum_{\ell=2}^{\lfloor \log_2 N \rfloor}\frac{1}{\sqrt{N}\ln N}\\
&\approx \frac{\ln N}{\ln2}\frac{1}{\sqrt{N}\ln N}\\
&=O\biggl(\frac{1}{\sqrt{N}}\biggr)
\end{align}
と上から抑えられる。ここで1行目から2行目の式変形では$\lfloor \log_2 N\rfloor \approx \log_2 N =\ln N/\ln2$を用いた。
最終的に$q(N)$の近似式は
\begin{align}
&\ln\ln N-\ln\ln2+O\biggl(\frac{1}{\ln N}\biggr)+O(1)+O\biggl(\frac{1}{\sqrt{N}}\biggr) \\
&=\ln\ln N+O(1)
\end{align}
となり、主項がメルテンスの第二定理と一致する。
したがって、$S_1(x)\approx1/\ln 2x$は命題3で示した$q(N)$の近似式の近似解である。
命題5の証明において、$q(N)$の近似式に$S_1(x)\approx 1/\ln 2x$を用いた。定数項は具体的に
\begin{align}
-\ln\ln2-\sum_{\ell=2}^{\lfloor \log_2{N} \rfloor}\frac{\mu(\ell)}{\ell}(\zeta(\ell)-1)(\ell-1)2^{\ell-1}\mathrm{li}\bigl(2^{1-\ell}\bigr)
\end{align}
と書ける。$\lfloor\log_2{N}\rfloor=100$まで数値計算した結果は$ 0.05210$であった。一方で、Meissel-Mertens定数$b$は$0.261497$なので定数項は一致しない。
$S_1(x)$の近似式が求まりました。 第4回 の命題5で考えていたものと一致しています。これにより素数計数関数$\pi(N)$の近似式を求められます。
それでは素数計数関数$\pi(N)$の近似式を導出します。
これまでの4つの近似の下で
\begin{align}
\pi(N)\approx\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)
\end{align}
が成り立つ。$\tilde{\pi}(N),\delta\tilde{\pi}(N)$はそれぞれ
\begin{align}
\tilde{\pi}(N)&= \sum_{\ell=1}^{\lfloor \log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{3}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr), \\
\pm\delta\tilde{\pi}(N)&=\pm\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\biggl(-\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{1}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr)
\end{align}
である。ここで$\mathrm{Li}(x)$は補正対数積分関数であり、
\begin{align}
\mathrm{Li}(x)&\coloneqq\mathrm{li}(x)-\mathrm{li}(2), \\
\mathrm{li}(x)&\coloneqq\int_{0}^{x}\frac{dt}{\ln t}
\end{align}
で定義される。
prin2の定理3の結果より$\tilde{\pi}_0(N),\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$はそれぞれ
\begin{align}
\tilde{\pi}_0(N)&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\biggl(S_\ell\biggl( \frac{N}{2^\ell}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{1}{\ell}\int_{1}^{ N/2^\ell }S_\ell(t)t^{1/\ell-1}dt\biggr), \\
\tilde{\pi}_1(N)&=\sum_{\ell=1}^{\lfloor \log_2{N} \rfloor}\frac{2}{\ell}\int_{1}^{ N/2^\ell }S_\ell(t)t^{1/\ell-1}dt, \\
\tilde{\pi}_{1/2}(N)&=\frac{\pi_0(N)+\pi_1(N)}{2}
\end{align}
と書ける。 具体例として$\tilde{\pi}_0(N)$について考える。
まず命題1の$S_\ell(x)\approx\mu(\ell)/\ell ・S_1\bigl(\sqrt[\ell]{x}\bigr)$の近似を用いて
\begin{align}
\tilde{\pi}_0(N)\approx\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\Biggl(S_1\biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\biggr)\sqrt[\ell]{\frac{N}{2^\ell}}+\frac{1}{\ell}\int_1^{N/2^\ell}S_1\bigl(\sqrt[\ell]{t}\bigr)t^{1/\ell-1}dt\Biggr)
\end{align}
となる。積分変数を$u=2\sqrt[\ell]{t}$と変数変換する。$t^{1/\ell-1}dt/\ell=du/2$の関係があり、$u$は2から$\sqrt[\ell]{N}$まで動く。つまり
\begin{align}
\tilde{\pi}_0(N)\approx\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\Biggl(S_1\biggl(\sqrt[\ell]{\frac{N}{2^\ell}}\biggr)\frac{\sqrt[\ell]{N}}{2}+\int_2^{\sqrt[\ell]{N}}S_1\Bigl(\frac{u}{2}\Bigr)\frac{du}{2} \Biggr)
\end{align}
と書ける。
次に命題5の近似式$S_1(x)\approx1/\ln 2x$を用いる。
\begin{align}
\tilde{\pi}_0(N)&\approx \sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{\ln \sqrt[\ell]{N}}\frac{\sqrt[\ell]{N}}{2}+\frac{1}{2}\int_2^\sqrt[\ell]{N} \frac{du}{\ln u}\biggr)\\
&=\sum_{\ell=1}^{\lfloor \log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{2}\frac{\sqrt[\ell]{N}}{\ln\sqrt[\ell]{N}}+\frac{1}{2}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr).
\end{align}
2行目の式変形では、補正対数積分$\mathrm{Li}(x)$を用いた。同様にして$\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$は
\begin{align}
\tilde{\pi}_1(N) &\approx\sum_{\ell=1}^{\lfloor \log_2 N \rfloor}\frac{\mu(\ell)}{\ell}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr), \\
\tilde{\pi}_{1/2}(N)&=\frac{\tilde{\pi}_0(N)+\tilde{\pi}_1(N)}{2}\approx\sum_{\ell=1}^{\lfloor\log_2 N\rfloor}\frac{\mu(\ell)}{\ell}\biggl(\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln \sqrt[\ell]{N}}+\frac{3}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr)
\end{align}
と近似できる。
図1では$\tilde{\pi}_{1/2}(N)$が$\pi(N)$をよく近似していたため、素数計数関数の近似式の基準とする。
\begin{align}
\tilde{\pi}(N)\coloneqq\tilde{\pi}_{1/2}(N).
\end{align}
そして$\tilde{\pi}(N)$と$\tilde{\pi}_1(N),\tilde{\pi}_0(N)$の差$\pm\delta\tilde{\pi}(N)$は
\begin{align}
\delta\tilde{\pi}(N)&\coloneqq\tilde{\pi}_1(N)-\tilde{\pi}(N) \\
&=\tilde{\pi}(N)-\tilde{\pi}_0(N)\\
&=\sum_{\ell=1}^{\lfloor \log_2 N \rfloor} \frac{\mu(\ell)}{\ell}\biggl(-\frac{1}{4}\frac{\sqrt[\ell]{N}}{\ln\sqrt[\ell]{N}}+\frac{1}{4}\mathrm{Li}\bigl(\sqrt[\ell]{N}\bigr)\biggr)
\end{align}
と書ける。以上より目的の近似式が得られた。
素数計数関数の恒等式の理論から$\pi(N)$の近似式が導けました。この近似式について、いくつかコメントしておきます。
まず素数定理との関係についてです。
第4回
の定理6として挙げた素数定理を再掲します。
素数定理
任意の関数$f(x), g(x)(\ne0)$に対して、記号「$\sim$」を
\begin{equation}
f(x) \sim g(x) \Leftrightarrow \lim_{x\to \infty}\frac{f(x)}{g(x)}=1
\end{equation}
と定義する。そのとき、primethmathやApostolのp.74より
\begin{align}
\pi(N)&\sim \frac{N}{\ln N}
\end{align}
が成り立つ。
$\mathrm{Li}(N)\sim N/\ln N$であるため、今回導出した近似式は素数定理と矛盾していません。
第4回
では素数定理から$S_1(x)\approx1/\ln 2x$の近似を仮定しましたが、今回は素数の逆数和$q(N)$の近似式から導いています。ただし、厳密に素数定理を証明するには、各近似についての正確な誤差評価が必要です。
それを踏まえ、導出で用いた近似をまとめます。
○近似1:素数計数関数の反転式で$\lfloor N \rfloor \to N,N-1,N-1/2$とする近似
結果:素数計数関数$\pi(N)$の近似式を、$S_\ell(x)\coloneqq \sum_{m=1}^{\lfloor x \rfloor}\mu_{\ell m}/\sqrt[\ell]{m}$を用いて表すことができます。
○近似2:1以上の実数$x$に対し$S_\ell(x)\approx \mu(\ell)/\ell \cdot S_1\bigl(\sqrt[\ell]{x}\bigr)$とする近似
結果:リーマンの近似公式の形から導かれ、$S_1(x)$のみを考えれば十分となります。
○近似3:素数の逆数和$q(N)$と自然数の$\ell$乗の逆数和$H_\ell(N)$の恒等式の反転式が、$H_\ell(N)$を近似式で置き換えても近似的に成り立つとする近似
結果:$q(N)$の近似式が得られます。
○近似4:1以上の実数$x$に対し$S_1(x)\approx1/\ln 2x$とする近似
結果:近似3による$q(N)$の近似式とメルテンスの定理から導かれるもので、これにより$\pi(N)$の近似式が得られます。
今回はこの4つの近似を用いて、素数計数関数$\pi(N)$の近似式を導出しました。近似2はリーマンの近似公式、近似4は図2から妥当と考えられますが、近似1と近似3は解析的に解くために用いたので検証が必要です。
これまでの「素数計数関数の恒等式」シリーズをおさらいし、特に
第4回
の結果から議論を始めました。本記事の近似式$\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$では、図3,4から$\tilde{\pi}_{1}(N)$が$\pi(N)$をよく近似していました。また近似式の精度は$N/\ln N$とリーマンの素数計数関数$R(N)$の中間に位置することが分かりました。
理論的には、素数計数関数の恒等式から始めて、近似を繰り返し近似式$\tilde{\pi}(N)\pm\delta\tilde{\pi}(N)$を導出しました。素数の逆数和$q(N)$と自然数の$\ell$乗の逆数和$H_\ell(N)$の関係式から$q(N)$の近似式を得て、$S_1(x)$の近似解を求めました。素数定理と無矛盾ですが、この方向で素数定理を厳密に証明するためには、各近似の誤差評価を要します。
本記事ではリーマンらの複素解析を用いず、素数計数関数の恒等式と解析的整数論の初等的な手法(メビウス反転、アーベルの総和公式)を用いました。今後の課題としては、近似の誤差評価、リーマンらの複素解析の結果との比較、リーマン予想との繋がりが挙げられます。
ひとまずは「素数計数関数の恒等式」シリーズの区切りとします。ここまでお付き合いいただき、ありがとうございました。今後も興味深い結果が出たら記事を書く予定です。