数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
查看: 6441|回复: 10

【转载】梅森数的一些性质

[复制链接]
发表于 2022-1-6 09:15 | 显示全部楼层 |阅读模式
梅森数的一些性质
阿谦  数论爱好者
知乎网  https://zhuanlan.zhihu.com/p/82356370

形如Mn=2^n-1的数称为梅森数,若Mp=2^p-1为素数,则称之为梅森素数。

性质1:若n为合数,则Mn为合数;若MN为素数,则p为素数。
证明:
由2^a-1, 2^b-1|2^ab-1知
若n=ab (a,b≠1)为合数,则Ma, Mb|Mn,所以Mn为合数
逆否命题为:若Mn为素数,则p为素数

性质2:n为大于1的奇数时,Mn的所有因子形式为8k-1或8k+1
证明:
令p为Mn的任一素因子(可能是Mn本身),则
(2/p)=([2+2*(2^n-1)]/p)=(2^(n+1)/p)=1 (括号为勒让德符号)
所以2为p的二次剩余,所有p≡±1(mod 8)

性质3:p为奇素数时,Mn的所有素因子可表示为2kp+1的形式
证明:
令q为Mp的任一素因子
由费马小定理得
q|2^(q-1)-1
∴q|(2^(q-1)-1, 2^p-1)=2^(q-1, p)-1
若(q-1, p)=1,则q|2^(q-1, p)-1=1,矛盾;
所以(q-1, p)=p
所以p|(q-1),
又因为q-1为偶数,所以存在整数k,使q=2kp+1
例如2^11-1=2047=23*89=(2*11+1)*(8*11+1)

性质4:若素数p≡3 (mod 4)
则2p+1是素数的充要条件是2^p≡1 (mod 2p+1)
证明:
必要性:
由p≡3 (mod 4)和2p+1≡7 (mod 8)
∴(2/(2p+1))=1 (括号为勒让德符号)
由欧拉判别条件得2^[((2p+1)-1)/2]≡1 (mod 2p+1)
即2^p≡1 (mod 2p+1)
由此我们可以推出若2p+1是形如8k-1素数,
则2p+1|Mp,
所以Mp是合数(p=3除外,因为恰好2p+1=Mp)
例如:23是形如8k-1的素数,所以23|M11,所以M11是合数
充分性:
首先易得2p+1没有大于等于p的素因子。
假设2p+1有素因子q
有q|2p+1|2^p-1
由费马小定理有q|2^(q-1)-1
所以有p|q-1,所以p<q,即2p+1有大于p的素因子,矛盾
所以2p+1是素数

性质5:正整数n 是一个偶完全数当且仅当n=2^(p-1)*Mp ,Mp 是梅森素数
证明:
充分性:
若 n=2^(p-1)*Mp ,
则 σ(n)=σ(2^(p-1)*Mp)=σ(2^(p-1))*σ(Mp)=(2^p-1)*2^p=2n
所以n 是偶完全数 ( σ(n) 为因子和函数,且为积性函数)
必要性:
设 n 是一个偶完全数。令 n=2^s*t ,其中s, t 是正整数且 t 是奇数
所以有σ(n)=σ(2^s*t)=σ(2^s)*σ(t)=(2^(s+1)-1) *σ(t)      (1)
又因为 n 是完全数,所以σ(n)=2n=1^(s+1)*t     (2)
由(1)(2)式得 (2^(s+1)-1)* σ(t)=2^(s+1)*t     (3)
又因为 (2^(s+1)-1, 2^(s+1))=1 ,所以有2^(s+1)|σ(t)
令 σ(t)=2^(s+1)*q ,代入(3)式,得
(2^(s+1)-1)*2^(s+1)*t
即 (2^(s+1)-1)*q=t     (4)
所以 q|t 且q≠t
(4)式等号两边加上 q ,得
T+q=(2^(s+1)-1)*q=2^(s+1)*q=σ(t)     (5)
下面我们证明 q=1
若 q≠1 ,说明 t 至少存在三个不同的正因子 1,q,t
所以σ(t) ≥1+q+t ,和(5)式矛盾,所以 q=1
由(4)式得 t=2^(s+1)-1
由(5)式得 σ(t)=1+t ,即 t 只有 1 和 t 两个正因子,所以 t 是素数
所以 n=2^s*(2^(s+1)-1) ,其中 2^(s-1)-1 是素数
由此可以得出,只要找到哪些数为梅森素数,就可以找出相应的偶完全数

接下来我们讨论如何高效地判别 Mp 是不是一个素数
Lucas-Lehmer判别法:
令 r1=4, rk=r<k-1>^2-2 (mod Mp)   (r后的1,k,<k-1>都是r的下标)
则 Mp 是素数当且仅当 r<p-1>≡0 (mod Mp)   (r后的<p-1>是r的下标)
Lucas-Lehmer判别法的证明可以见:
[1] H.W.Lenstra, Jr. ,"Primality testing," Studieweek Getaltheorie en Computers, 1-5 September 1980, Stichting Mathematisch Centrum, Amersterdam, Holland.
[2] W. Sierpinski, A Selection of Problems in the Theory of Numbers, Pergamon Press, New York, 1964.
2020.1.15更新

