Primality Certificate for (4466^2887-1)/4465

Andy Steward10,534 digits06 October 2005
Originally by A.A.D.Steward 2005

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.434092% factorization of N-1:

From Factorisation
44662 · 7 · 11 · 29
Φ23 · 1489
Φ337 · 61 · 8839
Φ63 · 43 · 154579
Φ1398047 · 33212520894697 · 19336850067427414874789737
Φ2618416506187 · 674103275167297 · 5069833662058052369
Φ37149 · c130
Φ3979 · 13213564027 · 1647557070337621056180841 · p52
Φ743257 · c128
Φ78p88
Φ11137 · 30315433 · 25734334573 · c244
Φ222223 · 443810723467 · c249
Φ4811045905120521 · c1565
Φ9622887 · 236653 · 6198167 · 66648916087663729 · c1545
Φ1443c3154
Φ2886p3154

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 :

3385 9126789333 0121324138 1052505207 0867432253 3739204134 0611385075 6644288256 3329959246 5753001349 8736141811 8955338030 1951235891 2607709951 9108624730 1299645773 1103877002 8830985828 7860788787 4204222991 5203595006 8388213376 2433413278 6917349507 6815538720 7695909780 1075693047 2086403601 4132546826 5567830395 9060356544 6806106393 0662352831 7291084096 7608831340 5114624130 7982639048 8735477480 7544146332 7716403220 0708276169 1778414948 2014239425 1043794289 0112287569 1297684747 9059187747 7404786419 6909745609 8186243527 6944881822 0556220063 3713535290 8446264082 8849327044 5972237627 9569711540 2817041424 7206746472 4499909099 8596116644 7151952878 0359782213 6877212361 0083765416 5304082467 7169800995 1528316639 2523010157 1796741673 4552152103 3992130942 5564853511 7384061159 2757237894 5770178134 9549150189 2470618396 8387364174 3884975718 0557148820 2852937282 1622555995 9386404275 1210288191 0465925614 7081836070 2062269860 4439978982 7418403953 6046014692 6614998273 6561002207 9163842330 9115778359 6746529505 1395274144 9788859048 9717429752 7343416809 9928412846 6289308384 8183069354 4522233196 0043211442 5247563669 5012406549 2648526049 4589089986 0685154917 3723910268 5734297797 7009525698 0320457285 1151756761 7764684751 6002817394 7259806925 7906321226 0843134301 4814219688 2705571061 3407038348 0148677393 9510839069 4095720140 5136989902 6952630451 0555109537 1958567937 8467306229 2747862470 4635594005 1728044974 6562394794 5445980512 9774731374 5006521320 7169837089 7051203433 1898829323 0578709203 0302651648 9655617495 5353247239 0888096400 2668206905 3607284175 4243078346 5745116401 0578216728 4620201669 4338166385 5322329193 1478756774 3729470349 5806597661 6962305823 5450739162 2786268289 8301392224 1820440713 8133682873 0425782366 6593326314 3618532949 2921114969 9925849754 4015230716 9487191038 8550173383 4726080420 2470922567 9500899048 0597865911 4205226152 3491805520 8606929914 9207827820 4338891392 8524405828 2829183865 4008435130 1913554729 8652582252 2367982119 3193832768 3116237884 7867192747 8267425075 5541970050 1177186704 2314117305 6901097508 5726671288 6921210622 3038231321 0237819923 7738111072 2285223097 4856268978 7690354536 0374296046 4040745970 5220006694 9239048042 6972615198 9036554713 1737528323 7813252663 8607321201 6466777531 6496470277 5542637717 3897925310 9479586692 4153893782 3257489669 1120312525 9911295244 3543011231 5081273371 9285148390 1337989954 2305946294 0016002964 4891747458 2804317342 1250810337 5854319155 2453249782 9858153262 1583129131 2085076342 6965720902 3790350167 1592455716 4756820962 3910842526 7704133828 5215151296 4859355042 1634596131 9813390896 7543953130 9140265872 8473539830 1652609860 8225644525 8975042295 1636964467 4691099631 0816402666 0673376764 9183170400 7537774082 8271341733 3831012124 7583340570 8689668907 8339318502 3319265575 6398085649 4331938361 9306288646 6759683545 1751279967 4297151551 7203359652 7118637648 5264232602 5929767933 4519566752 8962925223 0168803210 9440568642 7316331956 9323272988 9832792604 4203211402 4487114199 5549996745 2078779377 5930369967 4180665705 8700488658 4067482194 0687093532 0875529112 5747364316 6947569956 1429551999 7126936829 4294657178 1300456227 2427105946 9540751201 2372917722 1176144085 0590802017 3610402578 9000808356 2413486838 0716511062 3577975624 3152550300 4346396498 0330364590 7201831949 3953570086 9100073076 3229005044 6606597375 8020145681 9614156510 4763938131
39641171 7596065816 1731692044 8413617521 6375673110 0313058067 1731098933 5860770160 3554944731
23 0390444257 0980466689 7800550213 7059262585 8008787027
193368 5006742741 4874789737
16475 5707033762 1056180841
506983366 2058052369
6664891 6087663729
67410 3275167297
3321 2520894697
104 5905120521
44 3810723467
2 5734334573
1 8416506187
1 3213564027
30315433
6198167
236653
154579
98047
8839
3257
2887
1489
223
149
79
61

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

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