梅森素数--美丽的贝壳

一、价值五万美元的素数

   2000年4月6日,住在美国密歇根州普利茅茨的那扬·哈吉拉特瓦
拉(Nayan Hajratwala)先生得到了一笔五万美元的数学奖金,因为他
找到了迄今为止已知的最大素数,这是一个梅森素数:
                    2^6972593-1。
这也是我们知道的第一个位数超过一百万位的素数。精确地讲,如果
把这个素数写成我们熟悉的十进制形式的话,它共有两百零九万八千
九百六十位数字,如果把它以这个形式写下来,大约需要150到200篇
本文的篇幅。

   可是哈吉拉特瓦拉先生并不是一个数学家,他甚至很可能对寻找
素数的数学理论一无所知——虽然这使他赢得了这笔奖金。他所做的
一切,就是从互联网上下载了一个程序。这个程序在他不使用他的奔
腾II350型计算机时悄悄地运行。在经过111天的计算后,上面所说的
这个素数被发现了。

二、梅森素数

   我们把一个大于1的自然数叫作素数,如果只有1和它本身可以整
除它。如果一个比1大的自然数不是素数,我们就叫它合数。1既不是
素数,也不是合数。

   比如说,你很容易就可以验证7是一个素数;而15是一个合数,因
为除了1和15外,3和5都可以整除15。根据定义,2是一个素数,它是
唯一的偶素数。早在公元前三百年的古希腊时代,伟大的数学家欧几
里德就证明了存在着无穷多个素数。

   关于素数,有许多既简单又美丽,但是极为困难的,到现在还没
有答案的问题。其中有著名的哥德巴赫猜想,它是说任何一个大于6的
偶数,都能表示为两个奇素数之和。还有孪生素数问题。象5和7,41
和43这样相差2的素数对,被称为孪生素数。孪生素数问题是说:是不
是有无穷多对孪生素数?这里要顺便提一下的是,这些看起来很简单
的数学问题,它们的解决方法将一定是极其复杂的,需要最先进的数
学工具。如果你不是狂妄到认为几百甚至几千年来所有在这些问题上
耗费了无数聪明才智的数学家(有许多是非常伟大的)和数学爱好者
加起来都不如你聪明,就不要试图用初等方法去解决这些问题,徒费
时间和精力。

   古希腊人还对另一种数感兴趣。他们将它称为完美数。一个大于1
的自然数叫完美数,如果它的所有因子(包括1,但不包括本身)之和
等于它本身。比如说6=1+2+3就是最小的完美数,古希腊人把它看作维
纳斯也就是爱情的象征。28=1+2+4+7+14是另一个完美数。欧几里德证
明了:一个偶数是完美数,当且仅当它具有如下形式:
                   2^(p-1)(2^p-1)
其中2^p-1是素数。上面的6和28对应着p=2和3的情况。我们只要找到
了一个形如2^p-1的素数,也就知道了一个偶完美数;我们只要找到所
有形如2^p-1的素数,也就找到了所有偶完美数。所以哈吉拉特瓦拉先
生不但找到了世界上已知的最大的素数,还找到了世界上已知的最大
的偶完美数。嗯,你要问,关于奇完美数又是怎么样的情况?回答是:
我们现在连一个奇完美数也没有找到过,我们甚至根本不知道是不是
有奇完美数存在。我们只知道,要是有奇完美数存在的话,它一定是
非常非常大的!奇完美数是否存在这个问题,也是一个上面所说的既
简单又美丽,但是极为困难的著名数学问题。

   有很长一段时间人们以为对于所有素数p,
                       M_p=2^p-1  
都是素数(注意到要使2^p-1是一个素数,p本身必须是一个素数,想
一想为什么?)但是在1536年雷吉乌斯(Hudalricus Regius)指出,
M_11=2^11-1=2047=23*89不是素数。

   皮特罗·卡塔尔迪(Pietro Cataldi)首先对这类数进行了系统的
