级数与平方数

蔡家雄

兔子数列中的勾股数

1,1,2,3,5,8,13,21,34,55,89,144,233,377,610,

987,1597,2584,4181, ......



设兔子数列中的任意四个连续的兔子数：
第一个为a,第二个为b,第三个为c,第四个为d.

则 (ad)^2+(2bc)^2=(b^2+c^2)^2

蔡家雄

佩尔数列与再生数列的平方性质：

Pn = [(1+√2)^n - (1 - √2)^n]/√8
    = 1，2，5，12，29，70，169，408，......

Pn*P(n+1)*P(n+2)*P(n+3)+1 = 完全平方数。


Cn = [(1+√2)^n+(1 - √2)^n]/2
    = 1，3，7，17，41，99，239，577，......

Cn*C(n+1)*C(n+2)*C(n+3)+4 = 完全平方数。


兔子数的平方性质
fn = [((1+√5)/2)^n - ((1 - √5)/2)^n] /√5
    = 1，1，2，3，5，8，13，21，......

f(2n), f(2n+2), f(2n+4) 和 4*f(2n+1)*f(2n+2)*f(2n+3),
在这四个数中，任意两个的乘积，再+1，是一个完全平方数。
1*3+1=2^2
1*8+1=3^2
1*120+1=11^2
3*8+1=5^2
3*120+1=19^2
8*120+1=31^2


卢卡斯数的平方性质（讨论）
Ln = ((1+√5)/2)^n+((1 - √5)/2)^n
    = 1，3，4，7，11，18，29，47，......

蔡家雄

连续平方和趣题：

求出n+1个连续平方数之和等于n个连续平方数之和的通解公式。

3^2+4^2=5^2

10^2+11^2+12^2=13^2+14^2

21^2+22^2+23^2+24^2=25^2+26^2+27^2

36^2+37^2+38^2+39^2+40^2=41^2+42^2+43^2+44^2

等号左边第一个数是三角数w=m*(2m+1)

等号右边第一个数是       m^2+(m+1)^2


求解：毕氏方程
a^2+b^2 = c^4

7^2+24^2=5^4
119^2+120^2=13^4
527^2+336^2=25^4
1519^2+720^2=41^4
3479^2+1320^2=61^4
6887^2+2184^2=85^4

由我另类公式解：
a = (2k^2+2k -1)^2 -2,
b = 4k(k+1)(2k+1),
c = 2k^2+2k+1.

此时：
当 a < b 时，a为勾，b为股，
当 a > b 时，b为勾，a为股，
即 a 可为勾，可为股，b 亦如是。


罗士琳勾股数本原解公式

设 奇数Q=m+n,（m,n 互质 且 m>n, m,n 均为正整数）

则 [Q*(m-n)]^2+(2mn)^2=[m^2+n^2]^2 有 E/2组的本原勾股数。

其中，E 就是著名的 Euler 函数。但，不是朱火华的公式。


方程：\((b+1)^3+(b+2)^3+.....+(b+a)^3=p^3\)（b,a,p为自然数)有无穷解。

证：当a为立方数时，即\(a= i^3\)时，可求得\(b=( i^4-3i^3-2i^2-2)/6；p=(i^5+i^3-2i)/6\)。

易见，当i=3k时，方程无解；

当i=3k±1时，方程有解，从而原方程有无穷解。

例：当i=5,即a=125时，\(34^3+35^3+……+158^3=540^3\)

当i=10，即a=1000时，\(1134^3+1135^3……+2133^3=16830^3\)

当i=11,即a=1331时，\(1735^3+1736^3……+3065^3=27060^3\)

蔡家雄

由 \((y^2 - x^2)^2+(2xy)^2=(y^2+x^2)^2\) ,

设 \(x, y\) 为正整数，且 \(x < y\)，且 \(x与y\) 互素，

求 \(|y^2 - x^2 - 2xy| =2023\) 的 \(2^2\) 组 \(( x , y )\) 的通解公式，

即 两直角边相差 \(2023\) 的本原勾股方程 的通解公式。

\(x, y\) 是 \(A_{n}=\frac{(64 - 5\sqrt2) (1 +\sqrt2)^{n} + (64 + 5\sqrt2)(1 -\sqrt2)^{n}}{4}\) 中连续的两项，

