Primality Certificate for (16339^3613-1)/16338

Andy Steward15,219 digits20 January 2008
Originally by A.A.D.Steward 2008

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

From Factorisation
1633916339
Φ22 · 2 · 5 · 19 · 43
Φ33 · 1759 · 50593
Φ42 · 29 · 1429 · 3221
Φ6181 · 1474843
Φ77 · 2087 · 41229077 · 31590413914417
Φ1213 · 433 · 12661076731549
Φ1419025069730499281786622867
Φ21298579 · 15577087 · 164163920341 · 12125796545479 · 39097194578863
Φ2815569 · 33575775650475241 · p30
Φ426553 · 512277662921851 · p33
Φ431033 · 32046869 · 709618436439255848367485537 · p140
Φ84318263195881 · 791829475076751300837387469495141657 · p54
Φ8643 · 173 · 14621 · 189045733086155388529 · c149
Φ129492007 · 31832851749118532071 · 184561196311182232921 · c309
Φ1723613 · 1870682592277 · c339
Φ2581549 · 10837 · 33449701 · 1347559692092796509008484449 · c313
Φ3016083813 · 1624980407 · 905740901152697 · c1031
Φ51614449 · 17029 · 46957 · 1003684501 · 8494803569341 · c673
Φ602c1062
Φ90378443419373089 · c2110
Φ120454181 · 58997 · 2525813794324969 · 77512184961462965141 · c2079
Φ1806489427 · c2118
Φ36128275093 · 55164143581 · 333576390704518921 · p4212

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

56 0597924013 0417680088 2662865403 0721420910 9822443123 9353501164 9503307738 9624859156 4329892885 1037193508 4103734845 4261052941 3325223607 1862291593 3644192373 8869507720 3074897995 2608601987 0949944833 7008130080 3176457052 1839576300 0495454303 9152960525 8944395501 0882213796 3031809995 5517825457 1167342222 8240897371 3451515095 7688820931 4743464858 3330526034 1243090838 5592384700 2227415190 4992181943 8235345242 7243854015 9914684357 5276295356 7273021745 0318296305 4457263006 9779890047 7059208935 2021475590 8541725613 4282279446 8138956957 5587219624 2938049600 1014009465 4810480285 1905778416 9926573236 1535668023 6751347334 8631400889 3474387910 6391390989 5496314172 8617059934 4663999687 0902142027 4974898385 2569589727 0218825305 3845871649 8892738267 1703701854 9581123389 2247536851 9052674657 6026695030 4896929421 8614836153 3131235190 5110252594 0910432806 8655940812 2311750598 7575936875 8698481595 9329970899 1312286310 0675585524 9723531089 6683252714 9479394662 5211110029 5007719789 0902773079 0284757831 5832012235 3601563907 3869995470 4375528291 6215462704 0274914832 3745526888 0344571811 6114311943 1356345412 4075989448 1399664243 6670098735 5259881687 9161237794 5092798075 2680817875 7420549140 4513082374 8191476222 3690278155 0842241089 7575176825 9329035032 1303369581 6427432202 8691375175 6184036370 6232524211 0001954156 7682024841 9260959948 2880128658 1449580296 1582021640 0833895231 0132608812 3685457310 9108275477 4668310471 1526128737 3508121016 0954004188 9325862565 8838171742 2024549951 0238320841 1586208396 5729444648 2025636750 7513806484 3250426575 5049846673 6119637638 0458840011 4236538003 0757740246 4953284536 6809068798 9738235906 6234880620 9911641844 6480721465 5013651260 7050831624 7552403067 6080482346 1907552016 8719392460 2352693756 6881580021 7137425935 1424097454 8097745811 0437861066 6065377587 9652099458 9755913871 3628321025 7580153049 5084270154 0657475953 8450995317 2260884906 3968627832 4501755167 4872396375 1498331376 4795021576 9345760632 8382656856 5589512997 3986910380 1260285861 6233813908 3579597836 3167746143 8236511427 9136421168 5307704147 1611502069 6184679045 6944350671 0804771657 5173282866 4078583802 0312560275 4375303358 7516243110 2301293115 0542143828 8742337478 9254827641 8515599054 3158544007 2986475520 2362063950 5607822296 5366849193 6440974568 1417924227 8016072331 5647677622 7796983481 7573236288 7669084476 0090802371 7452608363 1155774265 9694434212 1002410886 7049873999 7558475469 0427103332 5504534828 5033952207 5533749316 0829356720 8701158758 6487534480 6924273552 0281731585 9978040308 1130889736 1202943559 2093900185 1646315313 7751463280 9347150069 7427498270 3458214933 6708446055 7779872667 9576383433 6923980321 6409963673 3289941014 9805302590 1562196584 4370210746 8997321421 9768106513 9231713705 6558902400 7030719473 9245977060 3360604218 2542436048 9529810073 6516827670 2497245739 4631952548 4552664927 5390107501 4781809308 4794955129 5120221566 4754556877 4481974486 9460770161 2973050467 7928480416 9483146842 3704198874 4509771176 1653097064 1872832159 0929039903 8203347479 5333517417 2537022392 3189441479 3694991466 2606305040 4912870527 4940918481 5346955521 1340710098 3717958735 7102084612 7663185956 0126912266 2016086237 0041766865 4070244074 2511474425 9636968494 2460869671 9466139514 0025851126 7172073790 4765461516 4879933061 0390699029 9426756850 4513107352 4206754281 1556321672 5581474159 4011043767 8040000750 6373166973 0580644680 3454108796 1031055451 5865385086 8899254820 6950808087 7974412130 8013304713 4651222877 8737632570 9008151018 0819724321 0794586134 3111150069 9084879771 2604627896 4814739392 8137720916 4083479513 3378821979 9084446586 2571211087 5610876152 8771859251 6807937867 0665940819 0723345395 9310811109 6050443678 5584669676 6279598998 4428867970 0018929050 9476380172 7122675364 6104850411 5334432609 7559914322 2910683363 7898143177 6427514121 5498192572 9230440712 5032593270 7865640557 7646467092 3281235322 8931389997 1938414729 3382701852 4409722799 7615347878 8976981479 5053841495 5896168520 3972046904 6406700085 5960857446 0843792847 1281196191 6906727050 1763249635 3486579420 3269415184 3996246257 4911869856 1418392644 6108518711 5476189340 2041804342 4427774500 8273745880 0413770462 4334647046 9546163338 4682247884 9274284730 7107750264 4945932183 7679956387 9924304115 3322275889 8895393318 5651695086 2903329651 9809364787 9127347984 7695664275 3609999661 6535914512 6874541422 4800091076 9672317596 2883343826 6659581321 4480277628 6020899752 0924711870 7570445517 8197555082 7137017751 2662134303 1026967846 1468374137
3842250012 3215424292 2656287969 5709963163 0830095355 2384221699 1236909352 9858220585 6173486683 2595679351 3069271477 6456912743 2141674891 3737409489
5199 8793230071 0736937294 3163588978 5305870014 1165139313
791829 4750767513 0083738746 9495141657
107 8416850522 8927002776 2416616787
6924985983 5089599142 8640527809
13475596 9209279650 9008484449
7096184 3643925584 8367485537
190250 6973049928 1786622867
1 8904573308 6155388529
1 8456119631 1182232921
7751218496 1462965141
3183285174 9118532071
33357639 0704518921
3357577 5650475241
252581 3794324969
90574 0901152697
51227 7662921851
7844 3419373089
3909 7194578863
3159 0413914417
1266 1076731549
1212 5796545479
849 4803569341
187 0682592277
31 8263195881
16 4163920341
5 5164143581
1624980407
1003684501
41229077
33449701
32046869
15577087
8275093
6083813
1474843
492007
489427
298579
58997
54181
50593
46957
17029
16339
15569
14621
14449
10837
6553
3613
3221
2087
1759
1549
1429
1033
433
181
173
432
29
19
13
7
5
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) = 33.115720%

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