研究。他在1603年宣布的结果中说,对于p=17,19,23,29,31和37,
2^p-1是素数。但是1640年费尔马使用著名的费尔马小定理(不要和那
个费尔马大定理混淆起来)证明了卡塔尔迪关于p=23和37的结果是错
误的,欧拉在1738年证明了p=29的结果也是错的,过后他又证明了关
于p=31的结论是正确的。值得指出的是,卡塔尔迪是用手工一个一个
验算取得他的结论的;而费尔马和欧拉则是使用了在他们那时最先进
的数学知识,避免了许多复杂的计算和因此可能造成的错误。

   法国神父梅森(Marin Mersenne)在1644年他发表了他的成果。他
宣称对于p=2,3,5,7,13,17,19,31,67,127和257,2^p-1都是
素数,而对于其它小于257的素数p,2^p-1都是合数。今天我们把形如
M_p=2^p-1的素数叫做梅森素数,M_p中的M就是梅森姓氏的第一个字母。

   用手工来判断一个很大的数是否素数是相当困难的,梅森神父自
己也承认他的计算并不一定准确。一直要等到一个世纪以后,在1750
年,欧拉宣布说找到了梅森神父的错误:M_41和M_47也是素数。可是
伟大如欧拉也会犯计算错误——事实上M_41和M_47都不是素数。不过
这可不是说梅森神父的结果就是对的。要等到1883年,也就是梅森神
父的结果宣布了两百多年后,第一个错误才被发现:M_61是一个素数。
然后其它四个错误也被找了出来:M_67和M_257不是素数,而M_89和
M_107是素数。直到1947年,对于p<=257的梅森素数M_p的正确结果才
被确定,也就是当p=2,3,5,7,13,17,19,31,61,89,107和
127时,M_p是素数。现在这个表已经被反复验证,一定不会有错误了。

   下面是我们现在知道的所有梅森素数的列表:(我们注意到梅森
神父的名字不在上面——这种素数已经由他的名字命名了,就把荣誉
分给最后确认者吧。)

序号   p           M_p的位数    相对应的   确认      确认人
                                  完美数的   年代
                                    位数
1            2            1            1      ----      ---- 
2            3            1            2      ----      ---- 
3            5            2            3      ----      ---- 
4            7            3            4      ----      ---- 
5           13            4            8      1456      佚名
6           17            6           10      1588      Cataldi 
7           19            6           12      1588      Cataldi 
8           31           10           19      1772      Euler 
9           61           19           37      1883      Pervushin 
10          89           27           54      1911      Powers 
11         107           33           65      1914      Powers 
12         127           39           77      1876      Lucas 
13         521          157          314      1952      Robinson 
14         607          183          366      1952      Robinson 
15        1279          386          770      1952      Robinson 
16        2203          664         1327      1952      Robinson 
17        2281          687         1373      1952      Robinson 
18        3217          969         1937      1957      Riesel 
19        4253         1281         2561      1961      Hurwitz 
20        4423         1332         2663      1961      Hurwitz 
21        9689         2917         5834      1963      Gillies 
22        9941         2993         5985      1963      Gillies 
23        11213        3376         6751      1963      Gillies 
24        19937        6002        12003      1971      Tuckerman
25        21701        6533        13066      1978      Noll & Nickel 
26        23209        6987        13973      1979      Noll
27        44497       13395        26790      1979      Nelson & Slowinski 
28        86243       25962        51924      1982      Slowinski 
29       110503       33265        66530      1988      Colquitt & Welsh
30       132049       39751        79502      1983      Slowinski 
31       216091       65050       130100      1985      Slowinski 
32       756839      227832       455663      1992      Slowinski & Gage 
33       859433      258716       517430      1994      Slowinski & Gage 
34      1257787      378632       757263      1996      Slowinski & Gage 
35      1398269      420921       841842      1996      GIMPS
36      2976221      895932      1791864      1997      GIMPS
37      3021377      909526      1819050      1998      GIMPS 
38     6972593     2098960      4197919      1999      GIMPS
39    13466917   4053947
40   20996011  6320431
41    24036583  7235734
42     25964951   7816230                           2005

   是不是有无穷多个梅森素数呢?数学家们目前还无法回答这个问
