Primality Certificate for (9751^2161-1)/9750

Andy Steward8,617 digits30 June 2001
Originally by David Broadhurst 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.515001% factorization of N-1:

From Factorisation
97517 · 7 · 199
Φ22 · 2 · 2 · 23 · 53
Φ33 · 4969 · 6379
Φ42 · 47541001
Φ55 · 11 · 701 · 3361 · 7901 · 8831
Φ695072251
Φ82 · 19409561 · 232890041
Φ93 · 19 · 127 · 1891189 · 3958687 · 15860989
Φ1061 · 151 · 164581 · 5963011
Φ121213 · 7453080642277
Φ158197551026101 · 9969298094170947901
Φ162 · 17 · 17 · 241 · 1409 · 416425559001428311638161
Φ1887013 · 1890523 · 5225511777949
Φ20p32
Φ2453353 · 926353 · 1653704448568165101289
Φ273 · 5347 · 187597 · 690607 · 23737516382719 · 6319399803477589 · 2037442693069561431732787297
Φ3031 · 2221 · 42335131 · 28043139173708977921
Φ3637 · 1297 · 14605318801 · 131324314169773 · 8027715903021612433
Φ4041 · 680054294341361 · 15594434721478201 · p32
Φ4516687987201 · 8821120327098325381 · p67
Φ48p64
Φ54271 · 757 · 2053 · 141697 · 426709 · 2608471 · 14982148430806907683 · 638231344970827973459572549
Φ601662781 · 58472095021 · p47
Φ7273 · 39097 · 31174501220470801 · 2721380621295698646591817 · p49
Φ801601 · 3277210143575041 · 624899502755752646854268801231629601 · p74
Φ9022051 · 156421 · 1229372191 · p78
Φ108109 · 15013 · 54217 · 2842346917 · 1001913658021 · 706951650490381981 · 12114690234177037095181 · p72
Φ1201201 · p125
Φ135811 · 2474551 · 93596851 · 264665341 · 577670403241 · 20705291062152933323069821 · c225
Φ14415990102224499370321 · 52000014773752256017 · c153
Φ180181 · 90703981 · c182
Φ216433 · 8236729 · 11209969 · 146068921 · 4985053417 · c253
Φ240105601 · 41895528226018219201 · p231
Φ270524610001 · 28134687961 · 1844063816103361 · p253
Φ360p383
Φ432515377 · 2945693311135818769 · c551
Φ540541 · 7523281 · 350705718918361 · c551
Φ7202161 · 178561 · 37890001 · 1640300780939761 · 623911371502385521 · 5588735368166552003041 · c695
Φ10807561 · 185193860152810884481 · c1125
Φ2160386641 · 38974584661446241 · 42641337696056178481 · c2256

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 :

