Primality Certificate for (5968^2341-1)/5967

Andy Steward8,836 digits06 February 2002
Originally by A.A.D.Steward 2002

This certificate uses a theorem of Konyagin and Pomerance 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 30.226524% factorization of N-1:

From Factorisation
59682 · 2 · 2 · 2 · 373
Φ247 · 127
Φ33 · 7 · 1696333
Φ45 · 5 · 1424681
Φ531 · 40928548278671
Φ637 · 962461
Φ93 · 73 · 517609 · 398590508578267
Φ1011 · 41 · 131 · 105251 · 203971
Φ121268572362999553
Φ1313 · 1409851 · 6686649919 · 16660705778209740375765900433
Φ1524181 · 66540105030674212537180621
Φ1819 · 109 · 268507 · 1753471 · 46338139
Φ205 · 439981 · 1448295481 · 505090733750441
Φ26521 · 13125061 · 9466739101 · 31530695249531336762077721
Φ30181 · 8892517023850222749894818221
Φ36268777 · 2061337357 · 2692958941 · 25596208801 · 53456340097
Φ392341 · 724777 · 893908003 · p73
Φ451990409851 · p82
Φ5253 · 157 · 7470592998025993681 · p68
Φ6061 · 541 · 2663711561941 · p44
Φ65285611 · 6384461201 · 1208308292939540122370326651 · p139
Φ7879 · 2670814361519038461822811 · 1327049706247824611722420520251 · p35
Φ909331502987881 · 72105083628931 · p64
Φ11720359 · 740177211899553712654267 · c244
Φ1301075452057551 · p170
Φ15617317 · 231447642661 · 4795666101898621 · 5441894660766380078361818161 · c123
Φ18093241 · 1011382959842713261 · c159
Φ1953121 · 55381 · 44133971406993781 · 4428261403282678651 · 46566062095787887065271 · 42681415320279838976920801 · c271
Φ234c272
Φ26077563721 · c355
Φ390131779117861 · 4154027531551 · c339
Φ46812637 · 47737 · 6095233 · 17901937 · 274264381 · 317232253735021 · c498
Φ5851171 · p1085
Φ780c725
Φ11703698371 · 3448769221 · 35434725824919331 · 4139983456795112249505781 · c1031
Φ23401495261 · 298743121 · 340273441 · c2152

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