题。

三、寻找更大的素数

   为什么要寻找梅森素数?为什么要打破已知最大素数的纪录?这
有什么用处呢?

   如果你所说的用处是指能够直接创造物质财富,那么我不得不告
诉你——梅森素数没有什么用处,多知道一个非常大的素数似乎也没
什么用处。即使我们知道了一个无比巨大的梅森素数,也不会使我们
的钱包增加一分钱(嗨等一等!如果你只对钱感兴趣的话,也请不要
立刻撇下我的文章。我其实是说,我上面说的话要排除我在这篇文章
题目中提到的那十万美元的奖金——你的钱包也许会因此鼓起来的。
所以请耐心一点)。

   但是人类并不只需要物质财富。博物馆里的钻石有什么用场呢?
为什么人类要收集它们?因为它们美丽而稀少。作为人类智慧的结晶,
素数、梅森素数和与它密切相关的完美数是非常美丽的。它们的定义
简单,却又如此神秘莫测,象欧几里德、笛卡尔、费尔马、莱布尼兹、
欧拉这样的伟大数学家都因为它们的美丽而对它作过大量研究;大家
也看到,两千多年来,经过无数代人的辛勤工作,我们一共只收集到
38个梅森素数,它们是非常稀少的。对于数学家来说,搜集素数、梅
森素数和完美数是和收集钻石一样富有乐趣的事情。

   人类还需要荣耀——也许更胜于财富。在体育运动中,能够跑得
更快一点,跳得更高一点,难道真的有实际物质方面的用途吗?不,
我们喜欢接受挑战,我们希望能赢。打破一个体育世界记录,攀登珠
穆朗玛峰,单身驾船横穿太平洋……,那是对人类体能极限的挑战;而寻
找更大的素数,则是一项对人类智慧的挑战。当我们完成了一项前所
未有的任务时,我们总会感到无比骄傲。1963年,当第23个梅森素数
被找到时,发现它的美国伊利诺斯大学数学系是如此地骄傲,以致于
把所有从系里发出的信件都敲上了“2^11213-1是个素数”的邮戳。

   在欧拉证明M_31是素数以后,下一个最大素数的记录由兰德里
(Landry)于1867年获得:M_59/179951=3203431780337。这不是一个梅
森素数。这个记录保持了九年。

   1876年爱德华·卢卡斯使用了一个比费尔马和欧拉的方法更先进
的手段,证明了M_127是一个素数。这个记录保持了七十五年。直到费
里叶(Ferrier)于1951年使用一部手摇计算机证明了(2^148+1)/17是一
个素数,它有41位数。

   借助手摇计算机的方法要算作手工计算方法还是要算做计算机方
法,大概是可以探讨的问题。不过技术的发展一下子把这种争论变得
毫无必要。值得指出的是,在人类寻找大素数的旅途中,数学理论的
改善要远远比具有强大坚韧的计算能力重要得多。卢卡斯的方法在
1930年被勒梅(Lehmer)简化后,卢卡斯-勒梅测试成为现在寻找梅森素
数的标准方法。

(卢卡斯-勒梅测试:对于所有大于1的奇数p,M_p是素数当且仅当M_p
整除S(p-1),其中S(n)由S(n+1)=S(n)^2-2,S(1)=4递归定义。
4  14 194  37634  1416317954 2005956546822746114
这个测
试尤其适合于计算机运算,因为除以M_p=2^p-1的运算在二进制下可以
简单地用计算机特别擅长的移位和加法操作来实现。判断一个梅森数
是素数的方法比判断一个差不多大小的其他类型数是素数的方法要简
单得多,所以在寻找最大素数的过程中,大部分纪录都是梅森素数。)

   在1951年米勒和维勒(Miller & Wheeler)借助于EDSAC计算机(这
种计算机还不如我们现在使用的一般计算器,它只有5K的内存)发现
了长达79位的素数180(M_127)^2+1。这个记录还是没能保持多久。次
年罗宾逊应用SWAC计算机,在1952年初发现了第13和第14号梅森素数:
M_521和M_607,后面连续三个梅森素数也在同一年被陆续发现:M_1279,
M_2203和M_2281。

   在那以后的年代里,为了打破巨大素数纪录而使用的计算机越来
