Primality Certificate for (13735^3433-1)/13734

Andy Steward14,202 digits09 August 2008
Originally by A.A.D.Steward 2008

This certificate uses a theorem of Coppersmith and Howgrave-Graham 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 29.303625% factorization of N-1:

From Factorisation
137355 · 41 · 67
Φ22 · 2 · 2 · 17 · 101
Φ33 · 127 · 495181
Φ42 · 94325113
Φ65479 · 34429
Φ82 · 233 · 7817 · 9769866433
Φ1123 · 11287 · 57795364309 · 15926446665640190523746629
Φ1213 · 37 · 6841 · 10815584281
Φ131219349 · 1934688596600378581 · 19108935190356825890062169
Φ22331 · 30713 · 1289971 · 1689698407 · 10782405089599933901
Φ24241 · p31
Φ26230907977 · 20687488873636034581 · 9435515507179638234593
Φ33137780985912445631041249188939827450866981 · p42
Φ391249 · 50758423 · 165317553103 · 266439382811036917812721 · p54
Φ442113 · 70313 · 77130461 · 92106169 · p59
Φ5253 · 459313817 · p89
Φ66765431353259112079 · p65
Φ7879 · 157 · 2887 · 79873201010227 · 2310556673778846714044593486471 · p48
Φ8889 · 353 · 68993 · 28456913 · 55928489 · 10917744902417 · 4544440229544073 · 114508835668425089139309601520590721 · p78
Φ10413313 · 127876394722649 · 274591034634460040907187441 · c154
Φ13237357 · 450696973 · p153
Φ1432861 · 63493 · 7263829 · 624181273 · p473
Φ15613 · 7177 · 12367837 · 123753241 · c179
Φ2645281 · 12409 · 73231707538783051056242961793 · p295
Φ286739622027 · 56880193543976207098116281 · c462
Φ312313 · 8737 · 2254590557724689483187409 · c367
Φ429c994
Φ57225275740050725089394041 · c971
Φ858859 · 125269 · 57790967521717 · c972
Φ11441124193929 · 14880490769 · 130317172767497 · p1953
Φ17163433 · 13729 · c1979
Φ34329825817 · c3966

We need the product F of all the prime factors from this partial factorization:

659 6270598207 8098643847 3970067034 0456460022 1669940634 2009761426 6532214468 8279660035 9588699341 8663941651 5554509903 1221033132 6396732435 7376160609 5946584612 3884500987 7684186079 3368755100 2116141935 0784053020 4445156761 9926494827 2599329376 2336700837 0090177205 8569168389 9307489498 8928946075 8130316690 0851291612 1894888911 9498424060 4632349899 1576449230 2327191170 0111526664 6487758419 3624023823 8219283415 3118215964 8687884178 6925041477 5528487977 7278318247 7844973856 6123505744 9446998634 8697238793 1172894629 4113635220 9336640524 8908509362 2742076918 9401003989 8519155119 0835139493 7938539964 1994863555 0062719255 0971272493 5376041531 1550446473 2214994306 6010326009 5011346972 3134980496 0291662079 3981603671 4213925228 9540944136 6364241860 4846883547 3538606893 9079760158 9869022105 9122546289 1074437152 9840349581 0030343727 1072593795 9891212525 4368976057 8350086242 0346866556 8687949221 3402653044 2018728731 6339851161 4393805394 6610760763 9464677106 3008364306 6847225756 9625165134 5732625487 6341077504 6747164709 8499535710 6840687656 1702958932 7647493223 2281017174 7391314997 0257770504 9264560880 9037025775 1464602010 3282087901 4812914255 9850312225 2393445311 0610034151 2366648803 3103060646 6701175038 5459483825 7322088749 4009028274 3914677024 2240721215 0755024654 9516377863 6171842651 5090616531 3798222357 2419453474 5599292160 8793299878 7851882920 9753801795 1293536762 6965048557 5409714505 9389048689 3774681798 2189904931 7651756864 6253823304 1535595468 5520187768 2596846191 7836462901 4703371700 3563960427 4728770376 1024851253 5323105845 3189653500 3858187608 2114936098 7574975776 3212986981 5935232152 7614420065 0643540848 3237452594 8844402089 0023022740 2812830735 4674519486 7843909885 9706577437 2639857251 3678153341 1710674090 6749243668 2405969097 1415821311 8101194820 2197542203 2471166272 7137817442 0124497807 5148553941 3044538121 5102075444 9015515030 0830836279 4426264628 5889037819 6317212463 9823412349 5906826120 0041477414 2545468361 7738740663 1823582097 2220502804 4429435370 7383696258 9238715485 1451316072 6230418745 5459893278 4463849633
420 4244630948 2385757850 4412434770 9915077064 9071888190 1312135175 3640272891 1834111710 5324419572 8900259329 6191641979 8458239555 8355931383 2204312647 3968720699 5157694885 5307313251 5367355396 1153982233 5835249077 1502368733 5851366476 6577819459 5947575975 0830263509 7241264682 8153311270 9863907759 6741905168 7793341982 1253877987 0881424975 9645945358 5799898451 5330853297 8190491655 2672876721 9002495868 6727116328 4962920265 3926908277 0430091113 4028712870 5004113811 1072917741 0120665687 0540949701
22138 1327737835 4078348222 1533705983 8996571405 3190358829 4344942407 0002827960 1577878336 2203688591 6943444225 4680289077 8782475769 5737835634 4940809174 3898751114 2553333191 4646268519 2760583433 1012573483 2731815168 0640444355 9298064564 1148494152 8285922666 6056775356 9374316280 6012244191 5495941601 0620189033
193 5931199218 2314887773 3808518705 5544478941 8151756043 9161230514 4776468098 1486222936 0670414567 0592437736 6801894900 6603075770 7143894239 4716311118 3806226841
834646445 6376704388 1969190360 9027020953 8683916752 6198426397 6455365906 2076435387 4897388101
16630434 6114631921 8978803200 9372671719 5376128054 5122134827 8540480434 4196149273
74593 0982156384 9819073009 8511169607 3175991808 2594643161 0680263359
540908120 7665024863 2215387315 0157925863 7160848928 0316709381
7275 5924378100 3609683248 3672595133 2926407277 4098619641
30748809 5012438511 9210742619 4709497176 1959282953
41 4335711207 7674956495 5886401600 9519016861
13 7780985912 4456310412 4918893982 7450866981
114508 8356684250 8913930960 1520590721
5 2554785452 0972817568 8115333361
2 3105566737 7884671404 4593486471
732317075 3878305105 6242961793
2745910 3463446004 0907187441
568801 9354397620 7098116281
191089 3519035682 5890062169
159264 4666564019 0523746629
22545 9055772468 9483187409
2664 3938281103 6917812721
252 7574005072 5089394041
94 3551550717 9638234593
2068748887 3636034581
1078240508 9599933901
193468859 6600378581
76543135 3259112079
454444 0229544073
13031 7172767497
12787 6394722649
7987 3201010227
5779 0967521717
1091 7744902417
16 5317553103
5 7795364309
1 4880490769
1 0815584281
9769866433
1689698407
1124193929
739622027
624181273
459313817
450696973
230907977
123753241
94325113
92106169
77130461
55928489
50758423
28456913
12367837
9825817
7263829
1289971
1219349
495181
125269
70313
68993
63493
37357
34429
30713
13729
13313
12409
11287
8737
7817
7177
6841
5479
5281
3433
2887
2861
2113
1249
859
353
331
313
241
233
157
127
101
89
79
67
53
41
37
23
17
132
5
3
25

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

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 = 77 suffices.