\(x, y\) 是 \(B_{n}=\frac{(88 - 43\sqrt2) (1 +\sqrt2)^{n} + (88 + 43\sqrt2)(1 -\sqrt2)^{n}}{4}\) 中连续的两项，

\(x, y\) 是 \(C_{n}=\frac{(64 + 5\sqrt2) (1 +\sqrt2)^{n} + (64 - 5\sqrt2)(1 -\sqrt2)^{n}}{4}\) 中连续的两项，

\(x, y\) 是 \(D_{n}=\frac{(88 + 43\sqrt2) (1 +\sqrt2)^{n} + (88 - 43\sqrt2)(1 -\sqrt2)^{n}}{4}\) 中连续的两项，

蔡家雄

x^2 - (4k+1)*y^2 = ±1 和 x^2 - (8k+2)*y^2 = ±1

存在通项公式的条件：

4k+1 或 8k+2 的(奇数)素因子均具有 4d+1 的形式。

A031396                Numbers k such that Pell equation x^2 - k*y^2 = -1 is soluble.        

2, 5, 10, 13, 17, 26, 29, 37, 41, 50, 53, 58, 61, 65, 73, 74, 82, 85, 89, 97,

101, 106, 109, 113, 122, 125, 130, 137, 145, 149, 157, 170, 173, 181, 185, 193, 197,

202, 218, 226, 229, 233, 241, 250, 257, 265, 269, 274, 277, 281, 290, 293, 298 , ...


A031397 数据 无4k+3 的素因子，此方程 存在 正整数解，

A031397                Nonsquarefree n such that Pell equation x^2 - n y^2 = -1 is soluble.               

50, 125, 250, 325, 338, 425, 845, 925, 1025, 1250, 1325, 1445, 1450, 1525, 1625, 1682, 1825, 1850,

2050, 2125, 2197, 2425, 2725, 2738, 2825, 2873, 2890, 3050, 3125, 3250, 3425, 3625, 3725, 3925,

4250, 4325, 4394, 4625, 4825, 4901, 4913 , ...


A031398 数据 无4k+3 的素因子，但方程 不存在 正整数解，

A031398                Squarefree n with no 4k+3 factors such that Pell equation x^2 - n y^2 = -1 is insoluble.               

34, 146, 178, 194, 205, 221, 305, 377, 386, 410, 466, 482, 505, 514, 545, 562, 674, 689,

706, 745, 793, 802, 866, 890, 898, 905, 1154, 1186, 1202, 1205, 1234, 1282, 1345, 1346,

1394, 1405, 1469, 1513, 1517, 1537, 1538, 1717, 1762, 1802, 1858 , ...


求 \(x^2 - ((2n+1)^2+4)*y^2= -1\) 的最小解

则 \(x=(2n+1)*((2n+1)^2+3)/2\) , \(y=((2n+1)^2+1)/2\) .

求 \(x^2 - ((2n+1)^2+4^k)*y^2= -1\) 的最小解，



求 \(x^2 - (n^2 -2)*y^2=1\) 的最小解

则 \(x=n^2 -1 , y=n\) .

求 \(x^2 - (n^2+2)*y^2=1\) 的最小解

则 \(x=n^2+1 , y=n\) .

求 \(x^2 - ((2n)^2 -4)*y^2=1\) 的最小解

则 \(x=2*n^2 -1 , y=n\) .

求 \(x^2 - ((2n)^2+4)*y^2=1\) 的最小解

则 \(x=2*n^2+1 , y=n\) .

求 \(x^2 - ((2n+1)^2 -4)*y^2=1\) 的最小解

则 \(x=(2n+1)*(2n*(n+1) -1) , y=2n*(n+1)\) .

求 \(x^2 - ((2n+1)^2+4)*y^2=1\) 的最小解

则 \(x=2*((2n+1)^2+4)*(n^2+(n+1)^2) -1\) ,

     \(y=(4n+2)*(n^2+(n+1)^2)*(n^2+(n+1)^2+1)\) .


若 5^(2k) -2 是素数，则 10 是素数 5^(2k) -2 的原根。

得 2k=2, 14, 26, 50, 126, 144, 260, 624, 1424, ......


若 5^(2k)+4 是素数，则 10 是素数 5^(2k)+4 的原根。

得 2k = 2，6，10，102，494，794，1326，5242, 5446, ......


用公式法求解特殊佩尔方程

设 \(p=4k+1\) 是素数，

求 \(x^2 - p*y^2=1\) 的最小解，