性质6:若 n>1 ,则 Mn 不可能为一个正整数的高次幂
证明:
令 m^k=2^n-1 ,由 n>1 知 m 为大于 1 的奇数
分两种情况讨论
若 k 为大于 0 的偶数,则 m^k-1=2^n-2
即 [m^(k/2)-1]*[m^(k/2)+1]=2*[2^(n-1)-1]
等式左边为 4 的倍数,而等式右边为 2 乘一个奇数,矛盾
若 k 为大于 1 的奇数,则 m^k+1=2^n
即 (m+1)*[m^(k-1)-m^(k-2)+…-m+1]=2^n
而括号里面有 k 项,所以m^(k-1)-m^(k-2)+…-m+1为奇数,矛盾
综上所述,若 n>1 ,则 Mn 不可能为一个正整数的高次幂
编辑于 2020-01-15 00:53
 楼主| 发表于 2022-1-6 09:20 | 显示全部楼层
梅森数的一个性质给出的启示
百度贴吧  哥德巴赫猜想吧
https://tieba.baidu.com/p/6050962709?red_tag=3075948998

liuluojieys
潘承洞潘承彪兄弟所著(初等数论),对梅森数的性质有若干描述,除了以前数学家们给出的若干相关定理以外,还证明了下面两条性质:(1)2^p-1的因子全形如2kp+1,(2)如果素数形如p=4n+3,则2p+1整除2^p-1的充要条件是,2p+1是素数。
不难看出,性质(2)还可以描述为:素数形如p=4n+3,若q=2p+1也是素数,则2^p-1的一个素因子是2kp+1,且k=1。
这就给了我们一个启示,存在两个等差数列,通项公式是p=4n+3,q=2p+1=8n+7,对任意给定的自然数n,如果(数对)(p,q)同时为素数,则2^p-1,必是一个梅森合数。这样的(素数对)(p,q)很容易对两个等差数列双筛得到。
通过研究,发现:对于任意的奇素数p,梅森数2^p-1,(1)有且仅有一个素因子q=2kp+1,k是奇数。
(2)若p=4n+3,则k=4l+1;若p=4n+1,则k=4l+3。(3)可表(2kp+1)(8mp+1)=2^p-1。称为梅森方程。
(4)2^p-1至多存在三个素因子。
于是猜想,当取奇数k>1时,素数q=2kp+1,在什么条件下对应的梅森数2^p-1一定是梅森合数呢?
得到了这个条件,追寻梅森素数就简单了。

雁荡山是我
根据素数作为因子数在2^N-1数列中的分布规律而编选梅森素数。一 在2^N-1平方根以内没有奇素数的,这个指数的数就是梅森素数,后面的指数能被这个梅森素数整除的都有这个梅森素数的因子,周期重复的都编上这个因子数。二 在2^N-1数列中N是合数的,2^N-1是由前面素因子和新增的因子数组成,只要除去前面的因子数就能得到新的因子数,新的因子数也以这个指数为周期重复出现,并在后面重复的数编上这个因子数。三 指数是素数的,剩下的因子素数根据减1能被指数整除来选合,选合后并在后面的周期重复的数编上这些因子数;剩下的因子素数在2^N-1平方根以内选完后,还有梅森数空着,这个梅森数就是梅森素数,并在后面周期重复的数编上为些因子数。

liuluojieys
三 指数是素数的,剩下的因子素数根据减1能被指数整除来选合,选合后并在后面的周期重复的数编上这些因子数;剩下的因子素数在2^N-1平方根以内选完后,还有梅森数空着,这个梅森数就是梅森素数,并在后面周期重复的数编上为些因子数。
这一条,请举例说明好吗?

雁荡山是我
指数11数值是2047,2047的平方根内还有43,41,37,29,23,19,13等奇素数(只要往后多编一些,留下的数就更少了),只要将这几个数都减去1,那一个数被11整除,只有23-1能被11整除,2047/23=89,23和89是它的因子数,两个因子数都是新增加的,因为它是梅森合数,凡是指数能被11整除的都有23和89的因子数,后面以11为周期编上因子数23和89;

liuluojieys
事实上,细心研究会发现:2^p-1的素因子
(1)p形如4n-1,存在于等差数列 2p+1,10p+1,18p+1,......,2(4k+1)p+1;之中。
(2)p形如4n+1,存在于等差数列6p+1,14p+1,22p+1,......,2(4k+3)p+1;之中。
式中k取自然数。
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-6 09:25 | 显示全部楼层
本帖最后由 yangchuanju 于 2022-1-6 09:30 编辑

梅森素数指数表
A000043
Mersenne exponents: primes p such that 2^p - 1 is prime. Then 2^p - 1 is called a Mersenne prime.