888 6260043106 9547535830 2698659710 8018539792 1651278324 7037839098 8000735081 7004261677 6262496426 2114724872 4999342379 9499861632 2581909247 2634228863 7486089756 7503034161 7124475760 7450748660 8690517029 3392563162 1114234609 5029492165 9639120345 4996278624 9070541134 6674005718 8407168450 0002117186 7263196942 9501634889 2377085829 7042357176 5471225727 2441658540 8631052763 5008278429 6304404659 5641968001
597 9772324991 1539668833 5006423706 5587609847 1217799841 4389610976 1178294102 7800875534 5007456220 8807984584 2639922682 9227341192 8227372682 3717587397 4443556649 0310369025 8038920176 3752448507 0964386735 2899261859 9009421846 0948464657 4587307254 1801922603 5397167481
4 5010134507 6578595676 6895497697 5532976856 3392367546 0874003103 2064679081 4195960832 1302818610 3331330603 4772382356 5626904005 9683283376 3604647555 3093552159 5498302331 4235336927 9132765404 5403642439 8033590092 2338524818 1075272527 5648907201
37156 0893547486 4331283002 1239779911 3285115668 5345381123 1784191089 3623154459 1751536720 7026058200 6518877253 1675254498 3062894801
12875764 4274386583 6021749745 9053629212 6099112805 5569374164 0957612324 6167566641
1361 0284534795 0295646711 0502503427 0677271393 6381346751 5758650092 7203715361
18 6436421688 0610331675 1141944722 0264774919 6154770392 7504100403 7804316217
3708959 1516100000 8673396942 3330361876 9319550849 1834052003 4347947821
6680 1544379642 2386894325 3052678264 2654135362 2046938804 4308328001
225488635 4889389826 4187789905 6093371627 7901960713
6870728 6910500029 2013136863 4154862828 3109574601
624899 5027557526 4685426880 1231629601
81 7322108929 5628901245 1086164001
15 3634737294 6012150251 9235294001
20374426 9306956143 1732787297
6382313 4497082797 3459572549
207052 9106215293 3323069821
27213 8062129569 8646591817
4164 2555900142 8311638161
121 1469023417 7037095181
55 8873536816 6552003041
16 5370444856 8165101289
1 8519386015 2810884481
5200001477 3752256017
4264133769 6056178481
4189552822 6018219201
2804313917 3708977921
1599010222 4499370321
1498214843 0806907683
996929809 4170947901
882112032 7098325381
802771590 3021612433
294569331 1135818769
70695165 0490381981
62391137 1502385521
3897458 4661446241
3117450 1220470801
1559443 4721478201
631939 9803477589
327721 0143575041
184406 3816103361
164030 0780939761
68005 4294341361
35070 5718918361
13132 4314169773
2373 7516382719
819 7551026101
745 3080642277
522 5511777949
100 1913658021
57 7670403241
5 8472095021
2 8134687961
1 6687987201
1 4605318801
4985053417
2842346917
1229372191
524610001
264665341
232890041
146068921
95072251
93596851
90703981
47541001
42335131
37890001
19409561
15860989
11209969
8236729
7523281
5963011
3958687
2608471
2474551
1891189
1890523
1662781
926353
690607
515377
426709
386641
187597
178561
164581
156421
141697
105601
87013
54217
53353
39097
22051
15013
8831
7901
7561
6379
5347
4969
3361
2221
2161
2053
1601
1409
1297
1213
1201
811
757
701
541
433
271
241
199
181
151
127
109
73
61
53
41
37
31
23
19
172
11
72
5
33
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.515001%

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 = 71 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 = 2 4164216012 8623747637 8281150737 2832135663 0076065267 3624566441 6161923578 5441807834 6437547506 9447237569 1907166980 0936825387 1148464014 4613448702 3004504413 0602205028 6837889464 0696713736 3467251714 8616060827 3892709280 2009046836 6908322212 9781256668 7550160638 0629972845 4176138349 8847243178 2810320303 6797564816 4548019259 1721740414 3247054552 6933208461 7155327816 3201641627 6600400781 6280091836 4628859133 2714635564 0581485558 3908003168 6578111947 2999283014 2834926527 4622511990 6803162282 3065597511 0106853449 8244841571 1732372135 2829704456 4497449081 9706458104 7758196257 3076887997 5437863393 9022866254 7164769557 4860172007 6920833973 8659517563 9498139446 5817319161 8028692882 6798086208 8911671745 2667669635 3051008182 7259327553 3207369932 1620732397 1675435186 1513322607 1644189687 2111643008 1864439131 3528202916 7745493917 2435737457 3402146301 2694161752 2364330559 2420461241 1972889336 3204745634 1910174026 9738359822 1473061396 7886780531 5794715505 3219269036 6627610239 0671119031 5754730654.