Given such a witness, Pocklington's Theorem shows that every prime factor of N ≡ 1 (mod F). As F4>N, N can have no more than three prime factors.

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 = 1 3436543078 4428103567 7721105972 8843064319 1651549393 5793878568 9930888629 2649111728 3333585796 9118184425 5422572287 3090099449 7549234726 2791549921 7908611926 6087328062 2543850523 9421215958 8658722308 6115492328 9238704095 2536787372 1360174146 5599036249 7696683911 5578018557 9558174686 3693182610 9297579450 3427658485 3863237405 6319257097 4564901286 3787990495 6113683975 4091088394 4595656808 8858505367 0946587590 5467345792 0541582526 5392801231 4576028306 1139901632 7642625592 7350845054 0992330194 5963261518 4650480976 7827794449 4362501568 4707644452 2346530828 3912796277 4599261589 7275917729 1366372483 1938916817 4989594126 4002993780 8117030884 5831578512 3179413724 9385195209 9068658184 2110077077 8891194439 0449068517 8433826826 8476970830 8582468241 4749286635 6461663044 3398272353 9830170363 3126780135 4026692849 7561635767 4975485325 6617460732 4875272077 9204269181 5695881965 8199193195 8831660432 4385228262 2891309929 1556742950 4809391122 3338028155 5136445164 4239302486 9310753161 7744891368 2747279222 4399945616 9431708970 8392446607 8461434271 4062834152 4919498339 6417643336 7140073610 7818917881 9863162562 8650336310 4448189783 2763393443 3244823126 8831295036 2976507126 8814459890 6547724554 7785324000 5704427831 7519099082 7995982532 3322105399 1437809486 2749266640 2037036179 1340147113 2678817714 5752534095 0485967665 0816320243 5093652439 8138853354 6249966623 0154621506 9896048217 5936438687 2568300351 6358213424 5427850480 3473669538 9700151747 2685671380 8095513719 8514424048 5300885459 2225292364 0161332456 8984915628 0590690370 8735093097 6428205565 7130705577 0593548466 3964653155 4894271639 1457047798 4872689238 3512272784 5962483888 0156892798 1015436801 5812133723 1499472211 6345043635 5382269870 5199460662 9539423271 6424316554 2368906334 6347722548 3736908649 0798145115 3349037408 5935191791 1081239911 6408829404 0160959525 8599704795 8788361069 6836358299 0559920785 4472590859 6798192862 8045618792 3175500051 7306058935 9950434595 7116818961 5698638075 6569129627 5377375233 0368106083 7173184470 1375846672 5983835482 5999813897 6240175002 4490195081 5034687811 2875019449 3000995691 4344100778 9969615932 0075384796 2204093000 6234478061 6233569705 2110705714 5272881938 7905722526 0383933394 9776742915 8557197850 9074904716 7949412163 4090725071 2399133390 5263410098 4196638652 7003971493 2770804015 0440333771 6677999372 5368162221 6555404704 3650737521 1507737541 0179116073 7878302468 2803393880 8569417816 7488213383 1601818560 9408996464 2428342210 6724051516 5898834133 3630785474 9271969644 6970673353 3607189441 2809283346 9020310796 6268542355 4482871060 1200171191 5119537364 0760615260 4326759383 6343143552 6625483594.

