Primality Certificate for (3164^3061-1)/3163

Andy Steward10,711 digits25 March 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.865746% factorization of N-1:

From Factorisation
31642 · 2 · 7 · 113
Φ23 · 5 · 211
Φ32917 · 3433
Φ413 · 770069
Φ531 · 1531 · 2112254761
Φ63 · 277 · 12043
Φ918517249 · 54180421887169
Φ105 · 11 · 1821570441011
Φ121157749 · 86562829
Φ151364598931 · 7357825588101955111
Φ17501841 · 31925434961753 · 3301838189862829 · 1907489389602414314393
Φ183 · 19 · 24247 · 725915351399167
Φ2061 · 164650070199599745449571781
Φ306211 · 2980231 · 90768481 · 5979742921
Φ34443 · 919 · 566529395835240589 · p33
Φ3637 · 73 · 3061 · 4357 · 133669 · 189073657 · 1105606290178674721
Φ45181 · 1645500151 · p73
Φ512857 · 21577693 · 40770217 · 1655953069 · 28598624857005794317 · 171530726416933013016361 · p42
Φ6057422075401 · 6376992211981 · p33
Φ681021 · 29173 · 235123261 · 21493450876357886715509053475333440217309 · p56
Φ851097062205726240321 · 130116060746784310913174166945461 · 430134683999497421859406028877881 · p142
Φ905666319021840452281 · 51067388823505531299355269871 · p37
Φ102103 · 44371 · 11161238089796083627561 · p84
Φ153307 · 428401 · 696457 · 573005279516944507507 · c302
Φ1701871 · 3177100081 · c212
Φ180541 · 1621 · 55980593821 · p152
Φ204409 · 5101 · 123217 · c213
Φ255861391 · 2520524743407455227852891 · c418
Φ30620809 · 20783315937408987591289 · c310
Φ34027541 · 248201 · 350248928041 · 10541427506321155321 · p408
Φ510c449
Φ612613 · 2986352339661253 · 34237422745253566217385913 · c629
Φ76510711 · 10426951 · 105804224641 · 8737117015012028641 · 53583120412505923573981 · p1281
Φ10204286879461 · c887
Φ153016831 · 126687061 · 253124611563506810311 · c1312
Φ30606121 · 813961 · c2679

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

