Primality Certificate for (9926^2521-1)/9925

Andy Steward10,072 digits10 July 2001
Originally by David Broadhurst & Sean Irvine 2001

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

From Factorisation
99262 · 7 · 709
Φ23 · 3 · 1103
Φ337 · 43 · 61933
Φ4353 · 279109
Φ55 · 11 · 71 · 2486106909971
Φ63 · 19 · 31 · 127 · 439
Φ744269 · 21606761038907777503
Φ817 · 113 · 5053237595537
Φ9199 · 4806097187778219964447
Φ1041 · 11964131 · 19787381
Φ1213 · 822793 · 907534489
Φ141179263 · 810944628735834377
Φ158650651 · 10891849209230322544147801
Φ183 · 73 · 2917 · 30097 · 2524807 · 19702153
Φ2014341 · 9165357241 · 716911140870115921
Φ21421 · 673 · 27608407 · 832304719 · 140483923916980904758314559
Φ2436697 · 5094601 · 504026566984153922833
Φ2829 · p47
Φ30181 · 331 · 331359241 · 4747142433340238401
Φ352381 · 5877395994943998907271 · 18692440451313813543101 · p49
Φ36p48
Φ40p64
Φ426763 · 11131 · p41
Φ45541 · 24159511 · 1303352374361581 · p71
Φ5611594521 · p89
Φ6061 · 61 · p61
Φ631009 · 1017715609 · 192556968460713896142223 · 13160701324088372631877591063702514094230340186379 · p60
Φ70631 · 3498504482761 · 286294112770513285047821 · p58
Φ721056956063617 · 2503632449809 · 6943339004708377 · 9644922583019106049 · p37
Φ841301497 · 1617757 · 440693388065411106277 · 160073073932606436308980193929 · p34
Φ901548541 · 100144711 · p82
Φ105211 · 2521 · 67164918451 · c176
Φ120241 · 4412993041 · 164378161801606958077431070441 · p87
Φ126883 · 8841096918495905156384610478787209 · p107
Φ14052081 · 922601 · 29341320208068640057861 · 2049877448116088359265702941 · c132
Φ168c192
Φ18018980353810981 · c179
Φ2101051 · c189
Φ252757 · 6301 · 1284101331500187649 · c263
Φ280281 · 36193361 · 122841041 · p366
Φ315687961 · 1455796402404894531451 · 2423864264756920741820110777280022061 · c513
Φ3606112753054473601 · c368
Φ42018061 · 31685221 · 4813927021 · 24607510296303901 · c346
Φ50412601 · p572
Φ630350281 · 25499819956546438651 · 2937984333367677126751 · 3416257792192746783729511 · 573516852561277743452556601 · c478
Φ84032317321 · c760
Φ12602436032686501 · 3663799169941 · 118516459680868316761 · 1350697619803053372481 · c1085
Φ2520334787041 · 4216028484338906761 · c2275

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 :