Given such a witness, Pocklington's Theorem shows that every prime factor of N ≡ 1 (mod F). As F4>N, N can have no more than three prime factors.

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

Brillhart, Lehmer and Selfridge's Theorem shows that N has exactly two prime factors if and only if c12-4·c2 is a perfect square.

Here, c12-4·c2 is ≡ 37 (mod 64) and therefore cannot be a square and this stage of the proof is passed.

Coppersmith and Howgrave-Graham

We are left with two possibilities for N: either it has exactly three prime factors or it is prime. The non-existence of exactly three factors is demonstrated by the Theorem of Coppersmith and Howgrave-Graham, here performed by a Pari/GP script written by John Renze and David Broadhurst. Here is the stdout:

realprecision = 5009 significant digits (5000 digits displayed)

Welcome to the CHG primality prover!
------------------------------------

Input file is:  IO\35A70D69.cin
Certificate file is:  IO\35A70D69.chg
Found values of n, F and G.
    Number to be tested has 14202 digits.
    Modulus has 4162 digits.
Modulus is 29.30362529% of n.

NOTICE: This program assumes that n has passed
    a BLS PRP-test with n, F, and G as given.  If
    not, then any results will be invalid!

Square test passed for F >> G.  Using modified right endpoint.

Search for factors congruent to 1.
    Running CHG with h = 6, u = 2. Right endpoint has 1718 digits.
        Done!  Time elapsed:  258718ms.
    Running CHG with h = 6, u = 2. Right endpoint has 1500 digits.
        Done!  Time elapsed:  263188ms.
    Running CHG with h = 5, u = 1. Right endpoint has 1182 digits.
        Done!  Time elapsed:  14484ms.
    Running CHG with h = 5, u = 1. Right endpoint has 977 digits.
        Done!  Time elapsed:  15157ms.
    Running CHG with h = 5, u = 1. Right endpoint has 566 digits.
        Done!  Time elapsed:  118671ms.
A certificate has been saved to the file:  IO\35A70D69.chg

Running David Broadhurst's verifier on the saved certificate...

Testing a PRP called "IO\35A70D69.cin".

Pol[1, 1] with [h, u]=[4, 1] has ratio=4.413418078 E-250 at X, ratio=1.237490834 E-815 at Y, witness=2.
Pol[2, 1] with [h, u]=[4, 1] has ratio=2.434797276 E-412 at X, ratio=1.265958758 E-411 at Y, witness=3.
Pol[3, 1] with [h, u]=[4, 1] has ratio=4.056777760 E-206 at X, ratio=6.38579335 E-206 at Y, witness=3.
Pol[4, 1] with [h, u]=[6, 2] has ratio=0.2998537915 at X, ratio=3.905988040 E-635 at Y, witness=7.
Pol[5, 1] with [h, u]=[6, 2] has ratio=7.77719552 E-219 at X, ratio=7.80611958 E-437 at Y, witness=2.

Validated in 3 sec.


Congratulations! n is prime!

The actual input file containing N and F and the output certificate are included in this file.