Primality Certificate for (2409^2521-1)/2408

Andy Steward8,523 digits20 November 2001
Originally by David Broadhurst & Sean Irvine 2001

This certificate uses a theorem of Konyagin and Pomerance to prove an integer N prime by making use of a partial prime factorization of N-1.

Factorizing N-1

As N is a Generalized Repunit, we make use of the algebraic factorization of N-1 to arrive at the following 32.295995% factorization of N-1:

From Factorisation
24093 · 11 · 73
Φ22 · 5 · 241
Φ35805691
Φ42 · 1453 · 1997
Φ51051 · 32057142031
Φ613 · 446221
Φ77 · 27932067147029528953
Φ82 · 1201 · 9137 · 1534513
Φ9184879 · 1057141729886149
Φ105 · 41 · 71 · 181 · 12778391
Φ1297 · 601 · 577698073
Φ14127 · 211 · 10039 · 726212846059
Φ15271 · 1171 · 506251 · 7057030632301471
Φ1819 · 37 · 199 · 1397051443931329
Φ201134212228063992156159856161
Φ21160704596288476927 · 237592640687909065816783
Φ24314718961 · 3603889704971511601
Φ2829 · 4201 · 1363086283357 · 230021005816152079492097
Φ3031 · 36602685366126790642480831
Φ35642087618645971 · 5087449822134492671088888899921 · p36
Φ362089 · 254506879117 · 71846172145122152107347757
Φ4020441 · 305133761 · 10549469681 · 142662422227441 · 137042785635579401
Φ42757 · 43597 · 408283 · 2836011974215418230001828563
Φ456301 · 1972891 · 11138761 · 752278309659434071 · p47
Φ56113 · 431369 · 4166563577417770147422755750124864833 · p37
Φ6061 · 2341 · 2161388521 · 4118269801 · 36128576401 · 28012969030435467601
Φ632808599347 · c113
Φ70102781120051 · p71
Φ7247881 · 7620459975058392433 · p58
Φ846882289 · 281596761013 · 43028878322977 · 1869645685766718057975421 · 9358438213537196568369649
Φ90991 · 8774789851 · p69
Φ1058592151 · 2204089441 · 19471893773633974051 · p127
Φ1201321 · 24547758297488955058239182908318530763681 · p65
Φ126379 · 55491024003369013041620575474357086932317 · p79
Φ140466560720781 · 109573924825273301 · c134
Φ168337 · 238615512245497297 · 1338355097163031921 · c125
Φ180121212361 · 385714171081 · 138630613335353448509269441 · p117
Φ21038983272931 · 4341542768617607083441 · 56820461446210442972677068841490181035881 · p90
Φ252c244
Φ280281 · 157081 · 64621201 · 127597961 · c302
Φ3152521 · 15121 · 497385045828564211093231 · c456
Φ3605926874584321 · c312
Φ420421 · 1393981 · 6512167848481 · 368226145908411951208322102270401 · c271
Φ5041009 · c484
Φ630631 · 507781 · 690284701 · 1691582131 · c461
Φ840p650
Φ126044101 · 116611984441 · c959
Φ2520383041 · 17866598508361 · 577318924700281 · c1915

From this partial factorization, we use sufficient of the largest prime factors of N-1 so that their product F is at least N 3/10 :

