Primality Certificate for (2302^2521-1)/2301

Andy Steward8,473 digits20 June 2001
Originally by David Broadhurst 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 31.907740% factorization of N-1:

From Factorisation
23022 · 1151
Φ27 · 7 · 47
Φ33 · 811 · 2179
Φ45 · 53 · 19997
Φ511 · 11 · 31 · 7489673981
Φ65296903
Φ771 · 127 · 617 · 7757 · 14057 · 245407
Φ889 · 315523180153
Φ93 · 49603310388729603091
Φ10821 · 26891 · 1271401
Φ1237 · 758961019849
Φ147 · 491 · 8807 · 4914006975277
Φ1542365221 · 18605629887167302831
Φ1819 · 433 · 18087994547440291
Φ205 · 10180601 · 15491699038220741041
Φ2143 · 463 · 436507 · p31
Φ24477146473 · 1652687857992236617
Φ2829 · 29 · 7962439213 · 3306904277648272723755673369
Φ3061 · 211 · 61294129656420890214961
Φ354492562551 · 3350959951258967245389605891 · p44
Φ36p41
Φ4041 · 177286245371021401 · p35
Φ4262581 · 508962843793260487 · 695542822257346477
Φ4523041 · 150211 · 10527481 · 470624136833881 · p50
Φ5617261696875186793 · 150297658265712554297 · p45
Φ60181 · 421 · 999434094040012021 · p31
Φ631559895333559583077 · 183955808847879105629257 · p80
Φ7030447971 · p74
Φ7273 · 577 · 685282921006969 · 25671748543489581332617 · p39
Φ846637 · 70981 · 20703481 · 38038126463653 · 3147987595792817404453 · p30
Φ9077041 · 24432153777427073288285818282636230571 · p39
Φ105c162
Φ120241 · 55441 · 3671401 · 10321950469801 · 728942791367665148372835271921 · p52
Φ12622807 · 39953065573 · 575045672041 · p95
Φ1403500075772398892022841 · p140
Φ1681009 · 36457 · 1756273 · 85990783369 · c137
Φ1801221095341 · c153
Φ210p162
Φ2522017 · 44101 · 5443404341331976842844025750373686101 · p198
Φ280281 · 499801 · 647081 · 549391081 · c301
Φ315372606817273934671 · c467
Φ3604109550481 · c314
Φ42045361 · c319
Φ504743401 · 1648081 · p473
Φ630631 · 284131 · c476
Φ840370441 · 34032601 · 396714361 · c624
Φ126030241 · 438054121 · 854172901 · 996360121 · 1436998501 · 7491378542761 · c916
Φ25202521 · 104317438681 · 14140459851739201 · c1907

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 :

