円周率

円周率とは、円の直径に対する円周の長さの比率のことをいい板倉聖宣、数学定数の一つである。通常、円周率はギリシア文字である π}}group="注" に代表される。円の直径から円周の長さや円の面積を求めるときに用いる板倉聖宣。また、数学をはじめ、物理学、工学といった科学の様々な理論の計算式にも出現し、最も重要な数学定数とも言われる。

円周率は無理数であり、その小数展開は循環しない。さらに、円周率は無理数であるのみならず、超越数でもある。

円周率の計算において功績のあったルドルフ・ファン・クーレンに因み、ルドルフ数とも呼ばれる。ルドルフは小数点以下35桁まで計算した。小数点以下35桁までの値は次の通りである。

1. 基礎
-----

1.1. 表記と呼び方
-----------

円周率を表すギリシア文字 π}} は、ギリシア語でいずれも周辺・円周・周を意味する el日本数学会あるいは elの頭文字から取られたgroup="注"。文字 をウィリアム・オートレッドは1631年に著した著書において半円周の長さを表す文字として用い、アイザック・バローは論文において半径 R}} の円周の長さとして用いた。ウィリアム・ジョーンズ  やレオンハルト・オイラーらにより円周の直径に対する比率を表す記号として用いられ、それが広まった日本数学会。日本では「パイ」と発音する。

数 を指す言葉には、日本・中国・韓国における「円周率」、ドイツの「Kreiszahl」の他、それを計算した人物の名前を取った「アルキメデス数」、「ルドルフ数」などがある。一般にドイツ語を除くヨーロッパの諸言語には「円周率」に対応する単語はない。

なお、「」の字体は、表示環境によってはキリル文字の ru に近い el などと表示されることがある。また、ギリシャ文字「」は、円周率とは無関係に、素数計数関数や、基本群・ホモトピー群、ある種の写像を表すのに用いられることもある。

1.2. 定義
-------

平面幾何学において、円周率 は、円の周長の直径に対する比率として定義される。すなわち、円の周長を C}}, 直径を d}} としたとき、

である。全ての円は互いに相似なので、この比率は円の大きさに依らず一定である。

ところが、この定義は円の周長を用いているため、曲線の長さを最初に定義していない解析学などの分野では、 が現れる際に問題となることがある。この場合、円の周長に言及せず、解析学などにおける性質の一つを の定義とすることが多いRudin。この際の の定義の一般なものとして、三角関数 cos //x//}} が 0}} を取るような //x// > 0}} の最小値の2倍とするもの、級数による定義、定積分による定義などがある。後述の#円周率に関する式も参照。

2. 歴史
-----

2.1. 古代
-------

円周の直径に対する比率は円の大きさに依らず一定であり、それは 3 より少し大きいことは古代エジプトやバビロニア、インド、ギリシアの幾何学者たちにはすでに知られていた。また、古代インドやギリシアの数学者たちの間では半径 r}} の円板の面積が //πr//sup であることも知られていた。さらに、アルキメデスは正96}}角形を用いて半径 r}} の球の体積が sfrac であることや、この球の表面積が 4//πr//supであることを導き出した。中国では 263年に魏の劉徽が3072角形を使用し3.14159と計算し、5世紀に祖沖之は祖沖之が十尺もの直径の円を使用して3.14159 26<<3.14159 27 と 求め、以後1000年これ以上正確な計算はなされなかった。祖の計算が正確であったことは、1300年頃に趙友欽が16384辺の内接多角形により確かめた。

2.2. 近代まで
---------

14世紀インドの数学者・天文学者であるサンガマグラーマのマーダヴァは次の の級数表示を見いだしている：

これは逆正接関数 Arctan //x//}} のテイラー展開の //x// =}} 1}} での表式になっている。マーダヴァはまた、

を用いて の値を小数点以下11桁まで求めている。

17世紀、ドイツのルドルフ・ファン・コーレンが正325}}億角形を使い、小数点以下第35位まで計算。1699年にエイブラハム・シャープが小数点以下第72～127位まで求めた。

18世紀フランスの数学者アブラーム・ド・モアブルは、ある定数 C}} を取ると、コインを 2//n//}}回投げて表が x}}回だけ出る確率は、n}} が十分大きいとき

で近似できることを、//n// =}} 900}} における数値計算により見いだした。この正規分布の概念は1738年に出版されたド・モアブルの『巡り合わせの理論』に現れている。ド・モアブルの友人のジェイムズ・スターリングは後に、であることを示した。

