今日は3月14日。そう、円周率の日です*1。
というわけで、今日は整数ではなく、円周率のお話をしましょう。
ラマヌジャン(Ramanujan)のMysteriousな公式
円周率に関する級数表示が多数知られており、
integers.hatenablog.com
でLeibnizの公式、Machinの公式、高野喜久雄の公式などを紹介しました。
今日は、2015年にイギリスで公開された映画"The Man Who Knew Infinity"でも取り上げられた"インドの魔術師"Ramanujanの発見した円周率の不思議な公式を紹介します。
Ramanujanについては"タクシー数"のエピソードが有名で
integers.hatenablog.com
で取り上げましたので是非ご覧ください。
次の公式がRamanujanが1914年に発表した円周率に関する驚くべき公式です*2:
私がこの公式に初めて出会ったのは高校生のときで、驚いて、興奮が止まりませんでした。
「タクシー数のような陳腐な話ではなく、この公式は本物だ。Ramanujanは本物の化け物である!」
と*3。
高校生のとき、誕生日プレゼントとして関数電卓を買い与えられた私は関数電卓で様々な計算を行うことが日課でした。Ramanujanの公式を知ったときも即座に関数電卓で試してみたのです。
その関数電卓は今も健在ですので、久しぶりに計算してみました:
公式の和においての部分のみを取った値からを計算すると
となって、小数点以下5位まで一致することがわかります(何も足さなくてもこの精度!!)。
次に、まで和をとって計算すると…
この関数電卓の精度ではと判断されてしまいました(笑)。
この公式は実際にの数値計算に用いられ、一時的に記録を保持したこともあるのですが、証明は中々付けられないままでした。
証明されたのはRamanujanが生まれてから丁度100年にあたる1987年で、の400ページを超える本
J. M. Borwein, P. B. Borwein, Pi and the AGM - A Studey in Analytic Number Theory and Computational Complexity, Wiley, N.Y., 1987.
に証明が掲載されています*4。
その後、も含めてRamanujan型の公式が多数発見されています。他にどのような公式があるかは、例えば、Wikipediaの記事
Ramanujan–Sato series - Wikipedia
を参照してみてください。
百年に一度の円周率の日から一年
円周率の日は毎年やってくるのですが、実は昨年の3月14日は『百年に一度の円周率の日』と言われていました。
理由はアメリカ式日付表記で3/14/15と表されるからというものです。この表記でいくと今年は3/14/16なので、近似値の観点から考えると今年の方が値が近かったりします。
さて、アメリカではパイを食べてお祝いしたりするそうなのですが、日本では『終わらない数の象徴』という意味合いが強い円周率にあやかってプロポーズや入籍、結婚式を挙げることなどが流行っているようです。円周率に永遠の愛を誓うということですね*5。
実は私の姉が昨年の円周率の日に入籍しました。
普段、数学と何らの関わりも持たない二人がこのような日を入籍記念日に選んでくれたことに感動した弟は「円周率は確かに終わらない数である!」ということの証明をプレゼントしました*6。
最も簡単な証明はNivenによるものであるが、私の最もお気に入りの証明はBeukersによる証明である。
と書いてありますが、Nivenの証明については既に記事にしています:
Beukersの証明というのは1979年の「Apéryによるの無理数性証明の別証明」に関する有名な論文に載っている証明です。の無理性証明がメインですが、同じ論法での無理性も証明できるため同時掲載されています。そうして、バーゼル問題
と合わせることによってが無理数であることが従い、特にも無理数であるという寸法です*7。なお、バーゼル問題についてはバーゼル問題の高校数学範囲内で分かる証明 - INTEGERSにおいて高校生でも分かる証明を書いておきました。
証明の解説
姉夫婦へのプレゼントは一ページにまとめたので当然省略があります。そこで、今回は(2人がこの記事を読むかは不明ですが)もう少し詳細の解説を行いたいと思います。これは、後にの無理性証明を紹介する際のウォーミングアップにもなります*8。
これは前回の記事数列lcm[1,2,…,n]のgrowthと素数定理 - INTEGERSで証明しました。
証明. のケースはを幾何級数展開して逐次積分すれば簡単に示される。の場合は望遠鏡和を作って証明するが、この補題を一般化したものを収束性なども厳密にみて別の記事で紹介するため、今回の記事では省略する*9。 Q.E.D.
証明. 定義式においてを二項展開してから階微分すればよい。 Q.E.D.
証明. 回部分積分を行えばよい。の定義においてLeibniz則を用いて微分を実行することにより、各部分積分において出てくる第一項は消えることがわかる。 Q.E.D.
証明. 大学一年生で数学を勉強した者は証明できる。後で時間が出来たら証明を簡単に書きます。 Q.E.D.
が無理数であることの証明
とする。このとき、補題2より整数列が存在してが成り立つ。一方、補題4を用いて積分変形を行えば、
と変形される。被積分関数を見れば明らかにであるので、補題5より
が得られる。よって、補題1から十分大きいに対してはと評価することによって
となる。が有理数、例えば、であったならば
となってしまう。すなわち、十分大きいを考えることによって我々はとの間にある整数を獲得してしまう。 Q.E.D.
*1:数学の日とも呼ばれています。英語では"Pi Day"です。ちなみにEinsteinの誕生日でもあります。更に、今年は一般相対論百周年!!
*2:S. Ramanujan, Modular equations and approximations to pi, Quart. J. Math. (Oxford) 45 (1914)で発表した公式のうちの一つです。
*3:陳腐と表現したのは、タクシー数の話は面白いけれども、私はそんなに凄いとは思わないという話です。
*4:Ramanujanの公式の証明を短くまとめた文献が中々見つからないので、いつか記事にしたいと思っています。
*5:終わらない数というのは数学的に言えば無理数ということです(とします)。進法表記に限定するならば無理数でなくとも終わらない数はあるのですが、任意ベースで終わらないという性質の方が良いと思われるので、無理数と考えることにしましょう。あるいは円周率が進法で巡回しない無限小数表示を持つ直接証明を探すくらいなら無理数性を示した方がよいという考え方も妥当です。さて、無理数は非加算無限に存在するわけですが、このような流行は「円周率が非常に慣れ親しまれた無理数である」ということを表しているのだと思われます。ここで、パラドキシカルなことは「円周率の無理数性」を知っている国民はどれぐらいいるのだろうか?ということです。だって、証明を知らなければ実は永遠の愛ではなかったとなって愛が消滅する可能性だってなきにしもあらず。。。ただ、証明が簡単なで類似を辿ろうとすると、1月41日がないことも影響していると考えられます。なんという無駄な考察!!
*6:二週間研究集会で出張している合間だったので作成が大変でした。英語で書いた理由は「見た目がかっこよいから」です。
*7:なぜこの証明法が一番好きであるかはここでは述べませんが個人的に私に会う機会がある方は聞いてください。