113 4917479097 5937747626 4428001622 5296522562 7699957459 5273233600 1592020448 0607450789 8082782423 1682047883 0766212817 8225735503 4392709324 8673927315 2296779605 8546467928 4494084099 0572993429 4622088255 9136123746 8730545289 5057196736 5018517458 4910548964 7089321770 2322311625 0891279977 0335026851 2560540920 3577990701 8465000526 8312332861 8773404122 2552280171 5158962838 0400749843 2941062653 9820982797 2661254298 3577675004 6523786826 8298246558 4107834951 6144094636 1038446422 0265036458 6979781641
24353313 4349168427 1822489347 3796546414 2288099353 5840705303 1423636809 2531726285 7033608894 6702882194 8736957537 6230897078 0649430144 0404471189 3390443056 1428473330 6806149052 9094392304 3516846248 3427963833
24 0362497245 6598590963 4233704519 3221650006 4226234632 7544915703 4387288735 3570640348 2179904689 0685726031 6483131400 9224094464 8288735032 7673529018 2770673783 7705711271
6870335754 9689184500 5782104683 5874385118 3205750567 8421525480 3076653526 0203524161 0788752342 7449787791 2078391887 1327108615 3620973349 8966688301
20723 8712168670 0057793756 3218005436 1446376295 5937604836 8642358590 0265489754 1466184543 0782667979
3784273611 9344572769 1758908612 6667590210 6217696316 6674968136 4880570146 5183278381
1611 2314278599 4402096746 8553802115 1090037267 2468140512 1649403178 1542016981
10 4769550722 8254134339 5011102018 5258509311 7308506961
2859735972 8806275199 6570133025 5169940067 4993482131
18901 3125278770 9684567830 2893841021 3067411601
3255 9339964621 1918849575 2079039316 4554520551
2 2144395610 0552559729 7805336852 6465847233
661764752 0698803365 2233391615 2110188537
260521757 0526158393 6416626835 0424481691
24432153 7774270732 8828581828 2636230571
5443404 3413319768 4284402575 0373686101
85551 2204616971 0402947889 2393393841
8 1652664142 2954995143 7134600361
2 5470324128 9663005509 8558452093
7289427913 6766514837 2835271921
4198733511 1316653089 4128196117
33509599 5125896724 5389605891
33069042 7764827272 3755673369
1839 5580884787 9105629257
612 9412965642 0890214961
256 7174854348 9581332617
35 0007577239 8892022841
31 4798759579 2817404453
1 5029765826 5712554297
4960331038 8729603091
1860562988 7167302831
1549169903 8220741041
165268785 7992236617
155989533 3559583077
99943409 4040012021
69554282 2257346477
50896284 3793260487
37260681 7273934671
17728624 5371021401
1808799 4547440291
1726169 6875186793
1414045 9851739201
68528 2921006969
47062 4136833881
3803 8126463653
1032 1950469801
749 1378542761
491 4006975277
75 8961019849
57 5045672041
31 5523180153
10 4317438681
8 5990783369
3 9953065573
7962439213
7489673981
4492562551
4109550481
1436998501
1221095341
996360121
854172901
549391081
477146473
438054121
396714361
42365221
34032601
30447971
20703481
10527481
10180601
5296903
3671401
1756273
1648081
1271401
743401
647081
499801
436507
370441
284131
245407
150211
77041
70981
62581
55441
45361
44101
36457
30241
26891
23041
22807
19997
14057
8807
7757
6637
2521
2179
2017
1151
1009
821
811
631
617
577
491
463
433
421
281
241
211
181
127
89
73
71
61
53
47
43
41
37
31
292
19
112
73
52
32
2

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) = 31.907740%

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 = 3 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 = 3 2905435593 4186080758 0653292535 8527401211 0756871044 3957256523 1085502474 3086200903 3202582938 2508695305 8633176495 1952381116 1235508123 4375735251 3602129542 9559663341 7347072274 7903517102 1831672922 8673652372 2061193453 0525641456 9977441927 2907982183 4639995171 3208996192 1783426965 8664368008 7385032511 4714756188 1130089901 0527147928 5643911334 7915560149 7069275068 5655205403 5582371434 9836520857 4988278012 9947613439 7084072840 8862847396 7792614589 6248412068 2049956407 6676154528 2915083737 9082427899 7961879225 4748445233 7661731390 3107214840 5973546840 7980809705 7393885548 2105267788 9608346479 7336982173 7567598505 8862616294 9634026515 8020802512 1780588929 5009556221 1413652955 6092122758 5761494861 1388316541 8160239222 9287602867 8657954588 8651748513 5470787572 9956432019 0594767816 0507074970 9817177067 5858739932 7521552871 7550245562 4662501532 1134705471 0642791170 6962314942 7102527974 0030980762 6397299261 6304369695 1628645568 6494587259 4970155133 6178288130 2833740675 8599546510 4831967556 0584863881 3172066691 3739041477 5008685317 2783505087 4191246221 8442008746 0590179258 5240655486 4112366393 1726236307 8972734166 1434408789 4224253514 6315464166 2726435157 2824659899 6533548354 2023706720 5277921111 2402629573 3719172313.

