0
高校数学解説
文献あり

素数計数関数の恒等式と素数計数関数

0
0
$$$$

はじめに

 本記事は「 素数計数関数の恒等式 」シリーズの第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}
と定義しています。

近似1の近似式の表記について

ここで定義した$\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種の近似式の比較 素数計数関数と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として再掲します。
!FORMULA[30][1731313890][0]の近似の検証 $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に示します。
素数計数関数!FORMULA[45][1667325170][0]との比較(!FORMULA[46][-118687712][0]まで) 素数計数関数$\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)$をよく近似していたため、相違が見られます。
素数計数関数!FORMULA[54][1667325170][0]との比較(!FORMULA[55][615581362][0]まで) 素数計数関数$\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$43.8220.596-4.439$\pm$14.904
$10^2$2523.8461.796-4.617$\pm$7.184
$10^3$1681618-4.276$\pm$4.481
$10^4$1,2291,18839-3.325$\pm$3.157×
$10^5$9,5929,354233-2.481$\pm$2.433×
$10^6$78,49876,9711,556-1.945$\pm$1.983
$10^7$664,579653,55311,114-1.659$\pm$1.672
$10^8$5,761,4555,678,19283,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)$と等しいため、それらの中間に位置づけられます。

$S_{\ell}(x)$の近似式

 ここからは$\pi(N)$の近似式を理論的に導出します。 第4回 では$\tilde{\pi}_0(N),\tilde{\pi}_1(N),\tilde{\pi}_{1/2}(N)$$S_\ell(x)$で表されていたので、その近似式から求めていきます。

$\ell\ge2$の場合

 $\ell\ge2$の場合はリーマンの近似公式を参考にして、$S_\ell(x)$$S_1(x)$で近似します。

$S_\ell(x)$の近似式 (近似2)

リーマンの近似公式から導かれる仮定の下、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)$のみを考えればよいことが分かります。

$\ell=1$の場合

 $S_1(x)$について調べるため、$N$以下の素数の逆数和$q(N)$と自然数の$\ell$乗の逆数和$H_\ell(N)$の恒等式の反転式を考えます。その準備として$H_\ell(N)$の近似式を求めます。

自然数の$\ell$乗の逆数和$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)$の近似式

素数の逆数和$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$の式と合わせて目的の式が得られる。

近似3の問題点

$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)$の近似式を導出します。

$S_1(x)$の近似式(近似4)

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の証明の定数項について

命題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)$の近似式を求められます。

近似式$\tilde{\pi}(N)\pm\delta\tilde{\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} 
と定義する。そのとき、primethmathApostolの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)$の近似解を求めました。素数定理と無矛盾ですが、この方向で素数定理を厳密に証明するためには、各近似の誤差評価を要します。
 本記事ではリーマンらの複素解析を用いず、素数計数関数の恒等式と解析的整数論の初等的な手法(メビウス反転アーベルの総和公式)を用いました。今後の課題としては、近似の誤差評価、リーマンらの複素解析の結果との比較、リーマン予想との繋がりが挙げられます。
 ひとまずは「素数計数関数の恒等式」シリーズの区切りとします。ここまでお付き合いいただき、ありがとうございました。今後も興味深い結果が出たら記事を書く予定です。

参考文献

投稿日:1時間前
数学の力で現場を変える アルゴリズムエンジニア募集 - Mathlog served by OptHub

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

LAJQ
10
597

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中