Primality Certificate for (4676^2161-1)/4675

Andy Steward7,927 digits19 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 30.471170% factorization of N-1:

From Factorisation
46762 · 2 · 7 · 167
Φ23 · 1559
Φ313 · 1682281
Φ4769 · 28433
Φ55 · 95635887595601
Φ63 · 31 · 235057
Φ81928401 · 247913777
Φ973 · 51517 · 2779544235599533
Φ106681971 · 71532031
Φ1237 · 61 · 373 · 567881341
Φ1513591 · 64081 · 68281 · 3842576893374451
Φ16245521 · 930909314133965966807537
Φ183 · 19 · 2719 · 99595261 · 677211427
Φ2041 · 214021 · 5274871141 · 4937918294801
Φ24601 · 4969 · 129230833 · 592224551136913
Φ27394363 · p61
Φ3012841 · 589291 · 12367501 · 2442745820371
Φ36397 · 1621 · 713175647721157 · 238080597085698833425789
Φ4012606641 · 58870198061808461758921 · 70387762066132286074419597641
Φ45271 · 91801 · p81
Φ481153 · 4657 · 369798053459713 · p38
Φ543 · 379 · 151471 · 346949538189407741998112581 · p32
Φ60181 · 241 · 4546864951261 · 330599250638123581 · 796677445630839294989941
Φ721297 · 1566842587030560313 · 164633667590254017768725350094833 · p35
Φ80291521 · p112
Φ90642151 · 483957897686033121646482179938110938011 · p44
Φ108109 · 541 · 2175661 · 109434673 · 407855521 · p105
Φ12035521 · 437641 · c108
Φ135105352381 · 100851538321 · 857830042652700692187198950941 · c216
Φ1442161 · 4553281 · 2758249801214497 · 20117173506396923809 · 37011146854165963055495791489 · p103
Φ1805749498395838723390945021 · 727931220375936443876133149680021 · p119
Φ21638449 · 169129 · 4459463641 · c245
Φ24010066769281 · 66745390081 · p215
Φ270850468681 · c256
Φ36015117481 · 150017486537161 · c331
Φ432433 · 1166833 · 7412942881 · c510
Φ540256786741 · 1652065201 · p511
Φ720c705
Φ10806481 · 317043064405056831117841 · c1030
Φ2160c2114

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 :

