有关梅森素数的研究.pdfVIP

第 26卷第 1期 佛山科学技术学院学报(自然科学版) Vo1．26 NO．1 2008年 1月 Journal of Foshan University(Natural Science Edition) Jan．2008 文章编号：1008—01 71(2008)01—0021—03 有关梅森素数的研究 王 洋 (成都理工大学 信息管理学院。四川 成都 610059) 摘要：介绍了梅森素数相关的定义、定理及算法。讨论了三个有待解决的数论问题。并将现已被发现的44个梅 森素数列举出来。 关键词：素数；梅森素数；数论问题 中图分类号：0156 文献标识码：A 在没有终点的自然数序列中，按性质分可将自然数分为素数、合数和“1”。其中素数引发了历史上许 许多多的数学家的浓厚兴趣，而梅森素数作为一类特殊的素数，更吸引着无数的数学家及数学爱好者对 其进行研究。目前人类只找到了44个梅森素数，相对于众多的素数而言可谓是凤毛麟角。目前，世界上 有150多个国家和地区近15万人参加了一个名为“因特网梅森素数大搜索”(GIMPS)的国际合作项 目，并动用了超过3O万台计算机联网来进行大规模的网络计算，以探寻新的梅森素数。本文主要从梅 森素数的有关概念、定理及算法设计；由梅森素数而引发的尚未解决的问题；梅森素数表三方面来介绍 梅森素数，从而使更多的人认识与了解梅森素数，并参与梅森素数的搜索工作。 l 梅森素数的有关概念、定理及算法设计 1．1 有关概念、定理 定义1 若2 一1为素数( 为素数)，则称2 一1为梅森素数，简记为M 。 如当户一2，3，5，7，13，19，31，61，89，107，127，521，607时，M 为梅森素数。 定义2 如果一个自然数的所有因子(除了它自身)之和等于它自身，则称这个自然数为完全数。 其中6就是最小的完全数。因6—1十2十3。在这里可以把6分裂为2*3。即6—2 *(2。一1)，而 2。一1是形式为2 一1的素数，即梅森素数。 1730年，欧拉(I uler)证明了定理“每一个偶完全数都是型如2 ·(2 一1)的自然数，其中 和2 一 1都是素数”。这是欧氏定理的逆定理。根据欧几里得与欧拉这两个互逆定理可知，找到一个型如2 一 1的素数，即找到一个偶完全数，而型如2一一1的素数恰为梅森素数。至此人们终于发现：偶完全数与梅 森素数是一一对应的。 定理1 “是偶完全数的充分必要条件是“一2 *(2 一1)且2 一1是素数。 定理2 若2 一1为素数(户1)，则户是素数。 证明 假设 不是素数，即 一是 ，则(2 一1)J(2 一1)即(2 一1)J(2 一1)，故2 一1亦为非素数，与 条件矛盾。 故命题得证。 定理3 若p=3(mod 4)为素数，则(2p4-1)lMp当且仅当2p+1为素数。 收稿日期：2007一l1一O7 作者简介：王 洋(1985一)．男．安徽淮南人．成都理工大学硕士研究生。 22 佛山科学技术学院学报(自然科学版) 第26卷 证明 令 一2p+1， 1 。由于2 ≠1(mod )，可知2 ≠1(mod )，2。 一1一(2 +1)Mp三O(mod )，由LUCaS检测可知 是素数。反之，若q一2p+1为素数。 由于q三7(mod 8)，(2 l q)一1。所以有整数m使得2=m。(mod q)，于是2p~2 ” --m 三1 (mod q)，即qlM 。 1．2 梅森素数判定法的算法设计 193 1)年，美国数学家雷默改进了卢卡斯的工作，给出一个针对Mp的新的素性测试方法，即卢卡斯一 雷默方法：对于所有大于1的奇数P， 是素数当且仅当M 整除 (户一1)，其中s(n)由s( +1)一 ( )。 一 2，s(1)=4递归定义。这个方法尤其适合于计算机运算，因为除以M 一2 一1的运算在二进制下可以 简单地用计算机特别擅长的移位和加法操作来实现。 以下是用c语言实现的可实际使用的卢卡斯一雷默判定算法： I UCRS—I emer(in