2025年探索素?cái)?shù)奧秘實(shí)驗(yàn)室揭秘之旅

素?cái)?shù)一.試驗(yàn)解讀素?cái)?shù)的問(wèn)題是數(shù)論里最富有魅力，最吸引人的問(wèn)題，在本次試驗(yàn)中我們將研究討論素?cái)?shù)的某些有關(guān)規(guī)律及性質(zhì)，在看似無(wú)規(guī)律的素?cái)?shù)中尋找規(guī)律。要想研究素?cái)?shù)，首先要理解素?cái)?shù)的定義，學(xué)會(huì)用定義去判斷一種數(shù)與否是素?cái)?shù)，考慮合數(shù)與素?cái)?shù)之間具有什么樣的關(guān)聯(lián)。通過(guò)某些有關(guān)結(jié)論猜測(cè)素?cái)?shù)的無(wú)窮性并給出證明，并分別用試除法和篩法列出一定范圍內(nèi)的素?cái)?shù)表，比較這兩種措施的有效性。在試驗(yàn)過(guò)程中，我們發(fā)現(xiàn)素?cái)?shù)并沒(méi)有一種簡(jiǎn)要的鑒別公式，那么就需要通過(guò)試驗(yàn)構(gòu)造出素?cái)?shù)的鑒別公式。在研究用整系數(shù)多項(xiàng)式來(lái)生成素?cái)?shù)時(shí)，最關(guān)鍵的是恰當(dāng)?shù)剡x擇多項(xiàng)式的次數(shù)與變量的個(gè)數(shù)。最終，通過(guò)研究一定范圍內(nèi)的素?cái)?shù)個(gè)數(shù)隨整數(shù)增長(zhǎng)而變化的關(guān)系，得出素?cái)?shù)的分布特性。觀測(cè)它的變化關(guān)系，并用函數(shù)將素?cái)?shù)的分布表達(dá)出來(lái)。1.素?cái)?shù)的鑒別與個(gè)數(shù)在不小于1的自然數(shù)中，只能被1和它自身整除的數(shù)稱為素在素?cái)?shù)研究中，一種最基本的問(wèn)題是素?cái)?shù)究竟有多少個(gè)，與否是無(wú)窮的。調(diào)大n至25,觀測(cè)并得出結(jié)論。再將n調(diào)至30,35……,結(jié)論與否發(fā)生了變化。根據(jù)以上的成果，猜測(cè)素?cái)?shù)與否有無(wú)窮多種的證明。2素?cái)?shù)表的構(gòu)造給出一種范圍，用Eratosthenes篩法和試除法列出所有的素?cái)?shù)，它們的原理為：Eratosthenes篩法的基本原理，將自然數(shù)列從2開(kāi)始按次序排列至某一整數(shù)N,首先，從上述數(shù)列中劃除所有2的倍數(shù)(不包括2),在剩余的數(shù)中，除2外最小的是3.接著，從數(shù)列中劃除所有3的倍數(shù)(不包括3),然后在剩余的數(shù)中，再劃去5的這個(gè)過(guò)程一直進(jìn)行下去，則最終剩余的數(shù)就是不超過(guò)N的所有素?cái)?shù)。試除法：假設(shè)我們已經(jīng)懂得前n個(gè)素?cái)?shù)p1=2,p2=3,..,pn,為找下個(gè)素?cái)?shù)，我們從pn+2開(kāi)始依次檢查每一種整數(shù)N,看與否能被某一種pi(i=1,2,….,n)整除，若N能被前面的某個(gè)素?cái)?shù)整除，則N為合數(shù)，否則N即為下一種素?cái)?shù)pn+1。為了提高效率我們只需要用不超過(guò)N^(1/2)的素?cái)?shù)清除就可以了。較少。將范圍調(diào)大，用這兩種措施列出10000,100000……以內(nèi)3素?cái)?shù)的鑒別公式得的余數(shù)。將m的值固定，變化n的值為2,3,……100取m=2,觀測(cè)2^(n-1)被n整除所得的余數(shù)取m=3,觀測(cè)3^(n-1)被n整除所得的余數(shù)取m=4,觀測(cè)4^(n-1)被n整除所得的余數(shù)取n=2,m=2,3,4,……,20,觀測(cè)m^(2-1)被2整除所得的取n=3,m=2,3,……,20,觀測(cè)m^(3-1)被3整除所得的余數(shù)取n=5,m=2,3,……,20,觀測(cè)m^(5-1)被5整除所得的余數(shù)形如2^n-1的數(shù)稱為Mersenne數(shù)，通過(guò)Mersenne數(shù)我們可以研究數(shù)論中的有關(guān)性質(zhì)。觀測(cè)并考慮Mersenne數(shù)與n的4.生成素?cái)?shù)的公式Fermat數(shù)：我們把形如F=22"+1表達(dá)出來(lái)的數(shù)稱為Fermat數(shù)。Fermat數(shù)與否都是素?cái)?shù)?在程序中增大n的值，很輕易懂得既然Fermat數(shù)不能作為素?cái)?shù)的生成公式，那么能不能尋求首先考慮一次函數(shù)，顯然是不行的。再考慮二次多項(xiàng)式，數(shù)，令f(n)=2n3+2n2+2n+1,f(n)=4n?+3n3+2n2+1,f(n)=n?+n3+n2+n令變量n的次數(shù)不停升高，觀測(cè)得出的成果有什么不一樣。若單變量整系數(shù)多項(xiàng)式不能生成所有的素?cái)?shù)，那么多變量整系數(shù)多項(xiàng)式呢?考慮兩個(gè)變量的函數(shù)，f(n,m)=n+m+4,將兩個(gè)變量的多項(xiàng)式的次數(shù)變?yōu)槎危頵(n,m)=n2+m判斷以上的f(n,m)與否生成的均是素?cái)?shù)，它們之間有什么規(guī)若還是無(wú)法找出這樣的兩個(gè)變量整系數(shù)多項(xiàng)式，再變化多項(xiàng)式的變量的個(gè)數(shù)和次數(shù)。得出概括性的結(jié)論。5.素?cái)?shù)的分布在上面的試驗(yàn)中我們已經(jīng)懂得了素?cái)?shù)是無(wú)窮多種的，并且素?cái)?shù)的生成公式并不是很明了，不過(guò)它的分布會(huì)不會(huì)具有什么樣的規(guī)律呢?π(100,200)、π(1000,1100)、π(10000,10100)、π(100000,100100)。從計(jì)算成果看，伴隨范圍的擴(kuò)大，素?cái)?shù)是越來(lái)越稀還是越來(lái)越密?深入，選用某些更長(zhǎng)的區(qū)間，做同樣的試驗(yàn)。將這些點(diǎn)畫(huà)在圖中，從圖中能更清晰的看出素?cái)?shù)的分布狀況。換一種角度考慮，從兩個(gè)相鄰素?cái)?shù)間距的大小同樣也可以看出素?cái)?shù)的分布，這時(shí)我們還可以發(fā)現(xiàn)某些更有趣的規(guī)律。先求出1000以內(nèi)的所有相鄰素?cái)?shù)的間距，并將點(diǎn)以(pn,d)的形式畫(huà)在直角坐標(biāo)系中，觀測(cè)圖像的特點(diǎn)；增大n的值，再在另一種圖中畫(huà)出，從這些點(diǎn)的分布可以看出素?cái)?shù)的間隔值的某些特性，以及它們的反復(fù)次數(shù)的多少，我們還發(fā)現(xiàn)：在增大N的值的同步，圖中的點(diǎn)也會(huì)隨之變高，也就是說(shuō)最大間隔值在變化，那么,存在最大間隔值嗎?給出結(jié)論及有關(guān)證明。用函數(shù)對(duì)素?cái)?shù)的個(gè)數(shù)進(jìn)行擬合。先進(jìn)行線性擬合，選用2到1000中所有的素?cái)?shù)進(jìn)行擬合，再變化擬合的多項(xiàng)式的次數(shù)，比較擬合效果。將點(diǎn)(n,π(n))標(biāo)在平面坐標(biāo)系中，并且用折線把這些點(diǎn)連接起來(lái)，觀測(cè)π(n)的變化趨勢(shì)，然后在程序中增大N的值，再觀測(cè)π(n)的變化趨勢(shì)，將π(n)的值與其他函數(shù)的值進(jìn)行比較，看能否找出最靠近π(n)的值的函數(shù)，即計(jì)算素?cái)?shù)個(gè)數(shù)的公式，注意此時(shí)n應(yīng)當(dāng)充足大。NumP[n_Integer]:=Num=Product[Prime[i],{i,1,n}]+1;Print[Num,"",PrimeQ[Num],"",FactorInte9699691FalseFalse{{33160490131True[211False{(61,1},(46491411False{{95False{{277,1},(34391False{(1063,1},{303049,1},{5988NumP[n_Integer]:=Num=Proct[Prime[i],{i,1,n}]+1;Print[n,"",Num,"",PrimeQ[Num20391False{(1063,1},{3030421492284554131False23False[{265739,1},(126924False{{131,1},(1039,1},{25551False{{2336993,1},{13848803,1},{,1}}再變化n的范圍為25到30,輸出成果為：25551False({2336993,1},{13848803,1},2739621128False{{149,1},(13203797,1},(63137,1},(83711293091False證明：反證法，若素?cái)?shù)的個(gè)數(shù)是有限的，且最大的素?cái)?shù)是pk。p\1,這與p是素?cái)?shù)相矛盾。故而假設(shè)不成立。得證：素?cái)?shù)的個(gè)數(shù)是無(wú)窮多種的。For[i=2,i<=n,i++,AppendToFor[i=1,Prime[i]<=Sqrt[n],i++,temp=t=Select[t,(#1==templ|Mod[#1,temp]!{2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,7,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,239,241,251,257,263,269,271,277,281,283,293,307,311,313,317,331,337,347,34967,373,379,383,389,397,401,409,419,421,431,433,439,443,449,457,461,463,467,47499,503,509,521,523,541,547,557,563,569,571,577,587,593,599,601,607,613,617,619,63,643,647,653,659,661,673,677,683,691,701,709,719,727,733,739,743,751,757,767,797,809,811,821,823,827,829,839,853,857,859,863,877,881,883,887,907,911,919,41,947,953,967,971,977,983,991j=1;divided=False;While[Prime[j]<=Sqrt[i]&&(!didivided=(Mod[i,temp]=If[!divided,AppendTo[t{2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,7,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,239,241,251,257,263,269,271,277,281,283,293,307,311,313,317,331,337,347,34967,373,379,383,389,397,401,409,419,421,431,433,439,443,449,457,461,463,467,47499,503,509,521,523,541,547,557,563,569,571,577,587,593,599,601,607,613,617,619,63,643,647,653,659,661,673,677,683,691,701,709,719,727,733,739,743,751,757,767,797,809,811,821,823,827,829,839,853,857,859,863,877,881,883,887,907,911,919,41,947,953,967,971,977,983,991For[i=2,i<=n,i++,AppendToFor[i=1,Prime[i]<=Sqrt[n],i++,temp=Pt=Select[t,(#1==templ|Mod[#1,temp]!=0)&]];]j=1;divided=False;While[Prime[j]<=Sqrt[i]&&(!didivided=(Mod[i,temp]=If[!divided,AppendTo[For[i=2,i<=n,i++,AppendToFor[i=1,Prime[i]<=Sqrt[n],i++,temp=Pt=Select[t,(#1==templ|Mod[#1,temp]!=0)&]];]j=1;divided=False;While[Prime[j]<=Sqrt[i]&&(!didivided=(Mod[i,temp]=If[!divided,AppendTo[對(duì)n=2,3,…,100,觀測(cè)m^(n-1)被n整除所得的余數(shù)當(dāng)m=2時(shí)運(yùn)行如下程序：M[n_Integer]:=Module[{y,k],m=2;k=m^(n-1);x=MPrint[n,"",PrimeQ[n],"23456789TrueTrueFalseFalse FalseFalse FalseFalseFalseTrue 0101210421812401184218281210429121212121212121212121212121212121212248118214814218124018411289121212121212121212121212121212121212121212121FalseFalseFalseFalseFalseFalseTrueFalseFalse2124184212121212121212121212運(yùn)行如下程序：M[n_Integer]:=Module[{y,k},m=3;k=m^(n-1);x=MPrint[n,"",PrimeQ[n],"1113113113111311311311131113113111311311131131113113113111391113111331139111301311311339111311139131311131從運(yùn)行成果可以發(fā)現(xiàn)：當(dāng)n(除3以外)為素?cái)?shù)時(shí)，3(n-1)被n整除所得的余數(shù)都是1。輸入程序如下：234214010516247118409172411441124114011211441241141624174811241140124141124141121141241142141121412411441124140121141211411241414211412411411214112114FalseFalseFalseFalse12141214414以上我們考慮的是m固定變動(dòng)n的狀況，接下來(lái)分析一下n固定，變動(dòng)m的狀況。M[m_Integer]:=Module[{y,k},n=2;k=m2(n-1);x=MoPrint[m,"",PrimeQ[m],"",GCD[m,23456789TrueFalseTrueFalseFalseFalseFalseFalseFalseFalseFalseFalse212121212121212121212121212121010101010101010101010101010112011201120112011201120112011201120112011201120M[m_Integer]:=Module[(y,k},n=3;k=m^(n-1);x=MoPrint[m,"",PrimeQ[m],"23456789FalseTrueFalseTrueFalseFalseTrueFalseFalseFalse1311311311311310110110110110111311131113113111311131131113111311131113113111311311131113113111311131113113111311131131113113113111311由以上那么多組數(shù)據(jù)，我們可以得出如下結(jié)論：若n為素?cái)?shù)，并且m與n互素，則m^(n-1).該命題的逆命題為：若m^(n-1)被n整除所得的余數(shù)都是1,并且m與n互素，則n為Print[n,"",FactorInteger[n],"",m,"",Prime由以上的成果可以看出，雖然有有m2(n-1)被n整除所得的余數(shù)都是1,并且m與n互素，不過(guò)n不一定就是素?cái)?shù)。如：Mersenne數(shù)的公式為Mn=2^n-1,通過(guò)該公式判斷與否所有Mersenne鑒別式得出的數(shù)都是素?cái)?shù)。Print[n,"",PrintQ[n],"",2^n-1,FactorInteger[216PrintQ[16]65535({3,1},{5,1},(17,1},(218PrintQ[18]262143{{3,3},{7,1},{19,1},{73,20PrintQ[20]1048575觀測(cè)以上程序輸出的成果，可以很輕易發(fā)現(xiàn)：當(dāng)n為合數(shù)時(shí)，它的Mersenne數(shù)也是合數(shù)。當(dāng)n為素?cái)?shù)時(shí)，Mersenne數(shù)就一定是素?cái)?shù)了嗎?觀測(cè)可得當(dāng)n=11時(shí)，M11是2047,是合數(shù)。數(shù)卻不一定是素?cái)?shù)。顯然2^n-1不是素?cái)?shù)。故得證原命題。其鑒別程序如下：If[PrimeQ[2^n-1]==True,a[t]=2~n-Print[b[i],"",a[i],"",PrimeQ[2~a[數(shù)的措施。Do[Fermat[n],{n,1,10}]提出問(wèn)題：素?cái)?shù)生成的公式究竟會(huì)是什么樣的呢?Module[y},y=n2t=0;Eu[n_Integer,m_Integer]:=Do[Eu[n,m],{n,1,10},{m,Module[{y},y=n+mt=0;Eu[n_Integer,m_Integer]:=Do[Eu[n,m],{n,1,10},{m,Module[{y},y=n^2+m*2+4;If[PrimeQt=0;Eu[n_Integer,m_Integer]:=Do[Eu[n,m],{n,1,10},{m,t=0;Eu[n_Integer,m_IntModule[(y},y=n"5+m^5+5;If[PrDo[Eu[n,m],{n,1,10],{m,t=0;Eu[n_Integer,m_IntegDo[Eu[n,m,k],{n,1,10假如多項(xiàng)式為二次的呢?令y=n2+m*2+k^2+5輸出成果為250把100到內(nèi)的素?cái)?shù)在坐標(biāo)軸上表達(dá)出來(lái)，程序如下：ListPlot[t,PlotStyle{RGBColor[0,0,1],PointSize[0.04]}]t:=Table[{i,PrimePi[i+100]-PrimePi[i]},{i,1ListPlot[t,PlotStyleRGBColor[0,0,1],PointJoinedTrue]14},{1000,16},{(1100,12},{1200,15},{1300,11),{1400,17),{1500,1t:=Table[{i,PrimePi[i+100]-PrimePi[i]},{i,1ListPlot[t,PlotStyleRGBColor[0,0,1],PlotJoinedT將i擴(kuò)大到1000到101000的范圍內(nèi)，輸出的圖形如下：ListPlot[t,PlotStyle{RGBColor[987連成折線如下圖：由上圖可以看出，素?cái)?shù)的分布是不規(guī)律的，就整體而言，素?cái)?shù)的分布伴隨整數(shù)的增大而越來(lái)越稀了。用將素?cái)?shù)按照大小次序排列pl=2,p2=3,……,用dn表達(dá)兩相鄰素?cái)?shù)間的間隔，dn=p(n+1)-pn。將1000以內(nèi)的素?cái)?shù)的間距求出的程序如下：Table[Prime[i+1]-Prime[i],t{1,2,2,4,2,4,2,4,6,2,6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2,4,,2,4,2,12,12,4,2,4,6,2,10,6,6,6,2,6,4,2,10,14,4,2,4,14,6,10,2,4,6,8,6,6,4,10,2,6,4,6,8,4,2,4,12,8,4,8,4,6,12,2,18,6,10,6,6,2,6,10,6,6,2,4,6,8,10,8,10,8,6,6,4,8,6,4,8,4,14,10,12,2,10,2,4,2,10,14,4,2,4,14,4,2,4,20,4,84,6,6,14,4,6,6,8,6,12,4,6,2,10,2,6,10,2,10,2,6,18,4,2,4,6,6,8,6,6,22,2,10,8,10,6,6,4,6,6,2,6,12,10,18,2,4,6,2,6,4,2,4,12,2,6,34,6,6,8,18,10,14,4,2,4,6,8,4,22,4,6,12,12,8,12,6,4,6,8,4,8,4,14,4,6,2,4,6,2,6,10,20,6,4,2,24,4,2,10,12,2,10,8,68,6,4,2,12,10,12,8,16,14,6,4,2,4,2,10,12,6,6,18,2,16,2,22,6,8,6,4,2,4,8,,6,12,2,4,2,10,12,2,16,2,6,4,2,10,8,18,24,4,6,8,16,2,4,8,16,2,4,8,6,6,4,12,6,14,6,4,2,6,4,6,12,6,6,14,4,6,12,8,6,4,26,18,10,8,4,6,2,6,22,12,2,16,8,4,8,6,6,4,2,4,6,8,4,2,6,10,2,10,8,4,14,10,12,2,6,4,2,16,14,4,6,8,6,4,18,8,1,14,4,6,6,2,28,2,10,8,4,14,4,8,12,6,12,4,6,20,10,2,16,26,4,2,12,6,4,12,6,8,12,28,2,6,6,6,4,6,2,12,4,12,2,10,2,16,2,16,6,20,16,8,4,2,4,2,22,8,12,6,10,2,12,10,2,10,14,6,4,6,8,6,6,16,12,2,4,14,6,4,8,10,8,6,6,22,6,2,10,14,4,6,2,10,14,4,8,18,4,6,2,4,6,2,12,4,20,22,12,2,4,6,6,2,6,22,2,6,16,6,12,2,6,12,16,2,44,2,18,24,10,6,2,10,2,10,2,10,6,2,10,2,10,6,8,30,10,2,10,8,6,10,18,6,12,12,2,18,6,,18,2,10,14,6,4,2,4,24,2,12,6,16,8,6,6,18,16,2,4,6,2,6,6,10,6,12,12,18,2,62,4,6,2,12,4,14,30,10,6,12,14,6,10,12,2,4,6,8,6,10,2,4,14,6,6,4,6,2,10,2,16,12,8,,12,2,6,6,6,28,6,14,4,8,10,8,12,18,4,2,4,24,12,6,2,16,6,6,14,10,14,4,30,6,2,6,4,2,6,22,6,2,4,18,2,4,12,2,6,4,26,6,6,4,8,10,32,16,2,6,4,2,4,2,10,14,6,4,8,104,2,6,30,4,8,10,6,6,8,6,12,4,6,2,6,4,6,2,10,2,16,6,20,4,12,14,28,6,20,4,18,8,6,4,6,6,10,2,10,12,8,10,2,10,8,12,10,24,2,4,8,6,4,8,18,10,6,6,2,6,10,12,2,10,6,2,6,6,6,10,8,24,6,22,2,18,4,8,10,30,8,18,4,2,10,6,2,6,4,18,8,12,18,16,6,2,12,6,102,6,10,14,4,24,2,16,2,10,2,10,20,4,2,4,8,16,6,6,2,12,16,8,4,6,30,2,10,2,6,4,6,6,82,6,8,12,4,14,12,10,24,6,12,6,2,22,8,18,10,6,14,4,2,6,10,8,6,4,6,30,14,10,2,12,102,18,24,18,6,16,18,6,2,18,4,6,2,10,8,10,6,6,8,4,6,2,10,2,12,4,6,6,2,12,4,14,18,4,,8,6,4,8,4,14,6,4,14,12,4,2,30,4,24,6,6,12,12,14,6,4,2Table[{Prime[i],Prime[i+1]-Prime[i]ListPlot[t,PlotStyleRGBColor[1,0,52000400060t:=Table[{Prime[i],Prime[i+1]-Prime[i]},{i,ListPlot[t,PlotStyleRGBColor[1,0,0]]200004000060000證明：假設(shè)相鄰的素?cái)?shù)的間隔值都為有限數(shù)，最大值設(shè)為N,不過(guò)可以構(gòu)造出N個(gè)持續(xù)的自然數(shù)，具它們都是合數(shù)，它們分別是(N+1)!+2,(N+1)!+3…(N+1)!+N+1,其中第一種能被2整除，其中第二個(gè)能彼3整除，直到最終一種能妝N+1整除，則此時(shí)存在兩相鄰素?cái)?shù)跟離至少為N+1,與原假設(shè)矛盾。即得證相鄰的章數(shù)，它們的間跟可以無(wú)限的大。b=Table[{n,PrimePi[n]},{n,2,g=Plot[ft,{x,2,10000},PlotStyle{RGBColor[1,0,0]},DisplayFunctionIdentity]f=ListPlot[b,PlotStyleRGBColor[0,1,0],DisplayFunctionIdentShow[g,f,DisplayFun60.9457+0.118897xb=Table[{n,PrimePi[n]