2 2280237647 7965808407 7660671952 5978297924 0186357581 7762729901 7281567281 0192255905 0074477300 2937390572 2767718060 3854221882 2956488124 3925429159 0267160613 5792499253 5461026423 2149796675 7284551199 5919818186 2635630608 0704606032 8701935960 5879440386 3317493974 8718778561 7507951612 9529899044 0216753862 1560354240 6457655620 9674467941 2784284983 2096606319 8325249881 9128023537 9967045591 7781664863 3377964079 1704708731 2117893342 1375798030 5442115224 7601455027 9699209557 2764925768 0651748511 6341735936 4026529413 0488580323 1455764888 1615236973 1957687150 1459299246 0813985144 7760569464 4046546062 1016098587 4961970866 3650669333 7028544418 9121169109 8707428742 7031354317 0400291744 8516230979 9281885151 1456022171 6874648843 3080509486 0356687981 4880755252 3705263733 1757045259 1223074975 6874785578 7314904182 3282591583 8639268399 1300445787 0045869015 5846007334 9996435192 9047711228 6177613161 0080308661 0730626674 5712751881 0237112822 0891075778 1826220309 5125688778 2964199452 0990428165 3218525913 1063310221 5273655869 0477594728 2685770251 5194610516 6037785726 1142380567 0298635409 8950178026 5720342055 8937011633 5428009506 8014347863 2402449319 8598646016 2624930901 7687329168 0298594560 9586997114 3316116479 8907509520 2679374867 4926076640 6386081977 1884014178 9877590735 7827198195 4643814944 9387909186 8196451300 2725597116 1143123459 4797403701
42482429 7262774373 3858000445 6124793005 0155282847 9001605985 1868661650 1607664324 6916653731 3815205692 9587286353 1001135487 9852133197 1412397551 4301919334 4801805740 2423566716 6517710196 7593002669 9105395489 7765207813 7066680717 4198300505 4281563825 9786902212 4860492567 8956702272 5647284517 0014478645 1056386565 6675258999 4496891018 1652437102 1623645815 4672842833 5855854164 8751268938 0344825521 5116024841 3335437983 6095934341
20 9089878709 3259155066 6151128221 2563040731 1659330336 2284754678 2642379736 7514110735 0397979282 7310878418 3160367336 7024192628 5631555742 9177290304 6048506261
16 8589484333 9092296827 4819626871 3310503411 7677801794 2160461242 2471790293 3285721292 2220082767 0403878710 2304173743 8849675883 2386518920 0521674361
1995 5178644508 0748369885 2724272675 4981783713 0342897895 3825886041 2733327754 1588726577
340 1722592994 6422035659 9140334171 9162985183 7201953222 3150207490 9384991691
676018 4848420939 8540688627 7075687109 7895314866 9641011433
49 8237642812 3214691783 6086148511 9594179561
2 1493450876 3578867155 0905347533 3440217309
3501311 2682540677 2115606772 0742944071
437 2246822489 5858393861 3369549229
430 1346839994 9742185940 6028877881
275 4792353776 0843186561 0766870701
130 1160607467 8431091317 4166945461
510673888 2350553129 9355269871
1646500 7019959974 5449571781
342374 2274525356 6217385913
25205 2474340745 5227852891
1715 3072641693 3013016361
535 8312041250 5923573981
207 8331593740 8987591289
111 6123808979 6083627561
19 0748938960 2414314393
5 7300527951 6944507507
2 5312461156 3506810311
2859862485 7005794317
1054142750 6321155321
873711701 5012028641
735782558 8101955111
566631902 1840452281
110560629 0178674721
109706220 5726240321
56652939 5835240589
330183 8189862829
298635 2339661253
72591 5351399167
5418 0421887169
3192 5434961753
637 6992211981
182 1570441011
35 0248928041
10 5804224641
5 7422075401
5 5980593821
5979742921
4286879461
3177100081
2112254761
1655953069
1645500151
1364598931
235123261
189073657
126687061
90768481
86562829
40770217
21577693
18517249
10426951
2980231
1157749
861391
813961
770069
696457
501841
428401
248201
133669
123217
44371
29173
27541
24247
20809
16831
12043
10711
6211
6121
5101
4357
3433
3061
2917
2857
1871
1621
1531
1021
919
613
541
443
409
307
277
211
181
113
103
73
61
37
31
19
13
11
7
52
33
22

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

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

Express N in base F

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 = 3369855 5970235720 1569808765 3481767072 0176068117 9186000908 4211357039 3956208205 3482714057 0115496216 7590091323 7549460165 4131635766 8979289496 9178509456 4843899433 9957783371 8033296273 6505507029 3431624186 4490979057 3218251718 5992847665 8956430947 6318374238 9921014663 6882123751 0008551392 1249974108 5455540654 7961332330 8943855941 8333771046 6314322905 0451346118 3069256756 6874942371 9087534732 6067890783 9519420560 9774086265 6781994886 6983175660 6065052201 6290287954 0599062729 3737758822 7775734104 3228318094 4187970385 7491640100 5084376997 1579429415 4016521224 2070600263 9890228424 6576154819 9945564367 7570185475 1474212059 0848152042 6391293125 1753971306 1195030329 0428695470 5568600767 1544375187 2157910518 3802274569 8364540254 6859082590 8393987754 6701406471 1556743446 5844051662 3290484062 3864522619 0249976558 1817935989 6115349662 2210560524 6732576029 3320959956 4806746568 0858574732 7677316549 9741915825 8474144969 7685307605 9964682784 9310782497 9195427775 9199189927 2235998754 0899762056 3779762479 5013759948 6612247454 8453238296 0247173836 7228690703 6300089092 9659331453 6879498691 1832818856 9414089430 1662093724 8322768033 6510269609 1096422434 0339190218 1157625551 8586583801 2069677074 6207373553 5337772446 6657886658 5461871145 3144232831 2893100003 7247604620 8379676180 9752391611 5480162892 5234735666 2964601329.

