Primality Certificate for (3863^3697-1)/3862

Andy Steward13,258 digits18 September 2005
Originally by A.A.D.Steward 2005

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

From Factorisation
38633863
Φ22 · 2 · 2 · 3 · 7 · 23
Φ343 · 347131
Φ42 · 5 · 17 · 41 · 2141
Φ63 · 4972969
Φ728603 · 68474701 · 1697145031
Φ82 · 457 · 1889 · 128979497
Φ11397 · 281227 · 28064015792539 · 236244759352933
Φ1213 · 349 · 757 · 2341 · 27697
Φ147 · 29 · 16365896539249983629
Φ162 · 193 · 257 · 401 · 6920065553 · 180145183937
Φ21p44
Φ229439 · 1011671 · 1894934273 · 40886057403674891
Φ24337 · 12889 · 7455193 · 26102641 · 58668529
Φ28197 · 9281285580793429477 · 6039793724927330028473
Φ331123 · 481207 · 3047617 · p57
Φ42p44
Φ4489 · c70
Φ48417553 · 2274264193 · 48833141357286274753 · 53030795959916104584433
Φ56122656543845706087169 · c66
Φ6667 · 199 · 150369350671 · 4454186086153 · p44
Φ77463 · 51437 · 2018149529260835587 · p190
Φ84673 · 8309533 · 234543877 · p68
Φ88881 · 59174003334061529 · c124
Φ112113 · p171
Φ132661 · 5281 · c137
Φ154c216
Φ1681009 · 15474649 · c162
Φ176c287
Φ231c431
Φ2643064179889 · 120545181241 · c267
Φ3085430349 · c424
Φ336p345
Φ4624621 · c427
Φ5283169 · 364321 · 4031125599601 · c553
Φ616617 · c859
Φ9247393 · 2776621 · 43209013 · p843
Φ1232565489 · p1716
Φ1848c1722
Φ36963697 · c3440

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

