Primality Certificate for (4411^2161-1)/4410

Andy Steward7,873 digits23 January 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.362390% factorization of N-1:

From Factorisation
441111 · 401
Φ22 · 2 · 1103
Φ33 · 151 · 42961
Φ42 · 1249 · 7789
Φ55 · 5381 · 14073875441
Φ613 · 1009 · 1483
Φ82 · 17 · 492257 · 22619209
Φ93 · 2455280371734634798831
Φ104051 · 93430256671
Φ1237 · 313 · 76081 · 429661
Φ1531 · 38821 · 158731 · 750080421110947561
Φ162 · 337 · 846161 · 1837393 · 136767057591041
Φ1819 · 2377 · 4987 · 32703948917551
Φ20603401 · 237514656604088574744041
Φ2480473 · 203449 · 8753680892343731233
Φ273 · 433 · 256771 · p58
Φ3061 · 181 · 12983342034135758874171601
Φ36369362629 · 305634412814977 · 480606332402925582517
Φ4041 · 9027525392801 · 4768733051544594330361 · 11636900936963005757281
Φ45758032565491 · p76
Φ48193 · 36721 · 2014752038298673 · 14728999605106561 · 97662145647019147009
Φ549446577790852621 · 85983669395248861 · p33
Φ60p59
Φ7273 · 570529 · 169373120402933857566671482329878521 · p45
Φ80241 · 343802561 · 1779062149963631201 · 289032496799527608267841 · p64
Φ908796421 · c81
Φ108109 · 2441449 · p123
Φ12029881 · 67801 · 44803982086199930783881 · c85
Φ135811 · 2161 · 46107333271 · 5160326699701 · 391120253426546101 · c216
Φ144577 · 1297 · 40609 · 1258011649 · 276883390657 · c144
Φ1801341166112091361 · 36141551010833279941 · 212802179243374854497708378641 · p111
Φ21638836153 · 13108300915707018801001 · p233
Φ2401201 · 7531518506401 · c218
Φ270271 · 2971 · 525099241 · 358239318554229398109779311 · p222
Φ360199081 · 77793532561 · 117398878858817880481 · c314
Φ432p525
Φ540541 · 30781 · 101026441 · 9626360221 · 349773803101 · 2637240057953461 · 3833961224754061 · c458
Φ72032401 · 3813545966569921 · 101682383559365777761 · 31455396336138941988001 · c638
Φ1080538921 · 159964776721 · 2182642674640314005492641 · c1009
Φ2160150682055761 · c2089

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 :

65063 3349426761 2254311886 1758705409 3453264036 7447284296 3078123589 4894749244 6763380950 8214593018 7085127003 9755247717 5708946013 5991798206 7711023490 2418993629 3724553069 6610554388 2562385897 3194006501 6397435059 9742914647 2120682539 7418628575 4794803295 1127926892 5665993163 9556004213 8041957237 0581455094 6976323525 1082227377 7336117690 6310359412 0936965675 8895000051 4428977329 0263555132 8025324824 1494520322 7239195242 2567557009 1743142495 0168614788 4036518301 3086446686 9684540372 3040550219 1626419615 2800171565 7260532945 8415712203 0097655521
501 0551774863 3700874251 8216574646 7600157400 7590486018 8218200520 9662316679 2201896935 2550421379 9439292810 1534489133 9943656575 4276312540 8752075700 5238281957 8021404985 8020628933 5560811459 3337036748 6464419569 0290850252 5753210858 7922321937
16 8415312978 0973234628 0518010796 3906226163 4822122901 6839277768 4583535624 1376713583 4431961428 5096516046 1702998166 1289630638 0597064661 7085848648 2227009798 1118884358 1764996907 1708106804 5058547755 5297793804 6362072359 1112017771
600 1501589341 7039470294 3332050222 9853212047 3811237948 5696221578 7967145033 5971847397 1692976137 8133038665 4057919044 6274931541
8 4006587753 9355286940 8673419727 5483000120 2783340555 9322476465 0207921202 5529506448 9074248577 9918748629 4697589501
388330 5189189668 2671005089 6447141268 3636750734 7788807246 7727462228 1097407531
9901 9767964523 8941786146 4110634862 3185174387 6282838105 3504452721
205396456 4520227022 5580164567 3245067638 6546619636 8324147681
11981519 9845490961 8185200782 4348102783 1651453150 9261692637
41729 6087957519 2169386661 2767176462 5364959153
169373 1204029338 5756667148 2329878521
492 0129795746 2907635851 5507105911
2128021792 4337485449 7708378641
3582393 1855422939 8109779311
129833 4203413575 8874171601
21826 4267464031 4005492641
2890 3249679952 7608267841
2375 1465660408 8574744041
448 0398208619 9930783881
314 5539633613 8941988001
131 0830091570 7018801001
116 3690093696 3005757281
47 6873305154 4594330361
24 5528037173 4634798831
4 8060633240 2925582517
1 1739887885 8817880481
1 0168238355 9365777761
9766214564 7019147009
3614155101 0833279941
875368089 2343731233
177906214 9963631201
75008042 1110947561
39112025 3426546101
8598366 9395248861
1472899 9605106561
944657 7790852621
383396 1224754061
381354 5966569921
263724 0057953461
201475 2038298673
134116 6112091361
30563 4412814977
13676 7057591041
3270 3948917551
902 7525392801
753 1518506401
516 0326699701
75 8032565491
34 9773803101
27 6883390657
15 9964776721
15 0682055761
9 3430256671
7 7793532561
4 6107333271
1 4073875441
9626360221
1258011649
525099241
369362629
343802561
101026441
38836153
22619209
8796421
2441449
1837393
846161
603401
570529
538921
492257
429661
256771
203449
199081
158731
80473
76081
67801
42961
40609
38821
36721
32401
30781
29881
7789
5381
4987
4051
2971
2377
2161
1483
1297
1249
1201
1103
1009
811
577
541
433
401
337
313
271
241
193
181
151
109
73
61
41
37
31
19
17
13
11
5
33

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.343270%

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 = 53 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 ≡ 57 (mod 63) and therefore cannot be a square and N is prime.