Primality Certificate for (10268^2521-1)/10267

Andy Steward10,109 digits10 July 2001
Originally by David Broadhurst & Bouk de Water 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.003359% factorization of N-1:

From Factorisation
102682 · 2 · 17 · 151
Φ23 · 3 · 7 · 163
Φ3181 · 439 · 1327
Φ45 · 5 · 113 · 37321
Φ511 · 41 · 61 · 1151 · 351078341
Φ63 · 19 · 127 · 14563
Φ71437409 · 815411999547636877
Φ86361 · 533633 · 3274729
Φ9487 · 4933 · 36433 · 39727 · 337050037
Φ10211 · 52676715845591
Φ1213 · 10616029 · 80544889
Φ147 · 5160737 · 32438674425192323
Φ1531 · 271 · 1840771 · 7989393898872489383431
Φ183 · 5023 · 83542429 · 930944153593
Φ205 · 84421 · 292729424752397034665957081
Φ21337 · p46
Φ24p33
Φ2829 · 23010886515121849253617 · 2058253996719560091684149
Φ305799072001 · 21309373074972685969501
Φ3534511 · 13224944531 · 12948481128046241 · 3840853858940700320602857612241 · p35
Φ3637 · 73 · 17700049 · p38
Φ40281 · 105935561 · 16539772008468921761338121 · 31009599022794954617064793241
Φ4243 · 19237 · 65111030522547423457 · 25504231643469945221227
Φ45631 · 1171 · p91
Φ56p97
Φ60661 · 555016741 · p53
Φ633529 · 10105518582601083363539413 · c116
Φ7071 · 117987089552365404412426775807233061914259081 · p51
Φ72756340035185291284726528035461376765747889 · p55
Φ84123733 · 1210021 · 8472393889 · 4011711817211719589809 · p54
Φ9043466005506345931 · 8534149827167366173441 · 102483856269510069040471 · p35
Φ10510096142019258485767921 · p171
Φ120123121 · 6342125161 · 42546019066147293432026445601 · p85
Φ126204751 · 648194415634459 · 886818637550051565134586493083031 · p92
Φ140421 · 238227522966217801661 · 815289212152904594875688101 · 20577640454128808069684788481 · c115
Φ16813314739109327229841 · 17733968981520143689 · p155
Φ1805581 · 9773782606036243261 · 1595680248178915687517581 · 7851965077308078827014927141 · c118
Φ2103431596498568917886761 · c172
Φ25215121 · c285
Φ2808681 · c382
Φ3158191 · 319411 · 698903731 · 1732554181 · 571908978632441041201 · 271246581635391893517211 · c506
Φ3602521 · 4252755684652201 · 609883420371658921101238288558561 · c334
Φ42016906850234297941 · c369
Φ5042017 · 1595751697 · c566
Φ63054181 · 10115269921 · 9340892016669157141 · p544
Φ84018481 · 44407743527281 · 597506276294994838357822321 · c726
Φ126052388093521 · 239971793041 · 3890718996973201 · 673973696515418501095432216536624361 · c1082
Φ2520c2311

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 :