23449 9364843328 0401326059 8563294959 1799791846 7828241287 7476698632 5162125788 7803781824 5649878318 3693945622 6331809378 8013941504 5365369717 7617158432 2786715449 4319123181 3099960793 3784099257 4726496604 4281403359 3836916248 0952404806 8266393037 8205818895 6833629917 5928096530 7165155817 2145313399 4922495152 8649381295 9398165106 9109821197 5414590249 6019077463 7181979592 0850272373 2846732097 8037519591 4200779372 1366291815 0012153043 9908958775 2689697250 6407368084 9201718668 7664017532 5648829220 1834073948 7424130813 9930871928 0023595108 9688029750 4974288477 4411761245 8058571956 5359315030 5502513028 7727595761 4717654904 8791173621 9671224260 2397463691 5073699786 0989296254 7022121378 6434597594 4665383584 8851426774 1600983330 3966778867 2509100283 7912378135 8416032504 1601556518 1360144384 5690809607 9019254342 5118812240 7121850190 5683131561 4108455564 1432334817 4236869805 0995799777 6894207082 0717274615 7185845630 0087963494 6061447316 1414764780 1958771296 2405428098 4881157234 3650395436 6043845161 1039135156 6844729884 0980723785 1914670858 0864152997 8879284212 0795950017 3918969887 3109662441 5914216935 4230570163 8783964287 1527352731
1615343185 6216889447 3172530299 3195222081 7253462909 2410917799 1720646481 8194239893 2101585307 8870146325 7663670146 6879758140 2514270015 4612031056 8556526834 2593173245 7466070751
788194462 9621792472 0716794797 3723834339 4912118387 1179161590 7201418364 6341592958 0329607484 1979157461 0641403366 4968506471 2293320387 4483453241
20 9386672215 1422208616 8360453130 1061349149 7772070435 2608411834 2015611207 7281001651
274 7388830857 4630075910 3037073711 2853786162 4605212438 1131662692 8338184231
67044157 3210338704 0033125071 6943516681 7939060305 6000896930 1482158513
6194 0408155542 0366874600 0860917358 3860053889 0936013778 2209535091
2946 0950047575 9817008173 9691946456 8289203061
14886 9704815295 6160888262 6063456279
1 3270497062 4782461172 2420520251
166607057 7820974037 5765900433
88925170 2385022274 9894818221
54418946 6076638007 8361818161
12083082 9293954012 2370326651
665401 0503067421 2537180621
426814 1532027983 8976920801
315306 9524953133 6762077721
41399 8345679511 2249505781
26708 1436151903 8461822811
7401 7721189955 3712654267
465 6606209578 7887065271
747059299 8025993681
442826140 3282678651
101138295 9842713261
4413397 1406993781
3543472 5824919331
479566 6101898621
126857 2362999553
50509 0733750441
39859 0508578267
31723 2253735021
7210 5083628931
4092 8548278671
933 1502987881
415 4027531551
266 3711561941
107 5452057551
23 1447642661
13 1779117861
5 3456340097
2 5596208801
9466739101
6686649919
6384461201
3448769221
2692958941
2061337357
1990409851
1448295481
893908003
340273441
298743121
274264381
77563721
46338139
17901937
13125061
6095233
3698371
1753471
1696333
1495261
1424681
1409851
962461
724777
517609
439981
285611
268777
268507
203971
105251
93241
55381
47737
24181
20359
17317
12637
3121
2341
1171
541
521
373
181
157
131
127
109
79
73
61
53
47
41
37
31
19
13
11
7
53
32
24

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

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

Express N in base F

As F3 < N < F4 and N ≡ 1 (mod F), we can let N = c3·F3 + c2·F2 + c1·F + 1. Let c4 = c3·F+c2.

Square Checks

For t = 0 to 5, we prove that Q(t) = (c1+t·F)2+4·t-4·c4 is not a perfect square. This is done by checking whether Q(t) is a quadratic residue modulo a variety of bases. If it happens to be a QR in all of the bases, we calculate s = floor(sqrt(Q(t))) and show that s2 < Q(t).

Continued Fraction

We approximate c1/F by a continued fraction u/v such that v is maximal while remaining less than F2 / N1/2 = 3739 3308421440 9630760241 2647548987 6539527231 7956172582 7782432766 2260499135 7526983995 9191146768 5159705455 0505089885 3568816853 2191715259 4436129195 1326998350 9474464143 5189514743 5630062972 8541518117 4990384503 5187242635 3726459573 2702908465 2589112481 8284474453 4273462504 6251678281 3695126973 1818369531 6263484485 8414033518 5379414069 5554553342 9025014844 3680712152 1415427240 6653780422 4714788443 0660349782 3897880448 3717354006 1103611899 7370398914 9400815859 2304229242 2882967406 6933790008 9706655160 9452526245 4907487289 1162590435 9125308064 1710177225 9173945622 2809837130 2157638654 4187067232 1630299384 9589358324 7992523808 2395225692 4373327153 8900818302 0870468067 8406621306 6700074381 0910224993 7607310178 4138876986 0087748821 0112128299 9322305388 5463860917 8188944328 6179686017 0900518149 3554040846 4687692579 1158046830 8304820773 7070930726 5237683839 7674136949 2869493310 8039246867 8751334724 7428874955 6915192037 1507566335 9694958873 1959766202 8943749737.

