Primality Certificate for (8917^2689-1)/8916

Andy Steward10,619 digits14 July 2001
Originally by A.A.D.Steward 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.377017% factorization of N-1:

From Factorisation
891737 · 241
Φ22 · 7 · 7 · 7 · 13
Φ33 · 181 · 146449
Φ42 · 5 · 61 · 130349
Φ62593 · 30661
Φ743 · 11692108883670269514769
Φ82 · 17 · 39113 · 4754173441
Φ124347613 · 1454200141
Φ147 · 49184059 · 1459961796933361
Φ162 · 769 · 766321 · 45064593617 · 752571229937
Φ215209 · 6133 · p40
Φ24673 · 54042841 · 77317417 · 14214106351441
Φ2829 · 230980251914370589 · 37727006394571316213590308257
Φ322 · 197497998433 · 11184384929246673341441 · p30
Φ42127 · 7477 · 11936511511 · p32
Φ48752833 · 1019571597086832093924496033 · p31
Φ5652678921 · 7575088781452421041 · 22546711132820654615210033 · p43
Φ642 · 193 · 6337 · c121
Φ8421974653 · 90405421 · p80
Φ9697 · 577 · c122
Φ112113 · 2129 · 11375747438016760177 · 150965744557828660769 · p145
Φ1282 · 182668289 · c245
Φ168337 · 384889 · c182
Φ1929601 · c249
Φ224449 · 44129 · 165089 · 10209473 · c360
Φ336c380
Φ384143617 · 7500991138499709313 · c482
Φ44812097 · 1364609 · 34774657 · 93318409409 · c730
Φ6727705472411617 · 49188455034313729 · p729
Φ8964151421569 · p1508
Φ13443784433451457 · c1505
Φ26882689 · 269329537 · c3022

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 :