设 \(x=2p*r^2 -1\) , 求 最小的 \(r=?\) ,

使 \(y=((2p*r^2)*(2p*r^2 -2)/p)^{1/2}\) 是整数。

推论：此时，

设 \(p=4k+1\) 是素数，

求 \(x^2 - p*y^2= -1\) 的最小解，

得 \(y=r\) ,  \(x=((2p*r^2)*(2p*r^2 -2)/p)^{1/2}/(2*r)\) .


设 \(d=8k+3\) 是素数，

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2 -1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2 -2)/d)^{1/2}\) 是整数。

设 \(d=8k+7\) 是素数，

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2+1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2+2)/d)^{1/2}\) 是整数。


设 p=4k+1 是质数，

则 x^2 - p*y^2=±p 都有解，并求出它的的最小解，

5 [5, 2] [20, 9]
13 [65, 18] [2340, 649]
17 [17, 4] [136, 33]
29 [377, 70] [52780, 9801]
37 [37, 6] [444, 73]
41 [205, 32] [13120, 2049]
53 [1325, 182] [482300, 66249]
61 [232105, 29718],[13795392780, 1766319049]
73 [9125, 1068] [19491000, 2281249]
89 [4717, 500] [4717000, 500001]
97 [55193, 5604],[618603144, 62809633]

由 两个 x 都是 p 的倍数，可以用公式法求解此方程，

由 x^2 - p*y^2=±1 可以推导出 x^2 - p*y^2=±p 的最小解。


求 \(x^2 - ((2n+1)^2+4)*y^2=1\) 的最小解

则 \(x=2*((2n+1)^2+4)*(n^2+(n+1)^2) -1\) ,

     \(y=(4n+2)*(n^2+(n+1)^2)*(n^2+(n+1)^2+1)\) .


求 \(x^2 - ((2n+1)^2+4)*y^2= -1\) 的最小解

则 \(x=(2n+1)*((2n+1)^2+3)/2\) , \(y=((2n+1)^2+1)/2\) .


求 \(x^2 - ((2n+1)^2+4)*y^2=±((2n+1)^2+4)\) 的最小解

\((2n+1)^2+4=5, 13, 29, 53, 85, 125, 173, 229, 293, 365, ..... \)

(5,[5,2],[20,9])
(13,[65,18],[2340,649])
(29,[377,70],[52780,9801])
(53,[1325,182],[482300,66249])
(85,[3485,378],[2634660,285769])
(125,[7625,682],[10400500,930249])
(173,[14705,1118],[32880380,2499849])
(229,[25877,1710],[88499340,5848201])
(293,[42485,2482],[210895540,12320649])
(365,[66065,3458],[456905540,23915529])

蔡家雄

【蔡氏完全循环节等差(9240)九生素数有无穷多组】

10是如下等差9生素数(公差9240)的原根，即：这些素数倒数具有最大循环节长！！！

1---(95339,104579,113819,123059,132299,141539,150779,160019,169259)

2---(7827167,7836407,7845647,7854887,7864127,7873367,7882607,7891847,7901087)

3---(9195167,9204407,9213647,9222887,9232127,9241367,9250607,9259847,9269087)

4---(32288903,32298143,32307383,32316623,32325863,32335103,32344343,32353583,32362823)

5---(59941697,59950937,59960177,59969417,59978657,59987897,59997137,60006377,60015617)

6---(72980177,72989417,72998657,73007897,73017137,73026377,73035617,73044857,73054097)

7---(77003567,77012807,77022047,77031287,77040527,77049767,77059007,77068247,77077487)

8---(121526753,121535993,121545233,121554473,121563713,121572953,121582193,121591433,121600673)

9---(121535993,121545233,121554473,121563713,121572953,121582193,121591433,121600673,121609913)

10---(171184589,171193829,171203069,171212309,171221549,171230789,171240029,171249269,171258509)


蔡氏完全循环节等差(120120)十生素数有无穷多组

蔡家雄

具有完全循环节的一条龙素数及其代码

设 n≥3 ,              
                                                         
若 (10^n - 1)÷9×2+1是素数，   
                                             
则 10是(10^n - 1)÷9×2+1的原根，
  
则 1/[(10^n-1)÷9×2+1] 具有最大的完全循环节长。                                          
                                                                              
有 n=3, 8, 11, 36, 95, 101, 128, 260, 351, 467, 645, 1011, 1178, 1217, 2442，......