2054308199 9670143874 8331094191 7641627426 9090588361 1866754656 7084687168 1535829939 2560525245 6247398451 4389428876 8277430798 9677245669 4992952996 9563424782 4883897278 2620306342 0520773989 0278971702 2069564451 6497976800 5082816071 2732517110 6826432061 7431200926 8425124198 2647040130 2077479702 2698866127 0786752632 9557518605 3808064196 5657711954 7570660640 9149732600 3392865594 2115092400 7794101874 7044761791 1538698989 2928583492 9352760211 6186625430 4704160684 4372043035 7784497439 7732791821 2566001520 6882538494 8019999107 4338119707 9776163016 5459978329 0749251066 8852876475 7524146488 3852865823 5840992598 6049773927 5954556730 0468575570 8777320228 9428873848 7450089838 7130243841
5775737 0207238107 5336874988 2249162825 2899322875 9786249527 2333201113 6083032459 0255672050 1205130500 2358578497 4368913554 6329755021
3284692 7324328379 7559393124 5883602865 5345559421 8287733545 8454118183 2546236620 8201481427 8526434019 0804674992 3794418481
2212885886 2655663326 2104509822 0082897291 8699474668 3496160497 8549318575 7692514224 1009068611
265010359 2695862436 1456583856 4870785025 6169779254 5142535973 6027019956 7674189247
1 4202019233 3368252501 3128796921 1138524568 2048234462 3535579889 3651899931
167792606 8027119465 2132108187 5949681892 1021071932 7340976560 5509704941
51034 3918473604 1838643570 8252512825 0972041228 8589485739 7601237241
39988839 5582322230 2045868333 4603375094 8814051011 4859555497
1400732 4927184912 0089571507 1995946721 1442804161
5 6820461446 2104429726 7706884149 0181035881
5 5491024003 3690130416 2057547435 7086932317
2 4547758297 4889550582 3918290831 8530763681
7184190 0356018910 8061120101 8969688681
4166563 5774177701 4742275575 0124864833
446486 6194582994 4430392381 2972199691
368 2261459084 1195120832 2102270401
5 0874498221 3449267108 8888899921
28360119 7421541823 0001828563
11342122 2806399215 6159856161
1386306 1333535344 8509269441
718461 7214512215 2107347757
366026 8536612679 0642480831
93584 3821353719 6568369649
18696 4568576671 8057975421
4973 8504582856 4211093231
2375 9264068790 9065816783
2300 2100581615 2079492097
43 4154276861 7607083441
2801296903 0435467601
2793206714 7029528953
1947189377 3633974051
762045997 5058392433
360388970 4971511601
133835509 7163031921
75227830 9659434071
23861551 2245497297
16070459 6288476927
13704278 5635579401
10957392 4825273301
705703 0632301471
139705 1443931329
105714 1729886149
64208 7618645971
57731 8924700281
14266 2422227441
4302 8878322977
1786 6598508361
651 2167848481
592 6874584321
136 3086283357
72 6212846059
46 6560720781
38 5714171081
28 1596761013
25 4506879117
11 6611984441
10 2781120051
3 8983272931
3 6128576401
3 2057142031
1 0549469681
8774789851
4118269801
2808599347
2204089441
2161388521
1691582131
690284701
577698073
314718961
305133761
127597961
121212361
64621201
12778391
11138761
8592151
6882289
5805691
1972891
1534513
1393981
507781
506251
446221
431369
408283
383041
184879
157081
47881
44101
43597
20441
15121
10039
9137
6301
4201
2521
2341
2089
1997
1453
1321
1201
1171
1051
1009
991
757
631
601
421
379
337
281
271
241
211
199
181
127
113
97
73
71
61
41
37
31
29
19
13
11
7
52
3
23

Note that all prime factors listed above have been proven. As primes of under 250 decimal digits can be verified in a few seconds, proof of their primality is not included here, in order to save space. Larger prime factors can take from hours to months to prove; certificates for all such factors have been PKZIPped into this file.

We set R = (N-1)/F. Note that GCD(F,R)=1 and Log(F)/Log(N) = 32.295995%

Finding a Witness to Primality

Next, we find an integer witness w such that for each prime factor p of N-1, w(N-1) ≡ 1 mod N and GCD(w(N-1)/p-1,N) = 1. In this case, w = 17 suffices.

Express N in base F

As F3 < N < F4 and N ≡ 1 (mod F), we can let N = c3·F3 + c2·F2 + c1·F + 1. Let c4 = c3·F+c2.

Square Checks

For t = 0 to 5, we prove that Q(t) = (c1+t·F)2+4·t-4·c4 is not a perfect square. This is done by checking whether Q(t) is a quadratic residue modulo a variety of bases. If it happens to be a QR in all of the bases, we calculate s = floor(sqrt(Q(t))) and show that s2 < Q(t).