With those constraints, the unique continued fraction is: {0, 4, 1, 1, 2, 2, 8, 1, 3, 82, 2, 6, 1, 1, 4, 2, 1, 13, 1, 1, 2, 3, 51, 2, 56, 1, 2, 1, 1, 41, 3, 1, 2, 2, 19, 1, 2, 1, 4, 7, 3, 3, 1, 3, 1, 3, 3, 1, 1, 3, 1, 11, 1, 2, 1, 1, 1, 1, 3, 1, 8, 4, 1, 1, 1, 1, 4, 3, 1, 5, 1, 48, 12, 1, 2, 1, 2, 2, 206, 51, 3, 1, 1, 1, 8, 2, 1, 2, 1, 14, 4, 1, 2, 3, 2, 597, 1, 4, 2, 1, 4, 2, 2, 2, 2, 3, 9, 5, 2, 1, 12, 2, 4, 1, 3, 3, 23, 13, 1, 1, 3, 2, 3, 1, 89, 2, 1, 2, 12, 4, 1, 1, 42, 2, 5, 2, 27, 1, 13, 3, 2, 2, 4, 1, 1, 1, 13, 1, 3, 3, 3, 1, 5, 1, 1, 1, 3, 1, 10, 2, 8, 6, 2, 1, 1, 1, 11, 44, 1, 19, 1, 1, 2, 2, 2, 1, 7, 1, 2, 2, 1, 2, 4, 2, 5, 1, 3, 2, 3, 3, 1, 12, 1, 1, 6, 1, 3, 1, 14, 1, 1, 4, 1, 2, 1, 4, 2, 2, 4, 17, 1, 1, 1, 1, 2, 6, 3, 2, 61, 1, 1, 3, 1, 1, 31, 2, 100, 2, 7, 2, 1, 3, 1, 1, 1, 1, 5, 3, 11, 1, 14, 1, 2, 1, 569, 1, 5, 5, 1, 310, 4, 1, 1, 5, 1, 6, 1, 3, 1, 2, 9, 1, 7, 1, 3, 10, 3, 1, 1, 2, 2, 4, 1, 2, 2, 1, 4, 7, 4, 17, 14, 1, 5, 2, 1, 18, 1, 3, 4, 8, 2, 1, 1, 2, 7, 11, 2, 1, 1, 1, 101, 3, 1, 1, 1, 1, 1, 2, 7, 1, 3, 1, 1, 1, 3, 2, 8, 2, 3, 5, 1, 1, 2, 29, 1, 4, 1, 2, 2, 1, 5, 1, 2, 1, 1, 1, 1, 61, 1, 1, 1, 42, 1, 1, 1, 3, 1, 16, 2, 1, 2, 4, 8, 1, 1, 1, 1, 1, 43, 1, 33, 4, 2, 1, 4, 226, 1, 3, 1, 1, 9, 1, 1, 1, 7, 2, 5, 14, 2, 2, 2, 1, 5, 1, 1, 19, 2, 7, 1, 38, 8, 1, 2, 10, 1, 2, 3, 4, 99, 1, 5, 3, 1, 5, 1, 2, 3, 1, 2, 28, 2, 2, 1, 1, 1, 2, 1, 1, 18, 5, 1, 2, 4, 1, 2, 1, 1, 45, 1, 3, 1, 1, 5, 1, 3, 5, 1, 5, 1, 2, 3, 1, 3, 1, 12, 1, 1, 1, 5, 2, 4, 1, 3, 7, 11, 1, 100, 30, 2, 1, 2, 3, 2, 2, 1, 1, 2, 12, 4, 2, 1, 2, 1, 2, 2, 1, 1, 1, 3, 96, 2, 46, 6, 5, 1, 5, 1, 1, 6, 1, 40, 1, 3, 3, 1, 22, 1, 1, 23, 1, 3, 1, 2, 1, 1, 1, 2, 4, 2, 2, 4, 6, 2, 2, 3, 5, 1, 5, 1, 1, 1, 1, 2, 3, 1, 2, 1, 1, 2, 2, 1, 1, 11, 1, 10, 2, 3, 2, 23, 1, 1, 4, 1, 1, 10, 8, 163, 1, 11, 1, 2, 1, 3, 1, 5, 1, 3, 2, 1, 1, 1, 1, 3, 9, 1, 3, 2, 2, 6, 2, 4, 33, 2, 2, 4, 5, 2, 11, 1, 2, 1, 6, 2, 2, 2, 1, 6, 5, 3, 4, 2, 1, 4, 1, 2, 1, 3, 2, 4, 1, 1, 2, 9, 5, 13, 1, 7, 1, 4, 5, 1, 2, 3, 13, 1, 1, 2, 1, 34, 1, 1, 1, 27, 6, 1, 7, 1, 9, 7, 6, 8, 11, 2, 1, 1, 2, 41, 5, 3, 1, 123, 1, 1, 1, 3, 1, 3, 1, 4, 4, 5, 1, 1, 1, 11, 1, 2, 1, 4, 2, 1, 6, 6, 6, 1, 2, 1, 3, 1, 2, 3, 4, 1, 5, 2, 20, 3, 1, 1, 26, 1, 20, 24, 1, 8, 1, 1, 2, 2, 10, 3, 2, 20, 1, 2, 3, 17, 32, 1, 3, 2, 1, 2, 4, 3, 1, 6, 1, 1, 1, 1, 1, 12, 1, 1, 3, 30, 3, 12, 2, 4, 1, 12, 2, 2, 3, 5, 1, 5, 3, 118, 10, 1, 2, 1, 1, 4, 1, 1, 1, 89, 2, 2, 10, 1, 1, 2, 195, 1, 1, 1, 5, 3, 1, 4, 1, 1, 6, 3, 2, 3, 2, 1, 65, 1, 5, 12, 1, 46, 3, 2, 3, 87, 8, 1, 10, 1, 10, 2, 275, 2, 2, 1, 12, 1, 1, 2, 1, 1, 1, 25, 1, 388, 2, 137, 1, 123, 4, 1, 4, 2, 1, 2, 1, 5, 1, 1, 1, 6, 1, 1, 3, 1, 8, 2, 1, 17, 7, 5, 3, 1, 2, 2, 1, 1, 1, 1, 5, 1, 4, 3, 2, 9, 1, 33, 4, 1, 8, 2, 2, 1, 2, 1, 4, 1, 2, 1, 4, 4, 1, 1, 2, 2, 2, 2, 3, 1, 1, 3, 1, 34, 2, 1, 1, 3, 1, 3, 5, 4, 1, 1, 4, 6, 2, 10, 6, 3, 13, 1, 2, 2, 2, 21, 1, 1, 1, 15, 1, 42, 1, 6, 1, 3, 14, 1, 1, 7, 2, 3, 2, 2, 1, 1, 3, 2, 4, 130, 1, 9, 1, 1, 13, 1, 1, 1, 30, 1, 55, 2, 71, 4, 1, 2, 4, 4, 1, 2, 1, 2, 1, 2, 1, 1, 41, 1, 1, 18, 5, 11, 1, 7, 1, 3, 1, 2, 6, 1, 2, 13, 3, 2, 18, 3, 8, 2, 7, 1, 1, 20, 2, 16, 14, 1, 4, 3, 65, 1, 1, 2, 1, 12, 26, 1, 1, 1, 7, 4, 3, 1, 2, 1, 3, 1, 1, 104, 6, 12, 1, 5, 1, 1, 3, 1, 7, 2, 2, 3, 1, 191, 4, 1, 1, 1, 4, 3, 1, 10, 49, 22, 1, 3, 1, 8, 2, 13, 1, 4, 3, 6, 2, 2, 2, 6, 1, 2, 1, 1, 105, 1, 50, 1, 6, 1, 6, 1, 2, 1, 14, 1, 10, 3, 4, 19, 1, 3, 1, 1, 3, 2, 1, 5, 2, 4, 3, 1, 21, 1, 2, 2, 1, 2, 1, 3, 12, 1, 1, 2, 26, 1, 4, 1, 1, 24, 10, 1, 4, 2, 7, 6, 4, 1, 1, 2, 1, 3, 3, 2, 3, 3, 6, 10, 1, 12, 1, 2, 1, 1, 10, 1, 1, 6, 2, 4, 1, 1, 1, 4, 1, 4, 1, 8, 2, 2, 1, 2, 6, 37, 1, 1, 1, 2, 3, 1, 1, 11, 1, 3, 1, 64, 12, 1, 21, 3, 1, 61, 1, 2, 4, 15, 1, 4, 4, 1, 3, 88, 1, 13, 3, 2, 1, 3, 2, 3, 6, 2, 2, 3, 1, 1, 2, 3, 3, 5, 1, 43, 1, 4, 576, 1, 3, 2, 1, 12, 5, 5, 3, 2, 1, 372, 1, 1, 1, 1, 37, 1, 4, 4, 1, 2, 1, 1, 1, 1, 1, 2, 24, 1, 41, 11, 2, 2, 1, 2, 2, 1, 3, 6, 12, 1, 1, 4, 25, 2, 4, 24, 2, 28, 2, 5, 1, 2, 4, 1, 6, 1, 1, 2, 1, 1, 2, 3, 1, 9, 8, 2, 3, 8, 2, 1, 67, 1, 2, 2, 2, 1, 4, 1, 1, 1, 16, 1, 1, 5, 3, 3, 5, 6, 36, 1, 8, 1, 4, 3, 5, 1, 2, 1, 47, 2, 13, 1, 1, 1, 2, 1, 9, 1, 1, 3, 6, 1, 2, 1, 2, 4, 1, 3, 1, 1, 14, 3, 2, 10, 1, 3, 14, 1, 6, 1, 17, 2, 4, 1, 1, 1, 2, 1, 1, 6, 2, 28, 2, 5, 3, 4, 4, 6, 1, 2, 1, 1, 2, 4, 1, 6, 2, 1, 16, 2, 1, 3, 1, 3, 2, 4, 1, 7, 1, 1, 2, 1, 6, 1, 1, 4, 3, 1, 12, 3, 2, 3, 1, 1, 1, 1, 4, 1, 3, 2, 1, 1, 1, 1, 1, 12, 4, 3, 4, 1, 2, 1, 1, 1, 1, 3, 1, 1, 2, 9, 2, 2, 2, 1, 1, 1, 15, 2, 3, 11, 5, 3, 1, 5, 1, 1, 10, 1, 19, 1, 1, 1, 188, 5, 1, 4, 3, 1, 7, 1, 763, 1, 1, 1, 4, 1, 2, 28, 4, 3, 1, 3, 3, 19, 6, 1, 1, 1, 1, 4, 1, 8, 1, 3, 5, 1, 3, 1, 4, 1, 1, 1, 4, 1, 80, 1, 10, 11, 1, 3, 1, 1, 2, 1, 1, 1, 5, 9, 1, 1, 43, 1, 3, 1, 5, 1, 1, 1, 3, 1, 4, 4, 1, 1, 12, 6, 1, 4, 1, 1, 5, 7, 2, 4, 1, 23, 2, 5, 4, 1, 11, 9, 1, 2, 15, 1, 1, 1, 3, 1, 2, 24, 1, 1, 1, 1, 8, 12, 1, 10, 1, 4, 1, 1, 29, 3, 1, 3, 18, 1, 1, 4, 4, 1, 13, 6, 1, 5, 2, 2, 8, 4, 1, 3, 20, 1, 3, 5, 2, 3, 1, 8, 1, 36, 1, 1, 12, 3, 1, 4, 1, 2, 2, 2, 8, 15, 1, 1, 1, 1, 2, 3, 2, 1, 4, 1, 30, 7, 7, 1, 84, 1, 166, 6, 1, 3, 1, 3, 1, 6, 1, 2, 1, 21, 2, 1, 89, 1, 1, 1, 2, 6, 1, 1, 2, 2, 1, 2, 6, 1, 14, 273, 13, 3, 2, 7, 3, 1, 4, 2, 3, 31, 1, 16, 1, 6, 1, 1, 1, 5, 5, 5, 4, 2, 1, 39, 81, 1, 4, 3, 1, 1, 12, 1, 6, 10, 6, 5, 35, 4, 6, 1, 68, 1, 2, 5, 2, 1, 2, 7, 2, 1, 1, 2, 1, 1, 1, 22, 7, 18, 2, 1, 2, 1, 1, 2, 4, 1, 3, 1, 1, 4, 1, 11, 6, 6, 2, 16, 3, 4, 1, 3, 2, 2, 1, 1, 1, 31, 4, 5, 1, 1, 2, 3, 3, 2, 5, 12, 1, 1, 1, 1, 1, 9, 2, 1, 1, 10, 1, 6, 3, 1, 1, 4, 1, 11, 1, 2, 1, 1, 1, 3, 1, 4, 2, 1, 11, 1, 1, 3, 1, 6, 4, 4, 5, 3, 108, 2, 5, 250, 1, 133, 1, 2, 2, 13, 44, 7, 3, 3, 1, 7, 59, 1, 1, 1, 1, 1, 1, 2, 8, 1, 1, 3, 3, 1, 3, 7, 5, 27, 2, 1, 13, 679, 5, 2, 6, 1, 2, 4, 1, 2, 7, 1, 71, 2, 3, 4, 1, 1, 19, 1, 1, 3, 3, 1, 1, 5, 3, 2, 80, 1, 3, 1, 1, 8, 1, 6, 1, 3, 8, 3, 2, 4, 1, 1, 2, 3, 1, 5, 20, 1, 1, 60, 4, 1, 116, 4, 1, 3, 1, 10, 4, 12, 1, 17, 1, 2, 1, 5, 1, 8, 3, 1, 3, 1, 1, 3, 1, 2, 5, 2, 2, 4, 2, 1, 1, 1, 1, 2, 2, 1, 3, 11, 2, 1, 7, 6, 4, 1, 5, 14, 1, 2, 2, 6, 17, 1, 2, 13, 1, 2, 3, 1, 11, 17, 1, 8, 8, 4, 34, 2, 3, 1, 2, 12, 2, 3, 1, 9, 1, 62, 1, 2, 1, 54, 1, 9, 8, 1, 3, 6, 5, 1, 4, 1, 1, 4, 23, 4, 1, 1, 1, 149, 1, 32, 1, 2, 2, 1, 1, 3, 3, 1, 2, 1, 1, 1, 3, 10, 1, 4, 1, 1, 2, 5, 28, 3, 2, 123, 3, 1, 41, 1, 98, 1, 6, 2, 1, 2, 3, 1, 5, 23, 4, 2, 28, 3, 2, 16, 1, 23, 2, 1, 1, 2, 1, 30, 1, 4, 3, 4, 1, 1, 1, 1, 2, 2, 1, 8, 2, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 4, 1, 2669625653 3115595543 9067257288 4680073400 2639482915 3027955229 5779662763 4479404575 9523096488 2466561092 5455227791 5472311732 7859415253 4958672501 0467559603 6069774285 0748178500 2262436499 7008705995 6617525504 7609414772 8924023208 4206618377 3013970207 4170498644 4831776004 8305891888 9198127338 0095140526 0541791935 9430325091 9780267391 8497076224 1798697758 2801733207 0656409381 9558529705 9702079481 8359193692, 1, 3, 1, 2, 43, 1, 3, 1, 5, 2, 1, 1, 2, 1, 1, 43, 1, 1, 8, 2, 17, 29, 1, 4, 2, 14, 2, 7, 2, 2, 5, 1, 1, 2, 1, 37, 1, 2, 1148, 1, 3, 3, 2, 4, 1, 3, 1, 3, 1, 14, 1, 8, 3, 2, 1, 2, 2, 46, 1, 5, 1, 2, 4, 1, 2, 1, 1, 1, 4, 2, 1, 3, 10, 24, 2, 1, 2, 7, 3, 1, 74, 11, 1, 2, 1, 5, 2, 1, 1, 2, 1, 3, 11, 5, 2, 1, 4, 2, 1, 56, 1, 1, 2, 7, 2, 5, 2, 13, 1, 3, 142, 3, 1, 3, 4, 4, 4, 1, 6, 2, 2, 1, 10, 1, 26, 9, 2, 6, 1, 1, 18, 2, 1, 2, 5, 1, 1, 8, 2, 2, 1, 3, 2, 41, 8, 1, 1, 1, 1, 5, 3, 1, 4, 1, 1, 1, 1, 3, 2, 1, 3, 1, 7, 1, 3, 6, 14, 1, 6, 2, 2, 136, 1, 1, 6, 1, 3, 10, 52, 1, 2, 2, 4, 5, 4, 3, 1, 7, 13, 1, 2, 12, 5, 1, 2, 28, 2, 1, 1, 8, 1, 2, 2, 8, 3, 1, 3, 1, 1, 2, 1, 1, 6, 21, 1, 3, 11, 6, 9, 2, 2, 7, 1, 1, 1, 1, 1, 4, 1, 6, 1, 4, 9, 1, 5, 58, 1, 1, 1, 9, 2, 1, 4, 16, 1, 1, 16, 1, 3, 6, 3, 62, 1, 1, 3, 3, 1, 127, 1, 16, 1, 25, 2, 1, 1, 4, 27, 5, 1, 1, 3, 1, 6, 1, 17, 11, 1, 2, 1, 429, 1, 17, 6, 5, 5, 5, 5, 4, 2, 1, 2, 1, 9, 32, 1, 1, 6, 2, 2, 1, 13, 1, 1, 2, 2, 2, 5, 2, 9, 2, 4, 1, 5, 5, 1, 2, 2, 6, 1, 1, 4, 9, 11, 1, 1, 2, 3, 5, 4, 1, 1, 2, 2, 1, 1, 6, 6, 1, 1, 1, 4, 5, 6, 6, 1, 14, 1, 2, 4, 32, 3, 13, 2, 1, 2, 114, 1, 1, 41, 1, 1, 2, 3, 4, 1, 2, 12, 1, 6, 1, 2, 2, 28, 1, 1, 1, 1, 3, 1, 4, 2, 1, 8, 17, 2, 5, 16, 1, 2, 317, 1, 4, 1, 1, 4, 1, 5, 1, 1, 10, 20, 3, 3, 1, 1, 1, 1, 1, 1, 21, 1, 1, 35, 1, 3, 1, 1, 14, 1, 2, 1, 2, 9, 2, 4, 4, 1, 1, 10, 3, 2, 7, 1, 2, 1, 8, 1, 3, 2, 5, 1, 2, 8, 12, 1, 2, 1, 8, 1, 6, 2, 2, 1, 1, 1, 1, 1, 8, 1, 3, 2, 2, 1, 2, 1, 1, 51, 1, 1, 2, 1, 12, 2, 2, 1, 2, 37, 2, 7, 1, 3, 1, 1, 3, 19, 1, 5, 1, 1, 13, 3, 1, 3, 3, 1, 1, 83, 2, 2, 7, 4, 12, 1, 1, 1, 6, 3, 4, 13, 32, 3, 1, 1, 1, 3, 1, 2, 4, 1, 1, 8, 3, 3, 1, 5, 2, 1, 3, 4, 5, 1, 13, 1, 2, 1, 2, 2, 15, 7, 2, 1, 6, 3, 1, 1, 2, 7, 3, 1, 1, 30, 1, 5, 1, 17, 6, 9, 7, 10, 1, 68, 1, 1, 2, 27, 4, 2, 4, 18, 1, 4, 1, 3, 2, 7, 1, 107, 4, 4, 15, 19, 2, 2, 4, 10, 21, 1, 2, 3, 1, 1, 3, 20, 3, 1, 8, 1, 5, 5, 2, 1, 619, 2, 4, 2, 1, 6, 11, 4, 1, 2, 29, 1, 1, 2, 3, 3, 1, 2, 1, 2, 1, 2, 1, 1, 3, 6, 2, 3, 1, 1, 10, 1, 1, 1, 19, 3, 2, 1, 6, 2, 54, 2, 64, 5, 1, 1, 1, 9, 596, 1, 1, 3, 1, 6, 1, 17508, 1, 1, 1, 14, 1, 48, 16, 1, 8, 1, 5, 1, 1, 18, 347, 2, 4, 1, 115, 2, 2, 1, 4, 103, 1, 45, 2, 8, 1, 16, 1, 8, 1, 1, 3, 1, 1, 1, 3, 2, 1, 2, 6, 1, 13, 1, 5, 1, 3, 2, 1, 16, 16, 1, 1, 1, 15, 15, 3, 1, 3, 10, 1, 13, 1, 3, 8, 1, 22, 1, 3, 24, 1, 29, 2, 1, 3, 1, 1, 1, 24, 4, 10, 2, 3, 1, 2, 3, 1, 1, 1, 1, 2, 3, 2, 1, 1, 4, 2, 7, 2, 1, 1, 1, 1, 1, 1, 3, 5, 3, 33, 1, 1, 4, 1, 45, 6, 1, 1, 3, 1, 6, 4, 2, 2, 2, 9, 1, 2, 1, 2, 6, 1, 2, 1, 7, 1, 3, 4, 1, 1, 1, 1, 1, 1, 1, 1, 3, 11, 2, 5, 7, 8, 1, 99, 1, 1, 1, 5, 6, 6, 1, 1, 11, 1, 1, 3, 1, 4, 1, 7, 7, 6, 2, 3, 3, 3, 4, 1, 5, 1, 5, 1, 34, 3, 3, 6, 286, 41, 1, 2, 1, 10, 1, 15, 4, 9, 1, 9, 1, 3, 2, 3, 2, 2, 1, 2, 3, 1, 1, 1, 2, 1, 3, 1, 1, 2, 2, 8, 3, 19, 1, 2, 3, 15, 3, 49, 2, 1, 1, 1, 13, 1, 1, 10, 2, 45, 7, 1, 12, 2, 7, 1, 5, 2, 278, 7, 1, 16, 2, 3, 1, 23, 1, 4, 1, 2, 8, 3, 2, 8, 3, 1, 89, 3, 11, 2, 4, 1, 1, 3, 1, 1, 2, 4, 1, 1, 3, 1, 1, 4, 1, 2, 2, 2, 1, 1, 3, 1, 2, 1, 5, 1, 2, 3, 2, 52, 3, 1, 2, 1, 1, 1, 5, 4, 1, 2, 5, 1, 3, 1, 1, 1, 4, 2, 5, 9, 4, 3, 1, 9, 7, 4, 1, 5, 2, 3, 2, 1, 1, 1, 1, 2, 1, 10, 5, 7, 2, 3, 1, 7, 11, 1, 1, 35, 4, 6, 5, 5, 1, 2, 4, 2, 5, 4, 1, 3, 2, 3, 2, 1, 6, 1, 1, 2, 4, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 2, 1, 1, 17, 1, 1, 1, 2, 1, 1, 1, 133, 3, 1, 1, 6, 2, 5, 14, 3, 1, 2, 1, 2, 2, 1, 2, 1, 7, 1, 3, 1, 2, 3, 10, 2, 7, 1, 24, 5, 12, 1, 2, 5, 12, 10, 1, 1, 15, 1, 6, 1, 1, 1, 2, 2, 1, 11, 1, 1, 3, 54, 1, 1, 5, 1, 2, 13, 5, 5, 6, 2, 2, 1, 8, 1, 1, 1, 1, 5, 5, 6, 3, 3, 3, 632, 27, 2, 16, 1, 3, 1, 9, 6, 1, 849, 3, 1, 58, 1, 2, 13, 2, 3, 25, 10, 2, 1, 1, 33, 1, 1, 1, 10, 1, 1, 1, 1, 7, 1, 2, 2, 2, 3, 2, 1, 2, 6, 1, 1, 13, 1, 2, 4, 8, 19, 1, 2, 1, 3, 51, 4, 1, 1, 1, 7, 4, 1, 4, 5, 2, 3, 2, 1, 2, 1, 1, 1, 1, 1, 9, 1, 2, 8, 2, 12, 1, 1, 4, 4, 1, 1, 1, 45, 1, 1, 1, 14, 1, 1, 12, 1, 1, 6, 13, 1, 3, 4, 20, 1, 8, 2, 1, 6, 2, 2, 2, 1, 8, 1, 2, 12, 1, 1, 1, 1, 4, 1, 2, 6, 8, 2, 87, 2, 53, 1, 1, 2, 2, 1, 10, 17, 1, 1, 3, 2, 1, 1, 1, 2, 1, 4, 1, 2, 617, 2, 4, 10, 6, 1, 2, 1, 9, 1, 2, 1, 1, 1, 1, 1, 6, 1, 4, 2, 1, 1, 1, 5, 3, 73, 3, 1, 1, 23, 1, 53, 7, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 2, 3, 1, 4, 2, 4, 1, 1, 6, 1, 6, 10, 1, 41, 1, 1, 2, 2, 1, 3, 2, 4, 2, 1, 2, 3, 1, 3, 5, 24, 1, 3, 2, 1, 6, 1, 3, 1, 3, 1, 3, 61, 1, 10, 1, 4, 12, 63, 1, 1, 27, 4, 1, 4, 5, 1, 1, 1, 4, 1, 1, 3, 4, 1, 2, 2, 1, 3, 26, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 5, 1, 1, 2, 1, 85, 1, 2, 78, 1, 1, 1, 9, 4, 2, 3, 75, 2, 31, 1, 1, 1, 4, 1, 36, 2, 9, 5, 2, 10, 1, 2, 1, 3, 4, 1, 1, 1, 2, 15, 1, 13, 3, 3, 1, 1, 2, 75, 6, 1, 2, 5, 2, 1, 39, 3, 5, 4, 1, 6, 1, 2, 1, 1, 1, 2, 2, 1, 14, 1, 1, 1, 4, 4, 1, 1, 1, 144, 3, 1, 1, 123, 1, 42, 1, 1, 2, 1, 168, 5, 2, 7, 1, 1, 3, 1, 2, 24, 1, 1, 15, 1, 9, 1, 14, 1, 1, 16, 1, 4, 9, 1, 19, 1, 7, 4, 2, 2, 2, 2, 13, 6, 1, 5, 4, 2, 4, 1, 1, 1, 2, 1, 11, 7, 2, 2, 4, 1, 4, 1, 1, 4, 1, 2, 1, 6, 2, 2, 4, 1, 1, 1, 16, 2, 5, 2, 1, 44, 1, 5, 2, 6, 1, 2, 1, 3, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 8, 2, 1, 4, 1, 2, 2, 8, 3, 1, 1, 2, 10, 5, 2, 1, 1, 3, 1, 7, 6, 2, 5, 23, 2, 5, 1, 39, 1, 15, 15, 1, 1, 19, 727, 1, 5, 4, 1, 2, 21, 4, 2, 4, 1, 10, 28, 1, 4, 2, 1, 1, 1, 3, 1, 3, 42, 9, 3, 1, 1, 12, 1, 10, 1, 1, 2, 1, 1, 1, 2, 6, 6, 4, 1, 1, 1, 2, 2, 1, 5, 1, 4, 2, 1, 50, 1, 1, 16, 1, 2, 1, 1, 4, 1, 2, 2, 2, 1, 2, 14, 1, 2, 9, 8, 1, 6, 1, 1, 9, 1, 4, 21, 2, 1, 12, 2, 3, 1, 7, 2, 1, 3, 3, 2, 2, 1, 3, 38, 1, 1, 2, 4, 1, 30, 2, 10, 1, 1, 4, 1, 8, 1, 9, 1, 1, 1, 2, 6, 2, 1, 2, 45, 1, 3, 6, 8, 8, 5, 1, 5, 4, 3, 3, 1, 3, 1, 23, 1, 3, 1, 1, 1, 1, 1, 5, 7, 2, 6, 1, 1, 3, 3, 1, 1, 2, 4, 12, 1, 2, 3, 1, 1, 4, 1, 2, 18, 2, 4, 1, 2, 5, 2, 7, 9, 1, 95, 1, 3, 6, 1, 5, 13, 9, 1, 8, 1, 13, 1, 1, 1, 11, 1, 2, 5, 1, 1, 3, 1, 2, 1, 9, 1, 2, 2, 2, 116, 7, 5, 2, 11, 1, 2, 4, 1, 13, 1, 1, 5, 2, 2, 1, 7, 5, 2, 2, 3, 2, 2, 8, 1, 8, 8, 1, 29, 1, 1, 2, 3, 1, 20, 3, 2, 1, 4, 1, 2, 6, 3, 9, 1, 11, 2, 4, 1, 4, 1, 32, 2, 1, 14, 4, 1, 6, 2, 23, 4, 1, 4, 1, 8, 1, 4, 2, 1, 8, 4, 1, 3, 5, 1, 3, 1, 2, 3, 2, 2, 2, 4, 3, 4, 1, 1, 1, 1, 4, 1, 1, 1, 1, 28, 3, 27, 3, 112, 1, 18, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 3, 19, 1, 2, 10, 1, 2, 3, 4, 5, 1, 1, 26, 1, 1, 1, 2, 2, 1, 121, 4, 1, 2, 4, 9, 4, 8, 24, 1, 1, 1, 4, 1, 2, 1, 1, 3, 1, 6, 1, 20, 1, 1, 2, 2, 28, 3, 1, 6, 1, 8, 1, 2, 1, 1, 1, 3, 17, 1, 2, 1, 1, 1, 36, 1, 4, 1, 3, 1, 1, 5, 5, 1, 14, 1, 2, 1, 2, 1, 2, 1, 2, 1, 32, 1, 1, 2, 1, 65, 1, 4, 4, 2, 2, 1, 1, 1, 1, 2, 1, 1}, giving these values for u and v:

We also need to calculate d = floor(c4·v/F + 0.5) = 831164880 3100592960 3109072920 4024184965 9877212354 6371435281 6869144619 2652114232 4232181777 6066770237 7857621518 6658143967 1519998726 7301458010 8213058249 2878057029 3589705367 8308447219 6330557267 9906502354 7547971564 6267595660 4147980244 1339397034 2733791650 2301963858 8299774565 7056040956 1066504187 7116268223 3789014197 8492224687 1372546397 3130530567 3679382089 7747212996 2769089249 4155195959 7124132785 4396876020 5527550308 4452769444 8132537201 8194858716 7197496298 2574410659 7185875511 1820317946 0464722644 6100073837 4501054627 1732885751 5045875059 6701582966 9723598179 2246435218 9720289653 6197268075 5912801538 2584944199 7183429110 2860501650 6965284593 1625377222 4532505505 0554707034 0828032163 0888785220 5136973203 0610856440 1125416639 7215636071 2172627807 7673378921 5771138584 0646990809 6939179895 4355365734 0308412160 9698840321 5420327729 9370924969 4867612178 2512140587 0223149508 3921408410 1168376317 9095702811 6589915044 6540693353 0227062062 7920611718 9145348180 0151870754 6783011230 6194684391 4180733395 8733382582 1638920781 8725244872 3690472440 3666500327 6132100930 4523263285 5336295448 8772605153 4341002702 2579300435 4708979645 7475065213 9205208067 0677013339 6936891809 0454390762 7567793555 2404237640 7308160428 2366710079 3104639613 1014409755 3536021088 8371222359 7315812013 4273264011 9708880546 6188198147 6343358275 4495400208 2003677551 5194289168 7049718706 3082795838 4764682806 6107274562 6756415215 4671812568 2774334005 1623407066 0038715374 0570265762 1291949255 0483304777 4906819653 1960414534 8403333759 9547298557 9005798186 8645822946 8937139016 2131521687 9003360722 4159384979 8114718317 9346219583 6754825469 9030870972 2261815574 8352080041 1406180155 1169286455 0103936127 1806928970 4057392515 4506367275 1607382937 4737201108 7436805940 1992275174 3765684124 7674019130 2937436337 2892484058 1180585297 6652731391 3570821980 5354494126 3560297770 5769149440 6754226049 1700943641 9972141067 1027766527 3644317966 0288746675 6317065596 5430664075 1204479222 1734716326 5649800729 5689081261 6781534546 5599965131 4163215186 2897119401 2651870598 3884832548 3117470614 0260964567 9841999219 8961457165 1057866245 1507108875 4112710818 0280586259 5587717851 4834991412 3889147536 5478552059 0279292980 4740503399 6866933874 8688347573 4869679020 5862146117 7292243898 0826767063 3776948402 2611165484 1427120393 4670816346 3965327293 2436056090 5820763779 7247005809 1293594046 8592201972 1554491334 8929808409 0567283120 3623522159 3784030457 0881372283 7033038544 3083284396 4674341090 1155126384 0946925289 8861509921 5718160680 7130428986 2408179324 6855055087 4100501071 5026474049 8530912547 7902000544 8210874941 7865554719 6704620067 4450594689 0772278265 8652697176 9946881384 2260023825 4514912382 5585830169 6658835786

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.