6 8285513514 0925463595 7529983899 0087060443 9587751802 0431827209 6730915405 3180960080 3270986534 0625562578 8106446059 7303355031 6280683044 9101949436 7925093810 6121620490 5943756193 3429965752 7649008664 9273128174 8770233018 7576626896 5913305865 6921586834 4088152690 4272330839 3325708772 0579365361 9044322077 3187378678 8716511984 9663674468 7570116872 5886564785 2762362993 2985549675 2330095043 0824181797 7090390886 3722877158 6603335859 2785504092 0794284128 4515425872 7319964324 2210518736 5667303587 3524455270 0384967078 1287287328 8035177261
11082 9855909037 3448625202 2904660653 6638747199 6355450285 0591981122 9277943105 5074602245 9986692778 5366690554 0573458437 2356002896 9678926979 6359560645 4784823802 9152132382 7835400707 9093269020 5159813769 0244565319 6375373441
340608762 2085176660 7573927049 0659817415 9777648473 2476986969 4795110748 1528738873 8984932707 1329180724 6474897390 4703783361
93 6082733672 8467464219 7906311739 4586731169 1741604209 7709237436 2171523221 3663922368 5636466110 2059042468 3408399681
22782 7324291695 8369790019 3792517639 1418094774 2434722725 2932217246 4689380024 9745563118 9198316555 2823560733
705 4460396680 1433091792 6967607290 6713831326 7736980260 2540277210 3351116369 3693527802 7546483469 8036874513
4 7992277137 4044473545 6123373547 9179138418 7594415848 1603544824 5608281614 5789415431
2 8963089286 4509515606 9998500569 0388343695 7980863333 2897495731
3841 8757484686 8510125032 5300113138 2379775541
483957897 6860331216 4648217993 8110938011
26308244 7679402464 9335289518 1204288337
35686 5088088352 4056251392 0172704777
727 9312203759 3644387613 3149680021
164 6336675902 5401776872 5350094833
19 1154608935 2612550546 8925372723
8578300426 5270069218 7198950941
703877620 6613228607 4419597641
370111468 5416596305 5495791489
3469495 3818940774 1998112581
57494 9839583872 3390945021
9309 0931413396 5966807537
7966 7744563083 9294989941
3170 4306440505 6831117841
2380 8059708569 8833425789
588 7019806180 8461758921
2011717350 6396923809
156684258 7030560313
33059925 0638123581
384257 6893374451
277954 4235599533
275824 9801214497
71317 5647721157
59222 4551136913
36979 8053459713
15001 7486537161
9563 5887595601
493 7918294801
454 6864951261
244 2745820371
10 0851538321
6 6745390081
1 0066769281
7412942881
5274871141
4459463641
1652065201
850468681
677211427
567881341
407855521
256786741
247913777
129230833
109434673
105352381
99595261
71532031
15117481
12606641
12367501
6681971
4553281
2175661
1928401
1682281
1166833
642151
589291
437641
394363
291521
245521
235057
214021
169129
151471
91801
68281
64081
51517
38449
35521
28433
13591
12841
6481
4969
4657
2719
2161
1621
1559
1297
1153
769
601
541
433
397
379
373
271
241
181
167
109
73
61
41
37
31
19
13
7
5
34
22

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

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 = 2 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 = 24634979 7527326637 1993469466 6535570919 5001035731 9119012635 1257139384 2859950235 8233061294 4332762903 2016881534 6599101750 7132073578 6603043752 7234106780 9688149287 1800002951 3232739445 4507895133 9341466855 4795559379 3702492985 4010007651 7841932925 7707107156 4234967346 0668830854 9608211485 0080414679 2996919340 6290140447 7264824823 6170513522 3556309953 4064037502 4579554734 9583074603 7582162036 6729011453 6005424127 5932890849 7790638715 3797632912 2196148813 8212483589 2528504164 6356144542 6206294957 2838864737 0251134983 6292150435 7207435194 6595727664 0554822322 6710444533 7542150220 2340996195 3728363149 9394616878 1062613309 1211099373 4635431065 0111990527 0655589153 4626738071 2817502016 1277588298 5000283783 4921361440 5454398487 4222410253 4010803525 0228977977 5500407412 6303847085 3092498986 7856952952 8792565455 0374600703 6066895044 2874088430 6778944909 3662706177 9776317629 2028849411 7707506020 4358910053.