Continued Fraction

We approximate c1/F by a continued fraction u/v such that v is maximal while remaining less than F2 / N1/2 = 3654 3746352781 6282635142 0953283702 4594346580 7488881061 6727754124 4142959394 4890864605 0269451128 5913636641 5262529471 7060222080 3830311021 4642797980 5284576171 8531434730 3878593211 8113184667 2726808198 0553752642 8272132266 3535523975 4038596053 2197201011 9503922395 8390816283 3309241291 2639435163 7645906777 2002610686 8529030404 0144145340 3523477038 6170498211 1696450398 2687398258 7191577835 4833128283 0805136647 6924234096 5336141665 2905919149 8007108625 7961482434 5735434999 8065078722 1781646311 2212678814 1334024606 6916704372 8695510040 1763705157 8497601212 7255113805 4024770837 1338966952 6772438295 0101743468 1973033860 5882407094 0811923474 5557213468 0324133312 3966640267 7639480593 3270679037 0331203541 5520136901 6171089073 4455852385 9871185669 6540464859 7646426930 2112899122 6837034402 6727617762 7347916339 6984904570 5956447855 9380666236 9384806033 4133975507 6691419062 2017326029 4310893972 2914396432 0545173560 3508287844 8874785457 1444292253 3877251342 9827984727 9927572205 1908243654 5431701367 6931045711 0474000949 8841798875 4431505701 1379400895 9193728508 2948945508 1520277501 9618263204 5020519022 0518468630 7447368917 3387799105 5601346664 7592954559 9652998486 3187865744 6454877285 1956909426 0271026987 8621628765 6235875178 8425595811 8377486632 2746489950 0208870097 6919601719 0840250882 2826219646.