With those constraints, the unique continued fraction is: {0, 1, 2, 2, 7, 1, 1, 4, 1, 1, 197, 47, 2, 17, 2, 5, 3, 2, 1, 4, 3, 1, 1, 1, 7, 2, 3, 1, 3, 2, 9, 1, 1, 23, 3, 12, 3, 2, 1, 3, 3, 5, 1, 1, 2, 1, 6, 1, 4, 1, 1, 1, 1, 2, 274, 1, 1, 6, 1, 5, 6, 2, 1, 27, 1, 17, 5, 16, 3, 3, 4, 1, 27, 1, 18, 1, 3, 1, 368, 13, 1, 1, 199, 1, 8, 2, 1, 2, 1, 1, 7, 1, 1, 1, 3, 3, 17, 4, 1, 7, 1, 1, 5, 1, 2, 1, 1, 1, 12, 9, 12, 2, 2, 1, 21, 1, 1, 1, 1, 2, 2, 3, 1, 2, 7, 2, 2, 28, 1, 8, 1, 1, 4, 3, 1, 4, 1, 9, 1, 5, 1, 1, 3, 1, 1, 2, 1, 39, 1, 4, 3, 3, 2, 6, 4, 1, 2, 2, 7, 1, 1, 1, 1, 4, 5, 1, 1, 3, 2, 3, 1, 11, 1, 1, 638, 1, 3, 1, 2, 14, 1, 10, 1, 9, 38, 1, 1, 1, 113, 2, 2, 9, 3, 2, 1, 1, 6, 2, 4, 7, 1, 18, 1, 1, 11, 9, 1, 20, 2, 3, 16, 1, 4, 1, 53, 4, 7, 2, 1, 6, 1, 1, 9, 108, 3, 1, 1, 3, 31, 15, 2, 16, 1, 2, 3, 1, 7, 6, 1, 1, 8, 4, 1, 1, 1, 1, 1, 1, 2, 91, 1, 2, 4, 2, 1, 1, 1, 1, 8, 20, 2, 1, 5, 1, 1, 2, 3, 1, 5, 4, 19, 1, 10, 2, 3, 1, 6, 1, 1, 1, 1, 6, 1, 2, 1, 1, 1, 1, 8, 388, 2, 7, 2, 1, 7, 1, 15, 1, 6, 1, 2, 2, 19, 6, 7, 1, 15, 2, 2, 12, 5, 2, 11, 6, 1, 1, 1, 120, 2, 2, 2, 17, 3, 1, 3, 2, 48, 1, 1, 1, 6, 1, 5, 1, 3, 8, 195, 1, 5, 1, 4, 1, 115, 4, 1, 1, 1, 5, 2, 1, 2, 2, 1, 1, 5, 3, 14, 1, 2, 4, 2, 1, 6, 12, 3, 2, 1, 2, 1, 1, 3, 3, 1, 2, 1, 2, 3, 2, 1, 1, 1, 2, 4, 1, 1, 12, 2, 1, 2, 2, 5, 3, 3, 1, 1, 5, 8, 1, 2, 1, 113, 1, 12, 1, 2, 2, 1, 1, 3, 4, 3, 1, 7, 1, 1, 2, 2, 1, 6, 1, 1, 2, 2, 10, 1, 1, 4, 4, 1, 1, 14, 20, 1, 2, 64, 7, 14, 1, 7, 1, 1, 6, 1, 1, 1, 2, 2, 1, 9, 2, 1, 4, 6, 3, 156, 1, 11, 1, 5, 1, 2, 2, 10, 8, 2, 1, 1, 3, 22, 7, 79, 4, 1, 2, 2, 1, 1, 1, 1, 3, 1, 48, 4, 1, 3, 1, 7, 4, 2, 2, 1, 1, 8, 1, 12, 11, 1, 1, 1, 12, 3, 1, 3, 14, 10, 1, 1, 1, 1, 1, 3, 1, 2, 7, 1, 1, 2, 1, 6, 26, 7, 8, 15, 2, 4, 3, 13, 38, 1, 2, 4, 3, 1, 1, 114, 1, 15, 28, 5, 5, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 3, 3, 25, 6, 89, 1, 16, 1, 4, 5, 1, 4, 5, 9, 2, 2, 22, 1, 12, 1, 72, 3, 1, 7, 1, 5, 1, 28, 1, 6, 1, 2, 1, 2, 10, 2, 1, 1, 2, 1, 1, 12, 1, 418, 1, 1, 7, 3, 5, 1, 2, 11, 4, 1, 2, 1, 9, 1, 3, 2, 258, 4, 2, 2, 2, 21, 4, 2, 35, 1, 2, 6, 1, 2, 133, 1, 74, 1, 4, 1, 1, 1, 17, 1, 10, 1, 7, 1, 1, 2, 1, 4, 1, 1, 3, 2, 7, 1, 2, 1, 4, 1, 348, 2, 16, 2, 9, 5, 1, 11, 1, 6, 1, 13, 1, 2, 2, 1, 66, 3, 1, 1, 7, 1, 1, 3, 15, 2, 2, 1, 1, 1, 5, 26, 3, 2, 4, 2, 9, 44, 1, 55, 1, 1, 1, 1, 3, 2, 2, 1, 1, 1, 3, 2, 1, 3, 5, 1, 2, 3, 3, 2, 6, 3, 2, 12, 1, 1, 1, 1715, 1, 4, 1, 1, 6, 1, 2, 1, 2, 2, 2, 1, 2, 2, 17, 5, 30, 3, 2, 2, 1, 1, 1, 23, 1, 2, 11, 2, 2, 2, 3, 1, 8, 15, 1, 2, 2, 3, 4, 1, 6, 38, 3, 2, 2, 2, 5, 1, 2, 2, 7, 3, 1, 13, 1, 2, 1, 1, 6, 1, 1, 1, 1, 1, 10, 2, 3, 1, 1, 149, 2, 7, 1, 1, 1, 10, 1, 1, 1, 1, 4, 1, 11, 1, 3, 1, 6, 1, 5, 1, 3, 2, 6, 1, 1, 9, 1, 4, 4, 3, 90, 2, 1, 2, 1, 2, 3, 4, 25, 1, 7, 1, 7, 9, 1, 10, 2, 2, 1, 2, 1, 7, 8, 4, 3, 1, 72, 112, 2, 17, 2, 2, 11, 1, 6, 1, 2, 1, 3, 7, 2, 1, 1, 5, 1, 1, 7, 4, 6, 1, 2, 1, 2, 3, 1, 1, 1, 1, 1, 22, 1, 35, 1, 5, 2, 3, 1, 1, 5, 2, 2, 2, 4, 13, 2, 2, 1, 1, 1, 1, 48, 10, 1, 1, 4, 6, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 16, 3, 1, 32, 25, 4, 1, 1, 4, 1, 10, 10, 17, 1, 1, 43, 1, 8, 311, 9, 1, 16, 1, 230, 1, 2, 2, 1, 3, 4, 1, 1, 3, 1, 1, 2, 1, 3, 1, 3, 15, 1, 5, 2, 5, 1, 28, 12, 2, 2, 4, 1, 4, 4, 10, 2, 1, 1, 1, 8, 25, 1, 3, 1, 9, 1, 5, 2, 5, 1, 3, 9, 4, 2, 46, 2, 1, 2, 2, 1, 2, 2, 2, 2, 1, 61, 4, 3, 4, 1, 3, 3, 12, 1, 6, 1, 6, 1, 3, 1, 1, 1, 12, 1, 3, 2, 1, 2, 1, 21, 1, 3, 1, 4, 1, 1, 4, 8, 10, 1, 6, 1, 19, 2, 3, 4, 1, 6, 287, 12, 2, 81, 138, 2, 11, 3, 1, 1, 1, 1, 1, 2, 1, 30, 1, 1, 1, 27, 1, 22, 31, 1, 1, 4, 3, 1, 10, 127, 15, 1, 23, 1, 3, 1, 25, 1, 8, 1, 1, 2, 2, 1, 1, 8, 33, 2, 19, 1, 1, 2, 10, 3, 3, 1, 2, 1, 4, 1, 13, 1, 1, 5, 19, 1, 1, 1, 1, 1, 3, 7, 3, 1, 3, 2, 2, 3, 1, 1, 1, 10, 1, 1, 101, 1, 2, 1, 1, 1, 4, 13, 2, 1, 33, 9, 2, 1, 1, 20, 1, 2, 4, 1, 88, 1, 2, 88, 1, 1, 104, 1, 1, 5, 1, 1, 1, 1, 2, 2, 1, 6, 4, 14, 5, 1, 1, 12, 1, 2, 1, 1, 2, 4, 17, 1, 1, 4, 2, 1, 2, 3, 3, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 40, 2, 5, 3, 15, 1, 2, 4, 1, 5, 1, 7, 1, 2, 4, 1, 4, 1, 22, 1, 24, 1, 9, 1, 4, 1, 6, 1, 4, 4, 2, 5, 7, 1, 4, 2, 1, 54, 3, 8, 2, 1, 4, 3, 3, 1, 287, 1, 6, 4, 1, 3, 1, 1, 5, 1, 25, 1, 3, 58, 3, 1, 2, 5, 1, 2, 2, 1, 81, 2, 1, 2, 1, 2, 7, 1, 10, 1, 5, 2, 1, 8, 17, 1, 6, 3, 3, 1, 2, 1, 1, 2, 1, 4, 1, 3, 6, 1, 1, 4, 2, 28, 1, 1, 3, 3, 1, 6, 6, 1, 3, 1, 2, 6, 1, 4, 11, 1, 1, 2, 1, 1, 67, 11, 4, 1, 109, 4, 1, 2, 22, 16, 2, 1, 8, 19, 2, 3, 15, 1, 64, 2, 1, 23, 1, 1, 2, 1, 1, 2, 5, 2, 1, 2, 1, 5, 1, 1, 1, 2, 2, 5, 2, 2, 1, 8, 1, 14, 2, 1, 1, 10, 12, 1, 51, 1, 1, 30, 2, 1, 9, 1, 3, 1, 2, 1, 10, 2, 9, 11, 9, 3, 1, 1, 2, 2, 1, 2, 1, 1, 1, 1, 3, 1, 28, 1, 2, 7, 1, 2, 2, 2, 9, 1, 4, 5, 1, 3, 1, 1, 19, 1, 1, 2, 1, 17, 3, 1, 1, 1, 1, 2, 1, 31, 3, 1, 1, 2, 1, 1, 1, 2, 29, 1, 7, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 1, 163, 2, 2, 8, 4, 1, 3, 1, 3, 8, 6, 4, 1, 2, 7, 1, 3, 2, 1, 1, 1, 5, 3, 2, 6, 1, 11, 2, 2, 8, 1, 2, 2, 1, 2, 1, 5, 1, 1, 3, 1, 6, 1, 2, 3, 2, 3, 1, 4, 4, 1, 6, 4, 2, 7, 2, 1, 1, 8, 1, 72, 4, 1, 15, 3, 2, 7, 3, 2, 2, 2, 2, 1, 1, 7, 2, 1, 12, 4, 2, 2, 1, 2, 1, 3, 4, 6, 2, 385, 8, 1, 3, 1, 1, 1, 16, 3, 3, 1, 1, 1, 1, 2, 1, 1, 5, 3, 1, 9, 2, 2, 28, 1, 15, 1, 3, 20, 1, 9, 1, 1974, 1, 15, 1, 1, 1, 1, 4, 1, 7, 5, 1, 2, 23, 2, 1, 1, 1, 1, 1, 2, 1, 6, 2, 7, 2, 9, 28, 1, 4, 1, 1, 2, 1, 3, 1, 4, 1, 1, 1, 4, 254, 6, 2, 1, 21, 6, 3, 1, 1, 4, 2, 11, 2, 1, 1, 1, 7, 1, 4, 8, 2, 1, 12, 1, 5, 14, 18, 1, 3, 1, 1, 1, 1, 6, 2, 3, 3, 3, 2, 1, 6, 1, 2, 4, 5, 8, 1, 1, 204, 11, 1, 1, 1, 1, 15, 1, 1, 28, 1, 1, 75, 25, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 3, 3, 2, 1, 3, 15, 1, 2, 5, 1, 7, 1, 9, 30, 2, 1, 5, 1, 14, 1, 1, 42, 4, 2, 1, 12, 2, 1, 2, 1, 6, 9, 1, 1, 8, 1, 7047, 1, 35, 1, 11, 1, 1, 13, 2, 1, 3, 2, 2, 1, 1, 1, 25, 1, 7, 1, 9, 2, 4, 9, 1, 3, 1, 14, 2, 1, 3, 1, 3, 2, 1, 2, 2, 3, 5, 1, 1, 1, 3, 2, 7, 3, 4, 1, 15, 3, 4, 221, 1, 5, 4, 1, 1, 1, 2, 5, 31, 1, 3, 1, 1, 7, 1, 3, 1, 3, 1, 1, 1, 8, 1, 8, 1, 1, 244, 2, 6, 5, 106, 3, 6, 8, 3, 1, 8}, giving these values for u and v:

… as taking one more term in the continued fraction would give a value of v of 6 1384803451 8375760084 9070876689 8738928277 1869098837 8128326311 3558142005 8380850809 6382758816 9496008617 1268487429 4412128531 5092445718 2363158773 4860718306 1275938336 6978080669 9594719100 5065640117 3015349970 4231548964 5055216889 4190214133 5983048172 8744372250 9773645818 3033283653 3568739940 7682955896 9327537384 1352799214 2778046214 7383071279 4963936882 1980433397 8158964056 7078686165 6470096372 9124464054 2530380893 2288444478 3388141711 0123192259 8536961659 9883507479 8124405990 6627908003 4791322733 8258010029 0000883590 5923608950 0350243993 1527235139 1817949205 8058852932 6381908832 0294353010 1544867265 3334866206 5222779200 0158516966 2387504574 1403650194 5837688768 0023575828 4576104305 5314389331 4240322394 9091000314 7791092211 4758444542 6248698922 5587468476 9517517620 6153583290 7145470469 5639459988 7683597266 0228268221 9269006974 0952202601 4368318782 3730585019 0746699366 9348187775 0550358981 1391353218 8831158270 4363091463 7648282533 8063642815 9172843473 4517242318 3988043855 3358527337, which is too large.

We also need to calculate d = floor(c4·v/F + 0.5) = 484811831 1545113846 1662452941 2845261775 2451046417 3152735585 7682888536 8947277034 5517017685 8522549376 6581317467 0184612559 0922489267 5985974641 8520449632 2248088201 7309532243 7368269484 6871839748 2985624293 8557793788 7211816830 2064177028 2617177366 6955903104 8046314972 7333039389 3521353627 6764654602 6988840725 1316923282 5622059522 6244569464 4296465917 9671648162 0606487425 6107484840 2388945808 4070598708 9185004613 7040203021 4861872617 7808848581 4881188468 7252455317 2147981598 7565888732 6677004302 1120453628 0906548033 1753275951 4490137203 0487661487 8999536512 6793310369 7610756347 7307130211 4094624984 7176219813 4158200845 5098226358 7946043292 3334754294 4777378082 3983411310 7991800642 4934994878 3082842553 6083917625 5207145696 5039798017 1243369214 6722114862 3420775333 3508798821 1446320839 2055578528 6749910810 6260828256 9749672714 6875740195 4705423812 6940815246 8745000163 3266474375 8464188132 2470302760 5997099324 9341071284 8766249034 8513347312 5127397660 3435840961 1947297629 8211780097 3410537638 0088230416 6710141553 0135835621 5008025925 7417192384 0718131273 5326515067 2627024768 6783634180 8505856300 2231756125 8918381945 4177598958 4881661507 4701737332 0945369258 2770529090 3988217504 3846557566 3134984864 2873862774 1763880802 1754346661 0307912063 2323135876 2778317188 2236160806 1101131473 2671967672 1955930759 8479315323 3043192632 8439390262 6667207754 6809746015 3259211737 7239434939 2497724791 8150799040 9205237601 3536965220 5908489819 3498111281 3544463404 6862763659 4188184926 1663327144 0822449525 0981249980 9390655493 3121228961 4024408944 3592221605 2126015145 9630596349 4872199591 3599983759 3156212885 6128801074 3892449646 3365261708 0990843743 0355069095 9670869727 3832493184 3657925791 4839734456 3230226380 9613748845 1798857380 6018944413 6139709510

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.