936344 6818014947 8605211211 1522544716 3198520921 3873726440 1307956987 8524885540 8122534370 7867294192 1469585908 1378274411 1712909937 1602824589 2646943142 2240010939 4149628998 6526746275 5240667258 6740368593 2574750110 6623656690 8779996650 5245236241 0790384320 5330187782 9797824887 0795612999 0431668485 1391015754 6836142864 3787677021 8175355742 5048696521 6103202439 2697567268 3237114643 8688168160 9280982488 1125647574 1488602255 5627519035 9649831730 2482184277 5260497756 6112420445 7847247933 7914660840 9203677925 4692183485 0065913543 7587278064 6911587856 0969541805 5382598750 3117528383 4676917747 8101252766 1431008910 6707040508 2787879123 4338459149 7252976432 8957399869 6564414458 2762065485 7786226661 8061089647 7935616953 4549314364 2456388952 0178135863 2691906494 1479278986 9935291798 2492312470 3645059309 7817606867 4369130785 0152520454 0924550023 5333137971 7243675548 8776990802 0535732124 4414497069 8643873266 5465246701 4969054853 9959050327 1208834499 3021658773 9972094005 6429699104 2291696802 4978729357 6063460340 9250790786 1621048374 9612987141 5724462439 1645595245 7450938921 0039879074 9559049219 2947570034 7797295336 5486145757 6043070274 2772937791 2157975593 1910547940 1806421493 2699376133 3453767142 0824205399 5951087212 1789978651 0659768471 8685562361 5240741995 0615918778 3791914278 0664279181 8508759879 0439242328 9092743583 3896955626 6367941369 7014574440 2795066517 9205715620 8959434158 4507847520 6147353791 3415841282 7997137640 5771982347 3437446144 3476184495 7377579291 0966594150 4634029169 4527727227 6951886930 4654344228 7955966408 0465776585 0490516575 7502481619 7689869557 3693400208 6356684695 2740589412 8554315329 5549849896 0242747511 9312204072 1239905113 5227563414 6171040860 9136197267 5292637467 2956151206 2084688809 1450657462 4295086341 1018905254 0364187226 5510827546 6366301502 1429969809
820 3860546874 3341852684 1407363179 8174641932 2183545556 0749016663 9256871410 5651661415 8249382952 0133090387 2418840541 1539023431 2549308027 6542192342 3476302622 7142802223 6077801472 2922338858 2053342085 3487158957 3301848318 1148852457 4002317595 1455737288 8563644835 1746598250 3027217605 1293275248 4832529088 0264123401 8223756135 0214104143 8494455736 5112918678 1093650006 0026864878 2774919423 3578786340 0276271190 5608902738 9832537894 9733801635 7051377476 0670180884 7374373963 9939090371 6687617102 1201566000 5687890530 8930703599 6247602880 9987462070 1207350834 4683772716 9706349075 6860787215 9098798979 3781253586 8160751667 4513204787 2910959166 1911700846 6715424333 1173006787 8848001390 0630245272 7425914032 6992968407 4073427819 5354922243 0864747579 1140451905 2331380921 4787726554 8403183636 1342080819 2756400231 9486285171 4650086347 1225261710 6034605921 9863769797 2060477756 8454763409
22119 3307109341 8852423618 7740608244 1037608660 9453420297 0901689045 6271158405 6827340141 2728905297 2643463870 6754924377 4298621101 6972361442 2224700572 9205911235 3758645832 5916142484 7890116744 6047105500 7199717805 6877329561 6260915190 8695553717 3752696959 5641048156 8103039780 3806753852 5699147624 5512629818 5432249649 9192680640 4901592305 8724338916 2882037121
3416332142 7055780140 2397194383 4757175659 7688033142 0843471463 8159264651 7250829265 2346407956 2734124003 2057221395 6900171535 3110069810 0928809364 0746876897 5200577672 6578096647 0995063392 9780516297
1 3161565440 5987196646 7670370994 4747089595 8725651217 6768178623 4368986760 2076736264 1566404287 0957320263 9155138093 2473890098 9845607385 7956119280 8607947090 7480167376 6329065297
92977367 0759795015 4162420743 1817235670 7108757202 9533686344 8494541617
3324372 2177573757 1812640737 9157866604 8451301954 9896993117
6134 1022651257 1093500418 3961590537 0436086411
1104 6098391746 9637343701 1427086975 6533584689
1104 0380949898 2469446016 7032049594 4615770929
530 3079595991 6104584433
60 3979372492 7330028473
1 2265654384 5706087169
4883314135 7286274753
1636589653 9249983629
928128558 0793429477
201814952 9260835587
5917400 3334061529
4088605 7403674891
23624 4759352933
2806 4015792539
445 4186086153
403 1125599601
18 0145183937
15 0369350671
12 0545181241
6920065553
3064179889
2274264193
1894934273
1697145031
234543877
128979497
68474701
58668529
43209013
26102641
15474649
8309533
7455193
5430349
4972969
3047617
2776621
1011671
565489
481207
417553
364321
347131
281227
51437
28603
27697
12889
9439
7393
5281
4621
3863
3697
3169
2341
2141
1889
1123
1009
881
757
673
661
617
463
457
401
397
349
337
257
199
197
193
113
89
67
43
41
29
23
17
13
72
5
32
26

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

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 = 6479 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 = 1828367 3545935651 0761657203 1085201756 1772817693 8724026083 4446547258 7855005933 8497090974 5043904563 9535246062 0665027722 4983061488 2313931817 9164117680 2413929133 3070861056 5206645112 1219025430 6556638021 0959012906 6306144558 7821115628 2210154292 9798724816 7989370537 4290004031 0581647756 7707867486 2221318938 5585524184 9590341792 2823262591 7173301435 0685717065 8299196257 5535420161 1563783585 8024292869 1797942308 5567124658 3834831029 0232158379 9026974520 1620957893 8363529711 7256546678 9824425365 0714719103 4373394711 6965074444 8091288473 1151118075 0238163412 1793477454 5230758174 6764856976 1245696510 0242507183 5441114618 0478798379 1526621546 9439594457 2275247172 4882334507 1857310952 0077386435 6979040685 0284457565 6570457863 5872890264 8811281399 9802192423 3475480021 4888618394 8287157679 1491033736 3701839755 2994012289 7590874963 0156110842 9560321232 2271364104 3736990230 5942556761 4054142510 8800456779 7674525277 6612726913 6718753978 2190541947 9551782097 7528679080 5026287507 9912690088 9798313097 2797688173 2931885883 4780517700 7421967418 2767359845 4643911137 1762455707 6426101150 5311482991 6047092800 3908745513 8738345929 8470051513 9407806435 2715819992 7772520785 1186605521 2515697597 5974136026 5274703272 8169153233 3765478204 6061862917 6647253541 8300874933 7146910075 7939439491 5128734712 7519861997 8804425205 5189171022 7486574882 0492996452 7795160696 6914910511 1151632844 0999400123 4658124277 1947298314 5211160576 1894935674 7801556616 7232456926 3861744315 6957167049 3379828969 2202461728 6627228992 9238294817 4465658275 1786273438 3706521102 0306871582 4287866074 1035509345 5008215951 6674733545.