8805 4386542656 0088514644 9870300100 9499974776 3093551982 1355606987 2264439071 2721903616 6579905874 8341854080 3829347319 7984544456 9144795167 5348346537 4102209301 7115963371 2940954548 9336728170 3223760936 1067734938 0392447309 3673762292 9121687137 0127176227 7569683283 1644546814 2924650414 0710203592 9334641886 9436308614 3844963707 6381602070 7916821840 3791084547 5755555940 9728333735 5098666360 5869764492 5969326971 2338828055 9974700674 9252998711 1525943426 8771443399 9317614396 7996225873 6185808919 4196945410 5829504644 7984759360 4665546156 0162515683 6523517013 5830931161
3 5253977818 5605608053 1522529072 8522085877 8814861624 9446917586 5601798978 0345907925 6167098411 5889583925 6408030648 3136399244 3766224636 0131961579 7721913480 2494204865 2837878781
15072 4091300058 2837735624 2839993925 1224554076 3075035704 1341235473 0175996680 9125381193 2569606177 8154298677 3302814296 1274464064 6771703041 8448297032 6731753449
1886516 6388791145 8453508204 0679370539 3035550387 5147386097 6349179545 2292089871 7511260697 9962412801
22 0153113364 0571928167 0921695756 8600644828 8216588488 9266881318 3837228179 5855149397 9969414731
2 5531385650 8161733815 6006375079 4312008999 4065876374 2988569335 4978644315 3722173021 3736663101
70165 1805153409 4067639966 9423819551 6214104372 2549585963 5951450297 6537934548 3848091721
24942 7050151714 9186169028 5883327247 6568221178 7257456209
3707 2026560210 0115545349 3552145741 3399773093 4089990457
416 1656410317 8113665177 4748938982 7286675825 7886659401
2 2522162976 7856861494 9844844118 8856115015 9018490851
407528 6268307431 8211879381 8097860358 2402457877
11798 7089552365 4044124267 7580723306 1914259081
75 6340035185 2912847265 2803546137 6765747889
28729711 6505354096 4888269753 7765234997
673973 6965154185 0109543221 6536624361
83103 6938479095 4955985235 1652996281
49624 3907017051 6539684389 1390922261
886 8186375500 5156513458 6493083031
609 8834203716 5892110123 8288558561
123 5625550070 7692607860 3002617601
3 8408538589 4070032060 2857612241
425460190 6614729343 2026445601
310095990 2279495461 7064793241
205776404 5412880806 9684788481
78519650 7730807882 7014927141
8152892 1215290459 4875688101
5975062 7629499483 8357822321
2927294 2475239703 4665957081
165397 7200846892 1761338121
101055 1858260108 3363539413
20582 5399671956 0091684149
15956 8024817891 5687517581
2712 4658163539 1893517211
1024 8385626951 0069040471
255 0423164346 9945221227
230 1088651512 1849253617
213 0937307497 2685969501
100 9614201925 8485767921
85 3414982716 7366173441
79 8939389887 2489383431
40 1171181721 1719589809
34 3159649856 8917886761
5 7190897863 2441041201
2 3822752296 6217801661
6511103052 2547423457
1773396898 1520143689
1331473910 9327229841
977378260 6036243261
934089201 6669157141
81541199 9547636877
4346600 5506345931
3243867 4425192323
1690685 0234297941
1294848 1128046241
425275 5684652201
389071 8996973201
64819 4415634459
5267 6715845591
4440 7743527281
93 0944153593
23 9971793041
5 2388093521
1 3224944531
1 0115269921
8472393889
6342125161
5799072001
1732554181
1595751697
698903731
555016741
351078341
337050037
105935561
83542429
80544889
17700049
10616029
5160737
3274729
1840771
1437409
1210021
533633
319411
204751
123733
123121
84421
54181
39727
37321
36433
34511
19237
18481
15121
14563
8681
8191
6361
5581
5023
4933
3529
2521
2017
1327
1171
1151
661
631
487
439
421
337
281
271
211
181
163
151
127
113
73
71
61
43
41
37
31
29
19
17
13
11
72
53
34
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.003359%

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 = 2 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 = 37 4670231288 4486392544 0573743050 8979732812 3071727531 2334481059 6379232239 3959839754 1265955645 7231550778 5200172564 1256536802 5564741918 8562388312 2284162192 3647989197 8089887758 0889450234 2018303703 6498776563 1482236876 8422569098 1045103760 1004579664 7422531035 3816472285 5692471152 0052705224 3215141272 8270633974 7744640810 6539084903 5827946376 3705154816 4524707251 6411719225 3742645901 8831875344 0449583088 6428481959 7375646996 6837526290 3302502571 6972345853 2568679604 6555924117 3458069819 3606920799 4086121820 3730061880 8675020300 3318364767 3575711109 1262080513 8957236515 0239102754 6890558842 0207517710 2710745862 8962623273 0533153495 7725558114 7023086449 0964380173 5707555408 9442130187 3585945323 9452005854 0901947395 9924147069 3311667954 4834601359 0137859940 9214649669 2693836586 8008717201 6505176061 4184313061 6126725865 0933258610 0358555747 7360124404 5538787763 2288294337 2794108809 0811814771 2097148623 9009130838 4347938855 1964988544 7216292371 3365673974 2327167274 1586420128 2564914535 7369677142 1244242382 9130107003 5501695649 6245468883 5579429200.

