Primality Certificate for (8854^2521-1)/8853

Andy Steward9,947 digits02 February 2001
Originally by David Broadhurst 2001

This certificate uses a theorem of Brillhart, Lehmer and Selfridge 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 33.352599% factorization of N-1:

From Factorisation
88542 · 19 · 233
Φ25 · 7 · 11 · 23
Φ33 · 26134057
Φ478393317
Φ561 · 53881 · 1870000151
Φ637 · 139 · 15241
Φ76553 · 918947 · 2760101 · 28988821
Φ86145511993475857
Φ93 · 134853931 · 1190836781977177
Φ105 · 2868011 · 428507281
Φ12415213 · 14800865857
Φ147 · 29 · 43 · 3697 · 14927054343244079
Φ1531 · 541 · 2251687561649759519352309361
Φ18481767063685648623358633
Φ2041 · 421 · 33845561 · 52724381 · 1226131292381
Φ21967 · 126211 · 295967911 · p31
Φ2473 · 71257 · 19468921 · 46977169 · 7938473929
Φ281709 · 1843256521 · 4497124881136049 · 16383678485047783561
Φ30331 · 2011 · 56744676526059918076163191
Φ3571 · 211 · 69162661 · 1139988294061 · 6362825972107921 · p55
Φ363287881 · 310659476797 · 4032651348541 · 56348562194726833
Φ40401 · 918495035561 · p49
Φ421597 · p45
Φ45205651 · 398837338124041 · 1165246949232271 · p60
Φ56617 · 12769330243093902692090272369 · p64
Φ60183067084801 · 2263899782281 · p40
Φ63379 · 4890614023 · c130
Φ70439868537921 · 241671302391168495818135521 · p57
Φ723801241 · p89
Φ8444764378014133 · 1002763675400087963393970184177 · p52
Φ90181 · 543061 · 671341002571 · 5967074360670950701 · 2429454286918527018391 · p35
Φ1052265360594791430557393571481 · p163
Φ12082235662969338481 · 449325025236672121 · p92
Φ126127 · 1009 · 30519721 · 95459869 · 782824046499566462444401917114967 · c89
Φ140281 · 837396550061 · p176
Φ168c190
Φ18036721 · 8380844057162341 · 150643812475098001 · 1417802894280598732830901 · 147461356689543777420522188446321 · c96
Φ2102377201 · 34993771 · 11844175351 · c166
Φ25212601 · 673093 · 2908333 · 150775129 · c260
Φ280106121 · 845639813201 · 820847964190623601 · c345
Φ315631 · 36541 · 119701 · 112416786005763961 · p539
Φ36077041 · 268363965241 · 976945153292894881 · c345
Φ4202521 · 2334781 · 2674561 · 1127007896042792497792084321 · c336
Φ5042017 · 244948096119889 · 366163289282037241 · 635231866838909779449662635369 · 98767045911566020041985979002537 · p472
Φ630258091301161459235611 · c548
Φ840c758
Φ1260541035628909826793121 · c1117
Φ25209791582415120192961 · c2255

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