18442065 1992924531 0681760493 7354263086 8476385486 2937164039 4252721634 4072298144 2120355779 9067308790 4513451461 5190752749 2179588770 6451131881 3551520502 5138623076 1241915174 8098577180 5802709154 8776287760 9699813608 5829761904 7627153381 4776195927 7393344870 3298419298 2452496280 4819083374 1900433714 9242829919 2868200321 7056389896 9591605921 0100652115 4306254414 2019666610 4089962170 1337447013 4746513178 7545409626 9650554576 8403057440 7459713392 9004743861 7063375843 0656995147 6842920702 4523416171 0825638353 0183903042 9606800533 1660534279 6202639758 8131381224 2961366030 7304721135 0634260109 5642315417 2519822341 5296864722 7867643799 1092628094 6701536439 8667099659 2460304158 2472561062 1034913144 6624410088 5292892790 4533750367 4872422534 5044895645 0285347959 2971863718 7834021003 2081270530 9000974244 6827525191 9788666643 4386851868 4447579671 5521093199 3801586678 1718843822 7226311582 7554704542 2430913337 7909782018 4760974765 6798259082 7956240991 1528217566 2640936595 4574152220 3069644245 7945559282 7054658784 3491136278 7888538569 4325750928 5753373938 2294652252 5379300848 7508177132 7134373334 5983768851 7118108547 2970160259 1897441530 7446329701 0089047938 2491413296 7852023548 3838524429 1765086425 8517587337 3140361218 2073520865 9989651682 2354061239 2963044324 3570211563 9793577373 8042277024 6264785720 1195660586 8493284992 9207580680 8131561582 3653763985 3580642278 3815597627 2054507224 3060370920 2659474975 9729053877 4262139238 5699706948 0789399345 0756821131 6358682021 6334705570 7142043454 6945521004 8818083069 3676937868 7941846557 7714441919 3080459116 9973799617 7920458369
730030060 4216477631 6539258225 3158305684 3025347872 2605875674 2909593142 6227870580 2864808298 3624093465 1131039102 0701624795 5873295464 4628148485 8526800143 5728308818 2856226719 2864941839 1245245360 1572058845 7201647596 1803833242 5055973861 8483175313 9293447669 0589531385 1209742796 6777559897 9232299599 4436215159 6385074699 2158845389 7857650055 5435087731 3527179350 3297474875 7638605798 6890199616 3993054885 0529276072 8512291132 6911611790 0473975538 1481218638 8793089525 4217623628 7198638332 0559941819 8966550945 3721902680 9544634314 2719477781 4401200670 8927904832 2633710408 3412998997 5441558955 4220866237 7643887132 5203971051 7719834245 1131141961 4815379498 1389179650 9985936745 8384754304 4978923337 3415262028 9956497403 2521772665 7173903407 3864657792 5273363617
98716 1795522954 5410670671 3635534949 5353562457 1552088403 1267264766 8025010152 2442132858 7303927501 4137703456 1456603038 8721738044 3889792396 3378010321
3214651559 4244158386 6553183058 1333965468 1953387032 9050427766 9913274229 8558401497
709 8107960307 2534470881 3982191890 9982200297
7909502556 3143476192 2284979758 9968611757
22 2979726421 7643637647 9286220061
2 0815358249 9127792380 4257895169
3616558266 3297635870 1222608897
377270063 9457131621 3590308257
10195715 9708683209 3924496033
225467 1113282065 4615210033
116 9210888367 0269514769
111 8438492924 6673341441
1 5096574455 7828660769
1137574743 8016760177
757508878 1452421041
750099113 8499709313
23098025 1914370589
4918845 5034313729
145996 1796933361
1421 4106351441
770 5472411617
378 4433451457
75 2571229937
19 7497998433
9 3318409409
4 5064593617
1 1936511511
4754173441
4151421569
1454200141
269329537
182668289
90405421
77317417
54042841
52678921
49184059
34774657
21974653
10209473
4347613
1364609
766321
752833
384889
165089
146449
143617
130349
44129
39113
30661
12097
9601
7477
6337
6133
5209
2689
2593
2129
769
673
577
449
337
241
193
181
127
113
97
61
43
37
29
17
13
74
5
3
27

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

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 = 41 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 = 76 4680914760 6413350407 9419466670 0153267689 2933474528 3772308978 8189183890 6616860715 9208703833 2524105271 1146090444 5137687527 5346784942 5106852789 2358207988 5583078919 2336434495 8521978719 4223736271 6246043847 2045948898 0061849234 4106709388 7575448803 1434777930 4988729122 3136318304 8306954112 7799433494 5486251060 6271138180 9730072682 3110396204 5287010139 7110856408 1787041414 4362206754 5428758003 1203817461 6570024490 7474139556 3279816533 3363619052 7971912095 7141607345 7085625367 2858139237 9880740826 9849968804 1368698006 0007763653 9808502308 7122654363 9766242466 0023949193 9550400841 5740108327 5883086312 5797241773 9273044091 9139596162 3624527934 9401953919 9365447928 2567640361 0083250050 1124777439 5095581705 7570186709 1665279161 9203980456 2852042159 5409620750 1014216843 3770648857 2374273747 8496974175 0688527958 0768876972 7820154391 2815969650 3394319422 2056985206 1502546560 9493506617 2405980626 0407422289 9947909820 9199069263 7885407859 0847200188 7364367433 0479573449 9486690286 1209320984 8110717534 8293726596 9248373875 5760996732 4943057989 2802527231 3840033467 8785322613 5230417661 3763231509 8336811823 5766698684 4534821435 9441136841 8457057632 4409583042 8725079677 5828646074 6535001795.

