关于π的趣事

哪里有数，哪里就有美———   Proclus Proclus (希腊，公元 410-485 年) 

传说华罗庚曾经讲过一个故事：从前有个教书先生，平日喜欢喝酒。一天，他把学生们关在教室里，要他们背诵圆周率，自己却提了一壶酒到山上的寺庙找老和尚对饮去了。先生回来后，考问学生有没有把圆周率背下来？有个学生就回应道：“山巅一寺一壶酒，尔乐苦煞吾，把酒吃，酒杀尔杀不死，乐尔乐”。先生开始很生气，以为学生是在讽刺他，但很快便意识到这其实是带浓厚乡音的普通话3.1415926535897932384626，顿时转怒为喜，把这个学生着实称赞了一番。  众所周知，圆周率是一个圆的周长与直径之比，而且不论圆的大小，这个比值是常数，永远不变。圆周率的近似值π ≈ 3.14 家喻户晓，关于它的文章和书籍更是多不胜数。由于希腊字母π 和英文单词 pie 同音，1988 年 3 月 14 日由美国旧金山的科学探索馆（San Francisco Exploratorium）提议，发起了每年这一天以吃 pie 来纪念和庆祝圆周率的校园活动。说来奇巧，3 月 14 日正好也是爱因斯坦的生日。2009 年 3 月 11 日，美国众议院还通过了一项决议，把 3 月 14 日正式确定为全国的“π 日”（National Pi Day）。更有趣的是，麻省理工学院的招生办公室每年都把新生录取通知书留到π 日才发送给新生。现在网络通信方便了，还有传闻说那里的工作人员安装了个计算机程序等到 3 月 14 日下午 1 时 59 分 26 秒才从网上把录取通知向全世界发放出去，对应于π ≈3.1415926。另外，因为 22/7 是比3.14 对π 更为精确的近似，每年 7 月 22 日也被公认为“π 近似日”。

                                                               （每年 7 月 22 日是“π 近似日”）

圆周率π 是一个不等于任何两个整数之比的实数，它在整数 3 之后再带上一个具有无限长度但又永远不循环的小数：π ≈ 3.1415926…。尽管人类对π 的认知可以追溯到远古，最先对π 值进行系统严格的估算者应当首推古希腊科学家阿基米德（公元前 287-212 年)，他得出不等式 3+10/71 <π < 3+1/7，平均值就是π ≈3.14…。在我国，三国后期魏国人刘徽（生于公元 250 年左右），留下了宝贵数学遗产《九章算术注》和《海岛算经》，他创始并使用的“割圆术”（即用圆的内接和外切正多边形来逼近圆的周长）为后来南北朝时期的数学家祖冲之（公元 429-500 年）的估算 3.1415926 <π < 3.1415927（密率π ≈355/113）提供了最基本的方法。前面提到的π ≈22/7，就是祖冲之给出的圆周率的约率。 

 