越强大,其中有著名的IBM360型计算机,和超级计算机Cray系列。大
家可以参看上面的梅森素数表来了解这个竞赛过程。在此其间只有一
次一个不是梅森素数的素数坐上过“已知最大素数”的宝座,它是
39158*2^216193-1,在1989年被发现。1996年发现的M_1257787是迄今
为止最后一个由超级计算机发现的梅森素数,数学家使用了Cray T94。

   然后,GIMPS的时代到来了。

四、GIMPS——互联网梅森素数大搜索

   1995年程序设计师乔治·沃特曼(George Woltman)开始收集整理
有关梅森素数计算的数据。他编制了一个梅森素数寻找程序并把它放
在网页上供数学爱好者免费使用。这就是“互联网梅森素数大搜索”
计划(GIMPS,the Great Internet Mersenne Prime Search)。在这个
计划中,十几位数学专家和几千名数学爱好者正在寻找下一个最大的
梅森素数,并且检查以前梅森素数纪录之间未被探索的空隙。比如上
面的梅森素数表中,最后那个素数的序号是未知的,我们不知道第37
号梅森素数和它之间是否还存在着其他未被发现的梅森素数。

1997年斯科特·库尔沃斯基(Scott Kurowski)和其他人建立了“素数
网”(PrimeNet),使分配搜索区间和向GIMPS发送报告自动化。现在只
要你去GIMPS的主页下载那个免费程序,你就可以立刻参加GIMPS计划
搜寻梅森素数。几乎所有的常用计算机平台都有可用的版本。程序以
最低的优先度在你的计算机上运行,所以对你平时正常地使用计算机
几乎没有影响。程序也可以随时被停止,下一次启动时它将从停止的
地方继续进行计算。

   从1996年到1998年,GIMPS计划发现了三个梅森素数:M_1398269、
M_2976221和M_3021377,都是使用奔腾型计算机得到的结果。

   1999年3月,在互联网上活动的一个协会“电子边界基金”(EFF,
Electronic Frontier Foundation)宣布了由一位匿名者资助的为寻找
巨大素数而设立的奖金。它规定向第一个找到超过一百万位的素数的
个人或机构颁发五万美元的奖金,这就是我们最一开始说到的哈吉拉
特瓦拉得到的奖金。后面的奖金依次为:超过一千万位,十万美元;
超过一亿位,十五万美元;超过十亿位,二十五万美元。

   搜寻结果的验证和奖金的颁发是非常严格的。比如说,得到的结
果必须是显式的——你不能宣称你的结果是一个有一百个方程组成的
方程组的解,却不把它解出来。结果必须由另一台计算机独立验证。
所有这些规则都在EFF网站上进行了解释。

   应该指出的是,通过参加GIMPS计划来获得奖金的希望是相当小的。
哈吉拉特瓦拉使用的计算机是当时21000台计算机中的一台。每一个参
与者都在验证分配给他的不同梅森数,当然其中绝大多数都不是素数
——他只有大约三万分之一的可能性碰到一个素数。

   下一个十万美元的奖金将被颁发给第一个找到超过一千万位的素
数的个人或机构。这一次的计算量将大约相当于上一次的125倍。现在
GIMPS得到的计算能力为每秒7000亿次浮点运算,和一台当今最先进的
超级矢量计算机,比如Cray T932的运行能力相当。但是如果GIMPS要
使用这样的超级计算机,一天就需要支付大约二十万美元。而现在他
们需要的费用,仅仅是支持网站运行的费用,和总共几十万美元的
奖金罢了。

五、网上分布式计算计划

   GIMPS只不过是互联网上众多的分布式计算计划中的一个,
GIMPS主页上就有这些计划的介绍。

   分布式计算是一门计算机学科,它研究如何把一个需要非常巨大