ForIf[n = 101;
PowerMod[10, ((10^n - 1)/9*2 + 1 - 1)/2, (10^n - 1)/9*2 + 1] == (10^n - 1)/9*2]


设 n≥3 ,         
                                                                                 
若 (10^n - 1)÷9×3+4是素数，  
                                                
则 10是(10^n - 1)÷9×3+4的原根，

则 1/[(10^n-1)÷9×3+4] 具有最大的完全循环节长。                                            
                                                                              
有 n=3, 6, 46, 394, 978, 2586, 2811, 2968，......

ForIf[n = 394;
PowerMod[10, ((10^n - 1)/9*3 + 4 - 1)/2, (10^n - 1)/9*3 + 4] == (10^n - 1)/9*3 + 3]


设 n≥3 ,   
                                                                 
若 (10^n - 1)÷9×4+3是素数，
                                                  
则 10是(10^n - 1)÷9×4+3的原根，

则 1/[(10^n-1)÷9×4+3] 具有最大的完全循环节长。                                               
                                                                                 
有 n=4, 10, 20, 26, 722, 1310，......

ForIf[n = 722;
PowerMod[10, ((10^n - 1)/9*4 + 3 - 1)/2, (10^n - 1)/9*4 + 3] == (10^n - 1)/9*4 + 2]      


设 n≥3 ,      
                                                                       
若 (10^n - 1)÷9×8-1是素数，
                                               
则 10是(10^n - 1)÷9×8 -1的原根，

则 1/[(10^n-1)÷9×8 -1] 具有最大的完全循环节长。                                         
                                                                           
有 n=3, 4, 6, 9, 12, 72, 118, 124, 190, 244, 304, 357, 1422, 2691，......  

ForIf[n = 118;
PowerMod[10, ((10^n - 1)/9*8 - 1 - 1)/2, (10^n - 1)/9*8 - 1] == (10^n - 1)/9*8 - 2]      


设 n≥3 ,     
                                                            
若 (10^n - 1)÷9×2+7是素数，
                                                   
则 10是(10^n - 1)÷9×2+7的原根，

则 1/[(10^n-1)÷9×2+7] 具有最大的完全循环节长。                                             
                                                                                       
有 n=3, 5, 14, 176, 416, 2505, 2759，.......

ForIf[n = 176;
PowerMod[10, ((10^n - 1)/9*2 + 7 - 1)/2, (10^n - 1)/9*2 + 7] == (10^n - 1)/9*2 + 6]


设 n≥3 ,
                                                                    
若 (10^n - 1)÷9×7+2是素数，
                                             
则 10是(10^n - 1)÷9×7+2的原根，

则 1/[(10^n-1)÷9×7+2] 具有最大的完全循环节长。                                          
                                                                           
有 n=66, 86, 90, 102, 386, 624，......

ForIf[n = 102;
PowerMod[10, ((10^n - 1)/9*7 + 2 - 1)/2, (10^n - 1)/9*7 + 2] == (10^n - 1)/9*7 + 1]

蔡家雄

用公式法求解特殊佩尔方程

设 \(d=2^{2n+1} -1 \) ，

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2+1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2+2)/d)^{1/2}\) 是整数。

1 (7, [8, 3])
2 (31, [1520, 273])
3 (127, [4730624, 419775])
4 (511, [4188548960, 185290497])
5 (2047, [5566204448, 123026943])
6 (8191, [137168708703221032895073625802589931406141894540466576152225731961377054107539335040, 1515606947046439809691055132298479939258772521204671748899078157656970867149893633])
7 (32767, [4989395453057620761318437024, 27563196825452342280188415])
8 (131071, [68012246112042192848000, 187859780807502143487])
9 (524287, [46997567320561778145191353610962878307987934593756793094427946763383336928816934242481086061476715595745859044117929495717576882495512013272009350122748657837297172805929900209856744164860003129102095470554027278444403614698310930195231667607295884200937188090279571190579061930767996574039047421403571939359840263187167295467057381580417022064357611621451352561562440161111807908201761856905236583371614657642824269432164359502751873003576866612568082872108067683541142293147523053373426326394532201231880879016720914034037527494035878025515276564347843206711912343432348946721305069707108143792557407643740316164224, 64906895008741122476222937583214609214243799343511042599770186693462505668024604821997293668498982744884496689072778902714164446881210975097103993081187378767551871768967055215221731508471912677599134669670778307686638503086861330068780326391434525717635311171812276918961490673285815482166096497112496202129039439516524661906563996866258177830913974661644392881180263757110495349973522174642669167999546080787182227944467163959051276031077066053223740563434912973906942490225545215546942556046567271856598842966734789024069142174111721341814056011491824374521545759295523871629793392166520592149699929568200351745])
10 (2097151, [94716390224614473408540822118113461048452313411820973215983824757335816500342324560789047374282847806009714778621972049086083922828444965267321229715897075070580677225533794783475554185165522956179202380273231115554458953955517542095864979978481990501635325801973427056692681449780979129744252153653028410826103368612341846012894446748561572667652160, 65404900180987526520505029585445844658185900613734410857615652737952867201530958096058535800003294047733713613634131099428721918752038293902093037577350108052529075142957670861394586353233358311356186212004814489737922917832914662488617368055701767664767158479023604548319172656813646225620266344943525094720499185141955249955405063452172806391807])


