実数$x$に対して, $\mathbb{P}$,$\mathbb{P}_x$をそれぞれ素数全体, $x$未満の素数全体とする. まず, 素数の逆数和について有限和で考えてみよう.
\begin{align}
\frac{1}{2}=0.5,\quad\frac{1}{2}+\frac{1}{3}=0.8\dot{3},\quad\frac{1}{2}+\frac{1}{3}+\frac{1}{5}=1.0\dot{3}
\end{align}
値が小さいし, 収束しそう...?そこで, pythonを使ってがばっと計算してみる(コード作りに慣れていないので見にくいかもしれませんがお許しください><).
素数
上のコードを使って, 整数$k$未満の素数の逆数和を求めると,
\begin{align}
\sum_{p\in\mathbb{P}_{100}}\frac{1}{p}=1.802\ldots,\quad \sum_{p\in\mathbb{P}_{1000}}\frac{1}{p}=2.198\ldots,\quad \sum_{p\in\mathbb{P}_{10000}}\frac{1}{p}=2.483\ldots,\quad \sum_{p\in\mathbb{P}_{100000}}\frac{1}{p}=2.705\ldots
\end{align}
となる. どれだけ$k$を大きくしようと, $10$すら超えそうにない. よって収束する, と安易に決めつけてはいけない. 実際,
\begin{align}
\sum_{p\in\mathbb{P}_x}\frac{1}{p}=\log\log x+O(1)\quad\pas{x\to\infty}
\end{align}
が成り立つのである($O(1)$は, $x\to\infty$で定数に収束する$x$の関数, と思ってくれればよい). 従って, $x>e^{(e^{10})}$ととれば$x$以下の素数の逆数和が$10$を越えるだろうと考えられる.
\begin{align}
9565<\log_{10} e^{(e^{10})}<9566
\end{align}
だから, $e^{(e^{10})}$は$9566$桁の数である. 上で求めた逆数和は$5$桁程度なので, いかに膨大な量を計算しなければ$10$に到達しないかが分かる. 話を戻すが, 理論上は素数の逆数和$\sum_{p\in\mathbb{P}}p^{-1}$は正の無限大に発散するのである. この級数は収束しそうで発散する代表的な例であり, 計算だけで判断するなという戒めにもなる.
そこで, この級数が発散することの証明を, Eulerの方法(1737?),Erdősの方法(1938),Clarksonの方法(1966), Vandeneyndenの方法(1980)の$4$パターン紹介する. 実はもう一つ簡単な示し方があるが, それは逆数和が発散することを示すのが目的ではないので, 次の機会にでも記事にしたいと思う. まず, 準備として次を示しておく.
$\sum_{n=1}^\infty n^{-\sigma}$は$\sigma>1$で収束し, $\sigma\leq1$で正の無限大に発散する.
$\sigma>1$とすると
\begin{align}
\sum_{n=1}^\infty\frac{1}{n^\sigma}\leq1+\int_1^\infty x^{-\sigma}dx=1+\frac{1}{\sigma-1}<\infty
\end{align}
より収束する. また$\sigma\leq1$の場合は, 対象が正項級数ゆえ$\sigma=1$のときのみ示せば十分である. そこで収束すると仮定すると,
\begin{align}
\sum_{n=1}^\infty\frac{1}{n}\geq\int_1^\infty x^{-1}dx=\log x\Big|_1^\infty=\infty.
\end{align}
これは矛盾である.
\begin{align} \sum_{n=1}^\infty\frac{1}{n^\sigma}=\prod_{p\in\mathbb{P}}\pas{1-p^{-\sigma}}^{-1}\quad (\sigma>1) \end{align}
以下を参照されたい.
私はゼータ関数のEuler積表示を無限積の絶対収束性を使わず示したい!
それでは, 本題に入っていく.
素数の逆数和$\sum_{p\in\mathbb{P}}p^{-1}$は正の無限大に発散する.
$\diamond$Eulerの方法・・・厳密には無限積の話や絶対収束の話が必要だが, 感覚的には命題$1,2$と対数関数のテイラー展開を知っていればすぐに理解できる証明.
$x>1$とする.
\begin{align}
\log\zeta(x)
&=\log\prod_{p\in\mathbb{P}}\frac{1}{1-p^{-x}}\\
&=-\sum_{p\in\mathbb{P}}\log\pas{1-p^{-x}}\\
&=\sum_{p\in\mathbb{P}}\sum_{m=1}^\infty\frac{1}{mp^{mx}}\quad\pas{\text{logのテイラー展開}}\\
&=\sum_{p\in\mathbb{P}}\frac{1}{p^x}+\sum_{p\in\mathbb{P}}\sum_{m=2}^\infty\frac{1}{mp^{mx}}
\end{align}
と変形できる. ここで
\begin{align}
\abs{\sum_{m=2}^\infty\frac{1}{mp^{mx}}}<\frac{1}{2}\sum_{m=2}^\infty\pas{\frac{1}{p^x}}^m=\frac{1}{2}\frac{\frac{1}{p^{2x}}}{1-\frac{1}{p^x}}\leq \frac{1}{2p^{2x}}\frac{1}{1-\frac{1}{2}}=\frac{1}{p^{2x}}
\end{align}
だから, 二項目は
\begin{align}
\abs{\sum_{p\in\mathbb{P}}\sum_{m=2}^\infty\frac{1}{mp^{mx}}}<\sum_{p\in\mathbb{P}}\frac{1}{p^{2x}}<\sum_{n=2}^\infty\frac{1}{n^{2x}}<\sum_{n=1}^\infty\frac{1}{n^2}.
\end{align}
と評価でき, 収束が確認できる. 従って, 素数の逆数和の収束発散は$\log\zeta(x)$の$x\to1+$での収束発散に置き換えられる. しかし, 命題$1$より$\zeta(x)$が$x\to1+$で正の無限大に発散するので, $\log\zeta(x)$も発散する. 以上より, 素数の逆数和は発散する.
$\diamond$Erdősの方法・・・非常にエレガントな証明. 証明を追ったすべての人を虜にしたに違いない. 自然数の性質をうまく利用している.
$n$番目の素数を$p_n$と書く. 素数の逆数和が収束すると仮定すると, ある自然数$k$が存在して
\begin{align}
\sum_{i=k+1}^\infty\frac{1}{p_i}<\frac{1}{2}
\end{align}
が成り立つ. $N$を$2^{2k+2}$以上の自然数としてとり, $N_1$を$p_k$より大きな素因子をもつ$N$以下の自然数の個数, $N_2$を$p_k$以下の素因子を持つ$N$以下の自然数の個数とする. 即ち, $\Z_{\geq0}$で非負整数全体を表わすとすれば,
\begin{align}
N_1&=\#\Pas{n\leq N : n=\prod_{i=1}^\infty p_i^{e_i}, \prod_{i=k+1}^\infty p_i^{e_i}>1\quad\pas{e_i\in\Z_{\geq0}}},\\
N_2&=\#\Pas{n\leq N : n=\prod_{i=1}^k p_i^{e_i}\quad\pas{e_i\in\Z_{\geq0}}}.
\end{align}
もちろん, $N=N_1+N_2$である. 一般に, 素数$p$の倍数である$N$以下の自然数の個数は$\floor{\frac{N}{p}}$だから, $N_1$は
\begin{align}
N_1\leq \sum_{i=k+1}^\infty\floor{\frac{N}{p_i}}\leq N\sum_{i=k+1}^\infty\frac{1}{p_i}<\frac{N}{2}
\end{align}
と評価できる. 次に, $p_1,p_2,\ldots,p_k$のみを素因子にもつ$N$以下の自然数$n$をひとつとり, $n=ab^2$のように平方因子を含まない$a\in\N$と平方数$b^2\in\N$の積で表わす. $a$の可能性は, 各$p_i\;\pas{1\leq i\leq k}$を含むか否かの$2^k$通りある. 一方,
\begin{align}
b^2=\frac{n}{a}\leq N
\end{align}
より$b\leq\sqrt{N}$だから, $b$の可能性は高々$\sqrt{N}$通りである. よって, $n$の可能性は高々$2^k\sqrt{N}$通りであり,
\begin{align}
N_2\leq2^k\sqrt{N}
\end{align}
と評価できる. 以上より,
\begin{align}
N=N_1+N_2<\frac{N}{2}+2^k\sqrt{N}\\
\end{align}
これを変形すると$\pas{\frac{N}{2}}^2<2^{2k}N$となり, $N$について解けば
\begin{align}
N<2^{2k+2}
\end{align}
が得られる. ところが, $N$は$2^{2k+2}$以上としてとっていたので矛盾が生じている. 従って, 素数の逆数和は発散する.
$\diamond$Clarksonの方法・・・素数の性質, 特に素因数分解の一意性をうまく利用した証明. 一見複雑だが, 難しい議論は必要ない. 有限和$\sum_{n=1}^r\pas{1+nQ}^{-1}$を$r$に依存しない形で上から抑えるという手法はとても面白い.
自然数$n$に対し$p_n$で$n$番目の素数を表わすとする.
収束すると仮定し矛盾を導く. もし収束するなら, ある自然数$k$が存在して
\begin{align}
\sum_{m=k+1}^\infty\frac{1}{p_m}<\frac{1}{2}\quad \quad\pas{1}
\end{align}
が成り立つ. $Q=p_1p_2\cdots p_k$とおくと,$1+nQ\;\pas{n\in\N}$は$p_1,p_2,\ldots,p_k$のどれでも割れないので,
\begin{align}
1+nQ=\prod_{m=k+1}^\infty p_m^{a_m(n)}\quad \pas{a_m(n)\geq0}
\end{align}
と因数分解できる. もちろん, $m$が十分大きければ, $a_m(n)=0$である. ここで, $\N$上の二項関係$\sim$を
\begin{align}
n\sim n^\prime\overset{\text{def}}{\iff}\sum_{m=k+1}^\infty a_m(n)=\sum_{m=k+1}^\infty a_m(n^\prime)
\end{align}
で定めれば, 素因数分解の一意性よりwell-definedである. 自然数$r$に対して, 商集合:
\begin{align}
A=\Pas{1,2,\ldots,r}/\sim
\end{align}
は$\N$の部分集合と同一視できるため, $A$に$\N$と同じ順序を入れ, $A$の元を小さい順に$x_1,x_2,\ldots,x_{r^\prime}$と書くことにする. このとき,
\begin{align}
\sum_{n=1}^r\frac{1}{1+nQ}
&=\sum_{n=1}^r\frac{1}{\prod_{m=k+1}^\infty p_m^{a_m(n)}}\\
&\leq \sum_{n=1}^{r^\prime}\pas{\sum_{m=k+1}^\infty\frac{1}{p_m}}^{\sum_{m=k+1}^\infty a_m(x_n)}\\
&\leq \sum_{l=1}^\infty\pas{\sum_{m=k+1}^\infty\frac{1}{p_m}}^l\\
&<\sum_{l=1}^\infty\pas{\frac{1}{2}}^l\quad\pas{(1)を適用した.}\\
&=1.
\end{align}
このことから,
\begin{align}
\sum_{n=1}^\infty\frac{1}{1+nQ}\leq1
\end{align}
が成り立つ. しかし,
\begin{align}
\sum_{n=1}^\infty\frac{1}{1+nQ}
&= \frac{1}{Q}\sum_{n=1}^\infty\frac{1}{n+Q^{-1}}\\
&\geq \frac{1}{Q}\sum_{n=1}^\infty\frac{1}{n+1}\\
&\geq \frac{1}{Q}\int_1^\infty\frac{dx}{x+1}\\
&=+\infty.
\end{align}
より矛盾が生じているので, 素数の逆数和は正の無限大に発散する. 但し, 二つ目の不等式にて
Euler-Maclaurinの定理
の定理$1$を用いた.
$\diamond$Vandeneyndenの方法・・・$1+p^{-1}$の無限積が発散することを用いた証明. よく知られた不等式$1+p^{-1}< e^{p^{-1}}$を使って素数の逆数和へ繫いでいる. 証明されたのが$40$年前と比較的最近なのも驚きである.
任意の素数$p$に対して
\begin{align}
\pas{1+\frac{1}{p}}\sum_{n=0}^\infty\frac{1}{p^{2n}}
=\sum_{n=0}\frac{1}{p^{2n}}+\sum_{n=0}^\infty\frac{1}{p^{2n+1}}=\sum_{n=0}^\infty\frac{1}{p^n}
\end{align}
が成り立つ. そこで, $x$以下の素数すべてに対して積をとる:
\begin{align}
\prod_{p\leq x}\pas{1+\frac{1}{p}}\sum_{n=0}^\infty\frac{1}{p^{2n}}=\prod_{p\leq x}\sum_{n=0}^\infty\frac{1}{p^n}.
\end{align}
$x$以下の素数の積で表せる自然数全体を$A_x$とおくと,
\begin{align}
\prod_{p\leq x}\pas{1+\frac{1}{p}}\times\sum_{n\in A_x\cup\Pas{1}}\frac{1}{n^2}=\sum_{n\in A_x\cup\Pas{1}}\frac{1}{n}
\end{align}
命題$1$より, $\sum_{n\in A_x\cup\Pas{1}}n^{-2}$は$x\to\infty$で収束し, $\sum_{n\in A_x\cup\Pas{1}}n^{-1}$は正の無限大に発散する. 即ち,
\begin{align}
\lim_{x\to\infty}\prod_{p\leq x}\pas{1+\frac{1}{p}}=+\infty.
\end{align}
ところで, 正実数$\alpha$に対して$1+\alpha< e^\alpha$を満たすから,
\begin{align}
\prod_{p\leq x}\pas{1+\frac{1}{p}}\leq \prod_{p\leq x}e^{\frac{1}{p}}=\exp\pas{\sum_{p\leq x}\frac{1}{p}}.
\end{align}
従って, 追い出しの原理より
\begin{align}
\lim_{x\to\infty}\sum_{p\leq x}\frac{1}{p}=+\infty
\end{align}
が成り立つ.
どの証明方法も美しく, かつ簡潔でいいですね. 私はErdősの証明が好きです. よくできすぎているところに神秘性を感じるのもありますが, 一番の理由は自然数の逆数和が発散することを用いずに証明できるからです. みなさんはどの証明方法がお好みですか?
$\;$[1] 素数とゼータ関数, 小山信也, 共立出版
$\quad$[2] Introduction to Analytic Number Theory, T.M.Apostol, Springer.
$\quad$[3] 数論2-類体論とは, 加藤和也・黒川信重・斎藤毅, 岩波書店