的计算能力才能解决的问题分成许多小的部分,然后把这些部分分配
给许多计算机进行处理,最后把这些计算结果综合起来得到最终的结
果。有时侯计算量是如此之大,需要全世界成千上万甚至更多台计算
机一起工作,才能在合乎情理的时间内得到结果。GIMPS计划就是在进
行这样的分布式计算。 

   但它并不是最著名的分布式计算计划。致力于寻找宇宙中智慧生
命的“搜寻地外文明计划”(SETI计划)中的SETI@HOME工程,已在全世
界招募了290万名(!)志愿者,利用屏幕保护程序来处理射电望远镜接
受到的大量的宇宙间传来的无线电信号。如果你参加这个计划,也许
有一天会在你的计算机上破译出外星人发来的问候呢。

   你也可以用你的计算机空余的计算能力为人类征服癌症作出贡献。
英国科学家设计了类似SETI@HOME工程的分布式计算屏保,它从有关网
站下载数据,分析化学物质分子的抗癌性能,然后将分析结果通过互
联网传回给研究人员,作为研制新型抗癌药物的参考。这项工程将于
2001年4月3日在美国加利福尼亚州正式启动。

   计算机硬件的更新令人目不暇接,上半年买的最新式的个人电脑,
在下半年就变成了大路货。三四年前的CPU,现在变得一钱不值——
也许不能这么说,你根本就买不到它们了——市面上最便宜的CPU也
要比它们强大得多。而一台普通的家用计算机连续运转五年也是没有
问题的。所以,对待计算机的最经济的态度就是:让它运转。

   而人类还有那么多的东西需要计算,还有那么多的问题需要找到
回答,还有那么多的难关需要克服。我们需要越来越巨大的计算能力,
我们也拥有这样的计算能力,只是太多太多被白白地闲置浪费掉了。
互联网已经使大规模的分布式计算计划成为可能。现在,我们唯一需
要的,就是这个网每一个结点上计算机用户的意愿和信心了。

时间: 2024-09-20 00:11:27

梅森素数--美丽的贝壳的相关文章

十万美元的悬赏 互联网梅森素数大搜索_应用技巧

一.价值五万美元的素数  2000年4月6日,住在美国密歇根州普利茅茨的那扬·哈吉拉特瓦拉(Nayan Hajratwala)先生得到了一笔五万美元的数学奖金,因为他找到了迄今为止已知的最大素数,这是一个梅森素数: 26972593-1. 这也是我们知道的第一个位数超过一百万位的素数.精确地讲,如果把这个素数写成我们熟悉的十进制形式的话,它共有两百零九万八千九百六十位数字,如果把它以这个形式写下来,大约需要150到200篇本文的篇幅. 可是哈吉拉特瓦拉先生并不是一个数学家,他甚至很可能对寻找素数

GIMPS 项目宣布发现已知最大素数

互联网梅森素数大搜索(GIMPS)项目宣布发现第49个梅森素数和已知最大的素数:2^74,207,281-1.发现者又是GIMPS志愿者.中密苏里大学数学计算机科学系的Curtis Cooper教授.最新发现的素数共有22,338,618位, 比之前的记录保持者长500万位.大素数对于加密具有重要价值,但最新发现的素数太长了,因此目前不具有实用价值.GIMPS是一个分布式计算项目,至今 已有20年历史,它利用志愿者的空闲CPU创建了一个遍布全球的超级计算机,它的 prime95软件此前发现了英特

如何测试电脑超频是否稳定?

1.计算高位圆周率的时间有多长 Superπ superπ是由日本东京大学金田研究室开发的一款用来计算圆周率的软件,设计者的初衷当初只是在HITAC S-3800/480超级计算机上使用,由于在计算π值时,考验到了CPU的多方面计算能力,因此后来被日本的超频爱好者移植到PC上使用,借助Superπ来测试超频后的性能,后来慢慢传入我国,许多硬件实验室也使用这款软件作为测试CPU稳定性的依据. 运行程序时,点击"运行计算"按钮,此时会弹出一个对话框请你选择要运行的圆周率位数,一般采用104