设 \(d=2^{2n+1}+11 \) ，

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2 -1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2 -2)/d)^{1/2}\) 是整数。

1 (19, [170, 39])
2 (43, [3482, 531])
3 (139, [77563250, 6578829])
4 (523, [81810300626, 3577314675])
5 (2059, [76292657571747809951, 1681336224678949560])
6 (8203, [69876915611822519318, 771520323981133629])
7 (32779, [80398835008645600324801235608771036666670182431880486400032257485064770658401842307691224942422999531536156600364499923056796970, 444070481702336563743442162036596677391329202863273716321349178397032836203330537445461630919283792936487491412555761435022659])
8 (131083, [712391055453831000489441596851496374586922211968388426389868928977886834600835988092416446450847, 1967638319811445164614618584089426421183233564951501613871659298300476433553720642927098589576])
9 (524299, [262370423419022052948006380362585230636886087574885576160106593835878424107999431011992139225761227061447388396667090, 362347576933161309620975692036326829940654492435422509507444016197507206937437226010673272525903902430923747403149])
10 (2097163, [66551049079488501983959284853373248448554106736699729803934265156536874676626434177518523856595115847326033553714770905049243078368745726019557433640728343803399, 45955639339182891136502078220721866137731621744981660957938911432411265785268208500581497233009913599223286111395352680221603536538078514618343979427535429580])


设 \(t>=0\) , \(d=4^{t+2}+3\) ,

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2 -1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2 -2)/d)^{1/2}\) 是整数。

0 (19, [170, 39])
1 (67, [48842, 5967])
2 (259, [847225, 52644])
3 (1027, [133150393, 4154868])
4 (4099, [6290679704195680536253612607008970662569749959960678667690, 98255894580947570793022586715791022877883393268276434001])
5 (16387, [11413027877434405409306870876631302, 89156118179665453435575511820337])
6 (65539, [6626987825251499940163644435438427657285605581778805650604671321439032704650227419569798702986515177884092626923423746897824217909664252366878028116392925544302544564813604764312894508009793862769489442946848410, 25886078713846848077540951727476971048660042171394841096343736446933897961080712295438431388466376567441157759365264502952910937004522407589434559131549013776596141391286134866897248681517810343464137085201679])
7 (262147, [461824715334079741401519303469504108330873412595592379829070562602875729178271042587960795189884479307292999341798093671040889345880714782458533465979666854073891767873038415959346557561109813822844640447229155501280451825605904774502577690098914740105088441699938619497758257893185226, 901996235887778937057743017583816500902733916131179182975046645499562912967305336150598729595983745922905890921818634606381383310249439070803040676565484234897617020025864775061633330791260216497072863048513377755273583285544543505939632797367764427395085022540331278012923129671185])
8 (1048579, [29776562622206178677730301142172546168291680332575443613253737607023929466971985325370114852698324417371451354347021086491467628966574696503298021764323853197195987312312976427633035743421980526040891707884185731169885041010222716536916751344923695073488237907927706327511181595881545468470585, 29078632838460032672666403050034204085460299673474260078426358718677635257573796496235301168927215703715536805459846463870144135282884727833097966936851168466377694310989456524746214821453380704124450764599593272341211324167432854315242930574639751471817538820290318871993958129089261314516])
9 (4194307, [25117668700757563653335457424693186499343253477468658168439843983648792643078346871742185090650829273128220192287509997202739528105193307898223534107084812545193085970193081959010079271378513867355783149228660891615201876345643100234555174885698, 12264482284171904013989712412289501754938887650714341377823346440598175567868634964373211558180631462153143335767287279311805837025346887031087964999355087525229310401332966841853143651995408496582378346852444671530903857372029317908160695423])
10 (16777219, [165538855338538541294520542671406644531477493498893623619963986289283279108371357103275177798948570988234668100614373256599076525178579541017528601789196799, 40414755990774785816699650975371229136677753557318861712986530079888102157828602884103566909987591921158096123624204149502584894033348180791836050488160])