With those constraints, the unique continued fraction is: {0, 32, 4, 2, 2, 1, 2, 2, 1, 1, 1, 3, 1, 2, 1, 1, 1, 2, 1, 1, 8, 1, 8, 1, 1, 1, 17, 1, 5, 2, 1, 2, 4, 1, 12, 59, 1, 22, 1, 4, 7, 2, 3, 7, 1, 1, 2, 4, 22, 3, 1, 4, 23, 1, 13, 2, 3, 1, 1, 1, 1, 1, 1, 37, 3, 2, 3, 4, 11, 1, 1, 1, 1, 1, 6, 3, 1, 1, 1, 1, 2, 117, 3, 1, 2, 2, 2, 1, 11, 1, 2, 37, 1, 2, 1, 1, 2, 4, 4, 1, 1, 1, 1, 2, 2, 3, 1, 4, 3, 1, 2, 2, 6, 9, 2, 1, 6, 1, 1, 9, 1, 1, 9, 1, 2, 2, 15, 1, 5, 3, 3, 2, 2, 2, 2, 2, 25, 41, 1, 22, 2, 3, 6, 2, 1, 1, 27, 3, 2, 2, 1, 1, 2, 1, 4, 3, 3, 1, 2, 1, 3, 1, 21, 1, 98, 2, 2, 9, 2, 2, 2, 1, 9, 1, 4, 7, 1, 2, 14, 16, 3, 2, 14, 1, 72, 5, 1, 4, 4, 1, 36, 1, 2, 3, 1, 6, 1, 2, 1, 4, 1, 17, 1, 2, 1, 1, 2, 2, 1, 2, 11, 1, 6, 6, 1, 1, 1, 1, 5, 24, 5, 3, 2, 7, 1, 37, 4, 2, 1, 4, 2, 2, 1, 19, 1, 1, 45, 1, 2, 1, 5, 4, 3, 2, 15, 130, 1, 15, 1, 2, 1, 2, 8, 2, 1, 1, 1, 1, 1, 2, 1, 4, 5, 1, 2, 10, 1, 36, 1, 2, 1, 36, 1, 5, 3, 1, 13, 1, 1, 1, 2, 3, 3, 2, 2, 7, 3, 1, 4, 2, 2, 12, 2, 3, 12, 1, 6, 1, 1, 5, 1, 74, 1, 20, 3, 1, 2, 8, 1, 1, 1, 1, 1, 8, 2, 9, 5, 1, 1, 10, 1, 34, 3, 1, 13, 2, 2, 1, 2, 1, 1, 1, 91, 3, 105, 4, 4, 2, 3, 2, 10, 7, 1, 6, 4, 2, 1, 16, 1, 10, 1, 1, 3, 7, 6, 5, 5, 1, 1, 7, 1, 1, 2, 5, 10, 1, 1, 4, 1, 10, 2, 46, 3, 3, 2, 213, 9, 2, 4, 1, 16, 5, 2, 40, 2, 1, 1, 1, 17, 2, 1, 1, 9, 3, 1, 1, 1, 3, 1, 6, 4, 1, 6, 1, 1, 1, 1, 1, 1, 3, 11, 2, 3, 1, 4, 3, 6, 1, 2, 1, 23, 2, 1, 2, 1, 1, 17, 1, 1, 1, 3, 5, 5, 45, 1, 1, 4, 1, 1, 7, 1, 2, 2, 3, 4, 1, 1, 3, 2, 1, 2, 2, 8, 6, 1, 2, 4, 1, 54, 2, 2, 1, 3, 6, 1, 4, 1, 1, 3, 5, 2, 1, 1, 1, 7, 2, 2, 9, 1, 1, 2, 1, 2, 19, 2, 1, 2, 1, 6, 1, 11, 20, 1, 1, 2, 1, 2, 1, 12, 25, 2, 5, 1, 1, 3, 72, 2, 2, 1, 2, 4, 1, 1, 2, 1, 1, 11, 1, 3, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 2, 5, 1, 2, 24, 4, 2, 1, 1, 1, 1, 1, 4, 1, 2, 2, 1, 6, 3, 1, 1, 1, 1, 1, 3, 2, 1, 1, 2, 1, 35, 3, 13, 3, 6, 1, 1, 3, 1, 2, 2, 6, 3, 1, 2, 6, 1, 3, 1, 6, 5, 2, 2, 8, 1, 1, 2, 7, 1, 1, 4, 1, 10, 25, 1, 15, 1, 1, 4, 2, 22, 2, 2, 2, 1, 1, 1, 1, 5, 2, 25, 7, 8, 2, 2, 1, 8, 11, 2, 3, 4, 1, 30, 1, 4, 10, 1, 5, 1, 6, 7, 1, 4, 1, 4, 2, 2, 1, 2, 1, 6, 5, 5, 1, 4, 2, 2, 1, 14, 5, 7, 8, 1, 1, 1, 1, 3, 1, 1, 2, 3, 1, 1, 1, 1, 4, 1, 1, 1, 1, 4, 1, 1, 1, 2, 13, 1, 2, 1, 1, 4, 6, 2, 2, 2, 1, 3, 1, 1, 7, 2, 1, 1, 9, 1, 1, 1, 2, 2, 1, 1, 2, 5, 1, 21, 1, 1, 2, 3, 3, 1, 1, 221, 7, 1, 1, 9, 1, 18, 1, 5, 2, 3, 3, 4, 1, 3, 1, 1, 18, 2, 3, 3, 1, 1, 1, 2, 1, 1, 19, 2, 2, 1, 1, 2, 1, 2, 2, 2, 15, 2, 1, 1, 1, 1, 2, 1, 2, 118, 5, 1, 2, 2, 23, 9, 2, 2, 2, 2, 1, 6, 1, 9, 2, 1, 1, 8, 15, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 28, 3, 6, 4, 1, 1, 2, 1, 16, 1, 9, 1, 3, 1, 1, 1, 8, 26, 2, 2, 28, 1, 4, 13, 7, 3, 5, 2, 1, 4, 1, 1, 5, 3, 293, 1, 1, 1, 3, 1, 4, 7, 21, 1, 4, 1, 11, 1, 19, 4, 1, 1, 1, 4, 1, 1, 12, 2, 3, 22, 4, 1, 4, 3, 1, 1, 148, 15, 3, 1, 8, 34, 1, 8, 1, 1, 1, 9, 1, 1, 1, 33, 2, 1, 183, 1, 8, 1, 3, 1, 4, 1, 1, 2, 1, 1, 1, 1, 5, 4, 1, 4, 53, 1, 20, 1, 4, 15, 1, 3, 1, 1, 2, 4, 1, 7, 112, 3, 2, 7, 1, 4, 38, 1, 17, 1, 1, 3, 1, 64, 3, 2, 2, 150, 1, 2, 1, 3, 1, 1, 2, 2, 1, 7, 10, 2, 13, 1, 3, 1, 5, 2, 1, 1, 1, 2, 2, 2, 1, 3, 2, 2, 1, 2, 2, 2, 2, 7, 8, 3, 4, 1, 5, 1, 2, 1, 6, 3, 4, 1, 1, 5, 1, 1, 1, 2, 2, 75, 4, 2, 7, 110, 1, 7, 2, 8, 2, 7, 1, 2, 1, 1, 2, 1, 2, 22, 1, 10, 2, 3, 1, 1, 1, 1, 1, 55, 1, 21, 12, 5, 9, 1, 10, 1, 2, 1, 7, 1, 1, 2, 5, 1, 6, 2, 86, 96, 1, 3, 72, 2, 23, 1, 2, 1, 116, 1, 15, 1, 5, 7, 2, 7, 8, 4, 4, 1, 35, 3, 3, 1, 10, 3, 1, 1, 2, 6, 3, 6, 2, 4, 1, 1, 69, 1, 2, 2, 5, 1, 15, 1, 3, 1, 1, 1, 2, 1, 6, 12, 1, 7, 135, 8, 1, 1, 2, 163, 1, 5, 1, 1, 3, 10, 8, 3, 2, 153, 3, 1, 23, 1, 9, 4, 19, 1, 8, 1, 1, 3, 12, 9, 1, 1, 4, 1, 3, 1, 1, 8, 4, 1, 1, 1, 1, 1, 2, 5, 30, 1, 1, 4, 4, 2, 2, 2, 1, 8, 1, 1, 1, 1, 1, 19, 4, 3, 4, 1, 1, 31, 1, 18, 2, 1, 8, 2, 2, 2, 1, 10, 23, 1, 13, 12, 5, 1, 1, 17, 1, 12, 3, 3, 1, 3, 1, 1, 19, 1, 1, 3, 4, 1, 1, 3, 2, 36, 1, 2, 1, 8, 2, 1, 2, 2, 2, 1, 9, 2, 2, 17, 1, 1, 1, 3, 1, 2, 25, 4, 3, 3, 1, 2, 2, 1, 6, 12, 48, 1, 22, 1, 3, 1, 3, 1, 1, 7, 1, 3, 1, 5, 107, 2, 1, 6, 1, 4, 3, 1, 22, 1, 1, 1, 3, 3, 1, 1, 12, 4, 1, 1, 4, 4, 2, 2, 8, 2, 15, 2, 72, 34, 1, 1, 3, 1, 1, 3, 1, 5, 1, 12, 1, 6, 12, 1, 4, 2, 2, 1, 42, 1, 1, 2, 4, 1, 349, 1, 1, 36, 2, 1, 3, 1, 1, 1, 12, 3, 1, 5, 2, 16, 4, 1, 1, 1, 9, 1, 1, 1, 7, 7, 1, 1, 3, 1, 2, 6, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 5, 3, 25, 1, 3, 1, 3, 1, 24, 1, 2, 3, 1, 1, 1, 1, 1, 4, 2, 1, 13, 3, 1, 1, 1, 2, 20, 1, 69, 3, 1, 7, 3, 1, 1, 2, 1, 1, 1, 5, 2, 10, 1, 2, 1, 1, 2, 3, 11, 1, 1, 18, 1, 6, 7, 1, 5, 26, 1, 4, 7, 60, 4, 1, 23, 1, 1, 7, 1, 1, 1, 15, 2, 1, 3, 2, 1, 1, 1, 69, 7, 1, 1, 1, 14, 4, 1, 1, 1, 1, 1, 1, 3, 2, 1, 3, 1, 1, 3, 3, 20, 1, 1, 1, 2, 1, 1, 1, 45, 2, 2, 1, 2, 3, 3, 1, 3, 2, 12, 248, 1, 2, 5, 2, 4, 268, 9, 1, 5, 2, 157, 4, 1, 1, 15, 1, 4, 2, 15, 1, 7, 1, 5, 2, 4, 5, 2, 4, 2, 1, 2, 1, 4, 15, 1, 18, 1, 1, 61, 2, 1, 1615, 1, 1, 1, 57, 1, 21, 2, 6, 1, 1, 1, 1, 2, 1, 10, 1, 4, 4, 1, 2, 1, 1, 2, 2, 2, 1, 16, 1, 3, 1, 18, 1, 22, 1, 1, 12, 19, 1, 3, 2, 1, 11, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 21, 3, 1, 2, 2, 1, 1, 1, 62, 67, 2, 1, 2, 8, 1, 5, 21, 1, 2, 2, 57, 3, 3, 8, 1, 4, 2, 3, 3, 4, 1, 2, 20, 2, 2, 3, 2, 1, 54, 3, 12, 2, 3, 4, 1, 6, 2, 1, 3, 16, 4, 5, 24, 2, 8, 13, 1, 4, 1, 73, 1, 2, 3, 1, 2, 1, 7, 1, 1, 1, 3, 4, 1, 1, 4, 2, 7, 23, 1, 1, 2, 1, 1, 1, 1, 2, 6, 14, 2, 7, 1, 3, 1, 1, 4, 1, 1, 3, 8, 2, 3, 2, 3, 4, 1, 2, 2, 1, 7, 2, 2, 4, 2, 1, 1, 3, 1, 2, 1, 1, 2, 1, 5, 2, 2, 12, 1, 1, 2, 16, 21, 1, 12, 1, 2, 1, 1, 1, 5, 8, 1, 1, 21, 7, 16, 7, 9, 1, 13, 1, 689, 1, 10, 17, 2, 2, 10, 1, 5, 1, 25, 1, 54, 2, 1, 13, 1, 1, 1, 3, 1, 10, 8, 1, 2, 1, 2, 2, 1, 1, 1, 5, 1, 2, 4, 1, 2, 2, 251, 7, 1, 16, 1, 4, 85, 1, 184, 1, 1, 6, 1, 1, 6, 2, 1, 6, 2, 1, 3, 45, 1, 1, 2, 1, 4, 1, 4, 2, 8, 7, 21, 8, 1, 8, 1, 3, 5, 3, 1, 2, 2, 1, 46, 3, 1, 6, 1, 4, 18, 96, 1, 11, 2, 53, 14, 14, 1, 2, 1, 2, 1, 3, 6, 1, 3, 3, 1, 3, 2, 1, 4, 1, 4, 7, 17, 2, 1, 1, 2, 2, 6, 115, 1, 2, 1, 1, 6, 1, 80, 2, 10, 1, 2, 2, 19, 10, 1, 1, 2, 1, 6, 1, 8, 1, 1, 7, 1, 1, 15, 60, 1, 2, 5, 3, 1, 2, 1, 1, 3, 2, 1, 2, 1, 19, 3, 1, 14, 1, 1, 8, 1, 2, 4, 3, 5, 18, 1, 2, 1, 9, 6, 1, 2, 1, 18, 1, 7, 8, 1, 812, 1, 1, 1, 2, 1, 1, 7, 1, 4, 17, 8, 1, 1, 1, 3, 3, 10, 1, 2, 3, 2, 2, 36, 18, 3, 3, 4, 1, 1, 2, 3, 1, 1, 2, 1, 8, 65, 1, 2, 2, 1, 24, 2, 1, 2, 1, 4, 7, 1, 2, 1, 2, 10, 502, 1, 5, 3, 3, 2, 5, 4, 3, 1, 4, 23, 2, 5, 3, 4, 2, 1, 110, 1, 12, 2, 4, 1, 4, 2, 11, 1, 1, 3, 5, 58, 1, 6, 2, 9, 4, 1, 7, 2, 1, 1, 6, 1, 13, 9, 1, 1, 3, 4, 5, 4, 33, 1, 3, 1, 4, 9, 1, 3, 1, 21, 1, 1, 1, 4, 4, 2, 2, 1, 2, 1, 1, 8, 1, 3, 2, 1, 3, 1, 1, 1, 1, 30, 2, 1, 1, 1, 1, 4, 5, 1, 1, 2, 2, 13, 1, 1, 1, 1, 1, 2, 2, 1, 4, 3, 47, 2, 9, 16, 2, 1, 1, 1, 6, 1, 2, 2, 2, 1, 1, 3, 57, 17, 1, 1, 8, 9, 1, 1, 2, 58, 2, 1, 2, 1, 1, 724, 11, 34, 1, 9, 2, 3, 1, 8, 1, 1, 2, 14, 1, 1, 1, 10, 1, 9, 1, 3, 1, 4, 4, 1, 1, 1, 1, 1, 1, 2, 1, 5, 3, 1, 3, 1, 13, 3, 5, 1, 8, 2, 9, 3, 5, 11, 4, 1, 1, 263, 1, 7, 3, 2, 83, 1, 1, 1, 6, 1, 1, 3, 1, 1, 3, 1, 1, 2, 1, 2, 3, 1, 2, 3, 2, 36, 1, 17, 1, 2, 7, 3, 2, 1, 1, 28, 1, 1, 1, 1, 6, 30, 10, 1, 1, 3, 3, 1, 2, 1, 6, 1, 22, 2, 8, 8, 11, 2, 1, 3, 1, 1, 1, 3, 1, 1, 50, 1, 122, 7, 1, 3, 1, 12, 2, 1, 4, 1, 6, 1, 1, 1025, 5, 9758, 3, 1, 2, 749, 2, 2, 1, 2, 2, 3, 2, 4, 11, 22, 1, 2, 2, 2, 1, 2, 1, 1, 2, 1, 2, 7, 1, 1, 2, 1, 5, 1, 51, 1, 30, 1, 1, 2, 1, 1, 1, 4, 2, 1, 5, 2, 1, 4, 2, 1, 10, 1, 1, 1, 1, 34, 1, 2, 3, 1, 1, 1, 4, 1, 3, 1, 2, 1, 7, 1, 27, 2, 19, 1, 2, 31243, 2, 5, 4, 1, 3, 2, 4, 1, 13, 1, 4214, 1, 2, 1, 1, 15, 2, 10, 16, 1, 1, 8, 4, 4, 7, 3, 2, 3, 1, 15, 2, 1, 2, 37, 5, 1, 1, 10, 1, 8, 1, 1, 9, 1, 1, 3, 1, 1, 1, 8, 2, 1, 4, 1, 3, 1, 3, 1, 1, 3, 2, 5, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 5, 1, 3, 5, 7, 1, 14, 108, 2, 6, 1, 1, 1, 7, 1, 1, 9, 1, 1, 1, 1, 7, 1, 2, 1, 1, 4, 1, 2, 2, 1, 2, 2, 1, 1, 1, 6, 1, 2, 1, 2, 5, 1, 1, 1, 4, 4, 5, 2, 4, 26, 8, 1, 20, 1, 18, 3, 4, 13, 1, 3, 5, 1, 1}, giving these values for u and v:

… as taking one more term in the continued fraction would give a value of v of 17222 4960821100 1875329345 1641060925 5238676166 9512624449 5654203258 3287058835 7993718036 5656368644 5686476330 0326016702 1302933042 8844789710 3082679240 8559800117 5351195496 8668806670 8774204444 7321917492 2446902879 3474742217 1594539982 2214873779 4106636561 5697766375 4025231410 2158196759 4745313407 7201960598 2541792079 5381405688 7768669934 7364278363 2117760753 9503531102 6748393652 2184283279 5644312578 3248323650 9930718972 7321636486 5725819612 5908273159 2155399804 8283990220 7119577473 6213780572 8584200745 1099505196 6727674488 0951113996 5852642238 7671856809 2883368498 9377841958 9377931486 5906561731 6428839056 8863547993 2899516244 9846462182 1356579297 6065172826 1987873104 9701731837 4328093148 1567820847 7265866925 3019613631 5310099719 6818744044 4249869285 6505483402 4236264536 4095036966 4147358558 4378073033 9101223957 8521955350 7628502734 0632769271 6804123463 8740212338 8094058005 5518594811 4324781410 5280547740 2094302104 6885974156 6237622304 7007225836 5145605824 1375684286 0181832811 8139285553 1306800693 0396250147 9369654988 8205298438 0541881622 9716269330 9652955810 6538874320 5737081869 7922848848 0194694985 8021547724 6505330625 7710838049 9341949188 2258951119 8875105723 7218138618 1287475791 6676255453 5936252898 0779653976 5900693133 8546808135 3382938041 5013257129 5542618613 8364420063 6991300828, which is too large.

