Primality Certificate for (3048^3121-1)/3047

Andy Steward10,871 digits14 February 2002
Originally by David Broadhurst, Bouk de Water and Andy Steward 2002

This certificate uses a theorem of Konyagin, Pomerance and Broadhurst to prove an integer N prime by making use of a partial prime factorization of N2-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 29.971874% factorization of N-1:

From Factorisation
30482 · 2 · 2 · 3 · 127
Φ23049
Φ39293353
Φ45 · 1858061
Φ586338074552361
Φ67 · 19 · 69829
Φ841 · 15073 · 139661369
Φ1071 · 2281 · 3691 · 144341
Φ1213 · 37 · 27337 · 6563929
Φ1379 · p40
Φ1531 · 571 · 57601 · 128971 · 56631400304911
Φ1617 · 17 · 25776376024271786271359713
Φ205 · 1489874373834166325460012941
Φ241321 · 105087121 · 53662066063585801
Φ26p42
Φ30181 · 541 · 811 · 380994961 · 246289692811
Φ39788659 · p78
Φ40881 · 16472837752801 · 29776658230218288601 · 128415993811951394201
Φ4811209729 · p49
Φ5253 · 157 · 542136026776678648133 · p59
Φ6061 · 61 · 15541 · 38679640022391184021 · 24809521912445161295574237601
Φ65131 · 600601 · 50052309470671 · 3236669770132241 · p131
Φ782693107 · p78
Φ80241 · p110
Φ104313 · 2129713 · 14060237754774433 · c143
Φ1208446249561 · 803656329361 · 66685164608334072481 · p70
Φ1303121 · 41341 · 241220564261 · c148
Φ15613 · 1196823111380209 · c152
Φ1952341 · 24811911541 · 1195644606391 · 330449604855155058565111 · c286
Φ2083329 · 182012273 · p323
Φ24011595165601 · 15440046371090401 · c197
Φ260524475901 · 158435164681 · 4213957752421 · 32164311518339128441 · 2205342979280185271280001 · c259
Φ312937 · 127593306361 · 3788963832963818747541558503807833 · c287
Φ3905714014044778998271 · 3502379080536953654125531 · c292
Φ520521 · 32143116721 · 309050199121 · 2481287326961 · c632
Φ624c669
Φ7803885181 · 41055388810021 · 1245648411025201 · 955620036715534313347381 · c610
Φ1040801841 · 607555841826124226881 · 1131791122008005231253995281 · c1285
Φ156092041 · p1333
Φ312037911121 · 2743734241 · c2659

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