27 2326555190 8078655699 5595382757 2090526061 9490622124 9996193905 0462264576 8162762307 5087610611 6132631186 9482029123 8805484101 0280529657 3493812215 3396154456 7061760422 7381219892 3680137928 9189926505 1504835836 9105248544 1402200407 4378730314 0827269650 8169468746 4460543178 7691176618 7539408514 7636978598 8371470707 6329570121 2384911999 6289288039 2118159368 9520546896 8858829656 7930608895 4921425477 0795882459 9842169277 8458908528 2606717669 8479774928 1533938754 4273922167 9288098769 7601112568 3078599741 8336364427 8957197478 2815759197 3834752620 9626027321 6618879319 4905762462 8917338997 6062197001
392329 7847959077 1623363960 9830333293 4692757097 3551506285 4972506155 7514356403 9776997188 5353731947 0783133796 4171748463 8542157507 2557686955 9144415976 2467372987 1272261691 2411828218 4557826982 3117035270 1895793805 7772494222 1948070362 2942543470 7980085429 0714604390 0900654842 1039136249 4189019035 4798419602 5909787801 8108746031 2688675068 6304435054 4264259521 7927952414 9318445721
9804081 1197149043 1362897448 8449378874 7503178686 1360669881 7153108782 6868303837 9327334315 6433718450 1121671083
721655106 6495850938 5932793970 7471711025 8887977699 5852289427 2696949149 3583636121 9900707481
4510059 4108070566 5783701486 9528922647 4131568022 3248006977 9159776410 8882596099 0092424681
53 9550096590 4036600239 7143162761 6684133963 7230945228 0650896771 6739017921 6017277451
4 9117386290 4116406545 2338240851 3169741381 4737359354 7582190173 1890064471
8879 4963649180 2049972450 2999522829 5295445510 8914671309 6816444401
2 3863199287 9384118057 5741512916 6959328889 7919512088 9019315581
2941154753 5746017702 1619140581 0032391397 2530571567 7795892813
13240428 4953649974 9241642168 9916825412 7583914414 3188356141
1316070132 4088372631 8775910637 0251409423 0340186379
319836770 3441420996 7724593055 7290274222 5272092501
91472647 7631748389 8525278195 1197803952 3476580801
3154229 2011985741 7682149945 2570940766 1278754769
1 2152382515 5027730533 5997954668 0926544967
4721585 7152098731 5328605490 7800271129
2423864 2647569207 4182011077 7280022061
8841 0969184959 0515638461 0478787209
5633 4082631922 8560940982 6545248693
1643781618 0160695807 7431070441
1600730739 3260643630 8980193929
20498774 4811608835 9265702941
5735168 5256127774 3452556601
1404839 2391698090 4758314559
108918 4920923032 2544147801
34162 5779219274 6783729511
2862 9411277051 3285047821
1925 5696846071 3896142223
293 4132020806 8640057861
186 9244045131 3813543101
58 7739599494 3998907271
48 0609718777 8219964447
29 3798433336 7677126751
14 5579640240 4894531451
13 5069761980 3053372481
5 0402656698 4153922833
4 4069338806 5411106277
1 1851645968 0868316761
2549981995 6546438651
2160676103 8907777503
964492258 3019106049
474714243 3340238401
421602848 4338906761
128410133 1500187649
81094462 8735834377
71691114 0870115921
2460751 0296303901
694333 9004708377
611275 3054473601
130335 2374361581
1898 0353810981
505 3237595537
366 3799169941
349 8504482761
250 3632449809
248 6106909971
243 6032686501
105 6956063617
6 7164918451
9165357241
4813927021
4412993041
1017715609
907534489
832304719
334787041
331359241
122841041
100144711
36193361
32317321
31685221
27608407
24159511
19787381
19702153
11964131
11594521
8650651
5094601
2524807
1617757
1548541
1301497
1179263
922601
822793
687961
350281
279109
61933
52081
44269
36697
30097
18061
14341
12601
11131
6763
6301
2917
2521
2381
1103
1051
1009
883
757
709
673
631
541
439
421
353
331
281
241
211
199
181
127
113
73
71
612
43
41
37
31
29
19
17
13
11
7
5
34
2

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

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 = 7 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 = 969 1560339749 6513724766 2342071877 0553042333 9160051009 3536820820 1356007672 7754071669 8989139068 7700073848 8054007439 2526005110 3089619044 7995707499 6537496274 2789905714 9807084650 5442003723 4937735651 6873316995 6998603548 0388567264 9254977266 9753543535 1054651022 9518351452 4483068024 0972117764 0281061331 8912212819 3631761184 5769017709 8153964133 7164431734 9716413486 2732920906 7653536323 5594662584 4254422016 7002237468 0391255043 4359117185 2546706421 5231624419 5273799822 2464159972 7875572363 8330145534 9266570250 3867968215 9017323151 6034394179 3545471818 7141380649 3649188741 3915250961 4193223634 9593675679 4135421228 8630301209 1346966138 0060411947 2517370425 5190065588 3058812505 4247362277 8903583555 0058581496 5425951699 8468632445 2405070004 7672284533 6334507653 7407050052 8287334555 7013009446 1625683230 7404619128 7520615445 8196657270 3278221214 6383742795 9114264971 5369699626 8623573181 0416292205 3202313075 1841507077 1140332960 8470050720 0419944361 4815228734 3353000439 8320612406 9405714972 4595441231 2341926327 3656220542 2716599229 3196149742 5771813487 3405982874 5311263804 5198470022 5851650961 5573592199 5662896142 9727768290 5226051124 8363149089 5626873303 7324372405 3423484370 2899778317 8932323123 5458859258 7795088565 6684183278 1575655997 7517806479.