With those constraints, the unique continued fraction is: {0, 1, 1, 6, 4, 2, 7, 8, 1, 1, 1, 7, 5, 2, 2, 284, 1, 2, 2, 24, 3, 17, 36, 4, 3, 6, 20, 1, 5, 1, 4, 7, 2, 1, 9, 1, 18, 4, 2, 47, 1, 1, 7, 1, 3, 1, 2, 1, 7, 1, 1, 4, 3, 1, 2, 6, 60, 7, 4, 17, 2, 1, 1, 1, 4, 4, 1, 2, 1, 8, 1, 1, 4, 1, 1, 1, 1, 45, 1, 34, 2, 1, 1, 1, 1, 12, 1, 8, 3, 1, 24, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 3, 44, 2, 1, 1, 1, 1, 408, 35, 1, 4, 2, 2, 8, 24, 1, 21, 4, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 27, 1, 1, 20, 6, 2, 3, 1, 1, 1, 4, 2, 1, 10, 26, 9, 1, 2, 16, 1, 4, 6, 2, 1, 3, 5, 30, 3, 11, 1, 46, 1, 27, 4, 3, 1, 1, 1, 1, 2, 2, 1, 3, 1, 38, 2, 2, 1, 1, 37, 3, 1, 1, 14, 1, 3, 1, 3, 1, 204, 21, 3, 1, 4, 2, 1, 3, 20, 1, 2, 3, 6, 2, 1, 2, 15, 1, 10, 1, 3, 1, 1, 29, 4, 2, 3, 4, 1, 4, 4, 3, 1, 3, 16, 1, 7, 6, 29, 1, 3, 1, 1, 9, 3, 5, 3, 1, 4, 2, 1, 18, 3, 13, 1, 4, 2, 1, 14, 1, 1, 13, 1, 4, 1, 1, 2, 3, 1, 1, 1, 1, 3, 3, 1, 6, 20, 1, 1, 1, 7, 1, 2, 9, 2, 2, 1, 19, 5, 7, 1, 16, 1, 2, 4, 2, 2, 1, 7, 1, 8, 15, 1, 1, 4, 31, 5, 122, 1, 1, 1, 3, 1, 1, 1, 1, 4, 1, 108, 1, 23, 2, 1, 2, 1, 2, 5, 4, 8, 1, 3, 1, 1, 1, 1, 1, 1, 1, 4, 1, 46, 4, 3, 1, 4, 1, 3, 2, 65, 1, 5, 1, 15, 2, 2, 5, 1, 1, 1, 26, 2, 1, 108, 1, 12, 2, 1, 1, 20, 1, 1, 2, 2, 1, 1, 2, 1, 8, 1, 1, 3, 1, 10, 1, 2, 1, 2, 18, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 15, 1, 3, 1, 6, 2, 3, 26, 1, 3, 2, 1, 2, 2, 1, 11, 3, 1, 1, 2, 1, 1, 1, 1, 3, 1, 21, 1, 2, 7, 4, 3, 1, 1, 20, 2, 1, 1, 2, 1, 1, 5, 4, 2, 3, 1, 10, 3, 1, 16, 1, 2, 1, 20, 3, 1, 4, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 12, 6, 2, 1, 275, 1, 2, 1, 5, 1, 1, 20, 18, 2, 7, 1, 9, 1, 7, 11, 1, 2, 1, 1, 3, 2, 626, 15, 9, 7, 17, 4, 1, 1, 7, 6, 1, 47, 3, 3, 1, 1, 1, 2, 1, 5, 3, 4, 1, 66, 1, 1, 5, 3, 3, 9, 1, 5, 1, 2, 1, 1, 3, 1, 5, 7, 1, 1, 27, 1, 3, 1, 24, 3, 1, 2, 1, 7, 2, 5, 10, 1, 1, 1, 5, 19, 1, 2, 3, 1, 7, 1, 3, 1, 3, 1, 1, 7, 1, 1, 1, 2, 1, 139, 4, 1, 1, 1, 3, 6, 1, 4, 1, 11, 1, 1, 1, 51, 2, 1, 4, 1, 1, 1, 2, 22, 1, 9, 1, 209, 1, 1, 4, 1, 5, 1, 61, 1, 1, 1, 6, 117, 1, 1, 4, 1, 1, 1, 3, 2, 52, 1, 2, 1, 10, 1, 1, 2, 3, 1, 4, 1, 53, 1, 12, 8, 2, 1, 19, 3, 10, 8, 2, 7, 10, 1, 1, 10, 1, 1, 1, 2, 2, 1, 1, 3, 3, 1, 1, 1, 16, 1, 65, 1, 2, 1, 1, 3, 1, 6, 2, 1, 1, 3, 4, 3, 4, 10, 3, 1, 1, 1, 3, 1, 2, 9, 1, 4, 1, 8, 1, 2, 3, 6, 4, 23, 1, 20, 9, 36, 24, 4, 2, 1, 41, 2, 38, 1, 1, 1, 3, 1, 1, 3, 1, 1, 28, 1, 11, 28, 2, 4, 4, 2, 1, 68, 23, 1, 1, 1, 3, 1, 3, 2, 1, 4, 4, 1, 1, 15, 1, 1, 15, 16, 1, 1, 2, 4, 1, 5, 1, 1, 8, 4, 1, 7, 1, 2, 1, 1, 3, 3, 3, 4, 1, 1, 5, 2, 1, 3, 2, 2, 2, 1, 4, 1, 114, 1, 67, 2, 1, 2, 1, 5, 2, 2, 1, 2, 8, 1, 1, 3, 5, 1, 6, 1, 10, 2, 1, 1, 3, 1, 12, 1, 1, 1, 1, 5, 1, 17, 2, 4, 6, 5, 1, 1, 1, 4, 1, 3, 22, 3, 2, 2, 2, 35, 1, 1, 1, 1, 6, 1, 5, 8, 5, 1, 2, 5, 26, 2, 3, 5, 3, 22, 21, 1, 13, 1, 1, 3, 1, 12, 1, 1, 1, 1, 14, 56, 3, 3, 8, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 3, 1, 10, 1, 4, 1, 2, 1, 1, 4, 2, 1, 1, 4, 4, 1, 1, 7, 8, 1, 2, 6, 1, 1, 1, 7, 2, 8, 1, 24, 3, 1, 30, 1, 2, 1, 2, 7, 5, 3, 2, 6, 1, 1, 2, 1, 1, 2, 1, 1, 5, 3, 1, 9, 1, 1, 7, 1, 5, 2, 1, 2, 3, 3, 2, 1, 1, 54, 1, 7, 1, 1, 2, 23, 6, 1, 11, 1, 2, 1, 23, 3, 1, 6, 2, 1, 1, 1, 28, 4, 2, 2, 4, 5, 3, 1, 23, 1, 3, 1, 8, 6, 7, 1, 4, 1, 3, 16, 9, 3, 1, 6, 1, 1, 3, 5, 6, 1, 4, 3, 5, 1, 6, 1, 1, 1, 1, 4, 8, 1, 3, 1, 1, 4, 1, 5, 1, 9, 1, 1, 13, 1, 1, 13, 7, 1, 4, 4, 22, 1, 3, 1, 1, 1, 10, 5, 2, 2, 50, 2, 2, 5, 2, 1, 3, 2, 1, 21, 1, 1, 4, 3, 4, 1, 5, 1, 5, 3, 1, 3, 3, 8, 2, 1, 1, 1, 1, 1, 2, 1, 1, 5, 393, 54, 1, 5, 1, 2, 12, 6, 1, 2, 53, 4, 6, 2, 2, 2, 2, 1, 1, 2, 4, 1, 2, 11, 1, 4, 16, 14, 1, 1, 1, 1, 1, 1, 28, 2, 5, 1, 1, 1, 2, 1, 15, 1, 59, 1, 7, 1, 7, 5, 1, 2, 2, 1, 28, 1, 2, 1, 4, 2, 17, 6, 1, 2, 1, 12, 1, 1, 309, 1, 5, 1, 3, 2, 1, 8, 45, 1, 3, 1, 1, 1, 1, 11, 17, 8, 11, 3, 1, 1, 3, 5, 3, 1, 2, 4, 1, 1, 1, 31, 1, 3, 1, 2, 1, 1, 1, 1, 1, 2, 12, 2, 1, 2, 21, 1, 48, 1, 1, 3, 13, 1, 367, 5, 103, 1, 94, 2, 5, 1, 2, 34, 6, 2, 1, 5, 2, 1, 30, 2, 34, 78, 4, 55, 2, 9, 1, 1, 1, 2, 1, 1, 15, 29, 2, 1, 13, 2, 7, 2, 1, 1, 3, 6, 1, 2, 2, 1, 8, 19, 2, 1, 3, 6, 5, 2, 3, 1, 5, 1, 1, 2, 1, 1, 4, 69, 1, 1, 3, 4, 2, 36, 2, 3, 3, 2, 2, 1, 5, 1, 1, 3, 1, 2, 2, 2, 1, 7, 1, 2, 1, 1, 2, 1, 14, 1, 1, 1, 7, 2, 7, 1, 1, 2, 4, 16, 1, 3, 1, 2, 1, 2, 11, 5, 2, 2, 3, 1, 1, 12, 1, 2, 1, 13, 1, 25, 1, 2, 23, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 9, 1, 2, 1, 1, 3, 1, 2, 1, 3, 1, 1, 3, 5, 6, 2, 39, 2, 3, 2, 140, 1, 76, 19, 5, 1, 10, 5, 3, 2, 1, 1, 1, 3, 1, 3, 2, 2, 2, 2, 1, 2, 1, 2, 15, 6, 2, 17, 7, 4, 64, 15, 28, 20, 3, 3, 2, 2, 36, 10, 2, 1, 1, 6, 2, 1, 1, 4, 3, 2, 30, 2, 1, 3, 1, 3, 2, 2, 10, 2, 3, 2, 2, 3, 1, 2, 1, 1, 2, 1, 556, 1, 16, 1, 1, 1, 1, 1, 1, 3, 3, 3, 2, 1, 3, 1, 8, 2, 1, 1, 3, 10, 1, 2, 20, 1, 14, 2, 2, 2, 1, 6, 22, 3, 1, 1, 1, 1, 2, 14, 2, 6, 1, 1, 18, 1, 2, 1, 2, 2, 1, 3, 5, 3, 1, 9, 1, 13, 2, 3, 2, 3, 3, 23, 15, 1, 133, 174, 1, 1, 1, 4, 3, 1, 1, 1, 2, 3, 4, 3, 1, 4, 2, 10, 10, 2, 1, 1, 1, 1, 5, 1, 1, 13, 145, 1, 1, 14, 158, 3, 1, 1, 1, 3, 7, 2, 1, 1, 21, 1, 3, 2, 1, 2, 2, 2, 3, 2, 1, 7, 11, 1, 34, 1, 5, 1, 7, 1, 4, 11, 7, 1, 75, 8, 1, 1, 1, 6, 1, 3, 1, 1, 1, 2, 2, 6, 10, 1, 2, 5, 2, 1, 45, 2, 2, 1, 1, 1, 1, 4, 1, 5, 1, 2, 4, 9, 1, 55, 4, 1, 4, 2, 2, 1, 38, 2, 84, 1, 3, 5, 2, 2, 2, 7, 2, 1, 2, 5, 4, 9, 4, 1, 2, 1, 1, 16, 1, 9, 1, 7, 2, 2, 1, 2, 10, 72, 3, 1, 2, 1, 1, 1, 1, 6, 1, 1, 1, 67, 1, 3, 2, 1, 1, 11, 11, 5, 61, 2, 1, 1, 2, 1, 5, 1, 2, 11, 19, 1, 31, 3, 6, 1, 2, 1, 5, 1, 2, 30, 1, 4, 3, 9, 1, 4, 1, 23, 1, 5, 2, 1, 1, 3, 1, 1, 5, 6, 1, 1, 7, 3, 1, 2, 6, 3, 6, 1, 5, 8, 1, 3, 2, 1, 11, 6, 2, 1, 33, 1, 8, 1, 1, 1, 3, 2, 1, 1, 20, 3, 1, 1, 1, 3, 2, 1, 41, 2, 16, 1, 3, 1, 4, 1, 2, 3, 1, 3, 1, 20, 3, 1, 1, 1, 2, 1, 1, 2, 3, 5, 1, 1, 1, 58, 2, 2, 76, 1, 1, 5, 1, 3, 3, 1, 2, 14, 3, 6, 1, 3, 1, 2, 2, 3, 2, 2, 2, 2, 1, 2, 11, 1, 3, 1, 6, 4, 1, 1, 2, 3, 1, 5, 1, 1, 4, 1, 1, 14, 5, 1, 20, 1, 10, 3, 27, 1, 1, 1, 1, 34, 1, 1, 2, 1, 1, 5, 12, 1, 7, 2, 7, 1, 13, 1, 4, 1, 2, 1, 2, 1, 5, 30, 1, 19, 1, 1, 3, 1, 2, 21, 1, 20, 12, 88, 11, 1, 27, 1, 4, 1, 1, 1, 1, 2, 27, 1, 4, 3, 1, 8, 1, 13, 3, 367, 1, 1, 4, 6, 3, 9, 40, 1, 3, 1, 4, 1, 6, 2, 1, 3, 1, 3, 1, 4, 1, 2, 5, 12, 5, 1, 1, 2, 1, 1, 1, 1, 1, 3, 56, 1, 7, 13, 3, 22, 1, 1, 1, 4, 6, 10, 36, 2, 1, 30, 1, 2, 1, 2, 2, 2, 4, 1, 2, 7, 1, 1, 6, 2, 1, 3, 1, 3, 2, 1, 1, 32, 1, 2, 2, 1, 5, 566, 1, 1, 17, 3, 1, 5, 1, 2, 1, 8, 7, 4, 4, 1, 1, 1, 63, 1, 1, 8, 1, 54, 3, 1, 2, 2, 5, 1, 1, 1, 4, 1, 12, 48, 1, 17, 2, 1, 4, 2, 3, 1, 1, 2, 2, 4, 1, 1, 2, 2, 1, 18, 8, 1, 10, 1, 1, 3, 5, 8, 1, 1, 1, 3, 6, 1, 1, 1, 5, 1, 1, 2, 1, 2, 9, 1, 6, 2, 1, 1, 3, 3, 1, 2, 1, 4, 1, 1, 1, 5, 2, 2, 2, 4, 8, 11, 1, 1, 3, 1, 8, 1, 2, 2, 48, 1, 4, 1, 1, 1, 1, 16, 1, 2, 2, 1, 1193, 1, 1, 7, 1, 4, 4, 1, 1, 2, 2, 1, 2, 16, 2, 1, 4, 669, 3, 1, 5, 2, 5, 1, 1, 2, 3, 2, 2, 2, 1, 2, 3, 1, 1, 14, 1, 4, 2, 20, 1, 8, 3, 1, 1, 1, 1, 1, 1, 4, 8, 18, 22, 1, 1, 3, 1, 2, 7, 5, 1, 1, 1, 1, 3, 6, 46, 1, 3, 2, 1, 1, 2, 1, 32, 1, 5, 1, 25, 16, 1, 2, 1, 1, 1, 1, 1, 14, 7, 2, 1, 1, 12, 2, 10, 1, 9, 1, 3, 7, 28, 1, 1, 3, 1, 1, 1, 1, 2, 124, 7, 1, 3, 148, 1, 2, 1, 2, 2, 1, 2, 5, 9, 1, 44, 2, 3, 28, 2, 3, 1, 1, 1, 5, 2, 9, 1, 2, 1, 7, 2, 5, 2, 4, 8, 2, 3, 8, 5, 1, 3, 1, 4, 1, 1, 25, 1, 48, 1, 6, 27, 1, 2, 431, 8, 18, 9, 2, 2, 2, 11, 15, 3, 26, 1, 1}, giving these values for u and v:

… as taking one more term in the continued fraction would give a value of v of 129 1266460734 8903046324 0013819120 4110292461 2540500375 4217972515 9830982712 5731863336 7873661559 1432038014 1987530865 9846739022 2001370804 7986505826 5386585910 8119074032 1270143749 8197576776 8295935083 8015450805 0878089353 1188380406 2948759490 2412168231 8683607978 2186748031 4594267571 5266549853 8459028055 9007577348 2735595014 5772039332 8964546680 8015291228 5255355615 6171571014 0985167866 0723151754 6780442311 6692120702 8857741745 1967964971 9317519122 1971722641 5891662437 3193281032 9359779183 8450830095 3447504982 6614196849 9440585137 5347780888 1813994602 1240997693 9976334677 5776169873 2207477566 0244434680 4076797098 4820665154 0053775019 0985568135 1808444884 0111442909 8438448112 4271993997 2368362423 3557912912 8186148402 4832692858 9003701653 6435276335 3356228165 2071579940 5348570960 5922864949 0792440505 8879473191 7878506875 4482843774 2405558239 3027709446 3781804503 3139342278 6267076405 6173974716 7369069140 4422681307 9698260744 1432605458 9962894895 4883520956 9724752536 0068900274 0435644583 2190981823 6935208375 8923781437 6401781739 2568782072 8135861952 7312572508 1504858292 6814070367 1558085866 8919486701 8128966128 0636545849 3011863146 8541586409 5173355054 4848924245 0859260846 2958114675, which is too large.