With those constraints, the unique continued fraction is: {0, 2, 2, 1, 1, 2, 1, 1, 22, 3, 2, 5, 1, 17, 1, 1, 11, 3, 2, 1, 2, 1, 2, 3, 1, 2, 4, 4, 1, 1, 1, 6, 5, 1, 5, 1, 23, 884, 1, 17, 2, 1, 1, 3, 36, 1, 1, 1, 26, 1, 16, 3, 1, 2, 1, 3, 1, 2, 2, 1, 2, 4, 1, 4, 2, 1, 1, 8, 42, 1, 4, 9, 2, 8, 24, 4, 1, 2, 4, 1, 3, 1, 5, 34, 2, 1, 1, 1, 7, 11, 1, 35, 11, 2, 2, 1, 25, 1, 1, 1, 2, 3, 31, 1, 1, 3, 9, 1, 2, 2, 2, 3, 1, 1, 1, 600, 1, 1, 2, 1, 2, 15, 5, 4, 43, 4, 7, 2, 1, 86, 2, 13, 7, 1, 3, 1, 3, 1, 9, 1, 57, 1, 1, 1, 1, 1, 1, 2, 1, 3, 4, 35, 1, 1, 1, 1, 9, 1, 21, 2, 1, 3, 1, 2, 2, 25, 1, 1, 8, 2, 19, 5, 1, 6, 17, 1, 1, 3, 3, 34, 1, 9, 1, 2, 2, 1, 3, 6, 7, 1, 39, 1, 13, 2, 3, 2, 8, 106, 2, 9, 1, 3, 1, 1, 1, 3, 3, 2, 1, 2, 1, 7, 1, 1, 2, 4, 2, 2, 11, 2, 1, 1, 1, 5, 1, 307, 1, 1, 1, 5, 1, 21, 1, 1, 8, 3, 1, 2, 2, 19, 4, 2, 1, 9, 9, 31, 1, 1, 1, 304, 2, 1, 4, 1, 7, 2, 1, 9, 1, 1, 3, 5, 2, 8, 1, 1, 3, 5, 2, 3, 1, 2, 3, 1, 1, 7, 19, 63, 1, 3, 12, 3, 13, 5, 1, 2, 1, 11, 2, 1, 2, 1, 17, 2, 3, 1, 5, 2, 1, 1, 1, 1, 1, 220, 1, 4, 1, 1, 1, 1, 1, 2, 1, 4, 1, 1, 2, 4, 4, 1, 2, 1, 8, 1, 1, 1, 4, 1, 1, 4, 1, 2, 1, 1, 3, 8, 1, 3, 23, 6, 1, 1, 8, 2, 7, 2, 10, 1, 1, 5, 4, 1, 1, 1, 11, 1, 1, 8, 33, 4, 1, 1, 1, 1, 27, 2, 37, 1, 14, 1, 17, 1, 4, 7, 5, 1, 4, 1, 1, 1, 2, 1, 1, 2, 1, 3, 1, 87, 1, 3, 3, 1, 13, 1, 7, 1, 1, 2, 2, 1, 2, 1, 1, 1, 1, 11, 1, 4, 8, 1, 1, 1, 2, 1, 6, 12, 3, 3, 9, 4, 4, 25, 4, 4, 1, 1, 1, 3, 2, 13, 3, 1, 2, 2, 3, 1, 1, 1, 10, 46, 3, 1, 16, 1, 4, 1, 1, 2, 3, 1, 17, 2, 1, 1, 1, 1, 6, 1, 1, 5, 4, 1, 3, 1, 1, 4, 2, 1, 1, 2, 3, 3, 2, 90, 24, 11, 4, 3, 4, 1, 4, 9, 2, 1, 6, 3, 27, 14, 3, 1, 6, 2, 18, 2, 1, 1, 3, 2, 3, 1, 12, 39, 2, 2, 1, 10, 1, 2, 1, 2, 1, 14, 14, 7, 1, 2, 1, 19, 1, 3, 43, 2, 2, 3, 3, 3, 4, 20, 1, 14, 2, 2, 19, 1, 2, 4, 30, 1, 16, 1, 1, 10, 7, 25, 1, 5, 1, 23, 5, 2, 2, 1, 4, 1, 2, 6, 3, 1, 2, 6, 4, 134, 1, 1, 1, 52, 25, 1, 11, 1, 1, 2, 1, 1, 11, 1, 1, 3, 4, 2, 4, 1, 3, 2, 4, 1, 13, 1, 8, 3, 1, 1, 1, 6, 1, 1, 1, 8, 15, 1, 1, 3, 26, 1, 11, 3, 1, 2, 34, 10, 5, 3, 1, 1, 2, 5, 25, 3, 15, 2, 2, 1, 1, 2, 1, 1, 3, 25, 9, 1, 1, 2, 1, 37, 1, 2, 2, 4, 2, 47, 2, 3, 2, 28, 1, 2, 2, 1, 2, 1, 1, 35, 2, 8, 1, 5, 2, 7, 5, 1, 1, 4, 3, 1, 1, 2, 2, 53, 1, 1, 2, 2, 8, 1, 1, 3, 3, 2, 1, 6, 2, 2, 2, 16, 1, 1, 1, 34, 10, 1, 24, 2, 4, 2, 3, 2, 3, 8, 1, 2, 3, 1, 1, 3, 12, 1, 5, 52, 1, 6, 2, 55, 2, 1, 51, 4, 1, 2, 17, 6, 1, 1, 1, 10, 10, 1, 1, 28, 1, 3, 110, 1, 3, 2, 1, 2, 18, 12, 3, 4, 1, 6, 4, 3, 27, 6, 1, 1, 2, 7, 1, 79, 2, 2, 176, 1, 90, 4, 2, 8, 1, 3, 12, 1, 1, 1, 25, 13, 1, 2, 26, 1, 5, 1, 13, 2, 4, 226, 4, 6, 1, 2, 9, 1, 4, 25, 1, 4, 1, 5, 3, 8, 1, 3, 51, 1, 5, 12, 3, 8, 1, 1, 1, 2, 2, 4, 6, 7, 6, 1, 4, 146, 4, 1, 1, 10, 5, 7, 28, 1, 2, 2, 156, 3, 3, 4, 2, 3, 3, 6, 1, 4, 1, 18, 1, 1, 1, 2, 1, 2, 1, 4, 1, 2, 1, 8, 93, 4, 1, 2, 1, 8, 1, 10, 1, 1, 4, 18, 38, 1, 1, 1, 2, 18, 3, 2, 2, 2, 1, 28, 13, 4, 13, 2, 3, 2, 1, 6, 2, 2, 9, 1, 6, 2, 1, 2, 1, 4, 1, 1, 5, 1, 1, 4, 2, 1, 3, 11, 3, 75, 1, 13, 1, 3, 7, 2, 1, 1, 1, 7, 4, 1, 24, 1, 1, 1, 1, 3, 4, 5, 3, 3, 1, 120, 2, 1, 3, 1, 4, 1, 1, 3, 4, 1, 6, 3, 2, 10, 1, 2, 1, 23, 2, 5, 71, 9, 1, 7, 1, 3, 1, 1, 1, 17, 2, 2, 4, 4, 3, 1, 2, 11, 1, 1, 1, 2, 1, 7, 2, 40, 3, 1, 2, 1, 1, 3, 1, 4, 1, 2, 52, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 6, 17, 2, 13, 1, 2, 4, 72, 2, 2, 3, 3, 7, 1, 1, 3, 1, 5, 1, 9, 3, 10, 1, 2, 1, 16, 28, 1, 5, 24, 1, 1, 4, 4, 1, 21, 3, 1, 1, 1, 19, 1, 1, 1, 12, 1, 81, 2, 2, 2, 1, 1, 6, 1, 392, 1, 10, 33, 2, 1, 2, 1, 24, 1, 2, 1, 3, 4, 1, 1, 4, 1, 1, 2, 2, 2, 1, 1, 1, 103, 2, 10, 5, 8, 1, 2, 4, 1, 4, 1, 1, 4, 9, 2, 60, 1, 3, 1, 1, 1, 1, 4, 7, 1, 1, 1, 1, 1, 1, 1, 3, 126, 1, 5, 1, 1, 14, 1, 7, 3, 2, 6, 1, 1, 1, 5, 4, 1, 1, 1, 8, 18, 1, 4, 1, 1, 7, 2, 1, 3, 7, 1, 1, 25, 1, 4, 1, 3, 4, 1, 1, 1, 4, 4, 2, 1, 3, 1, 1, 4, 2, 2, 1, 5, 2, 3, 90, 1, 7, 5, 11, 1, 3, 1, 1, 4, 1, 2, 1, 12, 1, 2, 2, 1, 12, 1, 1, 2, 23, 1, 12, 1, 2, 1, 14, 1, 1, 2, 6, 3, 1, 2, 26, 3, 2, 4, 2, 3, 1, 31, 2, 1, 1, 3, 1, 2, 7, 3, 1, 2, 3, 1, 2, 5, 25, 2, 1, 319, 2, 1, 9, 1, 5, 1, 2, 8, 1, 5, 1, 35, 1, 2, 3, 7, 1, 1, 2, 67, 1, 7, 5, 5, 2, 2, 3, 19, 1, 65, 6, 4, 2, 1, 112, 18, 1, 2, 2, 9, 6, 5, 2, 4, 1, 1, 1, 4, 32, 11, 3, 2, 10, 1, 5, 2, 4210, 1, 2, 1, 6, 1, 6, 3, 3, 2, 1, 1, 7, 1, 8, 1, 6, 1, 12, 2, 1, 22, 7, 10, 2, 1, 2, 7, 2, 3, 10, 4, 1, 1, 1, 1, 3, 3, 2, 7, 2, 2, 1, 3, 238, 2, 3, 6, 1, 14, 3, 1, 6, 1, 2, 1, 44, 1, 12, 12, 1, 1, 1, 2, 4, 1, 6, 1, 1, 2, 30, 4, 5, 2, 1, 2, 1, 4, 1, 38, 1, 4, 2, 2, 2, 1, 1, 4, 1, 7, 1, 9, 1, 3, 1, 5, 2, 3, 4, 1, 23, 1, 1, 2, 6, 2, 5, 5, 1, 25, 1, 1, 1, 1, 3, 1, 3, 4, 32, 4, 6, 1, 3, 2, 2, 1, 4, 1, 3, 1, 4, 2, 1, 1, 1, 3, 2, 665, 5, 3, 2, 36, 2, 2, 2, 1, 4, 1, 3, 2, 1, 2, 16, 1, 1, 3, 1, 3, 1, 21, 7, 1, 1, 35, 3, 1, 3, 1, 2, 1, 1, 1, 10, 1, 9, 1, 2, 1, 4, 3, 1, 13, 12, 1, 1, 25, 1, 2, 1, 3, 4, 2, 4, 1, 15, 1, 1, 1, 3, 1, 2, 9, 1, 2, 1, 1, 158, 2, 4, 5, 3, 2, 1, 17, 7, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 17, 19, 1, 1, 1, 14, 17, 2, 1, 23, 5, 7, 2, 1, 58, 3, 3, 1, 2, 1, 9, 3, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 8, 20, 1, 1, 1, 2, 73, 15, 2, 1, 3, 1, 11, 2, 4, 1, 15, 3, 59, 25, 5, 1, 2, 4, 2, 1, 35, 3, 7, 1, 3, 1, 2, 13, 7, 4, 1, 4, 2, 1, 6, 2, 2, 3, 1, 1, 1, 10, 1, 95, 4, 2, 2, 5, 4, 1, 1, 1, 2, 1, 1, 5, 2, 1, 1, 3, 1, 3, 1, 1, 1, 4, 2, 1, 2, 1, 1, 2, 1, 1, 5, 4, 1, 69, 1, 3, 1, 1, 1, 8, 5, 9, 2, 1, 2, 1, 2, 3, 1, 5, 3, 2, 1, 1, 1, 10, 66, 20, 2, 4, 1, 2, 2, 1, 2, 1, 14, 2, 1, 1, 7, 58, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 4, 1, 4, 2, 1, 2, 1, 42, 1, 7, 4, 1, 19, 1, 4, 5, 52, 4, 1, 1, 1, 1, 1, 1, 12, 1, 2, 8, 12, 1, 7, 1, 5, 1, 1, 5, 1, 1, 3, 4, 1, 5, 68, 1, 1, 21, 1, 2, 1, 1, 61, 1, 8, 1, 2, 3, 1, 3, 3, 43, 1, 1, 1, 15, 1, 12, 15, 3, 1, 1, 1, 25, 1, 1, 16, 1, 2, 2, 1, 1, 15, 13, 1, 2, 1, 1, 1, 1, 1, 3, 1, 11, 1, 1, 1, 3, 2, 2, 13, 1, 11, 5, 1, 7, 1, 1, 5, 3, 3, 1, 2, 5, 18, 1, 51, 1, 1, 1, 3, 2, 50, 5, 1, 1, 1, 2, 1, 7, 3, 1, 1, 1, 1, 3, 1, 1, 4, 2, 2, 2, 5, 10, 1, 163, 1, 2, 1, 2, 1, 3, 6, 2, 79, 1, 1, 5, 2, 1, 2, 3, 2, 5, 3, 2, 1, 14, 2, 1, 1, 1, 1, 8, 1, 71, 1, 1, 1, 1, 3, 1, 13, 3, 1, 2, 6, 1, 3, 8, 1, 5, 2, 3, 1, 2, 1, 3, 2, 3, 3, 1, 12, 2, 1, 4, 1, 4, 1, 1, 3, 12, 1, 1, 1, 1, 4, 2, 4, 1, 4, 2, 2, 4, 1, 5, 1, 6, 88, 1, 2, 17, 2, 8, 2, 13, 1, 11, 1, 1, 2, 1, 14, 2, 2, 1, 6, 2, 2, 3, 1, 3, 1, 2, 2, 1, 2, 2, 9, 2, 5, 1, 11, 4, 2, 1, 12, 8, 2, 2, 1, 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 123 4265247468 4959356886 9827741791 8990623858 0754283831 2005499426 3463932550 8352035470 4384392696 5041972204 6026752584 1783584994 4957820683 3067744809 2796448665 9077913686 1296574730 1985887066 7528762413 1778050421 8555964569 8683399349 8126360769 6274952754 7379430925 5354641896 7470461767 5734965516 8074597324 7479444093 5636533824 7973871621 1937286898 9948880721 4927196627 5302683518 8611235253 3125318356 1208359402 6189824267 2244092290 8009056365 0117525314 4576385420 4239077134 5278615976 4554003005 6637728955 5149605579 6727799056 4885240900 0556925075 5792912920 7966350220 9548936912 7310750285 6195450436 9095649567 2162402415 2127499502 5290028906 8693337464 3489228738 8214723895 4605802142 9251803149 5346645362 9797103584 3180716088 6303012458 0692164980 3989958148 0728052367 9702578808 2950607634 3872110840 6873023600 1889103945 7491400364 3754471545 0892074871 9121495537 0708768744 7259800309 4598686278 7906104140 8646263088 3603894152 5314044775 3703088589 9060582567 0884767037 0223491997 1825872584 1819863468 2797298384 2109158525 3448635439 3566984756 0043607607 3032886357, which is too large.

We also need to calculate d = floor(c4·v/F + 0.5) = 16 4600840672 8811358103 7362931580 2046773449 9915430907 2756307361 7710657530 1842489672 2018595810 6512524649 0570921843 0289956164 2429651170 7398520362 0882575457 6870783549 7004040052 6888819228 0311196201 2747682431 2128993133 3067130942 8919199336 8397727919 4640417892 8724259904 7931769014 1103663639 5021094262 0984627239 2478829176 2472429017 4486836827 1267020051 0604529082 0073956800 6691467700 4627929437 9311445523 5686886655 9912036952 6275014937 4202376402 7820488569 0874581454 3335572507 3278223953 9536132418 2766579413 5757546984 7520611583 2905128115 8548458122 4201381254 7449969648 4666113853 1963468119 8831054731 0946563431 7225060407 2973667929 1819034892 6681768146 8600975550 8375198843 5300401208 8963128358 6829533121 1367745932 1559490548 7664745225 8020989668 3421930815 0614191769 7653754557 7700280345 2717802469 0835659234 3921135836 1244066591 6323467872 0232078559 3758443698 6299443539 9771223070 3326174856 0753856400 4201669157 6312963065 9214540069 9280252189 6183216035 6576064500 9943569237 9685047799 9365611483 2212967672 8572874375 8854333274 8602467989 4492675889 4023570889 2616226882 5814806595 7275682360 8394412505 6795827065 9183366953 4700142147 8323410140 9862939448 3227823119 7861022416 2069068936 2940750642 1743058395 8382838598 1332450575 1315435459 0458265492 7408631126 1874220100 8610142419 0973078984 2894493737 6371396561 4178682342 6659854207 7454036603 6421104839 0161372649 0298354694 4271164249 3274983941 1916894438 1204548308 6234466953 9241589386 0447169594 8179318373 2088235592 0095807832 2346656971 8519758741 6412041231 6464359779 8690869341 2320303916 1327473109 6275103254 8047424891 8220588366 3005176016 6504498209 3727248817 2809375742 3769656980 1147707850 3698831487 5493506441 8611084636 5342627870 2051690098 5383678848 7090326689 9006443099 4942916083 9913836801 8846816673 7454245306 3038168499 8558466290 0410410141 8977419468 1136086620 0909753807 0897361292 2121073652 7914473577 4740765736 2583549873 8832950601 3689309255 2455322462 3643070393 0529348389 7189726589 4773608437 4144027694 1853922565 8723483095 2972402106 4999194349 1175762040 8528837667 4627409286 2289856228 2267728846 2905635112 4933088393 9986642608 7473091697

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.