云计算的未来畅想(二)

01本文摘要 ------ 人类未来最大的矛盾,是日益增长的数据处理与有限算力之间的矛盾,云计算的本质不就是计算资源与存储资源,如果最终用户那里也有可以使用的计算资源或者存储资源,为什么不能利用起来呢? 02作者简介 ------ 拖雷  袋鼠云CEO 国内最早一批Oracle ACE Director,前阿里巴巴研究员,先后担任嗨淘.无线事业部.数据事业部资深总监,生活服务事业部总经理,以及阿里云事业群总裁助理.2015年离职后,做了半年独立天使投资人,2015年底创办袋鼠云. 03正文 --

网格技术 网络时代新型应用领域

网格计算即http://www.aliyun.com/zixun/aggregation/13452.html">分布式计算 什么是分布式计算?所谓分布式计算是一门计算机科学,它研究如何把一个需要非常巨大的计算能力才能解决的问题分成许多小的部分,然后把这些部分分配给许多计算机进行处理,最后把这些计算结果综合起来得到最终的结果. 最近的分布式计算项目已经被用于使用世界各地成千上万位志愿者的计算机的闲置计算能力,通过因特网,您可以分析来自外太空的电讯号,寻找隐蔽的黑洞,并探索可能存在的外星智慧

走向主流的网格计算技术应该于服务器

&http://www.aliyun.com/zixun/aggregation/37954.html">nbsp;   分布式计算是利用互联网上的计算机的 CPU 的闲置处理能力来解决大型计算问题的一种计算科学.下面,我们看看它是怎么工作的: 首先, 要发现一个需要非常巨大的计算能力才能解决的问题.这类问题一般是跨学科的.极富挑战性的.人类急待解决的科研课题.其中较为著名的是: 1.解决较为复杂的数学问题,例如:GIMPS(寻找最大的梅森素数). 2.研究寻找最为安全的密码系统,

创造力的岁月里

 http://user.qzone.qq.com/648952046/blog/1258076026 发表于2008年10月19日       两年过去了 然而这两年我却没做出什么可以拿得出手的事来说我这两年的意义 只要有坚定的意志,不屈不饶的意志力和执着顽强的忍耐力和一颗伟大而崇高的心灵,这个世界终将听到你的声音. 记得高中时,留了好多奇怪的问题: 1.4个原子编号1,2,3,4,经过衰变之后,根据半衰期统计规律,应该是衰变掉一半,剩一半,比如衰变掉的一半是1,2,而剩一半是3,4,那么1,

牛逼顿

这斜乜的眼神,微微上扬的左嘴角,坚挺的鹰钩鼻,笑而不语-- 确实透着股秒杀春哥,视女人为粪土的轰轰牛气.   作者:singularitys   3月28日是牛顿的忌日,但是知道的人很少,我们毕竟更关心沈殿霞和张国荣.其实牛顿老师在科学圈里曾经很有权势,被女王封了爵位成了贵族,人称牛爵爷,官至皇家造币局局长兼皇家学会会长.如果阿尔伯特没有辞了以色列总统的话和他有一拼. 说他有权势并不仅是官大,主要是贡献大.如果17世纪就有诺贝尔奖的话,牛顿老师至少能连续垄断4届物理学奖(分光计.力学体系的构建.

Intel官方确认:Skylake处理器存死机Bug

Intel官方确认:Skylake处理器存死机Bug 虽然PC走势渐弱,Intel也把精力开始投入到了虚拟现实.物联网等领域,但谈及2015最激动的事,Intel CEO前不久亮相CES时还是把荣耀颁给了Skylake处理器. 作为第六代酷睿,Skylake惊喜并没有达到让所有人眼前一亮,性能就不谈了,上市后还被曝出像是i3 TDP弄错.处理器太薄遭散热器压坏等囧事,而现在又多了一桩.据外媒报道,Intel在今天承认,Skylake处理器中存在一项bug,可能会导致系统死机,目前正在寻找解决办法