1751年にヨハン・ハインリヒ・ランベルトは、x}} が 0}} でない有理数ならば正接関数 tan //x//}} の値は無理数であることを示し、その系として は無理数であることを導いた。さらに1882年にフェルディナント・フォン・リンデマンは が超越数であることを示し、円積問題は解くことができないことを導いた。

1873年、ウィリアム・シャンクスは彼自身の手で小数点以下第707位までを計算した。

2.3. 和算における円周率の取り扱い
-------------------

2.3.1. 江戸時代の初期の和算家の3.16
-----------------------

江戸時代初期の数学書である毛利重忠の『割算書』では円周率を3.16としている。その弟子の吉田光由の『塵劫記』でも3.16となっている板倉聖宣。しかし、当時の先進国中国の文献にはこの3.16という数値は見られず、中国の文献の数値を引き写したとは考えにくいという板倉聖宣。そのため、なぜ初期の和算家が円周率を3.16としたかの理由はよく分かっていない板倉聖宣。おそらく、毛利重忠とその弟子の吉田光由などの先駆者らは、円周率を実際に測定して3.14ないし3.16ほどの値を得たが、最後の桁の数字に確信が持てなかったため、「円のような美しい形を求める数値は、もっと美しい数値になっていいはずだ」と考え、「美しい理論」を求めた。その結果 sqrt が美しい数値として採用されたと推測されている板倉聖宣。その考えは日本で2番目に3.14の値を計算で求めた野沢定長の『算九回』の中にも見られ、その著書の中で「忽然として円算の妙を悟った」として「円周率の値は形＝経験によって求めれば3.14であるが、理＝思弁によって求めれば3.16である」として「両方とも捨てるべきでない」とした板倉聖宣。

2.3.2. 和算家が計算した3.14
-------------------

江戸初期、1600年代前半頃から、円を対象とした和算的研究である「円理」が始まる。その最初のテーマの一つが円周率を数学的に計算する努力であり中村邦光、1663年に日本で初めて村松茂清が『算爼』において「円の内接多角形の周の長さを計算する方法」で3.14…という値を算出した中村邦光。『算爼』では円に内接する正8角形から角数を順次2倍していき、内接2sup角形の周の長さで、

3.1415 9264 8777 6988 6924 8と小数点以下21桁まで算出している。これは実際の値と小数第7位まで一致している中村邦光。その後1680年代に入ると、円周率の値を3.16とする数学書はなくなり、3.14に統一された中村邦光。1681年頃には関孝和が内接2sup角形の計算を工夫し、小数第16位まで現代の値と同じ数値を算出した。この計算値は関の死後1712年に刊行された『括要算法』に記されている中村邦光。

日本の和算家に特徴的なのは、1663年に3.14が初めて導き出されても、その後1673年までの10年間に円周率の値を3.14とした算数書のいずれもが、先行者の円周率をそのまま引き継ぐことをせず、それぞれ独自の値を提出していたことである中村・板倉。この背景には当時の遺題継承group="注"運動に「他人の算法をうけつぐ」と共に「自己の算法を誇る」という性格があったためだという中村・板倉。そのため古い3.16の値が疑われてから、遺題継承の際に必ずといってよいほど円周率の値が変えられている中村・板倉。しかしながら江戸時代の3大和算書『塵劫記』『改算記』『算法闕疑抄』の増補改訂版では1680年代には3.14に統一された中村・板倉。

2.3.3. 3.14から3.16への逆行
---------------------

しかし、遺題継承運動は1641年に始まって1699年頃には終わってしまい板倉聖宣、いったん3.14に統一された円周率の値は江戸時代後半になると揺らぎ始め、古い3.16に逆行するという現象が生じた板倉・中村。文政年間に出版された算数書とソロバン書を悉皆調査した結果では、円周率の値を3.14とするものと、3.16とするものの2系統があることが明らかにされた板倉・中村。いくらか専門的な数学書では3.14とされているのに、大衆向けの小冊子の中では3.16の方が普通に用いられていた中村・板倉。

当時の識者である橘南谿は「いまに至り3.16あるいは3.14色々に論ずれども、なおきわめがたきところあり」と述べ、3.14はまだ確定していないとしている中村・板倉。儒学者の荻生徂徠も和算家の算出した3.14の根拠に納得しなかった中村・板倉。当時の和算家のほとんどは、円に内接する多角形の周を計算することで円周率を計算した。内接多角形の角数を増やすほど求まる円周率の桁は増えていくので、素人目にはその値が増大する一方に見える。「それがいくら増えても3.1416を超えない」ということを和算家たちはついに納得させることができなかったのである中村・板倉。