With those constraints, the unique continued fraction is: {0, 10, 14, 8, 1, 1, 10, 26, 1, 1, 5, 2, 15, 2, 2, 1, 22, 5, 2, 2, 1, 4, 4, 1, 10, 1, 6, 1, 5, 5, 11, 5, 1, 1, 4, 1, 876, 2, 5, 1, 3, 25, 1, 1, 1, 63, 1, 1, 1, 5, 1, 1, 3, 3, 1, 2, 3, 1, 5, 1, 1, 1, 1, 13, 1, 9, 1, 2, 1, 3, 2, 1, 1, 1, 2, 3, 1, 204, 1, 1, 5, 1, 2, 1, 1, 1, 7, 1, 1, 1, 1, 1, 2, 11, 1, 10, 2, 2, 2, 3, 1, 7, 2, 4, 17, 2, 1, 26, 1, 6, 3, 3, 1, 1, 8, 1, 2, 8, 1, 220, 1, 1, 9, 1, 4, 1, 1, 1, 61, 1, 21, 15, 13, 2, 6, 1, 3, 26, 6, 1, 4, 3, 4, 23, 1, 2, 2, 47, 2, 2, 98, 4, 1, 1, 5, 1, 1, 5, 1, 5, 2, 1, 2, 2, 2, 1, 1, 8, 2, 1, 1, 1, 3, 2, 1, 2, 2, 1, 4, 1, 13, 2, 2, 5, 12, 1, 2, 1, 1, 7, 1, 1, 1, 2, 1, 4, 1, 4, 5, 4, 2, 1, 1, 5, 1, 1, 2, 20, 1, 2, 1, 1, 1, 12, 1, 2, 2, 4, 3, 17, 6, 2, 1, 3, 42, 1, 7, 77, 1, 2, 5, 10, 5, 2, 1, 1, 3, 1, 1, 1, 5, 9, 1, 41, 11, 2, 1, 30, 1, 8, 6, 1, 1, 3, 1, 1, 2, 1, 1, 4, 2, 2, 1, 3, 3, 73, 1, 6, 7, 1, 1, 1, 1, 1, 1, 2, 2, 9, 3, 9, 1, 1, 1, 1, 5, 1, 96, 6, 1, 2, 3, 1, 3, 1, 2, 14, 2, 2, 1, 1, 1, 2, 2, 12, 1, 2, 14, 1, 1, 1, 1, 1, 3, 46, 3, 14, 3, 1, 1, 17, 1, 1, 1, 3, 2, 1, 2, 1, 1, 1, 1, 1, 64, 4, 23, 2, 2, 5, 5, 1, 2, 1, 1, 1, 1, 1, 2, 9, 24, 9, 1, 2, 1, 3, 2, 1, 11, 1, 4, 1, 1, 2, 3, 1, 12, 1, 1, 4, 1, 3, 2, 10, 2, 1, 1, 2, 1, 1, 2, 1, 6, 10, 1, 2, 1, 1, 2, 2, 1, 3, 50, 1, 3, 1, 9, 1, 1, 1, 8, 1, 1, 25, 1, 4, 29, 1, 2, 5, 13, 1, 4, 1, 4, 1, 10, 1, 1, 1, 1, 1, 6, 1, 1, 1, 3, 7, 3, 1, 1, 1, 2, 1, 2, 2, 1, 1, 3, 1, 1, 2, 2, 17, 1, 7, 1, 11, 1, 1, 1, 3, 1, 1, 2, 2, 6, 2, 18, 1, 11, 2, 1, 2, 76, 6, 2, 1, 1, 7, 56, 24, 1, 15, 1, 1, 2, 11, 1, 7, 6, 4, 1, 41, 1, 2, 798, 9, 9, 1, 6, 3, 1, 7, 1, 2, 1, 4, 16, 2, 3, 2, 3, 6, 1, 11, 1, 2, 1, 2, 2, 1, 1, 2, 2, 3, 4, 2, 2, 1, 3, 1, 7, 1, 12, 1, 5, 1, 4, 2, 1, 75, 4, 2, 1, 2, 1, 4, 1, 4, 1, 3, 9, 8, 2, 1, 3, 1, 4, 1, 1, 8, 1, 6, 1, 11, 1, 1, 2, 1, 1, 1, 6, 7, 3, 1, 7, 1, 6, 2, 1, 5, 3, 1, 16, 1, 2, 1, 39, 2, 1, 4, 1, 6, 9, 1, 7, 3, 1, 3, 1, 8, 5, 1, 3, 1, 1, 2, 1, 1, 5, 1, 1, 2, 1, 12, 1, 3, 5, 1, 1, 9, 3, 84, 4, 1, 36, 15, 4, 1, 3, 1, 4, 1, 41, 2, 2, 32, 2, 1, 2, 2, 1, 1, 3, 1, 1, 11, 1, 3, 1, 1, 1, 8, 2, 6, 2, 15, 1, 1, 1, 1, 69, 1, 2, 200, 1, 1, 2, 4, 8, 1, 1, 1, 2, 1, 1, 6, 5, 1, 4, 1, 6, 1, 9, 1, 4, 10, 41, 1, 2, 1, 5, 2, 1, 1, 59, 2, 5, 1, 3, 3, 11, 1, 1521, 4, 9, 44, 1, 22, 2, 1, 1, 2, 6, 2, 12, 2, 1, 1, 1, 6, 1, 6, 1, 1, 1, 4, 4, 1, 8, 1, 194, 1, 2, 1, 2, 3, 1, 2, 4, 1, 6, 1, 3, 2, 2, 3, 1, 5, 3, 2, 7, 4, 2, 2, 2, 3, 1, 1, 1, 1, 3, 1, 5, 1, 1, 3, 1, 5, 42, 5, 1, 1, 2, 1, 5, 3, 1, 12, 1, 27, 1, 3, 1, 4, 1, 68, 2, 10, 2, 1, 1, 1, 2, 3, 1, 29, 8, 67, 2, 2, 2, 4, 5, 1, 1, 2, 1, 6, 2, 6, 1, 5, 1, 1, 3, 1, 1, 4, 3, 7, 2, 1, 2, 1, 4, 1, 8, 1, 1, 1, 1, 27, 1, 1, 2, 4, 1, 1, 4, 1, 1, 3, 1, 2, 8, 1, 1, 1, 58, 1, 4, 1, 7, 1, 1, 1, 10, 1, 1, 3, 15, 1, 1, 8, 1, 48, 4, 8, 4, 3, 9, 1, 8, 1, 9, 1, 55, 1, 3, 1, 4, 3, 3, 13, 1, 9, 1, 3, 2, 6, 4, 13, 1, 17, 2, 4, 3, 3, 1, 1, 1, 52, 1, 14, 1, 2, 6, 2, 3, 1, 3, 2, 2, 1, 74, 20, 1, 5, 58, 2, 1, 1, 1, 4, 1, 11, 1, 47, 1, 1, 4, 2, 1, 3, 2, 2, 2, 3, 10, 1, 1, 1, 51, 1, 1, 2, 1, 17, 6, 1, 20, 2, 9, 1, 1, 1, 1, 21, 11, 26, 1, 77, 2, 45, 1, 12, 8, 2, 1, 1, 1, 1, 9, 4, 3, 1, 2, 7, 1, 2, 1, 1, 1, 2, 3, 2, 1, 1, 2, 10, 1, 12, 16, 2, 1, 1, 2, 1, 2, 5, 10, 1, 4, 4, 1, 1, 2, 1, 2, 3, 2, 3, 5, 1, 2, 2, 1, 2, 3, 2, 1, 2, 72, 1, 1, 11, 1, 2, 1, 18, 2, 1, 1, 2, 3, 4, 2, 104, 1, 3, 3, 3, 2, 1, 1, 1, 4, 3, 1, 3, 1, 2, 3, 1, 1, 1, 1, 1, 1, 3, 2, 3, 1, 1, 3, 18, 26, 1, 6, 1, 2, 1, 3, 2, 3, 2, 2, 3, 15, 1, 47, 1, 1, 1, 3, 1, 13, 1, 4, 9, 1, 3, 1, 10, 1, 4, 2, 7, 1, 15, 1, 1, 1, 1, 11, 2, 3, 6, 7, 1, 1, 1, 1, 3, 8, 1, 82, 1, 2, 1, 5, 19, 1, 4, 4, 26, 1, 4, 3, 56, 2, 1, 14, 3, 14, 5, 2, 1, 2, 2, 9, 10, 1, 2, 25, 2, 2, 2, 8, 2, 1, 1, 1, 1, 8, 2, 25, 1, 2, 2, 1, 2, 2, 3, 2, 7, 9, 2, 1, 15, 1, 3, 1, 2, 1, 8, 6, 6, 1, 2, 7, 1, 2, 1, 3, 3, 1, 1, 2, 2, 14, 2, 3, 2, 2, 1, 1, 15, 1, 44, 1, 2, 1, 1, 1, 5, 2, 1, 2, 1, 42, 1, 9, 1, 4, 25, 4, 2, 1, 4, 16, 6, 3, 1, 1, 4, 7, 75, 1, 2, 37, 1, 1, 2, 1, 13, 1, 1, 1, 3, 388, 1, 1, 5, 1, 4, 1, 3, 1, 2, 1, 64, 1, 2, 1, 5, 2, 2, 4, 9, 1, 4, 1, 29, 2, 1, 1, 6, 1, 5, 1, 1, 4, 1, 19, 2, 2, 3, 1, 5, 1, 5, 1, 3, 1, 1, 2, 1, 1, 17, 23, 4, 1, 2, 10, 9, 4, 3, 2, 2, 2, 5, 2, 1, 2, 1, 1, 5, 1, 2, 3, 3, 6, 1, 8, 1, 7, 1, 2, 4, 1, 6, 1, 118, 1, 21, 1, 16, 1, 1, 4, 1, 3, 5, 2, 2, 2, 1, 1, 2, 4, 1, 3, 3, 1, 7, 5, 2, 1, 2, 1, 103, 1, 9, 19, 7, 2, 1, 2, 3, 4, 5, 1, 4, 4, 2, 1, 1, 5, 1, 1, 1, 1, 7, 3, 4, 1, 3, 1, 30, 5, 6, 11, 1, 18, 1, 4, 2, 4, 1, 1, 3, 1, 3, 2, 33, 355, 1, 1, 2, 5, 1, 1, 2, 1, 19, 1, 2, 1, 7, 1, 2, 2, 2, 1, 5, 3, 1, 2, 1, 4, 1, 5, 1, 25, 17, 50, 5, 2, 3, 1, 27, 1, 1, 1, 11, 1, 5, 11, 1, 3, 1, 2, 5, 1, 5, 2, 2, 1, 1, 3, 1, 1, 2, 2, 2, 1, 5, 1, 1, 8, 2, 2, 8, 2, 2, 1, 18, 26, 19, 17, 1, 17, 2, 1, 1, 2, 1, 1, 4, 1, 2, 1, 3, 3, 1, 3, 1, 1, 3, 1, 2, 3, 3, 2, 2, 2, 1, 3, 13, 2, 3, 2, 7, 1, 2, 1, 1, 1, 2, 1, 3, 2, 1, 12, 1, 4, 2, 5, 1, 10, 7, 1, 2, 1, 1, 5, 1, 207, 1, 4, 2, 1, 1, 1, 2, 3, 3, 1, 8, 13, 1, 6, 1, 11, 6, 11, 20, 11, 1, 28, 3, 2, 1, 3, 1, 5, 3, 1, 1, 1, 11, 3, 3, 1, 4, 1, 8, 1, 24, 1, 5, 2, 1, 1, 5, 1, 8, 1, 2, 3, 3, 1, 3, 2, 3, 1, 5, 6, 2, 1, 24, 1, 2, 1, 1, 1, 5, 1, 16, 2, 12, 4, 1, 2, 1, 1, 2, 2, 1, 1, 2, 16, 1, 1, 1, 12, 1, 24, 12, 2, 1, 15, 1, 1, 3, 4, 6, 65, 1, 2, 6, 1, 25, 1, 38, 1, 1, 3, 1, 2, 2, 1082, 1, 35, 12, 2, 1, 9, 1, 1, 1, 5, 6, 1, 13, 18, 1, 1, 3, 1, 2, 9, 2, 1, 122, 2, 16, 1, 1, 3, 1, 1, 1, 1, 1, 5, 3, 2, 1, 1, 1, 1, 24, 1, 3, 2, 3, 1, 6, 1, 3, 1, 8, 3, 1, 2, 1, 1, 2, 1, 1, 5, 1, 2, 3, 1, 5, 1, 2, 172, 8, 1, 2, 7, 2, 1, 1, 2, 3, 12, 2, 1, 1, 21, 3, 10, 2, 1, 21, 1, 11, 2, 1, 1, 3, 2, 1, 1, 11, 5, 4, 6, 2, 1, 1, 18, 1, 1, 3, 7, 8, 1, 1, 3, 1, 3, 2, 2, 55, 1, 1, 8, 8, 2, 1, 2, 2, 1, 2, 1, 1, 1, 1, 3, 2, 1, 2, 1, 2, 5, 1, 24, 18, 3, 1, 8, 1, 2, 2, 9, 2, 1, 3, 1, 1, 1, 2, 1, 1, 8, 2, 5, 2, 3, 1, 7, 2, 10, 1, 3, 1, 5, 7, 1, 4, 1, 1, 1, 100, 1, 1, 5, 1, 2, 5, 3, 1, 154, 2, 2, 1, 1, 2, 1, 43, 1, 11, 1, 8, 5, 1, 2, 1, 1, 4, 32, 1, 10, 2, 5, 1, 1, 1, 2, 1, 31, 4, 9, 1, 16, 4091, 3, 2, 1, 2, 2, 2, 4, 3, 3, 5, 5, 10, 8, 2, 1, 1, 1, 5, 7, 1, 12, 3, 25, 1, 9, 46, 2, 8, 1, 7, 1, 7, 1, 1, 13, 1, 1, 1, 6, 1, 1, 2, 2, 2, 2, 2, 7, 1, 9, 4, 1, 1, 1, 1, 3, 3, 1, 9, 1, 8, 2, 6, 1, 6, 1, 19, 1, 1, 1, 19, 2, 1, 54, 1, 1, 2, 8, 1, 2, 2, 4, 16, 1, 1, 9, 36, 1, 10, 16, 1, 2, 1, 3, 25, 4, 1, 561, 6, 1, 1, 1, 1, 4, 1, 4, 4, 1, 1, 1, 3, 146, 2, 1, 1, 5, 2, 1, 1, 2, 1, 5, 3, 1, 3, 1, 4, 1, 11, 3, 2, 1, 5, 9, 4, 2, 15, 4, 1, 1, 1, 3, 1, 7, 1, 4, 1, 2, 9, 3, 3, 1, 1, 1, 1, 1, 11, 2, 26, 1, 2, 1, 5, 1, 1, 1, 1, 1, 1, 1, 6, 1, 8, 2, 1, 3, 3, 6, 1, 1, 1, 1, 1, 3, 5, 1, 9, 1, 3, 3, 1, 1, 2, 1, 271, 1, 4, 2, 1, 2, 3, 2, 2, 8, 20, 2, 1, 36, 1, 1, 74, 2, 3, 3, 4, 1, 3, 4, 2, 2, 8, 1, 8, 2, 18, 1, 6, 1, 284, 1, 21, 2, 1, 3, 2, 1, 2, 40, 1, 10, 2, 1, 2, 10, 1, 5, 1, 1, 2, 2, 1, 4, 5, 2, 16, 1, 2, 1, 5, 2, 1, 1, 51, 11, 2, 1, 13, 1, 2, 4, 3, 2, 1, 1, 1, 2, 2, 1, 1, 4, 3, 2, 4, 1, 1, 7, 13, 1, 1, 4, 3, 1, 1, 5, 3, 1, 2, 18, 1, 5, 2, 3, 3, 2, 3, 5, 1, 2, 5, 5, 1, 6, 1, 1, 6, 2, 2, 57, 1, 1, 1, 3, 1, 1, 1, 19, 2, 1, 3, 1, 6, 1, 11, 1, 1, 2, 1, 2, 1, 10, 19, 1, 6, 4, 1, 4, 1, 3, 4, 1, 1, 3, 1, 3, 1, 1, 1, 1, 4, 1, 93, 2, 2, 19, 31, 9, 1, 2, 3, 1, 2, 1, 15, 9, 3, 1, 3, 1, 19, 3, 1, 1, 6, 4, 22, 1, 1, 2, 1, 1, 4, 26, 1, 2, 1, 1, 1, 18, 2, 1, 1, 4, 1, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 9, 1, 1, 2, 3, 8, 1, 1, 1, 1, 1, 4, 1, 21, 1, 19, 2, 5, 5, 1, 25, 2, 1, 4, 1, 5, 1, 69, 1, 2, 1, 2, 12, 1, 3, 4, 16, 12, 2, 2, 2, 7, 1, 1, 3, 207, 2, 8, 2, 1, 9, 4, 1, 2, 1, 1, 1, 1, 1, 14, 28, 1, 3, 7, 10, 12, 2, 1, 1, 2, 7, 1, 3, 1, 10, 7, 1, 1, 3, 1, 1, 1, 3, 35, 6, 4, 1, 1, 2, 2, 1, 3, 3, 2, 2, 1, 2, 23, 4, 2, 367, 1, 3, 1, 1, 4, 1, 3, 2, 1, 1, 2, 4, 1, 1, 3, 2, 2, 1, 1, 2, 1, 3, 4, 2, 5, 2, 4724, 1, 1, 14, 2, 1, 1, 69, 2, 2, 3, 17, 6, 7, 1, 5, 6, 1, 3, 3, 1, 3, 8, 1, 2, 1, 1, 1, 2, 8, 2, 1, 1, 1, 2, 3, 1, 3, 8, 1, 3, 1, 2, 2, 3, 1, 1, 5, 1, 2, 2, 2, 2, 5, 1, 5, 1, 21, 4, 37, 1, 1, 5, 1, 1, 1, 1, 2, 1, 1, 32, 2, 1, 1, 2, 1, 25, 1, 1, 72, 1, 1, 1, 1, 1, 30, 1, 1, 13, 1, 7, 1, 1, 3, 1, 2, 9, 1, 2, 1, 1, 10, 1, 1, 2, 7, 30, 15, 1, 3, 3, 1, 8, 1, 3, 1, 3, 1, 1, 1, 99, 5, 1, 1, 2, 1}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 1538845004 3436163632 4468516765 6108183855 7170097791 0349106369 5569716898 8670514338 2074812426 8448305402 0043414223 9926887696 0976830742 1007858259 4579671951 7957780527 5982227396 4676463541 7970984459 7149467121 0943136744 5393349528 5039805923 6430095836 7909893348 9331343520 3562306541 9303236549 5142676925 1264056587 1075512960 4188243377 8729830266 3432021272 4948691138 0732230144 9213924934 5612492437 3942622151 5390519857 1533283577 9103550837 7193817072 2432073121 8826483019 8094042547 3933341588 7023026755 2704100769 7929516765 2660796516 4192861257 0069013815 0141373280 4590085716 9804356789 7971352431 6858949154 1344505541 6831181977 1926481868 9710941840 0031551463 1273521438 4392984323 5936743729 6539522769 6080741963 3019457684 1716649468 3305183444 2926518421 1303028589 8126875664 5438509069 9587312763 0472806353 5668808487 9492988637 3339553012 1357915221 7161654333 7575089263 6210846325 1217723858 7990318859 0545595161 8701218488 0138340050 8182926284 5787687370 0149049330 7253726295 2560304247 8732238194 4007909215 9499702996 2821934037 8386945063 1513737282 0962848180 8643261662 7040286429 7620089461 5811722287 9578153763 0912944455 5499900350 9130085077 0744632559 5198147302 4865166740 6251868248 0168290595 2961231498 7478643555 0557531472 4923195008 1545275674 8252079911 5826727360 3569414199 3610885011 9075218419 7706354887 7518905053 5378288608 5057954096 3759573502 0585249171 9301971371 2157860893 0747480981 4050019868 5514672763 4149983803 8594999569 4972020667 5259869906 2823692657 0764630793 8570004129 9704664568 1968675439 7859429539 6170475514 7301946518 9600294336 8101485863 0454097980 9891745249 6369419320 7186626680 3135153956 4169994706 9895523900 4766494960 5143115675 1944564327 8671708132 7012164019 9854817743 4121396503 0662231336 2439604089 1888733418 9534969674 7669309352 6530231716 7836502691 4129827370 6221761411 1275582425 2015684846 1896991102 6365355923 0129602838 6213592258 2178699488 3196376889 5348180479 5383538598 7304704775 2701828592 9850287046 9501974531 8033337020 2480097098 1407649430 3539784985 4883900620 7319168783 6095607605 4116079766 2890248660 6073859027 5353173269 6382454384 8822727427 1169807198 5037791783 5528496547 5207246160 4643768333 7320824440

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.