Primality Certificate for (4469^3257-1)/4468

Andy Steward11,886 digits18 September 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 35.004666% factorization of N-1:

From Factorisation
446941 · 109
Φ22 · 3 · 5 · 149
Φ42 · 2161 · 4621
Φ82 · 17 · 17 · 4657 · 148186057
Φ118950129 · p30
Φ2223 · p36
Φ37186481 · c127
Φ44618400421 · 22404642674041 · 994330488693600301 · p33
Φ7413913 · c128
Φ8889 · 617 · c142
Φ148c263
Φ29621313 · 76961 · c517
Φ4079769 · c1311
Φ8145124131 · 9328441 · p1301
Φ16281245421 · 6710617 · 16564901 · p2609
Φ32563257 · 146521 · 7238089 · c5241

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 :

102373581 0816788086 1869627830 3421525884 3823837598 7210386699 8461308345 7085660042 1711898016 5733815559 3759556784 6417980247 5007894466 8917108423 7510553853 3404410223 3185695993 9679551011 9076122486 1536884032 2681479553 6832252570 6376307059 1728265097 5100270295 5563289127 1383567392 7893622469 7483378264 8033586831 9927975239 0134602042 5912799419 2590316858 9653624719 6582449935 5308535627 7617727774 0549248943 2842649801 6724843698 6927600783 0123440936 6577489609 4428926789 1759997287 4412899554 6267049585 4255377682 6763808455 9355162158 0278956945 1347572575 7878327171 9943846732 2116531476 2284284706 0283722543 1546340990 5015965017 7995286742 9446559276 8139269111 2403108769 9660139347 0590692629 7441746666 0462894349 6212984103 5629686286 3409580998 9053947412 5273953731 7723793209 8058620773 9709850298 2966386395 3703277350 9818162912 8418272474 0049272145 8782971790 1539286959 4491644635 7352363753 8398476428 8701359660 2903621638 7712662298 0491214966 9067958934 7984239960 4283251941 8988984215 5157597263 0567812121 5501503724 9971266479 5992242299 6017745990 2224410752 9650188453 3009671552 8363382506 8464168814 0892441811 0330696877 6313373660 9508281365 7841563490 1315114679 3311290883 6893501072 7358510243 3775234016 0477141450 2102411453 9856764859 5039785024 9079506799 4486022422 2243985719 0395293988 5602915585 3157556175 9016781160 3214376390 1873769959 2257901041 4204911533 0360414060 9901646580 8359102078 6756396628 0682920936 7070301101 7797111410 3579761232 9758163442 7784471030 6942927285 8945225526 4181358883 8137252834 6450857015 6667936985 6937455928 1237472278 1261653499 1579726050 6740123977 7610068947 8327852321 9341591152 4915764500 8266984887 3440333695 3651710487 2787234905 1387523215 6887898405 6678834145 8997982526 3712950056 0145162115 7778988661 8476401898 3942579217 0552348033 7164480026 4186140332 6774754908 7018962035 3847321632 0079948788 5431158408 3072283455 0815731104 7957037435 1031842904 1280963953 0253335754 4478708267 1714814823 0518489667 9877424040 1885022535 9436304417 4800333763 4362623039 1808710814 0629108852 6596074647 5899330321 8818486202 4125269994 5035199950 7254761772 6172401807 1703098797 0957444419 7687327318 0289652220 9550256934 3302261314 1917399966 2916024056 9838210442 6270319594 1881058110 0254700164 1266962352 1949222750 3483003634 4181767858 2769646133 3182214061 4468188053 1427112714 5950315723 2984518451 5831345152 9826371759 6406517360 0474837026 5568466170 9223578691 0922192800 8474810363 8937927146 1348634493 5734287456 1561973619 5834633435 0539038571 7303104725 2370350294 1954350770 4261220963 2895340784 6536000532 8455653263 1280681702 2075546595 2039267390 6744885497 7572583510 2797010684 3333640969 5739119713 7431123914 4070619527 4767263146 7016919838 4763789080 9795413941 8878604197 2520761123 3477130441 8850900681 1996195593
2 4911251885 2725967424 8459331641 1548470775 4474502467 2418924143 7268709205 2745716466 1186532780 6394662884 2865464533 3742559870 8992932143 8146427267 6367657733 5063307703 3492268582 3529853366 5135126835 6413395380 5841675839 6819738434 0827365522 3693661933 1566378043 3367007405 3887888595 5878039308 4401588175 8114200723 6673327496 8533483957 6520382142 1524107517 4818636222 8514043107 8888384236 0916556948 0767186108 0825912712 1114321097 6223927381 7470648880 1336267658 2934688681 4105489734 9998534981 9908698059 6759711343 6424136768 2371588101 9957186224 6770909528 7009416127 4343543833 5465864773 7938413248 9932076728 2750538938 4239645193 1975628658 4431523882 6592814670 2376983293 0091186673 3624675890 2068683246 4966882021 7623854925 9880170166 0890181216 6403829232 2761030342 4129579568 6627455111 3442891729 9290949108 2829263378 6554942425 2737250997 7989296090 1597530661 8725807445 1976031849 9852945655 1942426819 5250467602 9810791908 8011029158 4129372912 4138022962 5528526032 6659045600 3799803848 1233791662 1182828705 5467001882 4529245861 6470420066 5524693596 0498923785 0966473238 5125327972 1870565836 9893208706 6366176272 3596030205 5262508144 5094062491 6428087225 3956021667 4796692700 4186196669 6833082295 7724333162 3114434932 6787306323 2719075853 8857706148 9971194142 9364286850 0666178517 9626042340 8033662604 6136408399 4099526023 2381942870 0261640832 7304667422 6085412531
138126 9881379987 0784292939 5800427507
732 9401569619 7566360499 8612599841

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

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