We also need to calculate d = floor(c4·v/F + 0.5) = 119463592 1167011388 6758372224 2746536979 1608618512 3076916301 3149476348 9188911521 9814829372 2758155269 5335247625 6623110634 8740076231 0210470076 7922996462 8084493261 2431962194 3047913621 4203986507 2787027915 3959648774 9649500812 2956181486 8280032546 4549063049 4842723313 7946252677 4426508163 2044675014 9035195925 7756961448 5235788073 7554025283 7108755474 2177447557 4193566210 7040047954 9638163676 7343432471 8312127866 6255670812 3540570700 6704362748 8804914772 6311135324 9117827073 8746464073 8070635013 3481879485 8932237053 6738087143 7555190315 8686035441 2475279168 6397856340 9381448683 5325110144 5186320902 8255690431 3942953630 3919668887 1018699854 0330510942 0950976625 4602293180 1387313179 5928784045 6410622833 7709187462 6942446434 8881442715 1147196299 2290861063 2967727500 2263500534 6957617792 8029470477 3475665347 0708831553 4578469620 5254872210 9160629277 1648056700 3555623911 7997957408 6639938192 0138903581 8835750488 3680184417 5591407564 7776770464 1067452914 8527392160 9991054492 5940762974 6401585887 3402632796 4678085962 3728701611 0631717309 8509337430 8058597682 3737768566 5653431037 4720537588 2076339261 1676465500 9041799294 3631364911 9329542468 2152251713 3757651702 5751960818 2052505111 9498495616 2176315686 5448199361 9611784085 6087408830 1327361712 8133264675 8746504415 8533456744 3582946546 8298227848 0849471831 0485614235 4142151398 7200986939 7223192937 9536133833 7759175676 1261693588 2218700110 3712504381 4626323165 7528312251 0658621434 6884714102 9537424040 8541413713 6954610094 0854926183 3541325171 2975804166 2108855916 3106581204 7050011327 2381135837 4344382206 0038957732

Cubic Polynomial

We now consider the cubic P(x)= v·x3 + (u·F-c1·v)·x2 + (c4·v-d·F+u)·x - d, which we express as: z1·x3 + z2·x2 + z3·x + z4, where:

We need to prove that this cubic has no integer roots r such that r·F+1 is a non-trivial factor of N. Clearly r (if it exists) must lie between 1 and R.

As z1 is > 0 and z4 is < 0, we know that P has at least one real root r1 > 0. This is easily found:

Note that the root is not an integer but actually lies in the interval (r1,r1+1).

To see if P has any more real roots, we examine the quadratic derivative P' of P = 3·z1·x2 + 2·z2·x + z3, which we express in the usual quadratic form as a·x2 + b·x + c

As b2-4·a·c is negative, P' has no real roots. Therefore, P has no turning points and is monotonic, implying that r1 is the only real root of P.

There are no integer roots of P in the interval (1,R), so the proof of primality is complete.