そのような和算家以外の素人たちを納得させるには、どうしても万人に納得させる「理」に基づいて計算してみせる他はない中村・板倉。それを行うには西洋で行われたように、「円を内接多角形と外接多角形ではさんで、円周率の上限と下限を示すこと」が必要であったが、和算家はついにその方法を取ることがなかった中村・板倉。

2.3.4. 宅間流和算の円周率
----------------

日本で唯一「円周を内接・外接多角形で挟み込んで円周率の上限と下限を示す」ことに成功したのは鎌田俊清が享保七年に著した『宅間流group="注"円理』である。その値は以下の通りである中村邦光。

内周：3.1415 9265 3589 7932 3846 2643 3665 8外周：3.1415 9265 3589 7932 3846 2643 4166 7鎌田は円周率の小数点以下24桁まで正しいと確信しうる円周率の値を算出することに成功していた中村・板倉。しかし、鎌田の方法は後継者を持たず、当時の識者に知られることがなかった中村・板倉。

2.3.5. 級数展開による算出
----------------

日本の数学史では級数による値の算出は広く一般的であった。円周率の級数による公式は多くの学者に研究されており、蜂谷定章、松永良弼、坂部広畔、川井久徳、長谷川寛らによるものがある。また、建部賢弘は円周率の二乗を求める日本初の公式を考案した。

2.3.6. 和算の限界
------------

日本の和算の弱点は単に理論面の弱さにとどまらず、万人が納得できる正しい円周率の教育・啓蒙への関心も失ったことであった中村邦光。そのため和算家たちがいくら円周率は3.14…と書いたところで、『塵劫記』の古い円周率3.16の値がそのまま残存する結果となった中村・板倉。『塵劫記』の重版などは古い円周率3.16のまま出版され続け、18世紀に大衆的な通俗算数書が大量に出版される際に、必ずというほど3.16という値を引き継ぐようになってしまった中村・板倉。

18世紀半ば以降の和算は数学的証明の概念の追求は無視され、せっかく宅間流の鎌田俊清がその独創的方法で正しい円周率を算出しても、全く継承されなかった中村邦光。江戸時代後半の和算家は家元制度的な秘密主義と保守主義と、権威主義が在野の独創性を無視し、結果として学問の進歩を妨げることとなった中村邦光。

2.4. コンピュータの利用
--------------

20世紀以降、計算機の発達により、計算された円周率の桁数は飛躍的に増大した。1949年に、電子計算機ENIACを使い72時間かけて、円周率は2037桁まで計算された。その後の数十年間、様々な計算機科学者や計算科学者など、あるいはコンピュータのアマチュアによって計算は進められ、1973年には100万桁を超えた。この進歩は、スーパーコンピュータの開発だけによるものではなく、効率の良いアルゴリズムが考案されたためである。そのうちの最も重要な発見の一つとして、1960年代の高速フーリエ変換がある。これにより、多倍長の演算が高速に実行できるようになった。

2022年6月9日に、Googleの技術者、岩尾エマはるかがGoogle Cloudで、チュドノフスキー級数を使い、157日23時間かけて100兆桁を計算したと発表。

3. 性質
-----

3.1. 無理性
--------

は無理数である。つまり、2つの整数の商で表すことはできず、小数展開は循環しない。このことは1761年にヨハン・ハインリヒ・ランベルトが証明したが、厳密性に欠けた部分があった。その部分は1806年にルジャンドルによって補われた。

したがって、円周率のコンピュータによる計算や暗唱、十進法表示での小数部分の各数字  の出現頻度は、人々の興味の対象となる。

3.2. 超越性
--------

さらに、 は超越数である。つまり、有理数係数の代数方程式の解にはならない。これは1882年にフェルディナント・フォン・リンデマンによって証明された。これより、整数から四則演算と冪根をとる操作だけを有限回組み合わせてもけっして の値をとることはできないことが分かる。

が超越数であることより、古代ギリシアの三大作図問題の内の一つである「円積問題」が不可能であることが従う。

3.3. ランダム性
----------

2022年10月の時点で、 は小数点以下100兆桁まで計算されている。そして、分かっている限りでは 0 から 9 までの数字がランダムに現れているようには見えるが、それが乱数列といえるかどうかははっきりとは分かっていない。たとえば が正規数であるかどうかも分かっていない。正規数であれば の10進表示において、各桁を順に取り出して得られる数列：