With those constraints, the unique continued fraction is: {0, 3, 4, 2, 4, 1, 3, 1, 1, 3, 1, 2, 1, 1, 1, 3, 3, 22, 3, 6, 2, 10, 15, 49, 5, 1, 1, 1, 1, 1, 40, 1, 5, 1, 2, 1, 1, 11, 14, 10, 1, 3, 1, 1, 1, 9, 1, 1, 5, 2, 2, 2, 3, 1, 2, 3, 2, 2, 1, 9, 7, 6, 2, 3, 1, 1, 13, 7, 1, 2, 1, 2, 2, 2, 9, 1, 4, 1, 1, 6, 7, 1, 1, 2, 1, 1, 1, 1, 5, 5, 1, 1, 1, 8, 1, 2, 1, 5, 4, 2, 2, 1, 3, 1, 5, 2, 45, 2, 1, 44, 1, 35, 1, 5, 5, 1, 1, 3, 2, 2, 1, 3, 1, 1, 1, 1, 2, 3, 1, 1, 85, 1, 22, 1, 4, 5, 6, 1, 28, 4, 4, 1, 2, 2, 1, 1, 5, 1, 3, 1, 11, 38, 1, 2, 1, 2, 46, 3, 1, 11, 1, 2, 3, 15, 4, 370, 9, 1, 2, 10, 10, 3, 12, 2, 5, 4, 1, 2, 13, 2, 4, 1, 142, 1, 1, 1, 4, 6, 1, 46, 1, 3, 2, 1, 2, 1, 4, 1545, 20, 5, 28, 1, 8, 1, 1, 4, 4, 4, 1, 1, 379, 3, 1, 10, 4, 26, 3, 16, 4, 2, 1, 5, 1, 1, 5, 6, 1, 9, 1, 7, 8, 1, 11, 6, 1, 23, 8, 1, 3, 1, 1, 2, 59, 3, 1, 4, 3, 1, 1, 1, 3, 6, 1, 10, 10, 1, 1, 1, 15, 1, 2, 7, 1, 18, 2, 1, 1, 1, 2, 5, 8, 1, 7, 2, 1, 8, 2, 5, 6, 7, 1, 5, 3, 1, 4, 5, 1, 51, 1, 15, 1, 557, 1, 2, 1, 268, 1, 5, 1, 29, 8, 1, 15, 45, 3, 60, 1, 1, 2, 1, 1, 7, 1, 1, 7, 2, 3, 1, 1, 1, 5, 5, 6, 1, 4, 3, 1, 4, 1, 210, 2, 10, 2, 1, 3, 2, 5, 9, 1, 6, 2, 1, 2, 1, 2, 1, 1, 2, 1, 1, 33, 1, 2, 2, 1, 1, 2, 1, 9, 1, 11, 1, 1, 17, 112, 2, 4, 106, 5, 1, 2, 6, 1, 6, 1, 7, 1, 23, 1, 11, 1, 1, 6, 1, 1, 26, 1, 15, 1, 2, 2, 7, 4, 2, 2, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 3, 2, 4, 13, 1, 3, 2, 2, 1, 2, 1, 5, 1, 1, 6, 1, 1, 4, 1, 1, 225, 1, 6, 5, 1, 2, 3, 61, 1, 3, 3, 1, 17, 1, 85, 5, 1, 1, 1, 123, 6, 1, 23, 19, 2, 1, 1, 4, 1, 2, 1, 13, 2, 1, 2, 7, 1, 11, 3, 1, 198, 1, 5, 2, 1, 1, 1, 13, 23, 1, 14, 5, 6, 12, 2, 2, 2, 3, 2, 1, 23, 1, 181, 1, 14, 2, 15, 2, 23, 1, 1, 1, 10, 1, 1, 1, 1, 1, 1, 1, 1, 3, 4, 4, 1, 16, 1, 1, 2, 1, 2, 1, 28, 8, 2, 1, 3, 1, 1, 3, 1, 5, 1, 3, 1, 39, 14, 7, 3, 5, 2, 7, 1, 6, 2, 49, 109, 8, 24, 1, 39, 1, 17, 5, 1, 2, 1, 2, 4, 9, 5, 1, 1, 1, 8, 33, 1, 8, 3, 1, 7, 1, 1, 1, 4, 1, 1, 1, 1, 1, 2, 1, 1, 3, 3, 13, 1, 3, 5, 2, 2, 1, 2, 1, 14, 1, 6, 1, 1, 2, 10, 1, 5, 1, 6, 1, 50, 2, 2, 1, 1, 4, 8, 1, 32, 1, 5, 1, 37, 3, 22, 1, 1, 2, 4, 7, 1, 2, 32, 1, 4, 2, 1, 1, 6, 5, 2, 1, 8, 2, 6, 20, 2, 2, 6, 13, 14, 2, 1, 1, 6, 1, 76, 1, 1, 2, 8, 3, 1, 2, 1, 1, 9, 1, 4, 2, 1, 1, 7, 1, 10, 13, 1, 1, 1, 6, 1, 9, 1, 1, 1, 11, 1, 5, 4, 5, 4, 1, 5, 3, 1, 1, 3, 130, 2, 1, 1, 3, 5, 19, 1, 2, 14, 1, 2, 2, 3, 1, 2, 1, 32, 2, 2, 8, 1, 2, 6, 2, 1, 1, 1, 3, 169, 1, 1, 27, 1, 15, 1, 308, 1, 9, 15, 1, 5, 87, 4, 2, 6, 4, 2, 3, 1, 2, 1, 30, 1, 15, 2, 1, 8, 1, 1, 2, 1, 3, 1, 1, 6, 2, 3, 2, 1, 1, 1, 3, 3, 1, 1, 1, 1, 2, 6, 1, 1, 37, 1, 14, 2, 2, 2, 8, 1, 4, 3, 8, 1, 2, 1, 6, 1, 15, 94, 1, 1, 2, 15, 1, 2, 1, 22, 2, 1, 1, 1, 7, 2, 1, 1, 1, 1, 1, 10, 2, 2, 11, 3, 2, 7, 38, 5, 11, 1, 2, 2, 1, 1, 1, 3, 1, 2, 1, 12, 1, 2, 2, 25, 1, 5, 2, 1, 1, 6, 2, 1, 4, 2, 1, 6, 1, 1, 5, 1, 6, 2, 1, 33, 4, 2, 1, 2, 1, 3, 3, 1, 2, 4, 1, 2, 1, 2, 12, 1, 1, 5, 1, 2, 3, 2, 16, 1, 1, 2, 3, 2, 1, 2, 1, 2, 2, 5, 7, 21, 26, 1, 4, 1, 2, 16, 1, 1, 6, 1, 10, 2, 2, 1, 1, 9, 1, 2, 3, 3, 5, 5, 4, 15, 1, 2, 8, 5, 8, 1, 4, 1, 10, 37, 11, 2, 10, 1, 1, 4, 1, 1, 6, 2, 1, 1, 14, 1, 24, 4, 1, 1, 5, 4, 1, 19, 34, 3, 2, 4, 5, 1, 1, 1, 2, 1, 1, 42, 1, 422, 3, 5, 2, 3, 1, 1, 1, 10, 1, 1, 1, 4, 12, 1, 3, 7, 1, 2, 1, 2, 3, 2, 4, 3, 2, 21, 23, 22, 4, 1, 3587, 2, 1, 1, 1, 1, 1, 6, 8, 3, 7, 6, 2, 6, 1, 6, 1, 3, 1, 1, 1, 75, 5, 1, 11, 1, 1, 5, 1, 16, 4, 2, 1, 7, 1, 5, 1, 1, 1, 2, 1, 1, 1, 2, 4, 2, 18, 2, 3, 4, 2, 2, 1, 7, 4, 2, 1, 8, 1, 3, 1, 3, 7, 2, 2, 3, 1, 2, 3, 7, 1, 30, 1, 2, 1, 14, 5, 3, 4, 4, 1, 1, 1, 1, 1, 1, 4, 1, 16, 1, 1, 7, 1, 1, 165, 5, 339, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 12, 1, 5, 1, 2, 4, 9, 1, 2, 1, 4, 4, 11, 7, 3, 2, 5, 1, 292, 1, 4, 1, 1, 3, 2, 13, 1, 5, 2, 1, 34, 2, 1, 5, 14, 8, 4, 1, 2, 11, 1, 1, 12, 1, 1, 1, 1, 2, 5, 2, 6, 1, 8, 1, 9, 1, 2, 2, 7, 1, 1, 5, 18, 3, 2, 1, 4, 8, 3, 1, 1, 1, 3, 3, 3, 1, 3, 2, 8, 1, 3, 8, 2, 2, 1, 3, 1, 1, 2, 4, 1, 32, 1, 9, 1, 8, 1, 1, 1, 3, 7, 23, 2, 1, 3, 1, 6, 18, 9, 1, 33, 3, 1, 2, 1, 2, 9, 5, 1, 1, 1, 5, 1, 2, 1, 2, 1, 10, 1, 1, 10, 3, 1, 3, 1, 4, 1, 1, 1, 2, 4, 1, 2, 1, 12, 6, 1, 1, 1, 1, 1, 2, 4, 3, 3, 2, 1, 20, 24, 2, 2, 2, 2, 7, 7, 1, 210, 2, 2, 2, 2, 1, 15, 1, 9, 11, 1, 10, 2, 12, 1, 4, 1, 3, 11, 5, 7, 1, 2, 4, 1, 1, 5, 3, 1, 2, 1, 13, 2, 1, 3, 6, 3, 14, 1, 4, 3, 1, 1, 11, 1, 4, 1, 1, 15, 1, 1, 1, 4, 1, 2, 7, 4, 1, 4, 2, 9, 18, 38, 3, 1, 1, 2, 2, 1, 16, 1, 2473, 1, 3, 1, 2, 2, 5, 2, 1, 1, 4, 8, 2, 1, 14, 1, 2, 8, 13, 1, 2, 1, 4, 1, 1, 2, 8, 2, 17, 1, 1, 1, 4, 1, 2, 1, 31, 1, 47, 8, 7, 2, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 3, 2, 2, 1, 1, 2, 2, 39, 1, 27, 15, 6, 1, 1, 2, 3, 1, 1, 2, 1, 3, 2, 10, 10, 2, 1, 1, 2, 1, 21, 8, 2, 1, 1, 1, 2, 2, 2, 2, 1, 4, 1001, 10, 1, 11, 57, 2, 5, 1, 40, 1, 1, 6, 1, 1, 1, 2, 5, 4, 4, 11, 3, 1, 11, 1, 3, 3, 1, 2, 8, 1, 10, 1, 4, 2, 44, 2, 1, 6, 35, 1, 9, 1, 2, 4, 1, 1, 1, 2, 2, 14, 2, 1, 11, 1, 2, 1, 3, 55, 1, 1, 3, 21, 3, 2, 2, 1, 4, 38, 6, 2, 1, 15, 1, 3, 183, 1, 2, 3, 1, 17, 1, 1, 7, 6, 1, 48, 5, 3, 1, 11, 2, 6, 1, 3, 1, 7, 41, 20, 22, 9, 17, 16, 2, 7, 17, 1, 2, 19, 2, 1, 3, 1, 62, 1, 2, 4, 1, 1, 1, 2, 6, 1, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 2, 1, 2, 4, 5, 1, 9, 1, 29, 21, 1, 1, 45, 1, 12, 1, 1, 3, 4, 25, 1, 2, 1, 1, 2, 4, 1, 11, 1, 3, 1, 2, 1, 1, 6, 1, 9, 1, 7, 2, 4, 10, 1, 1, 1, 1, 3, 3, 1, 45, 1, 10, 2, 1, 4, 1, 1, 2, 1, 4, 4, 3, 1, 2, 3, 1, 1, 6, 1, 16, 1, 1, 4, 2, 8, 8, 1, 1, 8, 1, 1, 1, 2, 1, 1, 2, 4, 3, 1, 1, 1, 4, 1, 1, 1, 81, 2, 1, 1, 2, 1, 2, 6, 1, 1, 1, 4, 10, 1, 2, 1, 2, 1, 1, 23, 2, 1, 1, 1, 1, 6, 1, 2, 2, 3, 1, 2, 2, 1, 5, 3, 5, 1, 1, 1, 1, 2, 3, 1, 2, 1, 28, 1, 14, 1, 5, 1, 9, 8, 2, 1, 4, 6, 1, 1, 1, 3, 2, 3, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 6, 1, 3, 2, 22, 64, 1, 1, 5, 5, 1, 118, 24, 1, 3, 1, 22, 1, 1, 6, 1, 2, 17, 3, 7, 1, 6, 1, 1, 1, 1, 1, 9, 1, 1, 11, 1, 3, 1, 2, 7, 770, 1, 6, 3, 4, 8, 6, 2, 1, 11, 32, 1, 1, 123, 3, 1, 5, 1, 1, 5, 1, 1, 36, 1, 638, 2, 1, 36, 3, 9, 1, 5, 17, 1, 2, 3, 28, 6, 1, 1, 5, 1, 1, 1, 3, 1, 11, 3, 3, 1, 1, 2, 2, 5, 1, 6, 8, 2, 1, 2, 1, 1, 2, 1, 33, 1, 4, 1, 1, 3, 5, 20, 2, 2, 2, 4, 3, 2, 3, 2, 1, 3, 2, 1, 3, 2, 3, 1, 2, 66, 1, 7, 10, 3, 4, 11, 17, 7, 1, 7, 1, 2, 1, 8, 1, 44, 1, 36, 2, 1, 2, 2, 1, 3, 1, 2, 7, 4, 1, 1, 2, 4, 1, 10, 2, 1, 1, 1, 378, 2, 6, 15, 1, 2, 27, 1, 2, 1, 1, 2, 9, 3, 1, 9, 25, 1, 1, 2, 11, 4, 1, 3, 5, 3, 6, 1, 3, 2, 3, 1, 1, 8, 4, 4, 1, 4, 1, 2, 1, 2, 1, 5, 2, 3, 2, 12, 1, 51, 1, 2, 1, 1949, 5, 1, 1, 2, 10, 2, 15, 7, 3, 1, 2, 1, 1, 1, 5, 1, 3, 2, 2, 1, 1, 15, 481, 1, 2, 25, 1, 2, 27, 1, 3, 2, 1, 7, 2, 9, 4, 1, 1, 4, 1, 3, 1, 1, 2, 1, 5, 4, 7, 1, 1, 28, 1, 3, 1, 3, 1, 1, 1, 1, 5, 2, 1, 8, 2, 1, 1, 1, 1, 4, 2, 7, 3, 1, 2, 14, 1, 28, 1, 5, 1, 2, 4, 3, 2, 2, 1, 1, 1, 1, 1, 6, 218, 4, 1, 3, 1, 1, 2, 1, 3, 1, 19, 1, 12, 1, 2, 1, 3, 6, 7, 5, 2, 8, 3, 2, 2, 3, 7, 1, 10, 1, 3, 2, 4, 3, 1, 1, 8, 3, 1, 5, 15, 1, 1, 1, 3, 13, 1, 2, 6, 27, 7, 1, 1, 1, 4, 4, 2, 1, 2, 8, 1, 2, 1, 5, 1, 1, 14, 3, 25, 1, 2, 3, 5, 2, 3, 8, 14, 2, 1, 3, 8, 2, 1, 11, 5, 3, 3, 21, 1, 6, 79, 1, 1, 1, 1, 5, 1, 2, 14, 1, 2, 1, 4, 6, 5, 6, 1, 4, 1, 1, 7, 1, 10, 11, 1, 4, 1, 1, 4, 11, 2, 7, 3, 1, 1, 3, 22, 1, 89, 3, 1, 1, 2, 3, 3, 9, 8, 1, 40, 5, 2, 93, 1, 13, 1, 1, 2, 1, 21, 8, 1, 213, 1, 2, 1, 2, 1, 1, 2, 11, 3, 10, 1, 2, 3, 1, 4, 1, 1, 1, 9, 2, 2, 2, 2, 1, 5, 1, 2, 1, 1, 1, 1, 11, 1, 16, 1, 3, 1, 3, 9, 3, 3, 1, 1, 6, 3, 13, 5, 1, 2, 1, 1, 7, 24, 1, 1, 2, 10, 1, 1, 1, 1, 1, 14, 3, 2, 2, 18, 1, 5, 18, 1, 4, 1, 1, 7, 6, 1, 5, 5, 9, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 11, 4, 1, 8, 2, 22, 1, 2, 2, 51, 1, 2, 1, 1, 1, 3, 3, 3, 2}, giving these values for u and v:

… as taking one more term in the continued fraction would give a value of v of 1478 3407908951 5906263107 1047181103 5759350740 4984322617 0084524972 1533265782 8833789491 0605106762 6031394336 8433596813 1275598979 7006068583 3936134516 9363763323 5170699851 2861420715 4774609196 7290070029 3849437398 5205027050 5070118381 7710272382 3577195191 3674311011 3513644643 0957686521 8948051102 6290798031 2876072287 8322333226 8645141314 6263143673 0355524350 7564411528 9002431117 0034617165 0740338988 2740895966 3617175436 2552256998 3096462411 1380752633 7228484723 2070185329 7961115714 6982400869 1144417269 6799131503 2352458410 4153697604 2580956868 8472430498 1688152642 7101782548 4484567876 3616885022 3247387158 3631797090 0802629065 3625967977 4011655706 0519474780 6182708508 5056597095 8374042799 2369274696 1922441115 2681720348 4271593848 6615022397 8538818495 7183470056 9884915742 3109567529 1318495796 5033258966 6779115292 1325060986 9397654389 6612833991 1275018319 0006745351 8784234458 6134861946 7097580815 7283208477 5831974945 3410141178 4277407416 9723966093 0991285933 2556455257 9989857774 9085463038 3049206675 4759856705 4040383862 5312316536 7189495490 6546805069 2544154521 5463055306 3385053376 8365957722 9470325411 1512427683 2947085519 4426346445 5277267538 6696685429 1849361132 7237148271 4901626515 4152248139 8208459311 9408117460 5856124363 1511215480 8966838031, which is too large.