790 2496930685 5181765756 4814327684 4582536136 6923202413 0897373301 6474512041 2273104144 9418084273 4587206556 5030071988 9014040395 2802263682 9612757532 8522441091 1643730937 2145683783 4905906844 8128147453 3233763766 9795998956 4291279690 9964368112 3574023232 1724124034 5845411843 3149056271 5471374468 5573671818 8319404329 3115087948 9672690347 9585377100 7081329374 2342595963 7413577804 4593975584 7567327320 3175671736 2168772787 8521878349 6850017911 7010355511 3599099245 7739563859 5769867228 6154159313 6600057943 0120233857 7459839495 1784423133 3338220335 8255096346 6734362468 3004880329 4298241123 5814192231 5197277498 7207316839 2606309110 9952662566 8130617920 9513821493 2746551817 7610761015 0001117395 3563733781 0758630431 3681376381 8060745860 6605369782 4342458674 3737405071 9912273369 7887139841 7355244719 0315858250 5125726700 6252953314 1846380757 6358316388 0012188749 9022574742 5913973300 2742799115 2165358090 0716649508 0921231722 8037031301 0429263885 8965738450 1105266446 8681276908 9373620294 8214700370 0103550444 0013356012 6203742858 4039090445 7206936124 4116537926 1852322659 8791031507 1679187963 5632992996 9591536279 4970071377 2341363304 0496262198 9264396577 8317021997 5101831908 2278164404 6626040068 7943719663 7246220053 7427873262 4318061538 9815147240 9979793778 8697066076 5573477153 8785861253 3979123279 7026235980 1412942262 3557418590 7107849811 8598922247 5667822962 4410002802 3360637160 0487722681
481 9721541575 8302517235 4022683947 4699313732 0721073453 4125763598 8699825481 5228970266 0069099934 5187409900 5241301842 0905469958 9307535031 7367482290 9783206129 1168128684 4085413718 9420834854 8980594673 1923284988 9651249511 9800865959 1473275152 0975721352 6526306515 2011623717 4492979447 1669784509 8659282046 6668610516 0881543268 3770280593
1 3402809359 0268767413 6031311084 5500815128 0841765791 1043856824 6007452906 3312559746 5143976895 1265260303 7841045456 9270309265 5528402581
1277796701 7231810354 4334617710 3652269021 1829137730 5799578396 6025416243 5429763273 4602082633 5587130712 8725499921
52399522 7922277577 1865372578 7949285280 4612346146 4904280545 6446312184 1091255299
15354934 1029065048 3781428160 7051793751 0362952596 9884710848 2692690090 9094699443
6803232278 0681502770 4005952005 7340878359 2979144925 9906092160 3417685561
916379292 1731603892 4262075091 9830697125 0912729768 7525453261
495044556 3104906002 0590499910 7243031192 7217189889
64 2742607506 6995949985 6921397803 2990026841
8141322696 6882846591 1952214221 5096760839
3788 9638329638 1874754155 8503807833
248095219 1244516129 5574237601
14898743 7383416632 5460012941
11317911 2200800523 1253995281
257763 7602427178 6271359713
35023 7908053695 3654125531
22053 4297928018 5271280001
9556 2003671553 4313347381
3304 4960485515 5058565111
6 0755584182 6124226881
5 4213602677 6678648133
1 2841599381 1951394201
6668516460 8334072481
3867964002 2391184021
3216431151 8339128441
2977665823 0218288601
571401404 4778998271
5366206 6063585801
1544004 6371090401
1406023 7754774433
323666 9770132241
124564 8411025201
119682 3111380209
8633 8074552361
5663 1400304911
5005 2309470671
4105 5388810021
1647 2837752801
421 3957752421
248 1287326961
119 5644606391
80 3656329361
30 9050199121
24 6289692811
24 1220564261
15 8435164681
12 7593306361
3 2143116721
2 4811911541
1 1595165601
8446249561
2743734241
524475901
380994961
182012273
139661369
105087121
37911121
11209729
9293353
6563929
3885181
2693107
2129713
1858061
801841
788659
600601
144341
128971
92041
69829
57601
41341
27337
15541
15073
3691
3329
3121
3049
2341
2281
1321
937
881
811
571
541
521
313
241
181
157
131
127
79
71
612
53
41
37
31
19
172
132
7
52
3
23

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) = 29.971874%

Factorizing N+1

We have the following 0.153629% factorization of N+1:

2
67
3863
5309
18221341

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