With those constraints, the unique continued fraction is: {0, 2, 4, 1, 1, 1, 2, 1, 305, 1, 2, 1, 10, 2, 1, 1, 1, 2, 2, 12, 2, 2, 4, 6, 2, 1, 12, 1, 14, 1, 1, 3, 2, 2, 1, 6, 1, 1, 2, 3, 3, 34, 2, 1, 1, 1, 1, 3, 1, 3, 5, 177, 1, 1, 1, 1, 3, 2, 1, 17, 1, 2, 10, 1, 1, 10, 5, 6, 2, 5, 17, 1, 1, 6, 2, 2, 2, 8, 1, 1, 1, 19, 41, 1, 1, 1, 1, 19, 3, 1, 3, 7, 1, 1, 2, 1, 9, 1, 2, 14, 3, 1, 1, 5, 2, 2, 11, 1, 131, 5, 4, 4, 1, 5, 3, 5, 2, 1, 1, 1, 3, 1, 4, 4, 4, 5, 3, 1, 1, 1, 1, 6, 1, 4, 1, 9, 1, 12, 1, 24, 1, 5, 1, 1, 1, 16, 3, 1, 1, 1, 6, 1, 5, 2, 1, 1, 1, 8, 1, 4, 3, 4, 1, 1, 1, 1, 2, 7, 1, 7, 1, 2, 1, 1, 26, 1, 30, 1, 9, 3, 1, 5, 3, 2, 1, 7, 2, 1, 2, 17, 1, 2, 1, 1, 6, 1, 1, 1, 26, 1, 16, 1, 1, 5, 6, 5, 16, 1, 4, 1, 1, 35, 4, 3, 3, 1, 1, 1, 5, 3, 9, 1, 1, 2, 1, 1, 4, 385, 5, 1, 2, 53, 17, 1, 1, 1, 1, 1, 1, 7, 1, 1, 3, 3, 3, 1, 1, 1, 6, 1, 1, 1, 1, 2, 10, 8, 3, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 1, 1, 1, 3, 9, 22, 1, 2, 5, 1, 18, 80, 5, 2, 4, 4, 109, 1, 6, 10, 3, 47, 2, 1, 1, 7, 1, 2, 1, 2, 3, 3, 26, 3, 11, 1, 3, 1, 4, 1, 1, 3, 468, 3, 1, 47, 1, 1, 8, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 4, 1, 3, 2, 84, 2, 1, 1, 1, 2, 7, 6, 3, 1, 1, 1, 2, 1, 2, 3, 2, 9, 2, 3, 1, 144, 2, 1, 2, 4, 3, 6, 2, 1, 14, 1, 1, 3, 1, 1, 1, 8, 1, 1, 1, 3, 2, 1, 1, 1, 1, 1, 5, 9, 1, 2, 2, 1, 6, 1, 15, 9, 3, 1, 2, 1, 2, 1, 1, 7, 4, 2, 1, 8, 5, 1, 1, 1, 7, 2, 1, 21, 1, 2, 2, 5, 4, 1, 1, 5, 1, 7, 1, 2, 1, 2, 1, 1, 9, 17, 1, 1, 1, 2, 1, 16, 1, 4, 50, 10, 1, 1, 1, 3, 2, 1, 1, 1, 2, 4, 1, 2, 1, 1, 2, 16, 712, 1, 1, 3, 1, 4, 2, 14, 1, 1, 32, 2, 14, 2, 3, 1, 2, 1, 14, 3, 5, 1, 103, 14, 1, 1, 2, 1, 4, 29, 1, 2, 2, 1, 6, 14, 1, 6, 1, 4, 6, 1, 14, 1, 2, 2, 8, 1, 4, 2, 13, 1, 10, 1, 2, 101, 4, 1, 5, 1, 2, 16, 1, 13, 2, 1, 1, 3, 2, 11, 2, 47, 5, 1, 20, 3, 1, 245, 1, 11, 1, 1, 1, 3, 2, 3, 16, 2, 1, 3, 1, 1, 2, 1, 1, 6, 2, 1, 6, 1, 34, 1, 10, 1, 1, 4, 3, 2, 1, 2, 16, 5, 1, 5, 1, 8, 1, 3, 4, 70, 2, 1, 2, 1, 2, 2, 4, 1, 1, 1, 1, 1, 4, 2, 2, 49, 25, 1, 38, 1, 3, 1, 1, 3, 3, 1, 1, 2, 1, 1, 2, 3, 2, 12, 1, 6, 1, 3, 1, 17, 13, 6, 1, 1, 1, 6, 2, 5, 1, 1, 1, 1, 2, 25, 1, 9, 1, 30, 4, 1, 16, 1, 2, 1, 5, 2, 4, 1, 2, 3, 1, 7, 1, 6, 2, 17, 2, 16, 1, 10, 9, 2, 5, 30, 4, 1, 35, 15, 4, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 3, 6, 2, 15, 4, 2, 3, 11, 2, 1, 2, 1, 56, 2, 1, 3, 1, 2, 3, 9, 48, 1, 9, 1, 1, 2, 58, 1, 7, 1, 1, 2, 3, 104, 14, 4, 1, 1, 1, 27, 1, 4, 4, 2, 1, 5, 1, 2, 3, 1, 4, 4, 4, 2, 24, 1, 2, 2, 1, 2, 1, 13, 1, 3, 3, 2, 122, 1, 2, 1, 3, 8, 1, 142, 3, 2, 2, 1, 3, 1, 2, 4, 1, 10, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 3, 1, 18, 7, 1, 5, 15, 1, 2, 1, 1, 1, 2, 4, 55, 5, 1, 31, 1, 18, 2, 8, 59, 1, 1, 1, 49, 1, 2, 1, 5, 32, 2, 1, 46, 1, 10, 194, 21, 309, 2, 11, 1, 12, 13, 1, 4, 61, 69, 2, 1, 1, 4, 1, 19, 13, 3, 1, 1, 6, 1, 23, 1, 2, 1, 1, 3, 1, 3, 5, 1, 4, 2, 1, 3, 4, 1, 1, 1, 1, 1, 5, 2, 2, 1, 1, 3, 1, 1, 1, 4, 3, 2, 1, 12, 2, 2, 118, 11, 1, 2, 3, 1, 2, 8, 2, 1, 1, 3, 2, 2, 3, 23, 1, 1, 1, 1, 1, 3, 2, 1, 174, 1, 13, 3, 1, 26, 1, 1, 4, 1, 1, 1, 8, 4, 1, 4, 1, 4, 1, 19, 2, 71, 1, 29, 17, 2, 50, 1, 9, 1, 1, 30, 1, 2, 1, 4, 1, 4, 23, 1, 1, 6, 1, 3, 1, 9, 55, 2, 1, 4, 1, 1, 1, 1, 4, 1, 2, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 3, 1, 3, 1, 8, 8, 6, 13, 7, 27, 3, 1, 3, 2, 42, 3, 4, 3, 1, 4, 2, 2, 3, 7, 1, 1, 1, 1, 5, 2, 3, 1, 2, 1, 22, 1, 2, 80, 3, 1, 8, 2, 35, 2, 1, 1, 1, 1, 1, 1, 12, 1, 1, 1, 1, 1, 2, 2, 3, 8, 3, 1, 1, 65, 28, 1, 2, 2, 34, 15, 7, 2, 46, 7, 4, 1, 1, 12, 14, 6, 1, 4, 1, 4, 1, 44, 1, 1, 1, 116, 11, 27, 1, 2, 1, 1, 11, 2, 4, 1, 1, 5, 5, 2, 3, 1, 1, 2, 2, 4, 28, 1, 30, 10, 4, 3, 1, 7, 2, 3, 1, 10, 5, 1, 2, 1, 11, 7, 3, 1, 9, 1, 49, 1, 4, 3, 2, 1, 7, 1, 1, 1, 6, 1, 2, 3, 3, 3, 2, 5, 2, 1, 1, 1, 13, 1, 1, 2, 3, 4, 1, 1, 2, 1, 3, 2, 3, 2, 1, 3, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 15, 1, 1, 1, 1, 5, 1, 1, 5, 3, 1, 3, 1, 3, 21, 1, 17, 1, 5, 4, 1, 1, 1, 1, 7, 5, 3, 18, 1, 3, 1, 2, 1, 1, 4, 2, 1, 3, 3, 11, 8, 3, 2, 2, 4, 3, 1, 2, 3, 12, 1, 26, 2, 1, 6, 1, 2, 2, 4, 1, 52, 1, 1, 10, 28, 1, 1, 1, 38, 86, 8, 1, 1, 404, 2, 2, 1, 5, 2, 1, 2, 7, 451, 7, 4, 28, 1, 1, 8, 2, 8, 2, 1, 2, 1, 1, 1, 4, 1, 37, 1, 162, 2, 1, 1, 8, 75, 3, 1, 12, 2, 7, 1, 28, 11, 1, 56, 25, 1, 2, 1, 3, 1, 1, 1, 243, 1, 7, 2, 2, 2, 1, 2, 4, 1, 1, 2, 1, 1, 1, 2, 1, 4, 2, 4, 1, 1, 13, 2, 4, 9, 2, 21, 31, 1, 2, 1, 1, 3, 2, 6, 3, 3, 2, 1, 31, 1, 4, 2, 3, 1, 1, 3, 2, 1, 2, 6, 2, 2, 1, 1, 4, 7, 3, 22, 1, 15, 1, 3, 2, 2, 2, 2, 1, 1, 1, 1, 2, 1, 5, 2, 1, 5, 3, 4, 161, 1, 1, 15, 1, 1, 1, 4, 1, 1, 1, 7, 2, 49, 1, 25, 1, 1, 1, 1, 1, 77, 1, 2, 2, 1, 1, 1, 14, 1, 4, 3, 1, 1, 3, 4, 1, 1, 3, 4, 1, 5, 14, 4, 1, 4, 4, 1, 1, 53, 4, 1, 1, 5, 1, 118, 4, 19, 1, 3, 1, 7, 1, 1, 1, 3, 1, 1, 1, 1, 41, 2, 5, 6, 3, 1, 1, 2, 3, 6, 2, 8, 1, 3, 3, 1, 1, 8, 4, 1, 2, 2, 1, 8, 2, 3, 1, 1, 44, 1, 3, 2, 1, 1, 3, 2, 6, 1, 1, 1, 4, 16, 79, 1, 1, 1, 3, 2, 1, 2, 1, 2, 4, 6, 2, 6, 1, 17, 2, 12, 4, 1, 1, 1, 4, 1, 9, 2, 3, 1, 1, 1, 1, 3, 1, 5, 2, 7, 1, 4, 5, 1, 1, 2, 1, 14, 6, 1, 1, 1, 6, 1, 1, 5, 1, 3, 4, 2, 2, 1, 4, 6, 7, 2, 2, 1, 1, 16, 1, 1, 1, 3, 3, 4, 2, 1, 3, 1, 14, 5, 3, 4, 1, 3, 2, 1, 4, 6, 2, 2, 1, 1, 15, 1, 1, 5, 11, 4, 1, 1, 1, 1, 2, 1, 1, 15, 2, 2, 1, 5, 1, 6, 15, 3, 4, 249, 5, 1, 1, 1, 1, 2, 1, 9, 1, 1, 2, 25, 2, 3, 3, 3, 2, 1, 1, 1, 18, 2, 1, 35, 4, 2, 5, 1, 3, 1, 6, 1, 2, 1, 2, 6, 4, 3, 4, 3, 44, 1, 3, 2, 1, 13, 2, 4, 1, 24, 3, 1, 3, 1, 1, 2, 3, 1, 2, 96, 1, 3, 7, 1, 14, 1, 4, 1, 1, 1, 3, 1, 1, 3, 28, 1, 17, 6, 18, 1, 5, 12, 5, 5, 5, 1, 1, 1, 1, 5, 2, 2, 8, 2, 92, 3, 6, 6, 1, 3, 5, 1, 5, 1, 8, 1, 4, 1, 11, 1, 2, 1, 4, 4, 5, 1, 1, 1, 4, 1, 4, 56, 2, 1, 4, 1, 1, 1, 4, 3, 3, 1, 118, 2, 8, 1, 8, 2, 1, 3, 20, 1, 5, 3, 2, 12, 1, 4, 1, 2, 2, 2, 5, 1, 4, 15, 1, 3, 4, 1, 10, 2, 2, 2, 3, 2, 7, 3, 4, 5, 2, 19, 2, 5, 32, 3, 2, 3, 3, 5, 1, 5, 2, 3, 7, 2, 1, 1, 48, 1, 1, 1, 2, 2, 5, 3, 1, 8, 1, 1, 24, 59, 12, 1, 1, 1, 3, 18, 1, 1, 3, 5, 1, 2, 3, 4, 1, 7, 3, 6, 1, 163, 1, 1, 53, 1, 6, 7, 2, 18, 1, 2, 1, 1, 1, 25, 1, 14, 11, 1, 3, 1, 13, 23, 5, 9, 1, 23, 1, 2, 12, 4, 1, 7, 1, 1, 1, 7, 1, 2, 1, 1, 1, 1, 1, 3, 3, 2, 1, 2, 1, 3, 1, 2, 2, 2, 642, 1, 1, 9, 10, 2, 5, 97, 3, 1, 1, 3, 1, 1, 12, 2, 9, 2, 4, 1, 5, 3, 2, 1, 4, 2, 2, 20, 2, 41, 5, 3, 1, 1, 2, 1, 3, 6, 22, 1, 1, 3, 15, 1, 5, 1, 1, 1, 1, 1, 4, 1, 1, 3, 2, 1, 1, 4, 10, 1, 1, 4, 1, 54, 3, 1, 2, 2, 6, 1, 41, 15, 21, 1, 2, 1, 1, 7, 1, 58, 1, 1, 1, 5, 1, 1, 3, 2, 2, 20, 1, 2, 3, 6, 17, 3, 147, 1, 1, 6, 1, 33, 1, 17, 20, 4, 1, 1, 2, 6, 2, 2, 5, 2, 1, 1, 1, 9, 1, 4, 4, 2, 13, 9, 1, 2, 9, 1, 3, 1, 2, 11, 1, 13, 1, 5, 3, 2, 2, 1, 6, 4, 17, 3, 1, 5, 5, 1, 1, 2, 1, 1, 11, 1, 3, 8, 1, 4, 2, 1, 10, 1, 9, 1, 1, 2, 3, 11, 1, 1, 5, 13, 7, 1, 4, 10, 1, 4, 3, 2, 1, 3, 1, 1, 1, 11, 3, 1, 4, 1, 1, 3, 1, 1, 6, 14, 2, 2, 1, 14, 1, 14, 1, 2, 3, 17, 1, 1, 3, 12, 1, 1, 1, 1, 18, 2, 5, 1, 25, 2, 1, 6, 1, 14, 1, 6, 9, 1, 6, 1, 3, 1, 10, 1, 252, 1, 34, 5, 2, 2, 3, 1, 1, 1, 2, 1, 1, 17, 2, 1, 19, 5, 22, 9, 3, 1, 6, 2, 4, 5, 4, 12, 9, 1, 3, 90, 1, 16, 1, 3, 40, 1, 5, 3, 5, 1, 1, 2, 2, 4, 1, 3, 1, 1, 7, 3, 2, 3, 5, 2, 1, 1, 2, 3, 2, 1, 3, 1, 2, 1, 1, 2, 6, 1, 1, 4, 1, 6, 9, 2, 13, 2, 2, 1, 1, 1, 1, 1, 2, 166, 84, 2, 1, 1, 3, 1, 3, 1, 1, 2, 2, 5, 7, 3, 1, 6, 11, 11, 2, 2, 2, 2, 18, 1, 14, 4, 4, 9, 1, 8, 3, 2, 3, 1, 4, 1, 14, 1, 1, 2, 2, 1, 6, 2, 2, 6, 2, 7, 1, 1, 1, 1, 23, 1, 2, 1, 47, 12, 1, 1, 1, 8, 5, 2, 1, 8, 4, 5, 1, 1, 2, 12, 1, 15, 1, 1, 1, 1, 1, 9, 1, 1, 3, 1, 4, 1, 2, 1, 1, 17, 5, 9, 5, 1, 6, 1, 26, 1, 1, 1, 4, 7, 1, 1}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 346 3355319241 6969596276 1202000662 5019118849 3612421468 7394416773 5849493030 4787786099 2033108738 4654202096 0379419678 4525036503 9617598287 5868934516 0109958443 8951607229 0375773865 8178675803 3976444506 7383398178 8710031126 6155507714 0076890016 2865688771 4879152603 6921062246 7535542554 1938443046 9179010861 4195324677 2876437039 8193860903 9216916355 0260700938 6362431204 8008206969 2653073441 7052345234 8249624415 7282998794 8724302344 8562884677 7885380079 5542191953 6744076746 5516462501 8808483793 6147643539 4620850521 6119858629 2341988015 6825553690 5889641014 1272956192 9429266477 9722594860 3618273236 0090587534 8559566271 8547132786 8196246983 1492350569 8523579421 9235501198 9687318775 7894184297 1149728692 2562023187 5233379392 3764139880 5031885831 9422758584 2498176191 1729268451 6705864667 3474148207 1752048521 8252708537 6051117141 3211448192 0693493152 6885096687 3254384632 4451289933 3505331819 6660556582 5077444646 6142427509 5282136802 7510146805 5842791971 3958312833 2093678158 2867544298 4772693964 3743565405 6697655300 1078830802 5433286853 1822758438 6084772468 2533409915 9530661986 0519142765 1612193967 6334744103 2906938129 5764607685 7926384599 1122627894 5277792031 1176779550 3089596291 7916346975 8468717319 8412853155 6758762631 8119319392 3861742572 0254220298 0339247902 2672999578 3474935048 4828400809 5050292091 2202413764 7681272667 9956367874 6939489517 0115360235 3654594683 7923146433 0849284512 9482283948 8726734727 4877159676 4294052312 1471629384 2100719257 5612803436 2295412043 8875423624 2071997867 7862511365 3225663730 5896236310 4795199446 9307629560 4383402197 1197332051 5864526679 5606514193 7992198621

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.

P has a single real root at:

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