With those constraints, the unique continued fraction is: {0, 6, 2, 1, 1, 2, 5, 1, 1, 1, 28, 1, 2, 1, 1, 1, 1, 99, 15, 2, 10, 3, 1, 39, 8, 3, 1, 8, 1, 2, 8, 4, 1, 111, 7, 1, 27, 1, 1, 2, 3, 398, 1, 17, 1, 6, 1, 2, 1, 1, 1, 1, 1, 1, 3, 8, 419, 8, 1, 1, 1, 6, 3, 69, 3, 1, 2, 1, 4, 1, 3, 1, 11, 2, 10, 3, 2, 3, 2, 3, 11, 5, 8, 12, 2, 1, 4, 2, 5, 1, 2, 3, 5, 7, 3, 5, 1, 2, 2, 1, 1, 5, 3, 1, 1, 397, 1, 1, 28, 12, 9, 78, 1, 1, 14, 1, 1, 2, 6, 2, 2, 6, 1, 10, 1, 1, 1, 7, 1, 7, 15, 1, 6, 2, 1, 1, 1, 6, 2, 1, 1, 1, 4, 4, 1, 1, 16, 1, 1, 1, 1, 3, 2, 38, 5, 2, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 5, 2, 3, 3, 2, 15, 1, 95, 4, 3, 6, 4, 2, 2, 3, 1, 1, 1, 10, 2, 4, 1, 1, 3, 1, 77, 1, 5, 1, 5, 15, 1, 5, 2, 1, 1, 1, 7, 27, 3, 2, 6, 1, 1, 1, 2, 1, 3, 2, 1, 3, 1, 1, 2, 4, 11, 1, 7, 2, 12, 1, 1, 5, 2, 1, 3, 2, 4, 1, 5, 2, 1, 1, 1, 14, 3, 1, 14, 1, 6, 8, 1, 1, 4, 2, 4, 2, 26, 2, 1, 1, 1, 16, 1, 1, 1, 4, 8, 1, 1, 8, 2, 3, 2, 1, 1, 2, 8, 1, 1, 1, 2, 1, 1139, 2, 9, 9, 1, 1, 6, 1, 5, 1, 23, 1, 8, 1, 1, 2, 1, 6, 2, 62, 1, 10, 3, 1, 1, 10, 2, 1, 1, 174, 11, 3, 2, 1, 16, 4, 2, 1, 7, 1, 3, 2, 6, 3, 1, 1, 6, 3, 1, 14, 1, 1, 1, 1, 90, 4, 1, 2, 1, 1, 1, 5, 7, 2, 7, 3, 1, 1, 4, 1, 1, 4, 1, 16, 2, 9, 2, 1, 2, 2, 1, 2, 1, 1, 28, 10, 8, 1, 1, 2, 26, 2, 9, 3, 1, 1, 1, 1, 1, 2, 10, 20, 3, 1, 1, 2, 6, 7, 1, 8, 1, 1, 3, 1, 1, 7, 4, 1, 1, 9, 9, 2, 2, 1, 3, 100, 11, 23, 1, 2, 2, 1, 4, 25, 1, 3, 1, 2, 6, 1, 1, 19, 1, 5, 3, 2, 2, 5, 2, 4, 2, 1, 2, 798, 1, 349, 12, 2, 2, 3, 2, 1, 1, 1, 2, 2, 11, 1, 11, 3, 3, 1, 2, 1, 1, 1, 1, 2, 1, 3, 4, 3, 2, 1, 4, 2, 1, 1, 3, 1, 1, 3, 3, 1, 1, 2, 1, 9, 1, 32, 1, 1, 5, 3, 2, 1, 1, 1, 2, 8, 2, 6, 2, 11, 2, 1, 25, 3, 1, 1, 4, 1, 1, 6, 2, 1, 4, 5, 2, 3, 2, 6, 3, 1, 15, 1, 27, 56, 1, 1, 2, 2, 4, 3, 1, 1, 2, 1, 2, 6, 1, 5, 1, 2, 2, 4, 1, 1, 2, 1, 22, 2, 2, 6, 1, 1, 3, 2, 27, 4, 1, 3, 1, 2, 8, 2, 9, 68, 7, 1, 11, 2, 1, 8, 8, 2, 1, 99, 1, 4, 1, 1, 3, 1, 3, 4, 1, 1, 2, 1, 4, 1, 14, 1, 4, 19, 1, 1, 1, 1, 2, 5, 4, 6, 2, 1, 1, 1, 20, 1, 1, 1, 1, 5, 1, 9, 2, 2, 1, 3, 1, 1, 7, 1, 1, 1, 2, 1, 3, 1, 23, 1, 2, 2, 3, 2, 3, 1, 1, 1, 3, 1, 16, 1, 1, 134, 4, 12, 1, 4, 1, 4, 2, 1, 1, 9, 1, 1, 57, 5, 1, 1, 2, 3, 2, 1, 8, 10, 7, 2, 1, 3, 2, 1, 1, 3, 1, 4, 1, 1, 1, 10, 1, 4, 5, 1, 6, 5, 1, 1, 1, 1, 1, 1, 1, 22, 2, 1, 1, 12, 1, 3, 2, 165, 2, 3, 1, 1, 1, 7, 1, 2, 1, 4, 1, 1, 3, 1, 2, 27, 2, 1, 6, 2, 1, 41, 1, 3, 1, 5, 1, 1, 1, 3, 3, 114, 1, 1, 2, 2, 2, 3, 1, 1, 1, 1, 4, 12, 1, 5, 2, 4, 1, 3, 1, 3, 2, 4, 2, 29, 3, 5, 11, 5, 7, 8, 1, 2, 66, 1, 3, 1, 2, 2, 4, 4, 3, 2, 1, 1, 2, 1, 1, 2, 1, 4, 1, 4, 1, 2, 2, 1, 4, 1, 1, 1, 4, 1, 1, 2, 1, 2, 4, 2, 1, 5, 1, 2, 1, 2, 1, 11, 1, 8, 1, 1, 1, 1, 1, 2, 2, 3, 3, 1, 17, 1, 223, 1, 1, 14, 1, 12, 4, 7, 1, 4, 1, 3, 3, 15, 2, 2, 1, 2, 1, 1, 11, 20, 1, 1, 1, 1, 1, 11, 1, 35, 1, 11, 1, 7, 1, 10, 8, 1, 53, 2, 5, 1, 7, 27, 3, 4, 1, 5, 2, 1, 11, 1, 1, 1, 2, 3, 24, 1, 11, 2, 17, 4, 1, 1, 1, 7, 3, 4, 3, 1, 8, 1, 1, 4, 1, 3, 1, 8, 1, 1, 2, 1, 1, 7, 1, 3, 30, 1, 3, 12, 1, 3, 1, 29, 2, 1, 3, 3, 57, 1, 3, 1, 3, 1, 2, 1, 2, 12, 8, 1, 1, 6, 16, 3, 1, 1, 2, 1, 3, 1, 1, 4, 1, 2, 1, 1, 1, 3, 1, 1, 3, 1, 10, 9, 2, 1, 1, 11, 5, 1, 2, 1, 1, 25, 2, 1, 2, 4, 7, 1, 9, 4, 1, 1, 3, 2, 29, 1, 2, 4, 2, 6, 5, 86, 2, 6, 3, 1, 1, 3, 1, 2, 1, 3, 4, 21, 1, 1, 6, 1, 1, 13, 10, 1, 2, 2, 22, 2, 5, 1, 1, 1, 5, 1, 1, 2, 1, 5, 3, 1, 11, 1, 2, 1, 3, 2, 30, 1, 776, 1, 1, 1, 2, 29, 19, 1, 1, 1, 2, 1, 6, 11, 8, 3, 2, 1, 1, 130, 1, 11, 2, 8, 3, 1, 1, 3, 1, 1, 1, 12, 1, 1, 2, 1, 1, 129, 2, 3, 7, 60, 17, 3, 5, 1, 5, 1, 1, 4, 1, 2, 3, 1, 1, 1, 66, 1, 33, 1, 2, 1, 1, 1, 4, 3, 3, 1, 1, 1, 12, 1, 6, 8, 3, 74, 1, 1, 2, 1, 2, 1, 1, 5, 1, 2, 5, 1, 2, 2, 2, 7, 3, 1, 5, 1, 5010, 1, 3, 5, 9, 1, 3, 2, 3, 4, 2, 21, 1, 1, 2, 3, 1, 108, 1, 8, 3, 1, 3, 2, 9, 1, 1, 1, 2, 2, 1, 9, 15, 2, 1, 5, 1, 3, 3, 2, 2, 1, 1, 1, 6, 1, 1, 18, 2, 2, 2, 3, 1, 1, 19, 1, 3, 2, 13, 1, 69, 4, 1, 2, 1, 7, 2, 1, 1, 1, 2, 103, 3, 1, 20, 1, 1, 1, 2, 2, 1, 4, 1, 1, 6, 2, 10, 2, 1, 15, 2, 12, 21, 1, 77, 2, 1, 1, 2, 5, 1, 14, 1, 2, 7, 1, 1, 5, 7, 3, 6, 3, 84, 3, 2, 3, 2, 2, 6, 15, 4, 2, 1, 2, 4, 1, 14, 6, 21, 1, 48, 5, 1, 1, 2, 7, 3, 1, 17, 2, 1, 1, 1, 7, 1, 2, 17, 3, 5, 1, 1, 3, 87, 4, 1, 2, 1, 19, 391, 2, 1, 3, 4, 2, 1, 2, 5, 2, 1, 1, 1, 2, 1, 2, 9, 2, 4, 2, 1, 1, 3, 3, 1, 2, 3, 1, 2, 1, 5, 17, 7, 1, 1, 1, 1, 6, 7, 1, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 1, 9, 1, 2, 12, 1, 38, 3, 1, 10, 1, 5, 18, 3, 1, 4, 2, 5, 5, 1, 2, 1, 1, 1, 1, 3, 4, 1, 35, 1, 4, 1, 4, 2, 1, 6, 9, 3, 39, 21, 1, 3, 1, 1, 1, 8, 1, 1, 9, 1, 6, 1, 14, 10, 9, 43, 2, 1, 3, 4, 2, 13, 1, 1, 4, 3, 16, 4, 3, 1, 4, 887, 2, 1, 12, 6, 1, 4, 1, 5, 1, 1, 1, 1, 1, 10, 2, 1, 5, 6, 1, 81, 1, 2, 5, 2, 3, 10, 1, 1, 2, 3, 81, 1, 14, 1, 1, 1, 1, 4, 2, 10, 6, 2, 3, 165, 3, 3, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 6, 2, 32, 1, 1, 2, 2, 9, 1, 3, 5, 1, 2, 3, 1, 4, 10, 1, 12, 1, 47, 1, 13, 3, 1, 6, 1, 3, 2, 8, 1, 3, 5, 1, 1, 1, 22, 5, 1, 1, 25, 1, 2, 1, 20, 6, 3, 1, 5, 1, 2, 6765, 1, 3, 3, 4, 4, 1, 1, 141, 1, 11, 9, 1, 3, 1, 3, 2, 5, 1, 9, 1, 3, 1, 4, 1, 9, 1, 3, 4, 4, 3, 1, 2, 10, 2, 2, 1, 2, 3, 1, 5, 1, 1, 2, 1, 2, 1, 18, 176, 3, 2, 5, 1, 229, 7, 11, 1, 1, 2, 2, 1, 4, 1, 4, 1, 12, 2, 35, 1, 17, 1, 4, 1, 1, 3, 2, 2, 1, 1, 1, 1, 1, 2, 2, 4, 2, 4, 1, 2, 1, 1, 2, 1, 1, 1, 5, 6, 1, 9, 1, 5, 1, 7, 3, 1, 2, 4, 1, 2, 8, 2, 2, 1, 1, 24, 5, 3, 1, 2, 1, 27, 1, 6, 1, 7, 2, 8, 2, 1, 1, 1, 20, 1, 3, 1, 1, 1, 1, 1, 1, 1, 5, 9, 2, 1, 11, 2, 1, 34, 2, 2, 20, 1, 2, 1, 2, 1, 2, 7, 25, 1, 2, 6, 4, 3, 15, 3, 2, 14, 4, 1, 3, 2, 1, 4, 1, 2, 188, 5, 2, 13, 1, 9, 7, 4, 1, 6, 3, 2, 1, 4, 1, 6, 1, 23, 1, 2, 3, 3, 4, 2, 1, 16, 8, 1, 1, 1, 1, 8, 2, 32, 13, 6, 2, 1, 2, 1, 2, 2, 4, 8, 4, 23, 14, 1, 1, 1, 1, 2, 18, 1, 27, 6, 2, 1, 1, 7, 1, 17, 3, 3, 1, 3, 10, 1, 2, 2, 1, 2, 2, 1, 4, 12, 1, 1, 6, 1, 3, 3, 1, 3, 28, 1, 2, 4, 69, 1, 4, 1, 2, 1, 2, 2, 2, 2, 2, 1, 2, 2, 10, 1, 228, 3, 2, 1, 1, 16, 3, 1, 5, 6, 1, 5, 1, 1, 5, 2, 3, 3, 1, 5, 3, 20, 1}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 9959369 3469416742 2561522715 4059207546 8105716008 8662254751 3414903535 4863516924 8995045894 1594911794 2999138412 7966693860 6630253803 3646411607 8617834136 4847788848 4355704682 1629294315 4364960631 5616047925 5682715434 9764070573 8859608437 2750351343 3243064233 0402825193 8609921005 0773943264 7732819778 5277068042 3604116315 4355837913 3860793200 3202270502 9158935952 9920208507 3279268940 3549957572 4473947042 3145368758 1429615858 9739860446 3603860071 3568622493 4758942754 2638571071 8905460987 9685608908 1637327580 3391329955 1405553385 1418018215 6217392258 2953179525 6263067941 1168220446 3668550217 2179732023 8229993735 6021930120 9547952677 3150985892 9082839723 5530657776 2738026629 9285857057 9784683834 8012828053 2546099739 9888657981 6688743405 9211039058 4229616854 8398909325 8819596839 5363030077 4132892385 4836048031 0026915326 5156801172 6566881404 9568011071 2285447257 2453997462 7474335145 0449601559 6679886909 9072544921 0420825378 6350295168 2060712126 1875785260 5072144044 1023609316 8813801868 3709468621 9956326457 0991929165 8158619999 0752173286 8751353826 1147805162 9582116797 6604724997 1438592425 3553119227 7801171592 0786429547 7578073167 3949060999 0841328981 4917678104 4098038827 1092594929 1894271923 6676100068 5587716686 5853868386 8008284084 4958748319 5713435038 9053198361 6168659692 0059730364 1487353517 5037138872 9278563996 3994056627 7860726841 7749847119 5518548971 4068403331 8672540703 5636450268 6798785113 1723111339 5838397473 5829174142 0388496198 0906158585 7625440814 3376141900 1235357726 5854364760 9852148607 2037145597 4118147939 4202547620 5509173402 6809457506 7994595323 8037784698 2922568994 2731289379 7664692238 4795140627 6169338174 4365096836 2376193712 7120088337 7981352252 2176208982 2933802684 9521770263 8278355410 9447057750 5171197748 9753362575 7107619242 4828795481 8741037087 5084703537 3350383870 4735845571 5275675698

Cubic Polynomial

We now consider the cubic P(x)= v·x3 + (u·F-c1·v)·x2 + (c4·v-d·F+u)·x - d, which we express as: z1·x3 + z2·x2 + z3·x + z4, where:

We need to prove that this cubic has no integer roots r such that r·F+1 is a non-trivial factor of N. Clearly r (if it exists) must lie between 1 and R.

The real roots of P are:

There are no integer roots of P in the interval (1,R), so the proof of primality is complete.