With those constraints, the unique continued fraction is: {0, 1, 1, 5, 1, 48, 1, 3, 2, 13, 2, 2, 4, 3, 2, 5, 3, 6, 2, 1, 1, 125, 1, 2, 2, 2, 1, 16, 1, 2, 3, 4, 2, 17, 4, 1, 1, 8, 1, 4, 2, 54, 2, 1, 2, 14, 3, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 8, 3, 1, 2, 1, 13, 1, 1, 5, 81, 1, 1, 2, 1, 1, 7, 1, 2, 2, 5, 2, 1, 2, 2, 1, 1, 6, 1, 1, 1, 5, 2, 1, 1, 1, 1, 2, 1, 2, 6, 1, 10, 2, 25, 1, 8, 17, 2, 4, 6, 3, 1, 1, 2, 1, 4, 6, 1, 4, 39, 1, 20, 1, 2, 3, 1, 9, 2, 1, 3, 1, 2, 1, 3, 1, 1, 3, 18, 3, 1, 5, 3, 9, 50, 3, 2, 4, 9, 4, 23, 1, 119, 2, 2, 2, 37, 1, 1, 1, 2, 1, 2, 3, 1, 2, 1, 2, 1, 2, 1, 10, 2, 2, 2, 3, 7, 1, 1, 19, 2, 22, 2, 1, 8, 3, 1, 3, 3, 2, 1, 12, 1, 1, 2, 3, 2, 1, 52, 1, 1, 3, 1, 2, 1, 72, 1, 2, 1, 2, 3, 29, 1, 1, 1, 2, 3, 3, 1, 1, 1, 1, 29, 3, 15, 1, 1, 15, 1, 1, 3, 1, 1, 1, 1, 2, 15, 2, 1, 1, 1, 1, 1, 5, 9, 4, 1, 3, 2, 54, 7, 1, 1, 16, 1, 2, 1, 2, 1, 18, 8, 3, 1, 1, 5, 1, 2, 10, 10, 3, 1, 1, 1, 1, 3, 3, 27, 1, 7, 1, 2, 1, 1, 1, 3, 3, 3, 5, 9, 1, 2, 9, 4, 1, 9, 13, 5, 10, 6, 12, 3, 1, 6, 9, 3, 2, 15, 30, 1, 2, 1, 10, 33, 2, 6, 1, 11, 2, 14, 2, 1, 4, 2, 4, 1, 5, 1, 9, 6, 3, 1, 3, 1, 4, 1, 4, 1, 2, 1, 1, 1, 17, 1, 1, 2, 6, 4, 1, 1, 3, 68, 13, 2, 24, 1, 19, 1, 16, 1, 2, 2, 2, 1, 3, 3, 2, 2, 25, 1, 4, 1, 90, 1, 1, 2, 2, 2, 1, 4, 1, 2, 6, 1, 3, 1, 1, 1, 4, 2, 1, 3, 1, 1, 3, 3, 1, 1, 2, 1, 2, 1, 2, 1, 6, 3, 1, 1, 3, 3, 2, 6, 2, 2, 5, 1, 1, 2, 1, 3, 3, 8, 1, 1, 1, 10, 41, 1, 3, 1, 2, 5, 3, 29, 1, 56, 6, 1, 4, 1, 2, 4, 1, 4, 1, 2, 2, 8, 1, 1, 56, 2, 1, 7, 1, 5, 1, 10, 2, 2, 7, 9, 3, 7, 5, 1, 2, 2, 2, 2, 5, 18, 1, 1, 13, 3, 12, 1, 13, 1, 10, 1, 1, 1, 2, 1, 1, 3, 2, 2, 2, 2, 3, 1, 2, 2, 1, 3, 1, 5, 1, 2, 1, 1, 1, 1, 3, 12, 1, 5, 1, 1, 5, 142, 1, 8, 3, 1, 5, 1, 4, 1, 1, 2, 3, 1, 2, 2, 5, 1, 2, 24, 3, 22, 1, 9, 1, 24, 1, 1, 2, 2, 46, 3, 4, 3, 4, 2, 4, 4, 1, 1, 3, 1, 1, 1, 1, 4, 1, 1, 2, 2, 1, 11, 2, 1, 1, 1, 10, 2, 1, 5, 1, 3, 3, 1, 1, 5, 1, 2, 1, 1, 1, 6, 4, 1, 1, 53, 1, 4, 2, 1, 1, 2, 1, 1, 2, 2, 1, 2, 2, 2, 1, 1, 10, 4, 2, 3, 3, 1, 3, 3, 6, 1, 6, 2, 1, 3, 15, 28, 1, 1, 4, 1, 2, 1, 3, 7, 1, 1, 1, 5, 1, 1, 5, 14, 2, 1, 15, 31, 1, 5, 1, 1, 1, 4, 40, 1, 3, 1, 2, 1, 3, 6, 1, 2, 3, 1, 1, 1, 1, 13, 1, 3, 1, 12, 2, 1, 25, 1, 6, 2, 2, 3, 1, 3, 12, 3, 6, 2, 17, 1, 8, 1, 290, 55, 1, 6, 2, 2, 1, 3, 1, 1, 1, 1, 7, 1, 4, 1, 1, 1, 3, 2, 1, 4, 2, 5, 1, 3, 1, 2, 3, 1, 13, 2, 1, 1, 1, 1, 2, 1, 2, 11, 2, 1, 1, 17, 1, 1, 8, 3, 6, 10, 1, 1, 1, 11, 1, 1, 2, 1, 1, 1, 1, 8, 3, 1, 15, 1, 1, 2, 3, 1, 2, 1, 1, 4, 27, 1, 1, 2, 31, 1, 4, 2, 21, 1, 1, 1, 2, 1, 9, 1, 1, 1, 4, 1, 11, 2, 1, 4, 3, 1, 4, 2, 1, 1, 2, 1, 1, 1, 1, 5, 4, 4, 1, 1, 2, 3, 9, 3, 1, 2, 18, 1, 2, 1, 6, 3, 7, 9, 8, 2, 1, 4, 3, 4, 2, 1, 1, 1, 1, 12, 2, 2, 1, 2, 1, 1, 1, 97, 1, 8, 6, 2, 3, 2, 17, 1, 1, 2, 3, 2, 4, 1, 2, 3, 1, 9, 1, 2, 1, 6, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 1, 28, 1, 8, 1, 4, 1, 2, 2, 1, 7, 1, 1, 1, 2, 1, 160, 1, 1, 2, 4, 1, 5, 1, 2, 3, 3, 3, 1, 3, 6, 2, 8, 1, 10, 1, 1, 1, 1, 1, 21, 4, 1, 1, 46, 1, 24, 20, 2, 53, 4, 1, 3, 2, 1, 7, 2, 24, 1, 3, 2, 1, 1, 1, 1, 1, 3, 8, 1, 2, 1, 1, 3, 1, 3, 89, 1, 1, 1, 9, 1, 2, 2, 2, 1, 6, 3, 1, 4, 1, 3, 6, 5, 2, 1, 3, 4, 1, 4, 3, 7, 1, 3, 6, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 3, 1, 64, 12, 2, 1, 2, 1, 8, 1, 2, 12, 1, 6, 1, 3, 3, 3, 1, 4, 9, 1, 7, 1, 1, 1, 3, 3, 3, 1, 1, 6, 20, 2, 4, 3, 2, 1, 3, 1, 1, 5, 5, 2, 2, 20, 3, 2, 2, 2, 5, 1, 1, 11, 1, 4, 7, 10, 2, 27, 4, 1, 3, 2, 5, 3, 2, 22, 1, 17, 2, 13, 120, 1, 1, 22, 3, 4, 1, 3, 2, 1, 8, 7, 4, 1, 2, 1, 3, 1, 1, 6, 2, 1, 1, 1, 2, 1, 73, 2, 3, 5, 2, 1, 2, 3, 2, 1, 1, 1, 2, 1, 2, 1, 4, 2, 1, 3, 4, 2, 1, 1, 1, 2, 1, 7, 1, 2, 11, 1, 1, 1, 1, 2, 1, 2, 1, 8, 1, 2, 1, 3, 1, 1, 1, 1, 7, 1, 1, 2, 2, 102, 3, 3, 1, 22, 2, 7, 1, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 1, 9, 10, 1, 2, 10, 3, 1, 1, 2, 1, 2, 1, 2, 2, 1, 3, 1, 1, 1, 1, 5, 1, 1, 1, 40, 1, 137, 6, 98, 1, 2, 1, 12, 1, 2, 2, 1, 1, 1, 1, 1, 2, 135, 4, 2, 14, 1, 2, 3, 2, 1, 1, 2, 1, 4, 2, 4, 1, 2, 2, 3, 1, 2, 8, 1, 19, 1, 1, 2, 1, 1, 5, 1, 1, 4, 1, 68, 3, 2, 2, 4, 1, 17, 1, 3, 3, 1, 3, 11, 25, 1, 1, 1, 4, 1, 1, 2, 1, 3, 2, 1, 4, 2, 3, 1, 2, 1, 1, 8, 1, 1, 2, 13, 1, 1, 1, 9, 3, 1, 13, 1, 6, 3, 3, 1, 23, 3, 17, 1, 9, 25, 2, 574, 4, 4, 6, 2, 5, 1, 1, 1, 15, 1, 29, 3, 2, 3, 1, 1, 1, 1, 8, 456, 2, 5, 1, 3, 1, 2, 3, 7, 1, 2, 1, 2, 11, 1, 27, 1, 1, 7, 1, 147, 1, 1, 5, 2, 113, 5, 5, 1, 4, 1, 4, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 19, 1, 1, 3, 19, 1, 1, 4, 2, 2, 5, 1, 18, 9, 3, 3, 1, 4, 5, 17, 3, 3, 1, 6, 2, 5, 1, 1, 2, 1, 1, 5, 1, 1, 1, 5, 1, 4, 2, 5, 2, 21, 1, 1, 3, 1, 1, 4, 18, 1, 1, 2, 1, 5, 3, 2, 2, 4, 2, 4, 2, 2, 1, 1, 2, 4, 14, 1, 5, 3, 2, 4, 1, 1, 1, 1, 2, 1, 6, 1, 2, 3, 3, 1, 2, 33, 1, 1, 1, 5, 1, 10, 4, 1, 1, 2, 1, 1, 1, 2, 2, 1, 4, 6, 1, 2, 2, 17, 1, 2, 1, 1, 2, 1, 2, 1, 5, 5, 1, 4, 1, 11, 1, 8, 2, 3, 3, 8, 5, 6, 1, 1, 9, 6, 2, 1, 1, 1, 3, 20, 1, 3, 451, 1, 5, 7, 4, 3, 83, 1, 5, 3, 4, 1, 18, 1, 1, 5, 1, 5, 1, 1, 1, 1, 3, 2, 2, 2, 5, 2, 2, 1, 5, 1, 9, 1, 3, 2, 3, 3, 1, 1, 11, 2, 3, 1, 1, 1, 9, 10, 3, 2, 4, 1, 2, 4, 1, 3, 1, 1, 1, 11, 1, 1, 1, 2, 5808, 1, 1, 1, 84, 4, 1, 5, 1, 1, 7, 1, 13, 1, 1, 1, 1, 1, 1, 1, 1, 47, 12, 3, 1, 1, 1, 7, 14, 2, 2, 2, 1, 1, 56, 1, 1, 20, 1, 5, 2, 3, 1, 1, 1, 1, 1, 1, 469, 1, 1, 1, 1, 1, 10, 6, 3, 1, 2, 11, 3, 2, 1, 8, 1, 1, 2, 1, 2, 21, 1, 3, 3, 3, 12, 1, 17, 4, 10, 1, 14, 2, 3, 1, 2, 2, 2, 15, 2, 3, 1, 8, 2, 4, 1, 1, 1, 41, 1, 1, 3, 1, 1, 2, 6, 1, 3, 1, 5, 1, 1, 1, 2, 10, 3, 3, 3, 21, 1, 1, 1, 2, 18, 10, 2, 4, 23, 2, 1, 4, 1, 54, 2, 93, 1, 11, 1, 1, 1, 1, 5, 6, 1, 1, 2, 3, 1, 2, 4, 1, 1, 7, 2, 67, 1, 3, 9, 1, 14, 1, 89, 29, 1, 1, 3, 1, 35, 1, 1, 2, 3, 1234, 1, 1, 29, 7, 1, 5, 2, 2, 3, 3, 1, 10, 2, 1, 128, 2, 1, 4, 2, 2, 1, 1, 1, 8, 5, 3, 1, 2, 7, 9, 4, 2, 1, 2, 292, 2, 1, 5, 4, 3, 2, 13, 27, 1, 1, 1, 3, 2, 7, 2, 1, 1, 1, 13, 1, 1, 2, 1, 2, 1, 8, 5, 9, 4, 1, 1, 6, 2, 1, 3, 1, 1, 2, 1, 2, 1, 1, 15, 2, 1, 9, 1, 40, 1, 1, 2, 2, 3, 4, 1, 39, 1, 2, 3, 1, 27, 4, 3, 2, 13, 1, 2, 1, 1, 4, 2, 143, 2, 1, 1, 2, 46, 1, 7, 1, 2, 3, 2, 2, 37, 35, 1, 2, 1, 1, 7, 3, 2, 10, 1, 1, 1, 9, 8, 1, 5, 3, 1, 4, 6, 3, 119, 1, 26, 1, 4, 1, 1, 1, 3, 14, 1, 3, 2, 2, 2, 5, 1, 8, 1, 4, 2, 4, 2, 6, 3, 1, 2, 4, 2, 23, 1, 1, 9, 1, 4, 2, 2, 1, 8, 4, 2, 9, 1, 87, 1, 30, 4, 5, 8, 1, 1, 2, 3, 10, 1, 33, 1, 1, 1, 1, 4, 2, 2, 1, 1, 1, 3, 1, 7, 1, 1, 6, 2, 1, 63, 3, 4, 2, 2, 1, 1, 4, 4, 5, 1, 1, 5, 1, 7, 262, 1, 3, 1, 2, 1, 3, 3, 5, 1, 1, 5, 1, 1, 14, 25, 1, 2, 1, 1, 3, 6, 2, 4, 1, 907, 116, 4, 1, 2, 7, 2, 1, 36, 2, 1, 1, 3, 1, 2, 146, 1, 106, 27, 1, 4, 3, 1, 21, 15, 2, 20, 1, 5, 22, 1, 6, 4, 1, 1, 4, 1, 6, 2, 1, 1, 1, 1, 1, 2, 1, 1, 6, 1, 9, 7, 1, 1, 1, 1, 4, 3, 4, 3, 6, 4, 13, 3, 25, 5, 2, 1, 4, 12, 1, 3, 692, 1, 19, 2, 1, 21, 1, 80, 1, 2, 2, 3, 15, 1, 9, 9, 1, 1, 104, 5, 2, 4, 1, 3, 1, 4, 1, 1, 1, 1, 13, 2, 1, 1, 95, 1, 5, 5, 7, 1, 44, 6, 1, 1, 101, 2, 19, 2, 1, 70, 4, 1, 2, 1, 6, 5, 2, 3, 2, 2, 1, 12, 47, 10, 3, 1, 2, 1, 1, 10, 2, 1, 1, 2, 1, 65, 4, 6, 1, 3, 1, 3, 4, 3, 1, 6, 1, 2, 8, 1, 5, 1, 1, 4, 2, 1, 1, 7, 4, 6, 1, 7, 3, 3, 6, 3, 2, 1, 3, 1, 3, 2, 1, 3, 4, 3, 9, 40, 1, 1, 3, 16, 1, 2, 2, 1, 1, 3, 2, 6, 1, 1, 37, 2, 1, 15, 2, 22, 1, 101, 1, 2, 11, 1, 5, 1, 11, 1, 3, 85, 3, 2, 1, 3, 2, 2, 1, 10, 1, 1, 3, 1, 1, 1, 3, 18, 2, 6, 1, 5, 1, 1, 4, 15, 2, 2, 1, 1, 2, 1, 5, 1, 31, 8, 1, 1, 9, 3, 1, 3, 6, 966, 1, 12, 5, 2, 10, 1, 1, 1, 2, 2, 1, 4, 1, 2, 2, 16, 4, 1, 9, 8, 1, 2, 1, 3, 1, 3, 2, 3, 1, 1, 1, 5, 1, 18, 1, 2, 47, 2, 156, 8, 1, 2, 4, 3, 1, 2, 4, 2, 2, 1, 5, 5, 4, 3, 2, 1, 2, 1, 8, 1, 3, 33, 7, 1, 11, 2, 2, 1, 4, 1, 6, 1, 1, 43, 3, 1, 1, 1, 2, 42, 1, 23, 5, 1, 1, 3, 1, 1, 63, 1, 12, 1, 2, 3, 6, 1, 24, 4, 1, 2, 4, 1, 1, 3, 2, 1, 10, 4, 1, 3, 5, 39, 1, 15, 6, 2, 1, 6, 1, 2, 12, 2, 8, 2, 2, 24, 33, 3, 6, 3, 3, 3, 1, 5, 1, 2, 1, 1, 1, 3, 1, 6, 1, 1, 2, 9, 1, 3, 1, 4, 6, 22, 1, 1, 6, 3, 3, 3, 1, 2, 4, 2, 1, 2, 2, 1, 1, 1, 2, 6, 33, 4, 6, 3, 1, 1, 1, 23, 4, 1, 7, 1, 3, 1, 23, 2, 8, 1, 1, 10, 1, 1, 1, 12, 3, 1, 11, 4, 8, 2, 7, 1, 1, 8, 3, 1, 4, 2, 4, 2, 2, 1, 1, 5, 1, 1, 1, 2, 1, 17, 1, 1, 6, 3, 8, 4, 2, 1, 1, 2, 1, 1, 2, 2, 1, 14, 1, 2, 9, 1, 6, 1, 12, 9, 2, 1, 1, 8, 1, 3, 1, 3, 6, 2, 1, 1, 13, 1, 5, 1, 9, 1, 3, 1, 2, 4, 1, 2, 1, 4, 235, 1, 1, 30, 2, 1, 105, 1, 5, 14, 5, 1, 2, 1, 1, 1, 2, 1, 16, 86, 3, 2, 1, 5, 2, 1, 7, 54, 1, 21, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 8, 1, 1, 1, 3, 4, 2, 37, 1, 5, 2, 5, 2, 75, 1, 4, 11, 1, 48, 24, 3, 5, 1, 3, 77, 2, 19, 189, 1, 9, 1, 11, 1, 2, 1, 2, 3, 1, 6, 1, 14, 1, 3, 2, 4, 1, 1, 5, 33, 1, 9, 1, 2, 1, 8, 1, 1, 4, 8, 3, 1, 1, 3, 2, 1, 1, 2, 1, 3, 1, 2, 1, 6, 13, 5, 5, 1, 7, 1, 9, 2, 2, 1, 6, 1, 47, 3, 3, 3, 6, 98, 1, 27, 4, 3, 2, 2, 1, 2, 1, 1, 3, 1, 2, 20, 9, 2, 1, 1, 1, 4, 5, 1, 8, 1, 11, 3, 1, 3, 2, 2, 37, 2, 90, 1, 1, 2, 1, 1, 9, 7, 1, 14, 3, 2, 4, 2, 1, 4, 2, 16, 1, 8, 1, 2, 4, 1, 2, 2, 1, 95, 1, 2, 1, 1, 1, 2, 1, 6, 4, 10, 2, 2, 1, 23, 1, 2, 14, 2, 1, 2, 2, 1, 3, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 2, 15, 3, 10, 2, 4, 14, 5, 3, 11, 1, 3, 1, 4, 1, 1, 1, 3, 1, 9, 17, 1, 1, 1, 1, 17, 5, 1, 6, 3, 8, 3, 1, 1, 3, 1, 2, 1, 13, 1, 48, 1, 2, 4, 1, 150, 1, 30, 10, 4, 6, 3, 2, 28, 2, 1, 5, 1, 4, 178, 11, 1, 1, 2, 4, 3, 1, 6, 3, 1, 10, 7, 4, 3, 10, 2, 3, 5, 3, 63, 1, 1, 1, 4, 5, 2, 49, 1, 2, 1, 1, 3, 1, 3, 2, 2, 2, 1, 1, 1, 7, 2, 50, 2, 1, 22, 1, 3, 1, 1, 3, 1, 1, 1, 5, 3, 1, 4, 1, 2, 1, 1, 26, 1, 5, 16, 1, 5, 19, 1, 76, 4, 2, 3, 4, 20, 31, 2, 6, 1, 1, 1, 9, 1, 1, 5, 5, 1, 2, 1, 1, 1, 4, 1, 6, 16, 2, 1, 3, 3, 1, 2, 1, 1, 57, 10, 5, 4, 1, 4, 4, 1, 1, 1, 12, 1, 6, 1, 3, 22, 1, 1, 1, 2, 14, 5, 4, 2, 41, 1, 10, 1, 3, 2, 2, 1, 2, 1, 10, 7, 1, 5, 1, 1, 1, 7, 5, 1, 5, 6, 2, 43, 1, 9, 1, 2, 1, 14, 1, 2, 12, 10, 16, 1, 2, 1, 1, 7, 33, 1, 1, 2, 1, 1, 2, 5, 1, 2, 23, 3, 3, 4, 4, 2, 1, 1, 2, 3, 2, 13, 1, 1, 2, 41, 2, 2, 8, 1, 28, 1, 2, 1, 1, 1, 7, 1, 1, 3, 1, 12, 4, 5, 2, 5, 7, 3, 1, 2, 5, 1, 14, 11, 1, 2, 4, 1}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 12 7521398749 3104256507 4501846520 9123342894 5065553901 4847883525 1055766994 9691948947 6363588403 3953669507 3416746141 3651582951 6605983245 9180142698 3254761548 5204658419 9775119320 7189624657 2190594591 6970126659 9581472677 3342225072 2115282677 4801794475 8244516035 7807190230 1363039955 9429599068 6638498378 9345133168 8794354435 9493596227 2004673618 1695815817 1716545512 3311044045 9526187382 0709094856 7260014765 4302562051 6952035154 9591671613 5826328994 7059373167 5577155758 4243182564 7864894511 0090256902 5211631741 0988534698 2870794556 1922512968 7389147254 9934958986 2298854983 6595237161 2884387020 4037148730 3869490451 0833145758 8370664951 3020582790 5667793370 5899218687 6063002627 5381483413 2670988162 0218877074 0916102895 3062887897 5203494922 2299721727 1156189055 7728877998 3945000028 7895392765 7714229490 8554056771 9847240899 2626195948 9557463746 2549310611 3833397414 0261803753 6224387279 9665999571 5250157189 2815101544 9632147434 9425270918 6920497716 7525745347 7760282201 5492893833 1115770821 5048783763 0134033659 6597277334 1266011047 3981145649 1727747927 9980155043 8541004129 9402410646 1615065904 4149278969 5848725280 6008018914 0009808771 3705309348 3432809657 6886882893 1885457978 8579062416 0311424057 5649262371 5546488183 5507986900 8927797318 6071657279 3342333208 6629903144 3218607818 2036210921 1578245521 0810494620 7762890571 4059906740 5831651596 3149627014 7980808489 0736356230 7386945669 2738971842 9927598029 3537882405 0994214164 7883816273 6153861478 4528728097 0087535255 3718985077 5099743242 9980598023 3680777429 5039518864 3318913826 8519503282 8437081221 6555526824 6720284872 1299826828 9010351440 7302309065 1856981471 9321823430 6210576653 5548497852 0193797058 8230546793 1621284735 6885879326 7112288569 0457534198 0003968866 3996849637 7009491174 1954913732 5674522948 4932339324 0393895003 0563434305 0210952247 1588643814 4553165749 3227056474 4327021889 1548755365 3945101565 6035627253 8964361638 3826697260 8418126337 2008561726 0627140266 4446589861 9441410970 4012222868 0900312538 7454322287 4800230515 1807494093 4879437822 0818750341 3346413218 1228438261 4163106074 2873218652 3321760354 1888462919 4477464416 9106790782 9422998365 5803646908 9211411870 2674371752 2184102478 8597956266 5495191596 4266970134 5915420699 7590364656 9316779789 7998323976 3952794732 3709548847 8198213119 7642175414 8954516323 3172261740 7472861218 6281116066 0476979194 8097809148 2755025504 3784321018 1782396653 1132125529 8491827929 1805752235 3358216780 8125531876 5349607396 8584017699 8872237790 6097307299 9474506623 8606320376 5369440769 2920483169 3731490975 1858027403 4670581607 3278334565 4388505092 9597904712 7316983049 5460104969 6465677158 7553989910 5032957421 7763197458 1123290783 7689598908 6850223442

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.