We set R' = (N+1)/F'. Note that GCD(F',R')=1 and Log(F')/Log(N) = 0.153629%

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

Finding a Lucas Series

Next, we find a discriminant Δ that is not a perfect square (mod N). We choose a pair of integers P and Q such that P2-4·Q = Δ. Let α and β be the roots of the quadratic f(x) = x2-P·x+Q. With α and β, we build Lucas series U and V. We then verify that UN+1≡0 (mod N) and, for each prime factor q of F', GCD(U(N+1)/q,N)=1. In this case, P=4, Q=-205375, Δ=236 suffice.

Ensure GCD(F,F')=1

As they stand, both F and F' have a common factor of 2. One of them has a single instance and the other has at least two. To ensure that GCD(F,F')=1 while making the minimum reduction in our factorizations, we divide out the single 2 from F' and multiply it into R' .

Ensure (F^(1/3))/6 > F'-1

log10(LHS)=1395.607318, log10(RHS)=16.398593

Ensure F'-1 > 3.N/(F^(10/3))

log10(LHS)=16.398593, log10(RHS)=0.477569

Express N in base F

Let N = c3·F3 + c2·F2 + c1·F + 1. Let c4 = c3·F+c2. Let t = c4 mod F'.

Square Check

We prove that (c1+t·F)2+4·t-4·c4 is not a perfect square. Here, it is not a perfect square: it is ≡ 3 (mod 63).

Continued Fraction

We approximate c1/F by a continued fraction u/v such that v is maximal while remaining less than F1/3 = 985325 0368231061 0478733402 2396561611 8562579592 2700199775 8880942628 3804878754 6081412970 5874867709 5755780649 8764821637 2839151383 7761749504 6426593385 5371697412 0115159317 5897230095 0327861761 8405305598 0423364911 6717055957 7121544389 6443067223 6575096664 9866220997 7407923331 0833267071 0404648015 1095283023 3143599429 2592976174 9724015619 0253315017 3303538608 7149834129 8660002279 5978216559 4222428090 1314851577 8144181898 4509035230 3175140680 7189789312 5188953910 2030834166 1117299199 2789271942 0353905140 6179707958 8575894376 7792207428 7072814172 6748390934 2471900170 5553561682 8891571604 3896502879 7238758601 9224124953 5982865729 2807033777 0079153172 6439645275 4484428263 9565900887 0411914242 0012689421 1804410392 9058757741 0068724467 7604467185 2325047790 2877506453 8565780993 0529191994 6146358278 6327445899 3356114328 3031792073 6047784075 5461584857 5621226138 6071234765 2948395237 3440964066 5290674981 6794938359 2864729646 5364745302 6806873058 2177179575 2605295647 5329788523 1122850491 7803398265 4514078958 1021906747 2471625881 1463433146 1611307596 8051294562 2016343138 0165004099 9578965366 2236285087 7493338176 1400510225 5395966855.

With those constraints, the unique continued fraction is: {0, 1, 2, 1, 1, 4, 1, 562, 1, 1, 7, 1, 4, 1, 25, 1, 1, 1, 6, 6, 1, 32, 1, 2, 1, 1, 8, 1, 3, 1, 42, 1, 7, 1, 1, 1, 11, 5, 3, 11, 1, 10, 2, 4, 2, 2, 1, 1, 8, 27, 2, 16, 3, 7, 3, 2, 1, 1, 3, 3, 1, 3, 3, 1, 3, 3, 9, 1, 1, 1, 10, 1, 4, 7, 2, 156, 1, 14, 1, 1, 2, 743, 1, 1, 1, 3, 2, 1, 9, 1, 14, 1, 1, 3, 86, 10, 1, 10, 1, 1, 2, 1, 2, 4, 1, 1, 4, 3, 1, 1, 3, 1, 2, 2, 1, 1, 1, 2, 1, 2, 3, 1, 4, 28, 1, 8, 3, 1, 2, 1, 2, 1, 9, 1, 4, 5, 15, 1, 1, 4, 1, 13, 5, 27, 5, 1, 6, 1, 7, 15, 2, 3, 4, 2, 1, 1, 1, 15, 4, 92, 3, 1, 89, 1, 1, 16, 1, 3, 1, 1, 119, 1, 1, 1, 1, 2, 1, 1, 2, 1, 4, 2, 1, 7, 1, 6, 2, 106, 40, 1, 4, 8, 2, 2, 1, 4, 2, 7, 5, 16, 5, 1, 2, 3, 4, 1, 1, 1, 1, 2, 1, 1, 2, 1, 2, 11, 2, 1, 4, 3, 4, 2, 2, 10, 1, 4, 7, 13, 4, 9, 1, 3, 17, 4, 1, 1, 146, 1, 3, 1, 11, 5, 1, 5, 1, 1, 3, 1, 28, 32, 1, 1, 23, 1, 1, 40, 12, 64, 5, 1, 1, 2, 29, 1, 1, 1, 3, 1, 1, 2, 1, 1, 33, 1, 8, 1, 1, 7, 1, 2, 1, 5, 3, 6, 3, 9, 2, 1, 8, 1, 1, 1, 1, 1, 3, 15, 9, 1, 2, 56, 1, 6, 13, 2, 1, 2, 1, 1, 5705, 1, 1, 2, 2, 11, 1, 1, 2, 8, 1, 279, 21, 1, 95, 2, 1, 1, 7, 15, 100, 1, 1, 1, 1, 8, 175, 2, 11, 2, 1, 1, 1, 1, 1, 1, 1, 11, 3, 1, 1, 33, 3, 1, 2, 4, 2, 1, 1, 7, 2, 1, 3, 1, 3, 2, 39, 1, 3, 2, 159, 1, 1, 17, 16, 3, 2, 3, 78, 4, 1, 2, 20, 1, 1, 5, 2, 1, 4, 1, 1, 8, 1, 9, 1, 36, 1, 1, 3, 11, 1, 3, 9, 2, 1, 3, 3, 5, 21, 1, 3, 3, 1, 101, 1, 10, 1, 9, 49, 1, 1, 3, 1, 4, 1, 2, 2, 2, 2, 1, 1, 2, 6, 11, 15, 1, 20, 2, 3, 8, 1, 9, 6, 1, 4, 2, 1, 290, 1, 12, 1, 2, 1, 2, 1, 2, 4, 12, 1, 2, 1, 1, 1, 4, 2, 3, 7, 2, 5, 1, 1, 4, 5, 2, 1, 1, 1, 2, 8, 3, 2, 1, 7, 1, 1, 2, 4, 2, 4, 2, 2, 10, 1, 14, 1, 5, 1, 4, 87, 8, 4, 17, 1, 9, 7, 3, 3, 2, 5, 1, 3, 3, 1, 2, 2, 6, 1, 1, 1, 1, 1, 2, 6, 2, 2, 1, 2, 6, 1, 1, 1, 3, 6, 2, 1, 2, 1, 2, 4, 2, 1, 1, 1, 2, 2, 15, 2, 2, 2, 4, 1, 1, 1, 1, 355, 1, 1, 1, 31, 3, 2, 1, 9, 2, 5, 1, 10, 1, 2, 5, 1, 6, 9, 1, 1, 1, 1, 1, 1, 51, 3, 1, 5, 1, 11, 21, 1, 3, 2, 2, 1, 1, 1, 2, 1, 600, 1, 1, 11, 5, 1, 4, 4, 3, 1, 2, 1, 32, 1, 1, 2, 2, 3, 1, 1, 2, 2, 1, 8, 1, 1, 2, 1, 7, 1, 10, 6, 2, 1, 1, 24, 33, 8, 3, 6, 4, 1, 2, 1, 1, 1, 8, 1, 3, 1, 1, 2, 1, 7, 1, 38, 372, 1, 1, 5, 2, 2, 1, 1, 1, 2, 2, 10, 1, 4, 4, 1, 10, 3, 4, 3, 2, 4, 1, 1, 1, 38, 1, 2, 1, 1, 1, 1, 17, 1, 1, 1, 2, 4, 2, 3, 1, 1, 2, 3, 102, 2, 3, 68, 1, 1, 12, 1, 8, 5, 1, 2, 1, 2, 1, 2, 2, 6, 1, 18, 16, 1, 2, 5, 2, 1, 1, 2, 9, 4, 1, 1, 1, 4, 1, 1, 1, 11, 1, 1, 2, 2, 10, 5, 2, 1, 6, 13, 1, 1, 1, 5, 1, 2, 2, 3, 2, 3, 2, 2, 1, 4, 1, 2, 1, 4, 1, 5, 3, 2, 1, 1, 30, 1, 2, 7, 5, 1, 5, 1, 1, 2, 1, 1, 11, 155, 1, 4, 7, 1, 4, 1, 2, 96, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 2, 1, 1, 1, 6, 2, 3, 5, 3, 1, 1, 2, 7, 1, 15, 1, 1, 14, 1, 32, 1, 4, 3, 1, 364, 2, 6, 1, 15, 1, 1, 7, 1, 13, 1, 1, 7, 1, 1, 2, 4, 1, 1, 5, 1, 14, 1, 2, 1, 1, 4, 2, 12, 5, 6, 11, 5, 1, 3, 1, 1, 1, 3, 3, 27, 1, 4, 1, 5, 1, 1, 2, 1, 1, 4, 3, 1, 7, 5, 1, 98, 46, 1, 2, 6, 2, 3, 1, 7, 1, 9, 1, 1, 3, 1, 2, 5, 1, 28, 3, 1, 1, 1, 4, 5, 4, 8, 1, 78, 2, 2, 1, 2, 11, 1, 6, 2, 2, 1, 2, 1, 3, 1, 1, 32, 2, 2, 5, 2, 3, 4, 2, 2, 1, 1, 6, 1, 1, 15, 1, 3, 2, 140, 2, 2, 1, 1, 5, 9, 1, 4, 11, 3, 2, 2, 3, 12, 1, 1, 2, 2, 7, 1, 6, 24, 1, 1, 7, 16, 1, 4, 1, 5, 1, 3, 110, 1, 3, 1, 4, 1, 7, 8, 3, 3, 4, 2, 1, 1, 1, 2, 1, 5, 1, 2, 1, 10, 1, 3, 5, 1, 2, 18, 1, 8, 1, 6, 1, 3, 12, 1, 2, 2, 1, 1, 6, 2, 1, 2, 3, 1, 4, 2, 1, 2, 1, 1, 1, 2, 1, 10, 13, 1, 1, 5, 1, 1, 4, 8, 3, 9, 7, 33, 3, 2, 10, 1, 34, 2, 1, 9, 1, 1, 1, 2, 12, 2, 1, 2, 1, 1, 2, 4, 1, 12, 2, 1, 2, 4, 2, 7, 8, 4, 2, 2, 1, 2, 1, 1, 3, 157, 2, 15, 1, 1, 43, 1, 1, 1, 1, 3, 2, 1, 1, 1, 1, 6, 1, 4, 3, 1, 1, 1, 4, 1, 2, 1, 1, 4, 1, 73, 3, 2, 2, 1, 1, 4, 6, 10, 4, 2, 2, 520, 1, 2, 1, 3, 19, 1, 6, 1, 1, 1, 3, 14, 47, 1, 1, 1, 1, 2, 23, 1, 1, 5, 1, 1, 7, 1, 1, 1, 1, 1, 4, 1, 2, 2, 2, 2, 2, 2, 1, 3, 4, 1, 2, 4, 1, 5, 2, 1, 5, 13, 1, 2, 4, 161, 2, 2, 1, 14, 4, 1, 1, 2, 16, 1, 20, 2, 2, 2, 1, 6, 1, 31, 2, 1, 9, 1, 14, 50, 1, 2, 1, 1, 191, 4, 1, 1, 1, 2, 1, 2, 2, 1, 1, 2, 3, 80, 1, 1, 1, 2, 4, 2, 6, 1, 1, 17, 1, 2, 2, 1, 2, 9, 1, 4, 4, 2, 397, 1, 1, 5, 5, 2, 5, 76, 1, 1, 3158, 1, 4, 5, 2, 1, 19, 1, 4, 1, 2, 1, 1, 3, 3, 2, 3, 1, 4, 4, 4, 5, 2, 7, 1, 20, 52, 1, 1, 1, 6, 1557, 1, 1, 5, 1, 15, 18, 3, 1, 1, 9, 1, 1, 1, 18, 9, 1, 1, 1, 1, 1, 1, 32, 1, 1, 5, 1, 1, 1, 1, 5, 1, 1, 5, 5, 1, 4, 6, 1, 1, 1, 1, 1, 1, 7, 2, 16, 1, 1, 38, 1, 5, 1, 14, 1, 7, 6, 1, 10, 1, 1, 3, 1, 1, 1, 3, 2, 7, 2, 6, 1, 1, 4, 1, 6, 1, 3, 78, 1, 1, 1, 55, 1, 1, 3, 1, 2, 12, 17, 1, 16, 2, 2, 1, 98, 3, 1, 8, 1, 4, 1, 10, 2, 367, 1, 3, 2, 1, 5, 1, 2, 5, 1, 16, 1, 31, 1, 1, 1, 2, 1, 2, 1, 2, 10, 2, 13, 1, 5, 1, 3, 3, 1, 4, 5, 1, 1, 2, 3, 8, 1, 2, 1, 5, 3, 2, 1, 1, 1, 7, 5, 13, 2, 1, 4, 1, 43, 1, 3, 3, 2, 2, 3, 4, 5, 1, 7, 1, 1, 86, 1, 1, 32, 2, 8, 1, 1, 2, 1, 1, 7, 5, 13, 8, 1, 34, 2, 2, 1, 148, 1, 1, 2, 1, 1, 1, 9, 3, 2, 1, 1, 5, 1, 2, 1, 3, 5, 4, 1, 7, 6, 1, 6, 19, 1, 1, 1, 1, 8, 5, 1, 3, 8, 1, 5, 1, 7, 1, 1, 3, 2, 4, 4, 1, 1, 4, 23, 1, 1, 3, 2, 1, 2, 1, 2, 2, 1, 1, 3, 2, 1, 75, 3, 1, 1, 8, 8, 7, 4, 7, 13, 2, 8, 1, 34, 65, 1, 3, 1, 1, 1, 1, 8, 1, 1, 38, 1, 2, 3, 1, 23, 1, 4, 8, 8, 1, 2, 1, 2, 1, 1, 6, 3, 6, 2, 1, 6, 12, 6, 5, 1, 1, 6, 4, 9, 2, 4, 13, 4, 10, 6, 3, 3, 17, 1, 3, 1, 6, 1, 1, 24, 1, 4, 7, 1, 3, 9, 3, 1, 1, 3, 23, 1, 2, 6, 1, 2, 4, 1, 2, 5, 1, 2, 16, 7, 1, 4, 1, 17, 9, 3, 2, 2, 1, 1, 6, 4, 6, 3, 1, 5, 4, 7, 1, 2, 13, 1, 2, 2, 3, 1, 1, 1, 1, 25, 3, 3, 13, 1, 1, 1, 14, 22, 1, 31, 3, 1, 7, 2, 3, 29, 3, 19, 5, 1, 3, 2, 3, 1, 1, 1, 1, 1, 2, 4, 7, 5, 2, 1, 1, 1, 1, 27, 2, 1, 5, 1, 6, 1, 32, 1, 2, 15, 2, 2, 17, 1, 1, 2, 1, 7, 5, 4, 1, 1, 10, 2, 20, 1, 1, 1, 1, 17, 7, 1, 1, 1, 13, 748, 1, 1, 4, 5, 1, 1, 6, 1, 2, 7, 24, 1, 1, 2, 8, 2, 3, 14, 1, 1, 3, 1, 9, 48, 3, 1, 1, 1, 1, 2, 1, 1, 1, 4, 3, 2, 4, 1, 1, 921, 1, 2, 3, 1, 18, 1, 1, 54, 1, 1, 1, 9, 17, 19, 1, 1, 2, 1, 9, 2, 4, 31, 1, 1, 3, 9, 1, 2, 4, 2, 2, 2, 1, 3, 26, 1, 1, 2, 1, 96, 2, 2, 7, 1, 1, 23, 2, 11, 1, 9, 1, 4, 2, 1, 3, 1, 6, 3, 2, 1, 4, 1, 6, 1, 2, 1, 2, 3, 1, 1, 7, 1, 1, 11, 1, 3, 4, 1, 1, 1, 4, 1, 5, 3, 17, 1, 2, 2, 7, 1, 10, 3, 1, 1, 4, 1, 4, 3, 1, 1, 3, 2, 1, 8, 1, 20, 1, 11, 2, 1, 1, 2, 1, 4, 1, 2, 3, 5, 2, 2, 1, 5, 62, 1, 1, 5, 1, 2, 1, 27, 3, 10, 4, 4, 2, 1, 1, 2, 3, 1, 2, 1, 2, 5, 2, 1, 1, 9, 5, 6, 2, 1, 4, 4, 1, 2, 1, 5, 2, 111, 4, 39, 1, 1, 1, 2, 20, 27, 8, 1, 1, 1, 1, 12, 93, 5, 43, 1, 2, 11, 1, 13, 3, 3, 2, 2, 21, 1, 3, 2, 1, 7, 4, 1, 2, 1, 3, 1, 3, 7, 1, 8, 1, 1, 6, 1, 1, 39, 5, 6, 3, 1, 1, 6, 2, 1, 1, 2, 1, 2, 1, 1, 8, 1, 3, 2, 3, 14, 1, 7, 137, 1, 4, 13, 1, 1, 2, 1, 28, 1, 4, 2, 1, 1, 5, 2, 4, 6, 20, 1, 1, 1, 4, 1, 4, 9, 1, 1, 10, 16, 1, 1, 3, 3, 3, 4, 2, 16, 10, 1, 1, 1, 10, 4, 2, 59, 1, 116, 3, 3, 6, 2, 6, 1, 2, 8, 1, 1, 1, 6, 1, 9, 125, 4, 5, 10, 1, 2, 1, 3, 3, 4, 1, 1, 5, 3, 2, 2, 5, 1, 55, 10, 1, 3, 1, 2, 57, 2, 1, 1, 6}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 5105310203 8379095492 8514111750 6433778506 6021227587 2535894454 2500362199 2742490467 8103053542 9351636168 8376678422 8194991617 1950356751 8186964781 1923362337 8231235124 5668349241 3105477968 9192956320 5545569701 3722520165 1213154677 8177617301 2497777102 8898700205 4464326849 6732153792 5888013095 5419640772 4132713657 2423745986 0661754329 3865208741 6661641566 4473539901 9269347575 4152555150 6836361972 5616813652 4496063900 7235910406 1835583468 1544475049 7449760565 6171237936 7824934022 9417980852 0584382572 1138597229 0468101678 4442664342 7665000114 1240075608 2324221791 5360634091 6260551142 5189714306 4927204813 7704131609 8835641898 8208167761 2353729562 1295566328 6184436455 3187570156 9966719921 0292299824 5554252758 2705860520 2051749018 6245540203 6481183023 1510117582 4875988384 0857285554 9547489992 5723804550 8896739525 0426503439 1516295593 0629258876 6195565687 3632178896 7721795070 3835484273 5247170431 2231184961 7298357569 4147696548 2918099766 7181338808 2478927150 3511668376 3709503603 6932764609 3752186645 2416831926 7171064464 8211269524 2332482607 9328430915 2069148100 2893290971 2049138738 0051962854 6899331573 1266354655 1462267985 2470720167 7421901520 3645497504 2781949567 5588896448 7978780437 1131748352 6497315925 9698399810 7568173835 2655576984 7922109094 7894415166 8322631612 8003228067 5280198322 1014968274 7099733809 3682390343 5911518729 3074409612 7265899610 9471169087 2986585341 0838735686 9669878227 9255868674 6640499159 4549658523 8394206128 5625823153 8494188239 9351411311 3025421713 3409225813 5063567099 8232938767 8260830694 8952539855 4563733434 2712118137 7261479716 6369991663 2844591104 8900962470 4893205501 0138481813 2024361388 4836160694 4086702978 0740289365 0468224958 6903232831 9796054595 6206287116 4702242950 4306118419 9986709393 0310028897 0633542940 5559164181 9722198048 7172142176 9092986451 2002162545 3744707568 2361545348 0249197264 0521581780 7926220635 9825714541 2116298424 6486980501 9310395913 5468383849 7344987156 9939447573 3394440929 1646589592 6564616721 8361153603 3865334106 9273407052 5860339691 3774669590 4774171951 4800666246 3970227500 2922321700 2151469802 2819496344 3732333884 4036279760 9867786514 7274560912 9252702124 5533588735 8330621519 9936499808 6956021077 1921795167 3766464895 6793125064 8729194178 1506281271 0251087547 3342367046 5707687589 3932928606 8330522125

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.