We also need to calculate d = floor(c4·v/F + 0.5) = 1869007 1858570715 4919268028 5074372554 5767859627 8716430943 0949145468 1157870663 0560595395 9909542930 3110658580 2598601171 6738042189 5993199919 8649167893 5982245271 4936247568 9485632218 7204238698 8669063959 7185902202 7476164921 3361266778 0661429826 7536199867 7858048463 7359040971 9317699132 1034007070 7771999057 2507802701 4625077977 9450395143 4425078377 7945874942 1525000596 3460727490 3230145204 4893914490 8303690769 9634171447 5801951449 1684114519 3736500559 5596385222 0441996273 2530338011 0737886541 1848415223 8842059334 5555614928 7816986365 8518190395 7066002175 2424068860 6236382690 2697486497 2903058075 2561762311 3891028408 1492853207 4515100759 1054173725 0485037876 3521550874 0183304289 0291815606 7728502787 1342190642 4067390117 0946689574 5571914505 8250132903 2648909730 7652761909 8146571044 3025410332 0176351631 1681164743 4586045422 2724690875 7354697264 3562360728 2884593048 7570716769 9724885654 7401539344 1194374287 8801678716 3648554518 1684488074 5737489114 2139538480 6549542893 8880991664 6300011985 2729561105 0879682687 6840910374 3410642783 3994786096 1066989857 5860112628 1466438324 1910420769 8063045667 2656385163 3973232084 2551160369 4077918545 1116085637 6370259924 7329822529 0923085893 1838101323 4032675024 3134666579 7343386857 7342748764 4820816742 1283567163 4885967056 4229648022 5517132164 9430169215 1890959550 3584104920 1783805722 1825742201 2292859992 6248317319 5486393607 4203338124 3581523061 8564680097 2704270958 6637355206 4715962292 4115694500 6212274023 3380180779 9498976893 8012224571 5732719215 1828376213 6880424531 8006623021 2740021750 7920287536 8735818206 4829875272 8982569004 3976547052 7501435156 5583461178 1906963650 0304588046 1689123199 4720147074 6749391321 7976445844 5596246316 6950979782 3630988009 9364077937 3993961106 6772900046 5081882141 3366785110 3015669228 6945644499 9177375152 8754574011 0520070036 0845309088 2976674461 9480074308 2472633379 8109663160 3971666354 4448439522 2887093700 7484819049 8131371955 5786369536 0982517432 9838342403 7982345651 2993033363 8825835066 6931531743 7090373801

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.

P has a single real root at:

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