787683611 7083581737 6178598070 2500145359 9762500909 5504321180 6032870606 2305568900 9505410934 8675229641 0759080107 4935395061 1007155230 2480233946 6603529167 3722094040 2853277370 0034491293 4575140006 1786695274 4408327794 4931908038 5404726917 2762835767 3241181105 5553607736 8419327010 2867114514 7701062340 2161536218 5089052555 6516326394 8568432517 7801282017 9852626731 2282741735 2247363548 1718437274 0528937563 0625243781 9167956753 9925491161 1765066859 8455178205 1927319887 5376000117 4582173591 3043100150 4212855451 6062352615 6705293442 7793834626 9933414928 2606401511
21 5323500939 4792284476 6141665034 3467338149 7510664677 7533558421 3882733052 7959932387 6167340694 6887976966 5024166892 4458667268 1756955860 5681734543 9560270736 0480157738 5755362596 2709675721 6379793407 1803394807 3452685303 2376691216 4048244037 9901000569 7575818953 5263472852 8865111168 4380065526 4934031543 8546552365 9969959165 2698379094 8485983923 2344801215 4913555842 1016129329 4126418560 6727140188 3634137925 0861305991 0055734955 4650651334 8251519604 2501651880 3138193588 5532171225 6683755209
123327 3401944841 1554615956 3836873955 2016286916 6270544435 9138811761 3457942830 8385593718 6740595130 5201000686 4712556409 8580399777 9949203856 6385086509 6745924012 1525891680 3312759041
128 1175303583 7467821574 1116644445 2506897438 7098239639 2932553905 4278768350 8864564678 9695640471 0696468153 4042890653 4337171341 8455622665 5256632167 9736261549 9824642451
55 0609574738 9229757321 3582542207 1927947048 6954722170 0068607736 2434687208 3190757077 9689167721
141717348 6131257423 2497767953 5754410740 8049454330 6827899904 3314266450 9328939247 5550041001
6837 4653811823 9775722331 6915688107 8035315035 2643075909 3310672217
5636426059 2025965990 6250436537 0357209360 3028655515 2061152981
5068150 3387631637 5498663204 3858460800 0738197914 1593337751
71669 8399867053 8788825905 0467792172 2604673436 3498386271
12 0009957237 1342217256 5464253708 5264298731 3912238401
387267539 6074406806 0448480671 5458884761 1632804601
14535 1106911126 2148091025 4082789194 6640314843
3441635116 1771190809 4629908047 6402846781
56312 8981681942 7230944751 1768588041
782 8240464995 6646244440 1917114967
147 4613566895 4377742052 2188446321
98 7670459115 6602004198 5979002537
6 4247591159 6055183982 2714668353
1 0027636754 0008796339 3970184177
6352318668 3890977944 9662635369
127693302 4309390269 2090272369
22653605 9479143055 7393571481
22516875 6164975951 9352309361
11270078 9604279249 7792084321
2416713 0239116849 5818135521
567446 7652605991 8076163191
14178 0289428059 8732830901
4817 6706368564 8623358633
24 2945428691 8527018391
5 4103562890 9826793121
2 5809130116 1459235611
1638367848 5047783561
979158241 5120192961
596707436 0670950701
97694515 3292894881
82084796 4190623601
44932502 5236672121
36616328 9282037241
15064381 2475098001
11241678 6005763961
8223566 2969338481
5634856 2194726833
1492705 4343244079
838084 4057162341
636282 5972107921
614551 1993475857
449712 4881136049
119083 6781977177
116524 6949232271
39883 7338124041
24494 8096119889
4476 4378014133
403 2651348541
226 3899782281
122 6131292381
113 9988294061
91 8495035561
84 5639813201
83 7396550061
67 1341002571
43 9868537921
31 0659476797
26 8363965241
18 3067084801
1 4800865857
1 1844175351
7938473929
4890614023
1870000151
1843256521
428507281
295967911
150775129
134853931
95459869
78393317
69162661
52724381
46977169
34993771
33845561
30519721
28988821
26134057
19468921
3801241
3287881
2908333
2868011
2760101
2674561
2377201
2334781
918947
673093
543061
415213
205651
126211
119701
106121
77041
71257
53881
36721
36541
15241
12601
6553
3697
2521
2017
2011
1709
1597
1009
967
631
617
541
421
401
379
331
281
233
211
181
139
127
73
71
61
43
41
37
31
29
23
19
11
72
52

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

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 F2 < N < F3 and N ≡ 1 (mod F), we can let N = c2·F2 + c1·F + 1.

Brillhart, Lehmer and Selfridge's Theorem shows that N is prime if and only if c12-4·c2 is not a square.

Here, c12-4·c2 is ≡ 29 (mod 64) and therefore cannot be a square and N is prime.