设 \(t>=0\) , \(d=4^{t+2}+7\) ,

求 \(x^2 - d*y^2=1\) 的最小解，

设 \(x=d*r^2+1\) , 求 最小的 \(r= ?\)

使 \(y=((d*r^2)*(d*r^2+2)/d)^{1/2}\) 是整数。

0 (23, [24, 5])
1 (71, [3480, 413])
2 (263, [139128, 8579])
3 (1031, [651737448664200, 20297537082877])
4 (4103, [36076932, 563221])
5 (16391, [52472886125273249168576935, 409856877328651972146708])
6 (65543, [41780761229710759094849505962993865917706803726191499359145061011994559849082511321952117034371527885710237688, 163197383105723735856182252498452836185524714027012978030385027507904070465599631685229749329039036679421949])
7 (262151, [72919501603120923013270672024141768578706836731811057353841647799969090234358400005317356344197149147239301670533135933749679112145897442497905988746497586725820191966341208118863314358760, 142419000082493600157839034265172632792305052672094541862089637471232086383865308126948109745798108066784404636850506703442647442173796351018334532948853648798894885826723428723794562307])
8 (1048583, [418229444023502271361812404617873926746992317871256949587464210764596500805651130752157377718599287500969695548929536492357103050025530995451192612782974155158251130179318873881058564997395623111811664356164877077981455623488750429751206076422161607460692253488914657675151863967623958352674686098942388140029808502995477518988282157821379091543923746773312815140290495995102220328, 408425828163197962003832720560210417423586434034709599430402144473753992372140369540717657424064448265632501677112941950805810763684256393223756058397515042741908003309627714276604060934296253320289390884127450710539654042120642168558451580131953471841986579370509209368405079585980274711779811151691110790037070509897410512742552098517092232184815889988689044506228933231818499])
9 (4194311, [51177054942435006277229234536620832079444803721125541132021470455820858871814804794857324858228402746532292339295515713235199790870528664021032810391685783026504588734199209057085306058647981957842238967078368667128165377355648672803222147567944156720419449185256655797769215055178564765985227094358168236975681995839991327487079917829920872375810321019832137169884560607208961828411383883914582356557322270066580657182774406297808991338679202909994511114671232104409331590335330835903373122527635746167780, 24988775506360319141826550074426223136987102897607547179104530942622750976201631646786329811140689212765079836744732351129725319420842450572656257702441606875729687408639638340611805165566960030674676979535174702602425251628674213251143161381403282727563043351599695274933985418257973274185624894371898795501862187860896736106867593870932741384823998126520338238497521194509632197513334484152395781003888156585302734563446428778824202320580122873260868579459893393948197917193337027663923215160591349397])
10 (16777223, [3352794142606776554737283396545455877813836162704255416005798405010451840621496841467507112294049620163554858705779009781925131361764415613114860399109051987646983019114153032902195590259008902005124786195991201925356749188180604386663455448808967965603181185620401171730273085591330497829237152, 818553086708894498533150539407674059525330169184097833631896269142138216535942857956080806385820036837989137822853140264221488514005189740716793850834153613502827027947737843684168670016020634664318345945443367495486986629629040276212942493832107589655351652267594369708052478558593705632669])

蔡家雄

判定梅森质数的卢卡斯序列

卢卡斯级数的通项公式

Ln=Round[((1+√3)/√2)^(2^n )/2]

L1=2,
L2=7,
L3=97,

L4=18817=(2^5 -1)(2^5*19 -1)=31*607 = 两个梅森质数的乘积，

并且：2^31 -1 与 2^607 -1 同为素数。

即有：2^30*(2^31 -1) 与 2^606*(2^607 -1) 都是 完全数。

L5=708158977,

L6=1002978273411373057

   =(2^7 -1)(2^7*61698958748239 -1)

   =127*7897466719774591 = 两个梅森质数的乘积，