We also need to calculate d = floor(c4·v/F + 0.5) = 266 7595554073 7510049497 9755821450 3529597781 5155264291 3827321750 5873912258 3040915052 4979991730 2612659613 3868965536 1673185605 8310428104 5007897057 6228750898 4624536900 9047744525 2932631019 0861092721 7166173102 3719889015 5902283576 8406470266 5088793459 8626670169 0772115516 6421064420 8427781913 6441378122 5577058753 6864890517 8841825890 8867896773 3828759911 8991850202 6467191675 7098641538 1355522876 3663377996 4620377059 5475055997 5864206396 5502166141 7839122856 9502249485 0270736815 2837949496 1207152068 4007080749 8321881968 0045607156 1071234560 1196211338 4065448625 0622988829 8994776312 1026422805 6835372744 1079306809 0195411591 3619011973 6947959883 0294488125 0097682620 6969458342 1618379581 5578757311 4778628754 0101400277 0233444001 3075427979 7786277517 5185541722 5562006518 0296935393 0722693835 5887115479 0500002625 0475883620 6640989506 8572972857 0949759293 9476403000 8980953814 6032095420 3499839040 2372419518 0976550876 3718565507 6761682802 6882167820 0765496846 4608212527 3718159140 9947353371 9647906350 2890973156 7789035087 3286263288 4065618192 0480206873 8340045996 2848561076 6763735909 6005337349 1758447152 0971150363 2298187894 2088556690 2781597327 5559990556 7337593149 9672711649 4197140420 0576028860 2301450050 4489761566 4971980202 5868106525 1020316853 6328434727 4723743748 5367123104 7920544530 2738413601 4749803246 8814395999 8443599262 6937165733 3683860007 2579912205 2546683095 2959734828 9383473919 6966554961 9229409580 9695200613 6287479535 7335077885 6021359566 1648148228 9054752088 3039148618 7042597377 5287457717 7030931615 4815900280 8937168480 8796578252 2467008135 1389685426 9438752685 6098535104 3595951306 0590969442 8333719228 3966025575 7492943050 5677858679 0525155151 7775167934 1209127652 2218408423 1478748697 6686480161 6212548578 7702300441 8442128386 3160712377 7265551932 3382452495 6571200519 8172493235 6204153855 3949429033 3782181372 7974633337 9171587341 2529980259 7776839095 2677260647 3903506769 0905189611 9619702997 9416304028 8983555286 7842105231 8467192026 1240100308 4980849592 4839510167 3987218968 7334759503 7173437985 0574874351 6813859156 5190840201 2461416585 5165969006 3931040035 5662835648 7538407534 6990776236 0347424742 4988953094 6345933115 6682062379 8186007182

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.