With those constraints, the unique continued fraction is: {0, 1, 5, 1, 2, 3, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 2, 1, 19, 1, 12, 1, 6, 1, 17, 3, 2, 1, 1, 3, 1, 1, 2, 1, 1, 2, 8, 1, 19, 1, 10, 3, 8, 1, 12, 1, 1, 1, 11, 2, 2, 3, 1, 3, 1, 10, 1, 1, 5, 20, 1, 3, 3, 10, 7, 2, 1, 1, 1, 2, 1, 1, 20, 2, 1, 1, 2, 4, 2, 2, 5, 2, 3, 3, 5, 1, 2, 1, 1, 1, 12, 244, 5, 4, 10, 2, 2, 1, 13, 1, 1, 1, 1, 3, 13, 1, 13, 1, 1, 6, 1, 2, 1, 5, 2, 3, 4, 2, 9, 1, 1, 4, 1, 1, 7, 1, 35, 2, 3, 2, 1, 1, 5, 6, 3, 1, 122, 1, 5, 3, 3, 1, 4, 1, 29, 2, 2, 3, 6, 1, 13, 2, 1, 2, 5, 1, 1, 13, 1, 2, 1, 3, 4, 2, 3, 3, 1, 42, 12, 3, 19, 2, 1, 24, 1, 6, 1, 2, 1, 4, 2, 1, 3, 2, 2, 1, 6, 6, 2, 1, 31, 1, 33, 1, 2, 1, 2, 1, 11, 4, 3, 1, 4, 3, 5398, 1, 2, 1, 33, 20, 2, 1, 3, 1, 195, 1, 3, 4, 1, 1, 2, 5, 9, 16, 1, 5, 1, 819, 2, 1, 2, 2, 7, 3, 2, 7, 2, 31, 3, 4, 2, 4, 3, 1, 3, 1, 8, 1, 1, 1, 1, 1, 1, 18, 15, 1, 3, 1, 1, 1, 2, 2, 7, 1, 1, 1, 7, 2, 7, 6, 1, 4, 1, 1, 5, 1, 1, 2, 7, 1, 8, 1, 2, 14, 1, 10, 3, 1, 1, 7, 1, 4, 1, 1, 3, 4, 1, 1, 1, 3, 5, 1, 2, 2, 3, 1, 1, 8, 5, 2, 6, 3, 24, 1, 7, 19, 2, 9, 1, 2, 1, 3, 5, 1, 7, 2, 1, 14, 1, 16, 2, 1, 3, 1, 3, 9, 1, 14, 3, 1, 1, 1, 1, 1, 7, 1, 15, 2, 2, 6, 30, 1, 1, 1, 69, 1, 1, 32, 2, 1, 2, 2, 1, 2, 1, 12, 1, 1, 7, 1, 3, 8, 2, 1, 62, 2, 1, 1, 2, 4, 1, 42, 3, 2, 8, 3, 1, 5, 1, 8, 2, 2, 1, 1, 1, 1, 5, 6, 11, 1, 7, 1, 1, 47, 1, 2, 1, 1, 3, 1, 13, 3, 1, 2, 1, 2, 1, 7, 1, 83, 1, 1, 3, 1, 1, 1, 5, 8, 1, 1, 4, 1, 1, 1, 1, 2, 209, 1, 6, 35, 3, 5, 1, 2, 2, 3, 1, 1, 2, 6, 8, 73, 1, 6, 1, 2, 2, 593, 33, 1, 3, 1, 4, 1, 1, 1, 1, 2, 1, 3, 1, 1, 6, 5, 1, 28, 1, 5, 1, 1, 8, 1, 2, 2, 1, 2, 7, 5, 1, 8, 5, 6, 1, 10, 1, 11, 1, 8, 7, 4, 3, 1, 1, 8, 1, 2, 14, 1, 3, 1, 4, 2, 7, 2, 3, 4, 1, 2, 4, 11, 1, 1, 1, 2, 3, 1, 4, 1, 45, 4, 1, 9, 5, 1, 6, 3, 4, 3, 2, 2, 1, 6, 1, 1, 4, 4, 1, 4, 2, 1, 4, 75, 1, 54, 1, 10, 1, 3, 2, 2, 3, 2, 3, 28, 1, 7, 1, 1, 31, 9, 2, 2, 3, 80, 1, 31, 1, 5, 2, 30, 3, 1, 38, 1, 13, 1, 1, 2, 1, 6, 2, 2, 3, 2, 28, 2, 6, 1, 1, 1, 2, 1, 2, 5, 9, 1, 2, 6, 1, 1, 5, 2, 1, 4, 1, 8, 2, 1, 6, 1, 5, 1, 1, 26, 1, 1, 6, 1, 2, 1, 5, 1, 2, 2, 2, 5, 1, 3, 1, 4, 2, 3, 13, 7, 12, 10, 1, 1, 13, 3, 4, 1, 3, 1, 15, 1, 1, 3, 6, 1, 1, 5, 35, 310, 1, 10, 1, 1, 2, 1, 1, 1, 5, 1, 3, 6, 14, 5, 2, 8, 3, 4, 1, 1, 13, 17, 1, 1, 1038, 1, 5, 1, 2, 1, 4, 4, 524, 1, 1, 5, 3, 1, 5, 2, 1, 1, 1, 3, 633, 1, 5, 8, 1, 3, 1, 2, 1, 10, 1, 1, 2, 2, 5, 5, 14, 19, 1, 1, 1, 1, 1, 1, 1, 3, 36, 1, 2, 1, 4, 1, 3, 3, 1, 18, 2, 2, 5, 22, 1, 2, 1, 1, 2, 4, 1, 2, 1, 39, 1, 2, 4, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 6, 4, 16, 8, 4, 2, 4, 45, 2, 2, 1, 4, 9, 3, 2, 1, 1, 5, 4, 1, 7, 70, 1, 2, 2, 4, 1, 4, 2, 1, 1, 10, 1, 6, 1, 8, 1, 1, 6, 2, 5, 1, 1, 4, 1, 1, 3, 1, 2, 3, 2, 7, 1, 1, 2, 1, 3, 3, 18, 1, 4, 3, 1, 4, 1, 3, 1, 1, 4, 2, 65, 1, 34, 7, 8, 5, 2, 1, 4, 1, 1, 1, 15, 1, 61, 1, 2, 1, 1, 18, 11, 2, 1, 1, 2, 1, 10, 1, 2, 2, 1, 1, 1, 2, 1, 1, 6, 1, 3, 1, 1, 2, 3, 2, 2, 6, 1, 1, 2, 5, 1, 2, 91, 1, 1, 3, 5, 2, 1, 1, 4, 2, 2, 2, 1, 5, 1, 53, 1, 25, 20, 1, 8, 1, 3, 1, 245, 36, 1, 10, 4, 1, 1, 1, 3, 1, 2, 2, 13, 1, 2, 10, 2, 1, 4, 1, 5, 1, 1, 1, 1, 1, 6, 1, 1, 9, 8, 2, 1, 2, 225, 1, 2, 8, 5, 1, 1, 5, 2, 80, 1, 1, 6, 2, 5, 2, 2, 9, 6, 16, 1, 1, 1, 15, 3, 3, 1, 5, 1, 1, 1, 2, 14, 399, 2, 1, 1, 2, 4, 1, 4, 2, 2, 3, 9, 7, 1, 2, 30, 1, 2, 1, 10, 1, 30, 17, 1, 1, 1, 45, 11, 1, 15, 1, 1, 1, 4, 6, 1, 2, 1, 3, 1, 2, 5, 1, 2, 1, 1, 1, 354, 2, 2, 4, 1, 2, 18, 8, 3, 71, 1, 4, 3, 1, 6, 2, 3, 10, 2, 1, 4, 1, 1, 1, 1, 2, 2, 1, 1, 2, 4, 1, 12, 4, 7, 1, 49, 2, 1, 1, 1, 1, 1, 2, 8, 1, 316, 1, 33, 2, 5, 11, 1, 1, 1, 3, 1, 1, 89, 1, 1, 7, 1, 3, 2, 3, 3, 1, 1, 3, 1, 1, 1, 1, 6, 1, 3, 2, 1, 2, 1, 1, 1, 2, 1, 8, 7, 2, 1, 2, 1, 2, 2, 3, 1, 1, 3, 2, 156, 1, 6, 3, 25, 31, 2, 2, 2, 4, 1, 2, 1, 2, 3, 1, 6, 6, 4, 1, 3, 2, 1, 8, 1, 5, 1, 2, 1, 4, 1, 2, 1, 3, 8, 1, 2, 51, 4, 3, 2, 1, 1, 25, 2, 1, 1, 1, 11, 4, 3, 21, 1, 1, 1, 1, 3, 1, 1, 11, 1, 1, 1, 59, 3, 3, 2, 3, 2, 1, 39, 1, 1, 3, 12, 64, 2, 1, 18, 1, 9, 2, 4, 1, 2, 11, 1, 3360, 1, 4, 1, 1, 1, 2, 4, 1, 12, 1, 6, 1, 3, 1, 8, 1, 1, 7, 2, 5, 3, 1, 5, 2, 1, 3, 1, 2, 1, 2, 1, 1, 394, 30, 1, 2, 8, 35, 4, 56, 1, 1, 13, 1, 5, 4, 2, 3, 1, 3, 4, 48, 1, 6, 1, 2, 1, 5, 5, 8, 2, 2, 3, 2, 1, 1, 3, 1, 11, 2, 8, 3, 5, 4, 2, 2, 4, 1, 2, 2, 1, 1, 30, 3, 1, 7, 1, 2, 130, 2, 3, 10, 28, 3, 1, 7, 3, 2, 5, 1, 1, 1, 4, 2, 5, 2, 4, 1, 3, 1, 1, 4, 1, 1, 1, 1, 3, 38, 2, 2, 1, 4, 6, 1, 2, 1, 2, 15, 1, 1, 20, 1, 5, 6, 1, 3, 2, 12, 2, 1, 14, 7, 1, 2, 7, 2, 2, 1, 3, 5, 7, 1, 3, 1, 40, 1, 3, 1, 2, 3, 7, 7, 3, 3, 1, 2, 1, 11, 1, 1, 1, 3, 2, 2, 9, 1, 2, 1, 4, 1, 1, 1, 3, 1, 3, 1, 1, 3, 1, 4, 1, 3, 8, 2, 2, 1, 1, 1, 2, 12, 4, 6, 5, 7, 7, 33, 18, 2, 1, 2, 2, 6, 3, 3, 17, 2, 2, 2, 11, 322, 1, 1, 2, 1, 2, 3, 2, 44, 1, 10, 18, 2, 47, 3, 3, 23, 3, 1, 2, 5, 1, 13, 4, 1, 1, 3, 11, 1, 3, 2, 1, 1, 7, 1, 15, 4, 52, 2, 11, 1, 13, 7, 4, 2, 1, 2, 1, 5, 14, 1, 2, 5, 6, 1, 2, 5, 3, 2, 1, 1, 1, 5, 3, 1, 1, 2, 1, 10, 3, 7, 30, 3, 1, 1, 1, 7, 8, 1, 2, 3, 1, 1, 21, 10, 1, 1, 1, 4, 5, 1, 2, 1, 1, 1, 2, 1, 2, 2, 2, 6, 1, 5, 1, 2, 5, 1, 3, 1, 205, 15, 4, 1, 2, 1, 2, 1, 1, 1, 2, 6, 7, 2, 3, 45, 9, 5, 6, 1, 1, 2, 1, 1, 4, 4, 4, 1, 1, 4, 2, 49, 3, 2, 3, 4, 1, 1, 2, 6, 1, 1, 2, 1, 1, 1, 6, 4, 2, 7, 1, 2, 1, 1, 4, 2, 16, 2, 1, 1, 3, 12, 2, 1, 1, 5, 2, 1, 2, 4, 6, 10, 1, 3, 1, 1, 1, 18, 1, 5, 3, 1, 1, 1, 220, 3, 1, 1, 1, 1, 5, 1, 2, 5, 1, 24, 7, 3, 5, 35, 1, 1, 38, 4, 30, 1, 1, 5, 1, 1, 3, 1, 5, 1, 2, 5, 1, 1, 5, 1, 54, 2, 2, 8, 1, 1, 1, 4, 4, 1, 12, 2, 2, 6, 1, 1, 13, 1, 2, 2}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 12343589 0427661894 2150569501 1063830550 6140378776 4681494642 1299405361 2353021788 1377880298 2299292189 2848537997 9652665360 1749216227 0711101794 2335388556 4134737917 5127385087 5220008258 4270586922 3480969648 9825517302 2630285624 8319594214 2178913472 6627048648 7673783218 3174322433 9573035949 5066034659 7236044064 5598262686 4330085375 2997974896 1037359399 3377391041 9089757559 1045069606 1574184551 7491402050 3718468616 0325405361 0295666286 7373524261 5876811260 0936474252 8010415343 1211657625 9556133221 8904653551 9414126621 5948852269 3291278433 1031722713 5482081180 5916162584 1172089073 6524031518 8072105291 2922265901 2560218924 0816163386 9827918700 3652684673 5372960843 3290228453 0086019828 5972767438 7822203920 2240871241 4054002832 5352098976 2254517747 3702172072 4767437059 7180239355 9955635985 0806048899 4746132739 2941367055 1121823928 8732079851 2210808111 6251044314 2318537641 3229978262 7526550176 6681480242 7095532088 9191753111 6118567910 6228986022 3824885923 7159105140 4140973784 4139070130 9777097470 2654592723 2175253212 3178874041 0287208961 4841597880 7831062613 5297259156 6806870474 4019907882 3918847103 4169902063 7672891107 1743606226 0305590817 9519589322 6637768434 8754325427 2704817635 5457319440 8380836565 9222720835 4829379959 6221440964 1134314335 6522228252 9309150214 0380595119 2487562586 8862500911 3970330359 2821980018 4330262050 0279226073 7747221295 5735286841 5071591701 2480334926 4165718861 3070989832 9608849542 7435458707 9846632865 5902316000 3639779041 7328775648 2986826849 5042170291 4539208596 8512874728 6459325987 8052524933 6176884665 2396359274 4312766286 3823376013 1314827967 9057101078 3217518079 5188164720

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.

The real roots of P are:

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