已证：127 和 7897466719774591 都是素数，

已证：2^127 -1 是素数，据此 7897466719774591 是梅森素数，

即有：2^126*(2^127 -1) 是完全数，

以及：2^7897466719774590*(2^7897466719774591 -1) 是 完全数。

蔡家雄

素数倒数最大循环节长定理

设 k 为非负整数，

若 30k+7 和 120k+29 都是素数，

则 1/(120k+29) 具有最大循环节长d= 120k+28.


蔡氏完全循环节问题

设 n>=3,

设 P 和 2^n*P+1 都是素数，

且 10^(2^n) -1 不能被 2^n*p+1 整除，

若 2^n*P+1 ≡ 17或33（mod  40），

则 10 是 2^n*P+1 的原根，

则 1/(2^n*p+1) 具有最大循环节长d= 2^n*p .


设 k 为正整数，t 为非负整数，

若 30k+17 和 2^(4t+3)*(30k+17)+1 都是素数，

则 10 是素数 2^(4t+3)*(30k+17)+1 的原根。

若 30k+29 和 2^(4t+3)*(30k+29)+1 都是素数，

则 10 是素数 2^(4t+3)*(30k+29)+1 的原根。


设 k 为正整数，t 为非负整数，

若 30k+1 和 2^(4t+4)*(30k+1)+1 都是素数，

则 10 是素数 2^(4t+4)*(30k+1)+1 的原根。

若 30k+7 和 2^(4t+4)*(30k+7)+1 都是素数，

则 10 是素数 2^(4t+4)*(30k+7)+1 的原根。


设 k 为正整数，t 为非负整数，

若 30k+11 和 2^(4t+5)*(30k+11)+1 都是素数，

则 10 是素数 2^(4t+5)*(30k+11)+1 的原根。

若 30k+23 和 2^(4t+5)*(30k+23)+1 都是素数，

则 10 是素数 2^(4t+5)*(30k+23)+1 的原根。


设 k 为正整数，t 为非负整数，

若 30k+13 和 2^(4t+6)*(30k+13)+1 都是素数，

则 10 是素数 2^(4t+6)*(30k+13)+1 的原根。

若 30k+19 和 2^(4t+6)*(30k+19)+1 都是素数，

则 10 是素数 2^(4t+6)*(30k+19)+1 的原根。



梅森素数与蔡氏完全循环节问题

设 素数 p >=7，

且 (2^p -1)=30k+1 或 (2^p -1) =30k+7 也是素数，

若 (2^p -1)*16^m+1=(2^p -1)*2^(4m)+1 是素数，

则 10 是素数 (2^p -1)*16^m+1 的原根，

则 1/((2^p -1)*16^m+1) 具有最大循环节长d= (2^p -1)*16^m .



蔡氏完全循环节问题

若 3^(2n)+2^(2n+1) 是素数，

则 10 是素数 3^(2n)+2^(2n+1) 的原根。

10 是素数 3^2+2^3=17 的原根，
10 是素数 3^4+2^5=113 的原根，
10 是素数 3^6+2^7=857 的原根，
10 是素数 3^12+2^13=539633 的原根，
10 是素数 3^22+2^23=31389448217 的原根，
10 是素数 3^32+2^33=1853028778786433 的原根，
10 是素数 3^36+2^37=150094772735952593 的原根，
10 是素数 3^46+2^47=8862938260389989451257 的原根，
10 是素数 3^80+2^81=147808829414348341167722439464732709953 的原根，
10 是素数 3^154+2^155=29969067287845284806900763424259354345695037325432711901413781193689500137 的原根，
10 是素数 3^236+2^237=39867234790105605031052158475473603885214702979674478224207279045242447503005351752119734009677347787327998900593 的原根，
10 是素数
3^250+2^251=190683748116796615589766511371277507701260429967651126103174761897479526555470283571068048447508517209471538025816027497 的原根，
10 是素数 3^992+2^993=201504468751837621839727977404685926835150439376946832443975059933977838339689818555216935505351023009283611903196552713908086356058621995407002396777618301177811343575418503951721787383304872793005232962627245632571105557936833267194815304366778547765493764740044293673929368039251196446446723938278591809721982615857996114286344721172111746288859733499249290109426007919969472996555103682274449696455944323255567670689270574629521811105681761209036785000603422679143272833 的原根，