说不尽的圆周率

前言

2019年11月26日，联合国教科文组织第四十届大会宣布，将每年的3月14日定为国际数学日 (IDM, International Day of Mathematics)，即俗称的$\pi$日。有人还挖出3月14日15时16分9秒，称为所谓的$\pi$时（$\pi$秒）。可见，圆周率是多么的迷人。作为中国人，最引以为豪的就是祖冲之运用割圆术取得的成就：3.1415926，及密率：$\frac{355}{113}$，和约率：$\frac{22}{7}$ 了。不过，采用割圆术计算圆周率，尤其是在现代的计算机出现之前，靠人工计算，计算量实在太大了，难有突破。圆周率计算的真正突破，是有了无穷级数的算法之后，特别是现代计算机及计算技术的发展，使得圆周率的计算如井喷一般飞速发展。截至2021年6月14日，谷歌宣布圆周率已经计算到小数点后31.4万亿位。本文介绍几种常见的计算圆周率的算法，献给$\pi$日。

运用割圆术计算圆周率的迭代公式

最早用割圆术计算圆周率的是阿基米德（Archimedes）。他用单位圆的外切正多边形和内接正多边形的周长来逐步逼近单位圆的周长（$2\pi$），从而计算出圆周率。令$a_i$为单位圆的外切正多边形的周长，$b_i$为单位圆的内接正多边形的周长，从正六边形开始计算，每次倍增。即$i=0$对应正6边形，$i=1$对应正12边形，$i=2$对应正24边形，以此类推。迭代公式如下：
$$\{\begin{array}{l}{a}_0=4\sqrt{3},b_0=6\\ \\ a_{n+1}=\end{array}$$