古人计算圆周率，一般都是用正多边形来逼近圆的。阿基米德用正 96 边形逼近圆而得到π 小数点后 3 位的精度，刘徽用正 3072 边形逼近圆也只能得到π 小数点后5 位的精度。由于这些计算方法效率极低，估算π 的历史进程十分缓慢。直到微积分问世以后，情况才大为改观。1706 年，英国天文学家 John Machin 发现了一个简单的解析公式并用它来计算π 达到了小数点后 100 位的精确度： π−=arctgarctg)239/1(4)5/1(16  微积分的出现，不但给出了许多关于π 的解析估计，更大大地加快了其数值计算。且不说优雅漂亮的 Gregory-Leibniz 和差公式 

 

后来发现能够用来对π 作快速近似计算的公式还有很多。1914 年，自学成才的传奇印度数学家 Srinivasa Ramanujan 就曾经写下过 14 个关于π 的无穷级数展开公式。电子计算机出现以后，人类大规模高精度的计算能力得到了无以伦比的飞跃，使得圆周率计算的进程突飞猛进，实现了计算π 的（二进制）数字长度的一次又一次的重大突破。1985 年，美国数学家 Bill Gosper 用 Ramanujan 的一个公式计算圆周率精确到了小数点后 17,500,000 位。1994 年，苏联裔美国数学家 David 和Gregory Chudnovsky 兄弟在同一个公式的基础上计算π 到了 4,044,000,000 位。1995 以后，日本数学家金田康正（Yasumasa Kanada）及其团队更在世界上遥遥领先，在 2003 年把对π 的计算提高到了 1.24 万亿位。这个纪录在今年 1 月 8 日由法国一位程序员 Fabrice Bellard 打破。他改进了 Chudnovsky 兄弟的公式，用了131 天在一部台式计算机上成功地把π 计算到 2.7 万亿位。这个数字有多长呢？如果您平均用一秒钟来报出一个数字的话，那么您需要 8.5 万年才能读完它的所有数字。厉害吧？可是，这个世界纪录到 8 月份就被打破了：日本工程师近藤茂（Shigeru Kondo）与美国西北大学计算机系的香港学生余智恒合作，把π 值算到了 5 万亿位。谁知道，当欢呼声犹在耳际，9 月 17 日英国广播公司 BBC 又作出了惊人的报道，说 Yahoo 科技公司的研究员、原香港科技大学毕业生施子和（Nicholas Tse-Wo Sze）采用“云计算”技术，利用 1000 台计算机同时计算，历时 23 天，将π 计算到小数点后 2 千万亿位，为 8 个月前法国人 Bellard 记录的近一千倍。呵呵，这可不是只增加了大约一千个数位，而是由一千个、每个都有 2.7万亿位长的数字连接起来。No kidding! 施子和兴奋无比地向全世界宣布了他的计算结果，还说他发现了“The two quadrillionth bit of π  is 0!” 

您或许不以为然，说用计算机来计算不足为奇。那么，用人脑来进行计算和记忆就不容易了吧？如果我告诉您，说有人能正确无误地记住π 的数值从 3.14 开始一直到小数点后几万位数字，您会觉得匪夷所思吧？其实，多少年来一直有着凭人脑记忆来背诵π  的数值的吉尼斯世界纪录（GuinnessWorld Records），目前领先的是中国人吕超。2005 年 11 月 20 日，当时是西北农林科技大学研究生的吕超用了 24 小时零 4 分的时间，连续背诵出π  的值到小数点后 67890 位。 

 

读到这里，我知道您如果不是数学粉丝的话其实并不会特别激动：知道π ≈ 3.14 就已经足够了，谁去关心它长长尾巴上的无穷多位数字呢？  说来也是。造物主有这么多东西可以创造，为什么偏偏造出这个神奇的π 来，让许许多多的数学家和粉丝们两千多年来呕心沥血地去追寻不休呢？不过，我到最近才开始明白，原来造物主并不是像一些人所想象的那么无聊。 

你可能早已听说过，π 与地球上的许多河流有关：河流弯曲河道的曲线长度与河道首尾直线距离之比通常都接近于 3.14 —— 河道越是蜿蜒曲折，这个近似值就越好—— 亚马逊河便是一个例子。最令我惊叹的是，今年 11 月《科学》杂志上报道了德国哥廷根大学马克-普朗克动力学与自组织科学研究所和伯恩斯坦计算神经科学研究中心的科学家 Matthias Kaschube 及其团队的一项研究成果，说他们在相同哺乳类中的三种远亲动物嬰猴(galago)、樹鼩(tree shrew)和矇眼貂(ferret)的大脑中发现，它们的视觉皮层的螺旋结构（不妨类比于太空中银河系的螺旋结构）平均分布密度都是一样的，而且都落在与π 的值相同即 3.14 的 %95 置信区间之中。哇噻，这不会是巧合吧？虽然这个 3.14 不能与用计算机得到的 2 千万亿位数字相提并论，但对于理论生物学家来说，π ±%2 已经是一个令他们难以置信的超级精确度了。 

 

科学家们相信，今天的生物大脑都是经过漫长时间通过自组织演化过程而逐渐形成的，并且都具有某种最优的结构和功能。如果是这样的话，那些哺乳类动物大脑视觉皮层中基本螺旋结构分布都不约而同地呈现出以 3.14 为平均密度之特征，是不是自然进化的最优结果呢？不会是造物主在创世开始时，不经意地在它们的大脑中写下了那个神奇的数字π 吧？