2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011, 24036583, 25964951, 30402457, 32582657, 37156667, 42643801, 43112609
48-51#素数:57885161,74207281,77232917,82589933
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-6 09:32 | 显示全部楼层
梅森素数表(1-18#)
A000668
Mersenne primes (primes of the form 2^n - 1).

1 3
2 7
3 31
4 127
5 8191
6 131071
7 524287
8 2147483647
9 2305843009213693951
10 618970019642690137449562111
11 162259276829213363391578010288127
12 170141183460469231731687303715884105727
13 6864797660130609714981900799081393217269435300143305409394463459185543183397656052122559640661454554977296311391480858037121987999716643812574028291115057151
14 531137992816767098689588206552468627329593117727031923199444138200403559860852242739162502265229285668889329486246501015346579337652707239409519978766587351943831270835393219031728127
15 10407932194664399081925240327364085538615262247266704805319112350403608059673360298012239441732324184842421613954281007791383566248323464908139906605677320762924129509389220345773183349661583550472959420547689811211693677147548478866962501384438260291732348885311160828538416585028255604666224831890918801847068222203140521026698435488732958028878050869736186900714720710555703168729087
16 1475979915214180235084898622737381736312066145333169775147771216478570297878078949377407337049389289382748507531496480477281264838760259191814463365330269540496961201113430156902396093989090226259326935025281409614983499388222831448598601834318536230923772641390209490231836446899608210795482963763094236630945410832793769905399982457186322944729636418890623372171723742105636440368218459649632948538696905872650486914434637457507280441823676813517852099348660847172579408422316678097670224011990280170474894487426924742108823536808485072502240519452587542875349976558572670229633962575212637477897785501552646522609988869914013540483809865681250419497686697771007
17 446087557183758429571151706402101809886208632412859901111991219963404685792820473369112545269003989026153245931124316702395758705693679364790903497461147071065254193353938124978226307947312410798874869040070279328428810311754844108094878252494866760969586998128982645877596028979171536962503068429617331702184750324583009171832104916050157628886606372145501702225925125224076829605427173573964812995250569412480720738476855293681666712844831190877620606786663862190240118570736831901886479225810414714078935386562497968178729127629594924411960961386713946279899275006954917139758796061223803393537381034666494402951052059047968693255388647930440925104186817009640171764133172418132836351
18 259117086013202627776246767922441530941818887553125427303974923161874019266586362086201209516800483406550695241733194177441689509238807017410377709597512042313066624082916353517952311186154862265604547691127595848775610568757931191017711408826252153849035830401185072116424747461823031471398340229288074545677907941037288235820705892351068433882986888616658650280927692080339605869308790500409503709875902119018371991620994002568935113136548829739112656797303241986517250116412703509705427773477972349821676443446668383119322540099648994051790241624056519054483690809616061625743042361721863339415852426431208737266591962061753535748892894599629195183082621860853400937932839420261866586142503251450773096274235376822938649407127700846077124211823080804139298087057504713825264571448379371125032081826126566649084251699453951887789613650248405739378594599444335231188280123660406262468609212150349937584782292237144339628858485938215738821232393687046160677362909315071
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-6 09:42 | 显示全部楼层
本帖最后由 yangchuanju 于 2022-1-18 08:01 编辑

梅森数表(2^2-1至2^163-1)
n        digits        number
2        1        3=3
3        1        7=7
5        2        31=31
7        3        127=127
11        4        2047=23·89
13        4        8191=8191
17        6        131071=131071
19        6        524287=524287
23        7        8388607=47·178481
29        9        536870911=233·1103·2089
31        10        2147483647<10>=2147483647<10>
37        12        137438953471<12>=223·616318177
41        13        2199023255551<13>=13367·164511353
43        13        8796093022207<13>=431·9719·2099863
47        15        140737488355327<15>=2351·4513·13264529
53        16        9007199254740991<16>=6361·69431·20394401
59        18        576460752303423487<18>=179951·3203431780337<13>
61        19        2305843009213693951<19>=2305843009213693951<19>
67        21        147573952589676412927<21>=193707721·761838257287<12>
71        22        2361183241434822606847<22>=228479·48544121·212885833
79        24        604462909807314587353087<24>=2687·202029703·1113491139767<13>
83        25        9671406556917033397649407<25>=167·57912614113275649087721<23>
89        27        618970019642690137449562111<27>=618970019642690137449562111<27>
97        30        158456325028528675187087900671<30>=11447·13842607235828485645766393<26>
101        31        2535301200456458802993406410751<31>=7432339208719<13>·341117531003194129<18>
103        32        10141204801825835211973625643007<32>=2550183799<10>·3976656429941438590393<22>
107        33        162259276829213363391578010288127<33>=162259276829213363391578010288127<33>
109        33        649037107316853453566312041152511<33>=745988807·870035986098720987332873<24>
113        35        103845...191<35>=3391·23279·65993·1868569·1066818132868207<16>
127        39        170141...727<39>=170141183460469231731687303715884105727<39>
131        40        272225...647<40>=263·10350794431055162386718619237468234569<38>
137        42        174224...471<42>=32032215596496435569<20>·5439042183600204290159<22>
139        42        696898...887<42>=5625767248687<13>·123876132205208335762278423601<30>
149        45        713623...311<45>=86656268566282183151<20>·8235109336690846723986161<25>
151        46        285449...247<46>=18121·55871·165799·2332951·7289088383388253664437433<25>
157        48        182687...871<48>=852133201·60726444167<11>·1654058017289<13>·2134387368610417<16>
163        50        116920...807<50>=150287·704161·110211473·27669118297<11>·36230454570129675721<20>


回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-18 08:04 | 显示全部楼层
p        位数        分解式
167        51        2^167-1<51>=2349023·79638304766856507377778616296087448490695649<44>
173        53        2^173-1<53>=730753·1505447·70084436712553223<17>·155285743288572277679887<24>
179        54        2^179-1<54>=359·1433·1489459109360039866456940197095433721664951999121<49>
181        55        2^181-1<55>=43441·1164193·7648337·7923871097285295625344647665764672671<37>
191        58        2^191-1<58>=383·7068569257<10>·39940132241<11>·332584516519201<15>·87274497124602996457<20>
193        59        2^193-1<59>=13821503·61654440233248340616559<23>·14732265321145317331353282383<29>
197        60        2^197-1<60>=7487·2682880399...33<56>
199        60        2^199-1<60>=164504919713<12>·4884164093883941177660049098586324302977543600799<49>
211        64        2^211-1<64>=15193·60272956433838849161<20>·3593875704495823757388199894268773153439<40>
223        68        2^223-1<68>=18287·196687·1466449·2916841·1469495262398780123809<22>·596242599987116128415063<24>
227        69        2^227-1<69>=26986333437777017<17>·7992177738...31<52>
229        69        2^229-1<69>=1504073·20492753·59833457464970183<17>·467795120187583723534280000348743236593<39>
233        71        2^233-1<71>=1399·135607·622577·1168681298...31<57>
239        72        2^239-1<72>=479·1913·5737·176383·134000609·7110008717824458123105014279253754096863768062879<49>
241        73        2^241-1<73>=22000409·1606194743...39<66>
251        76        2^251-1<76>=503·54217·178230287214063289511<21>·61676882198695257501367<23>·12070396178249893039969681<26>
257        78        2^257-1<78>=535006138814359<15>·1155685395246619182673033<25>·374550598501810936581776630096313181393<39>
263        80        2^263-1<80>=23671·13572264529177<14>·120226360536848498024035943<27>·383725126655170964501315730676446647<36>
269        81        2^269-1<81>=13822297·6862598850...63<74>
271        82        2^271-1<82>=15242475217<11>·2489277578...91<72>
277        84        2^277-1<84>=1121297·31133636305610209482201109050392404721<38>·6955979459776540052280934851589652278783<40>
281        85        2^281-1<85>=80929·4800921529...19<80>
283        86        2^283-1<86>=9623·68492481833<11>·2357954301...73<71>
293        89        2^293-1<89>=40122362455616221971122353<26>·3966452270...47<63>
307        93        2^307-1<93>=14608903·85798519·23487583303<11>·78952752017<11>·1121774764...61<57>
311        94        2^311-1<94>=5344847·2647649373910205158468946067671<31>·2948036813...31<57>
313        95        2^313-1<95>=10960009·14787970697180273<17>·3857194764289141165278097<25>·26693012026551688286164949958620483258358551879<47>
317        96        2^317-1<96>=9511·587492521482839879<18>·4868122671322098041565641<25>·9815639231755686605031317440031161584572466128599<49>
331        100        2^331-1<100>=16937389168607<14>·865118802936559<15>·2985426249...19<72>
337        102        2^337-1<102>=18199·2806537·95763203297<11>·726584894969<12>·7877804732...69<68>
347        105        2^347-1<105>=14143189112952632419639<23>·2027034530...93<83>
349        106        2^349-1<106>=1779973928671<13>·34720396273212657799920861294559<32>·1855539364...99<62>
353        107        2^353-1<107>=931921·2927455476800301964116805545194017<34>·6725414756...63<67>
359        109        2^359-1<109>=719·855857·778165529·65877330027880703<17>·370906580744492785430299503112990447<36>·100361196281293745682520861860411315001<39>
367        111        2^367-1<111>=12479·51791041·78138581882953<14>·301311116540899114446723859201<30>·1975574008...81<56>
373        113        2^373-1<113>=25569151·7524403464...41<105>
379        115        2^379-1<115>=180818808679<12>·6809649408...53<103>
383        116        2^383-1<116>=1440847·7435494593<10>·503823044204581129045587727<27>·15174923558680812616818436353130417<35>·240522700235167893496900256599634325263<39>
389        118        2^389-1<118>=56478911·4765678679<10>·4684435266...19<100>
397        120        2^397-1<120>=2383·6353·50023·53993·202471·5877983·814132872808522587940886856743<30>·1234904213576000272542841146073<31>·6597485910270326519900042655193<31>
401        121        2^401-1<121>=856971565399<12>·2136958965524920285681<22>·594538100848945223169882301931953<33>·4743358775...93<55>
409        124        2^409-1<124>=4480666067023<13>·76025626689833<14>·3881196575...29<97>
419        127        2^419-1<127>=839·903780021613921<15>·5800422716722833271214743<25>·10287968884341772230096159036619433593<38>·29919490848598531825060153417921002916701815927<47>
421        127        2^421-1<127>=614002928307599<15>·8819779591...49<112>
431        130        2^431-1<130>=863·3449·36238481·76859369·558062249·4642152737<10>·142850312799017452169<21>·1807482391...57<70>
433        131        2^433-1<131>=22086765417396827057<20>·737748363812546584876297<24>·5028667832511715101284999<25>·2707040788...21<63>
439        133        2^439-1<133>=104110607·127321491658223<15>·122551752733003055543<21>·8738801468...69<90>
443        134        2^443-1<134>=887·207818990653657<15>·1232194392...73<117>
449        136        2^449-1<136>=1256303·6871197486841<13>·3578620616468306981503<22>·526385733768051189441947168504909159<36>·8939669339...41<59>
457        138        2^457-1<138>=150327409·2475539419...19<130>
461        139        2^461-1<139>=2767·358228856441770927<18>·7099353734763245383<19>·8461346092...33<99>
463        140        2^463-1<140>=11113·3407681·448747600991881<15>·239932071009857681156251129<27>·385606580062688087218266143<27>·1514826404...17<62>
467        141        2^467-1<141>=121606801·1148984537...47<58>·2727318928...41<76>
479        145        2^479-1<145>=33385343·6293443049<10>·6834814451...61<57>·1086918660...81<71>
487        147        2^487-1<147>=4871·8203321996...37<143>
491        148        2^491-1<148>=983·7707719·110097436327057<15>·6976447052525718623<19>·19970905118623195851890562673<29>·3717542676439779473786876643915388439<37>·14797326616665978116353515926860025681383<41>
499        151        2^499-1<151>=20959·1998447222711143545931606352264121<34>·3907550462...33<113>
503        152        2^503-1<152>=3213684984979279<16>·12158987054135300783<20>·1873030665061080894263<22>·3578015615...77<96>
509        154        2^509-1<154>=12619129·19089479845124902223<20>·647125715643884876759057<24>·1075116828...69<104>
521        157        2^521-1<157>=6864797660...51<157>
523        158        2^523-1<158>=1601887783...63<69>·1714176918...89<90>
541        163        2^541-1<163>=4312790327<10>·6115209994009<13>·77146448294831869472151022106713<32>·904106555381898703253733254107256591<36>·3913112911...79<73>
547        165        2^547-1<165>=5471·172720604638150729<18>·157552196128597543991<21>·3094370432...83<124>
557        168        2^557-1<168>=3343·21993703·4565508951414453792719<22>·22053236096920219270632521586535711<35>·6372521752...11<101>
563        170        2^563-1<170>=2815747080256641401887817<25>·13299213974872825932452460424235657<35>·8062466865...03<111>
569        172        2^569-1<172>=15854617·55470673·182602768015690099110572536951<30>·110582329556343704552404016904356632250881<42>·1088066029...41<87>
571        172        2^571-1<172>=5711·27409·6969336604...27<73>·7084851186...39<91>
577        174        2^577-1<174>=3463·132305774316967<15>·1079633141...51<157>
587        177        2^587-1<177>=554129·2926783·39483330766889<14>·73208283304744901303<20>·3728298863422039632638351<25>·2898153126...33<107>
593        179        2^593-1<179>=104369·3061144307110551703729<22>·4027016481761119763553952463965455549881<40>·2519699187...11<113>
599        181        2^599-1<181>=16659379034607403556537<23>·148296291984475077955727317447564721950969097<45>·8398047009...83<114>
601        181        2^601-1<181>=3607·64863527·6433907608...93<60>·5513226964...63<110>
607        183        2^607-1<183>=5311379928...27<183>
613        185        2^613-1<185>=44599476833089207<17>·332817722770314187794325446534549089<36>·2290082526...17<133>
617        186        2^617-1<186>=59233·68954123297<11>·1577519781...01<51>·8441284558...71<120>
619        187        2^619-1<187>=110183·710820995447<12>·1093786816...39<66>·2539567680...33<105>
631        190        2^631-1<190>=333628015107245479<18>·474640860193534882628078580680807822523991<42>·5627294314...23<131>
641        193        2^641-1<193>=35897·49999·1173835097<10>·2401258891949526685926151441<28>·7452763007...99<69>·2420161564...79<79>
643        194        2^643-1<194>=3189281·2253242905...53<71>·5079095912...99<117>
647        195        2^647-1<195>=303303806129303896428103<24>·11502383762931955619183430785983<32>·1673946143...23<141>
653        197        2^653-1<197>=78557207·289837969·16008527538753578495897<23>·3044938125828889184039273431<28>·3786707123184812723340717607<28>·8893122618...73<103>
659        199        2^659-1<199>=1319·11527429277532648241<20>·626564962613678012662146877852049<33>·2510867382...97<144>
661        199        2^661-1<199>=1330270433<10>·3370159489168519<16>·1009157848082361225065617<25>·118420287267066844820208926433723871<36>·1785876939...59<116>
673        203        2^673-1<203>=581163767·41283139633378645724930694480520226273492263<44>·6839676957...19<59>·2388252511...09<92>
677        204        2^677-1<204>=1943118631<10>·531132717139346021081<21>·978146583988637765536217<24>·5362511269...67<53>·1158330612...99<98>
683        206        2^683-1<206>=1367·4348364991...69<78>·6751379697...09<125>
691        209        2^691-1<209>=906642603313<12>·10488160032325844521<20>·2833637724427940664433391497<28>·16636201944470267267182958207<29>·2291886288...41<122>
701        212        2^701-1<212>=796337·2983457·28812503·1073825104511<13>·9983923992673<13>·15865578195367<14>·40686928318417<14>·2220668173...89<140>
709        214        2^709-1<214>=216868921·1391646221255097953<19>·8923613716...47<187>
719        217        2^719-1<217>=1439·772207·7375728433...69<51>·3364863383...51<157>
727        219        2^727-1<219>=1760629171...27<98>·4009949972...01<122>
733        221        2^733-1<221>=694653525743<12>·14399141148866077141941966959399761<35>·13948432601001273785699253868326270649<38>·3238594546...33<138>
739        223        2^739-1<223>=184603056517613273120809<24>·48050683584092004380805463790111<32>·2561503388...63<82>·1272721565...51<87>
743        224        2^743-1<224>=1487·1219280833<10>·14904366017<11>·118722715461092305629361<24>·4721525455401597740684262559<28>·4977047949106985392753791512048265888682003683833<49>·6137204207...03<102>
751        227        2^751-1<227>=2276402451...87<66>·6493500319...97<67>·8013068084...73<94>
757        228        2^757-1<228>=9815263·561595591·5722137022...93<79>·2403382164...59<134>
761        230        2^761-1<230>=4567·6089·738686421813192728921171408273447<33>·2802302669...07<69>·2107048624...13<121>
769        232        2^769-1<232>=1591805393<10>·6123566623856435977170641<25>·3185460664...47<198>
773        233        2^773-1<233>=6864241·9461521·1165626988664823792549971253972697<34>·101625278121114813541046756104385520302531369<45>·3102804258...07<52>·2081217395...81<91>
787        237        2^787-1<237>=9951597611230279<16>·96512008100928793<17>·1711247935...03<63>·4952443062...47<142>
797        240        2^797-1<240>=2006858753<10>·54573369937<11>·104757762864135516671<21>·747609928190821086986322983<27>·21023711323746974956423747989180911<35>·2896793131...57<55>·1595588224...01<85>
809        244        2^809-1<244>=4148386731...37<61>·8229761617...03<183>
811        245        2^811-1<245>=326023·1205618453...21<66>·2311981238...63<83>·1502740709...43<92>
821        248        2^821-1<248>=419273207·109840721427977<15>·286121480219517473<18>·2534833058...23<62>·4186644906...71<146>
823        248        2^823-1<248>=1460915248436556406607<22>·1534086200463688788034864584049<31>·1431846377...47<95>·1743071575...67<103>
827        249        2^827-1<249>=66161·1637241673<10>·1656489332706171370072190422631<31>·2550707823005230793212015514781257<34>·1339459522...99<80>·1459864713...23<93>
829        250        2^829-1<250>=72953·89654239955972104838901500177<29>·5473339866...31<216>
839        253        2^839-1<253>=26849·138561000316919<15>·377801626929390823<18>·113910681722635191781067775764311<33>·2289649272...09<185>
853        257        2^853-1<257>=2065711807<10>·513740645819473<15>·4727997228805279065233<22>·1844945010...89<80>·6488017349...13<132>
857        258        2^857-1<258>=6857·2893063747...27<60>·6625204725...91<75>·7311640476...79<120>
859        259        2^859-1<259>=7215601·13183242101527799620391251382520507315529510625017<50>·5860989871...37<62>·6894464928...03<141>
863        260        2^863-1<260>=8258911·169382737·175642891399<12>·2928599996...01<97>·8546802933...99<137>
877        265        2^877-1<265>=35081·1436527·1839699736969340949026772071177<31>·4478527345975241075738524755872079528958635991<46>·9138345796...23<71>·2655649011...53<107>
881        266        2^881-1<266>=26431·6099757718...21<261>
883        266        2^883-1<266>=8831·63577·258777491057348926546569104663<30>·4438638313...47<228>
887        268        2^887-1<268>=16173559·139714099995626040378049<24>·2382154297313396649423300881<28>·1958905756...79<73>·9785361819...03<137>
907        274        2^907-1<274>=1170031·3256645177<10>·131895016205024767<18>·2152830659...63<241>
911        275        2^911-1<275>=1823·26129303·201955048939840841121786425435009<33>·1799520676...07<232>
919        277        2^919-1<277>=33554520197234177<17>·1526510740...21<126>·8651972246...11<135>
929        280        2^929-1<280>=13007·340388595097<12>·9123887298...77<51>·1213545181...77<61>·9257164238...21<153>
937        283        2^937-1<283>=28111·2419437071<10>·3481316282...43<57>·4665472298...69<59>·1051662080...73<154>
941        284        2^941-1<284>=7529·2919837555...37<69>·8455317738...87<211>
947        286        2^947-1<286>=295130657·25749931927<11>·4621646208862937<16>·11892980076500863942962129100776937<35>·4247385134...73<104>·6705118430...89<113>
953        287        2^953-1<287>=343081·562070136841<12>·40496650491842199555367<23>·824501439026637805612489876810823<33>·2106041424...21<60>·5614631499...11<155>
967        292        2^967-1<292>=23209·549257·1438495818172960049<19>·73214830052905408299354903681341531907328727<44>·9291078206...73<219>
971        293        2^971-1<293>=2391710497...13<53>·8344823973...19<240>
977        295        2^977-1<295>=867577·1813313·2069655374719577273<19>·49858990580788843054012690078841<32>·2383040575...51<62>·3301809434...97<171>
983        296        2^983-1<296>=1808226257914551209964473260866417929207023<43>·4520983959...09<254>
991        299        2^991-1<299>=8218291649<10>·41473350001<11>·231620367206687<15>·6721885469...59<127>·3943730076...91<137>
997        301        2^997-1<301>=1675608165...07<57>·7993430605...53<244>
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-18 08:19 | 显示全部楼层
p        位数        分解式
1009        304        2^1009-1<304>=3454817·198582684439<12>·20649907789079<14>·21624641697047<14>·30850253615723594284324529<26>·1134327302421596486779379019599<31>·3280162939...47<89>·1560053963...37<116>
1013        305        2^1013-1<305>=6079·71486759597273<14>·4112091256...73<59>·4912076350...01<229>
1019        307        2^1019-1<307>=2039·75407·8243323067...89<53>·1140356877...41<76>·3886812454...31<171>
1021        308        2^1021-1<308>=40841·795808241·51731968308257<14>·1866013803867738316856548085894075729799351216503<49>·7162207820...01<232>
1031        311        2^1031-1<311>=2063·435502649·8588431101...31<74>·2982096305...51<225>
1033        311        2^1033-1<311>=196271·36913223·88910338973359434076776038119074340716281<41>·163771629900828769671363963472861890953631838273<48>·8724820123...79<210>
1039        313        2^1039-1<313>=5080711·5585366661...53<80>·2075818194...89<227>
1049        316        2^1049-1<316>=33569·459463·1059099980653317121<19>·3692661882...53<288>
1051        317        2^1051-1<317>=3575503·2157380128...49<63>·3050179060...99<69>·1025501468...99<180>
1061        320        2^1061-1<320>=4681722635...33<143>·5277396428...47<177>
1063        320        2^1063-1<320>=1485761479<10>·6651755792...33<311>
1069        322        2^1069-1<322>=17481727674576239<17>·37473613084215372416028665821312617743022228898297<50>·5557036167...21<70>·1737449627...77<187>
1087        328        2^1087-1<328>=10722169·2144921195591<13>·33990250708335239<17>·2105485879475562449<19>·421733524843190780249388261839156579767<39>·23864222009193938317456687011256092421563991<44>·2664797814...17<61>·3756253784...67<131>
1091        329        2^1091-1<329>=87281·551978359·791343585595817<15>·97535130364653673535422113471207017<35>·496474000537225566899310934393330825097522839<45>·1437013471...83<221>
1093        330        2^1093-1<330>=43721·111487·42565880168703357316051447<26>·4611633294...49<122>·1109058770...11<173>
1097        331        2^1097-1<331>=980719·4666639·22926719175799<14>·1734492765338661710171833<25>·36079584237927411511333522577<29>·205295973898620029528272537868263<33>·1259505554...43<220>
1103        333        2^1103-1<333>=2207·4126110275598714647074087<25>·1193279152...23<305>
1109        334        2^1109-1<334>=30963501968569<14>·85608965982066833903<20>·2461601921...19<146>·1065805713...67<156>
1117        337        2^1117-1<337>=53617·5209739697...09<78>·6373626673...07<254>
1123        339        2^1123-1<339>=777288435261989969<18>·4321749374...87<86>·3391902395...69<235>
1129        340        2^1129-1<340>=33871·833798113·2682863551...11<139>·9624487049...87<188>
1151        347        2^1151-1<347>=284278475807<12>·81344898763887260484533897<26>·658622107254745112979229913<27>·39146265717288735265327499826695585290297<41>·8311919431...47<103>·6171954834...79<139>
1153        348        2^1153-1<348>=267497·84755607199<11>·1012236096...93<89>·5331087677...29<243>
1163        351        2^1163-1<351>=848181715001<12>·337097300570078978047<21>·1042816042...87<73>·4201777709...63<246>
1171        353        2^1171-1<353>=153606920351<12>·1234867096...13<128>·1690806550...69<215>
1181        356        2^1181-1<356>=4742897·196834168670106625096599487<27>·199668911054466914791190422986407<33>·1808422353...87<73>·9742499217...01<218>
1187        358        2^1187-1<358>=256393·113603023·1761533117965696212626325049982897<34>·766302834531296176581826196660125454171041129<45>·1533278332...69<63>·3486555298...69<204>
1193        360        2^1193-1<360>=121687·8522732620...29<104>·1297065115...17<251>
1201        362        2^1201-1<362>=57649·1967239·8510287·2830858618432184648159211485423<31>·9303484663...59<61>·1354777128...99<253>
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-19 13:11 | 显示全部楼层
预测52号梅森素数的指数在91552704上下

梅森素数指数表:
1 2,  2 3,  3 5,  4 7,  5 13,  6 17,  7 19,  8 31,  9 61,  10 89,  11 107,  12 127,  13 521,  14 607,  15 1279,  16 2203,  17 2281,  18 3217,  19 4253,  20 4423,  21 9689,  22 9941,  23 11213,  24 19937,  25 21701,  26 23209,  27 44497,  28 86243,  29 110503,  30 132049,  31 216091,  32 756839,  33 859433,  34 1257787,  35 1398269,  36 2976221,  37 3021377,  38 6972593,  39 13466917,  40 20996011,  41 24036583,  42 25964951,  43 30402457,  44 32582657,  45 37156667,  46 42643801,  47 43112609,  48 57885161,  
48-51#素数:57885161,74207281,77232917,82589933

第51号梅森素数指数是82589933,4811740号素数;第41号梅森素数指数是24036583,1509263号素数;第31号梅森素数指数是216091,19292号素数;第21号梅森素数指数是9689,1196号素数;第11号梅森素数指数是107,28号素数;

梅素号码        指数        素数号        素数差        间距
21        9689        1196        9582        958.2
31        21609        19292        11920        1192
41        24036583        1509263        24014974        2401497
51        82589933        4811740        58553350        5855335
在41-51号梅森素数之间,平均5855万左右就有一个梅森素数;
用31-51间梅森素数的间距回归,可得回归式:y=5267.7*x^2-139241*x-744558
号数x        y        实际间距        指数
31        962231         1192         21609
41        4082565         2401497         24036583
51        8456439         5855335         82589933
52        8962771         预测指数        91552704
53        9481638         预测指数        101034342
54        10013041         预测指数        111047383
预测52号梅森素数指数在91552704上下,53号梅森素数指数在101034342上下,54号梅森素数指数在111047383上下,……
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-19 14:46 | 显示全部楼层
Lncas-Lehmer梅森素数判定法

采用太阳试除法虽然能判定一个梅森数2^p-1是不是梅森素数,该方法正确有效,但它不是一种简捷的判定方法。
若要判定移动梅森数2^p-1是不是梅森素数,通常采用Lncas-Lehmer梅森素数判定法。

Lncas-Lehmer梅森素数判定法
设p是素数,设第p个梅森数为Mp=2^p-1,定义r1=4,对于k≥2,利用rk≡<rk-1>^2 – 2  (mod Mp),0≤rk<Mp;可以递归得到一整数序列,那么Mp是素数当且仅当<rp-1>≡0 (mod Mp)。

【附注】文中尖括号<rk-1>和<rp-1>中的k-1、p-1都是下标,<rk-1>和<rp-1>都是整体符号,尖括号是笔者外加的;Mp、rk中的p和k也都是下标。

例:考虑梅森数M5=2^5-1=31。那么r1=4、r2≡4^2-2≡14 (mod 31)、r3≡14^2-2=194≡8 (mod31)和r4≡8^2-2=62≡0 (mod 31)。因为r4≡0 (mod 31),故可知M5=2^5-1=31是素数。
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-1-20 08:08 | 显示全部楼层
如果p模4余3,p和2p+1都是素数,则2p+1能够整除梅森数2^p-1;
如23|M11,47|M23,167|M83,263|M131,359|M179,383|M191,479|M239,503|M251,……
M11、M23、M83、M131、M179、M191、M239、M251等都不是素数。

如果p模4余1,p和2p+1都是素数,则2p+1不能整除梅森数2^p-1,但能够整除2^p+1。
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-6-21 06:11 , Processed in 0.109275 second(s), 15 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表