3, 1, 4, 1, 5, 9, 2, 6, 5, 3, 5, …には、0 から 9 が均等に現れるはずだが分かっておらず、それどころか、0 から 9 がそれぞれ無数に現れるのかどうかすら分かっていない。もし仮に正規数でないとすれば、乱数列でもないということになる。

5兆桁までの数字の出現回数は以下の通りである。全てほぼ等しく、最も多いのは 8 で、最も少ないのは 6 である。

0：4999億9897万6328回 1：4999億9996万6055回 2：5000億0070万5108回 3：5000億0015万1332回 4：5000億0026万8680回 5：4999億9949万4448回 6：4999億9893万6471回 7：5000億0000万4756回 8：5000億0121万8003回 9：5000億0027万8819回

3.4. 連分数
--------

分母を整数と分数の和で表すことを続けていった表示を連分数という。「整数」を最大にしていくと、分子を全て 1 にできる：

は無理数であるから、円周率 の連分数展開は有限項では終わらず無限項の連分数となる：

上記の正則連分数展開を途中で打ち切ると、 の良い有理数近似が得られる。その最初の4つは 3, sfrac である。これらは古くからよく知られ使用されてきた近似値である。これらはそれぞれ分母が大きくないどの分数よりも に近く、 の最良有理数近似である。

さらに、 は超越数であることが知られている。一般に、正則連分数の分母に現れる整数部が循環するのは二次無理数に限られ、 は二次無理数でないため循環連分数として表せない。加えて の正則連分数は規則性を示さないが、 の一般化連分数では以下の規則をもつものが知られている：

4. 未解決問題
--------

次の数は超越数か？
は正規数か？

5. 円周率に関する式
-----------

についての式は非常に多い。ここではその一部を紹介する。数式によってはそれ自体が の定義になり得るし、 の近似値の計算などにも使われてきた。

5.1. 幾何
-------

半径 r}} の円の周長：
半径 r}} の円の面積：
半径 r}} の球の体積：
半径 r}} の球の表面積：
長半径 , 短半径 の楕円の面積：
rad

5.2. 解析
-------

ライプニッツの公式、#2千年紀、逆正接関数も参照
マーダヴァの公式、逆正接関数も参照
ウォリスの公式
ビエトの公式
逆三角関数は主値を取るものとすると
逆三角関数の公式より
逆三角関数の公式より
逆正接関数のテイラー展開による：
オイラーによる：
双曲線関数の公式より
ニュートンの無平方根公式 この式は、 のマクローリン級数となっている。
マチンの公式 と書かれることもある。 4}} と sfrac が二進法と相性が良く、収束も早いため、コンピュータでの円周率計算によく使われる公式の一つである。
sfrac の連分数表示
ガウス＝ルジャンドルのアルゴリズム 初期値の設定： 反復式：//asub が希望する桁数になるまで以下の計算を繰り返す。小数第n}}位まで求めるとき logsub回程度の反復でよい。//////// // の算出：円周率 は、//asub を用いて以下のように近似される。//////// // 非常に収束が早く、金田康正が1995年に42億桁、2002年に1.24兆桁を計算したスーパー に使われていた。

David Bailey, Peter Borwein, およびサイモン・プラウフによるもの、Adamchik と Wagon によるもの、Fabrice Bellard によるもの等については、あまりに高度なため割愛する。

5.3. 複素解析
---------

後者はオイラーの等式の一般化であり、1}} の n}}乗根の総和は 0}} になることを示している。//n// =}} 2}} とするとオイラーの等式になる。

5.4. 特殊関数
---------

5.5. 数論
-------

整数全体から無作為に2つ取り出すとき、その2つが互いに素である確率は sfrac である。

5.6. 力学系・エルゴード理論
----------------

ロジスティック写像 1=//x//sub により帰納的に定まる数列 mset を考える。初期値 //x//sub を 0}} 以上 1}} 以下に取るとき、そのほとんど全てで、次が成り立つ。

5.7. 統計
-------

幅 1}} の無数の平行線の上から長さ sfrac の針を落とすとき、その針が直線と共有点を持つ確率は sfrac である。

5.8. その他
--------

河川の長さの水源・河口間の直線距離に対する比率は、平均すると円周率に近い。

6. 暗唱
-----

6.1. 語呂合わせ
----------

の桁を記憶術に頼らずに暗記する方法が各種存在している。

日本語では、語呂合わせにより、長い桁を暗記するのも比較的簡単である。有名なものとして、以下がある。