せきゅーんです。「素数と量子計算」という大それたタイトルですが、これらが関係あるのかは私は知りません。というか、量子計算の勉強は私はしたことがありません。
今回の記事の目的は、次の不等式を証明することです:
この問題は物理学を学んでいる後輩から質問を受けたものなのですが、どうやら量子計算・量子コンピュータに関する論文中に出てきたとのこと。該当の論文は
R. Cleve, A. Ekert, C. Macchiavello, Quantam Algorithms Revisited, http://arxiv.org/pdf/quant-ph/9708016.pdf
で式(6.3)に先ほどの不等式が現れています。
ちゃんと読んではいないですが、二つの整数が互いに素であることが効いてくる場面があるらしく、確率計算の過程で用いるようです。
とにかく、後輩の質問に答えましょう。
解答1
示すべきことはです。ちなみに、この和は
なる値に収束しますので、確かに所望の不等式が成り立っていることが分かります。
これでは当然不満が残るでしょう。2003年の東大の入試問題
の解答として、「なので。」と書いてあっても、それではペケを付けざるを得ない採点者の気持ちです。もう少し論証っぽく見える証明を模索しましょう(何らかのレベルで計算に頼ることになりますが。この記事では「円周率」および「自然対数」の数値計算は満足にできるものと仮定します)。
解答2
いわゆるバーゼル問題に頼った素朴な解答を紹介します。バーゼル問題の証明についてはバーゼル問題の高校数学範囲内で分かる証明 - INTEGERSを参照してください。
任意の自然数に対して
が成り立つため、所望の不等式が得られるまでを大きくしていって素朴に計算するという証明法です。この方法を手持ちの関数電卓で試したところ、
と比較的少量の計算で証明できました。
解答3
もう一つの方法を紹介しましょう。あまり日本語の記事をみかけない気がするのですが、素数ゼータ関数と呼ばれる関数があります。とりあえずなる複素数に対しては絶対収束級数
と定義されます。つまり、素数ゼータ関数のにおける特殊値こそが今計算したい値なのです。
素数ゼータ関数に関する基本的性質として
が成り立つことをメビウス関数 - インテジャーズ
で証明しました。ここで、は同じ記事で定義したMöbius関数です。この公式より、
が成り立つため、Riemannゼータの偶数での値を計算できればの値も計算できることが分かります。ところが、リーマンゼータ関数 - INTEGERS
で紹介したEulerの公式によって円周率(と関-Bernoulli数)を使ってを計算することができます。この方針でを証明してみましょう。
Möbius関数の定義より無平方であるようなのみを考えれば十分であることに注意しておきます。とりあえず、
を用いて
を計算するととの真の値にそれなりに近い値が得られます(解答2の方法より優秀)。ここまで計算しておけば、かなりいい加減な評価でで押さえられることを論証できます(いい加減にせず同じ計算をで続ければどんどん真の値に近づきます)。
を用いると、
が得られます。ではとなるだけをピックアップして、までは全て上記値で計算することにより上から評価しましょう(は単調減少であることに注意)。なお、からまでにとなるようなものは
の個あります。については次の極めて粗い評価を行います(の次にとなるのが)。Riemannゼータの記事で証明した不等式により、
が成り立つので、に対して成立する不等式を用いると、
と評価できます。以上により、
を得ます。
追記
②の最初の二項の和は
ですが、
という変形に気づけば①の証明で使った式
より、
が論証できます。この解法はTwitterで知りました(解答もつぶやかれています):
— ずきの (@w7t8h6) 2017年2月5日