Primality Certificate for (6588^2647-1)/6587

Andy Steward10,105 digits18 September 2005
Originally by David Broadhurst & Bouk de Water 2003

This certificate uses a theorem of Brillhart, Lehmer and Selfridge 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 35.463945% factorization of N-1:

From Factorisation
65882 · 2 · 3 · 3 · 3 · 61
Φ211 · 599
Φ3331 · 131143
Φ613 · 3338089
Φ77 · 1163 · 2319353 · 4330548063241
Φ912007 · 16163299 · 421266597949
Φ1429 · 127 · 22194936512292310439
Φ1819 · 73 · 199 · 379 · 19927 · 160579 · 244243
Φ2143 · 617972713 · p36
Φ27109 · 1346491 · 1966718426833 · p49
Φ42967 · 2437 · 36814471 · 609974711356921 · 126327122913591601
Φ497 · 119953 · 784323812287 · 63477872214611 · c129
Φ5491369 · 535486141 · 32311439700193 · p42
Φ63631 · 883 · 2295873469 · c123
Φ989475947517 · c151
Φ12645732961 · p130
Φ14783113801 · 55813065501901 · c300
Φ189757 · 1166131 · 16319392804747560403 · c385
Φ2942684809 · 44321414047 · c304
Φ378828577 · 1742486749669 · c395
Φ441172445939317 · 215596876447 · c940
Φ88236530677 · 682000327 · 139410082477 · 467337762551899 · c921
Φ13232647 · 219964627 · p2876
Φ2646550135153 · c2879

From this partial factorization, we use sufficient of the largest prime factors of N-1 so that their product F is at least N 1/3 :

163157 7202364273 8382575038 3213541209 0031788394 4304560313 4393499147 1946415546 6370386433 6379152680 0695314628 3102086182 7341245697 5261478873 3555910718 7020333158 5015446016 2897005320 6267543601 1781598988 3494692580 2325182449 0125578842 0422235725 4173430097 9681351488 7056526208 8918250351 1475742605 5853799365 0145159118 4347577970 0295462937 9788274560 2620748339 2028095055 6466865245 5435160367 3344814643 0681671452 6076725426 5351068676 2334168188 9338556720 3184482735 9265471979 5536953368 6982081224 0765314130 6386645236 7388664250 2422776639 6267833622 0775155519 0516442792 1483151006 5212257958 5697913546 7564359449 5117554449 6830029113 8355851055 1631862319 8328140742 3697064387 6502135468 9080753673 8711038402 5314429723 1951391630 4723511649 5038317173 1562838658 1264762776 8944604203 1158976269 2498099199 9385903053 1024905676 4370176241 4311010219 6288340931 8135088427 3299973624 1314039400 3374889110 8178740819 8962999463 2780933839 7868632657 0564772874 9821685906 4578329548 3351697615 4048557270 4455118156 6225244001 5302799796 2898619215 4936270685 4142815641 9252629757 1060991584 2212228061 1366058606 3083107360 9945235684 1539362560 2232034360 8801151317 8639825754 6581656107 7887966302 4404323199 4851921892 5378552567 8497147103 6238570920 2114103608 5821760460 5877111639 5678910149 7436796449 3634625507 8466943098 6087118403 5643825535 9024674482 4344711052 1007487471 3220631911 5122020208 7121599223 4776375144 5718593595 8475099343 9879831895 2406007142 1671014849 9916125291 3225908081 2967303335 0016025103 7684027443 9901248158 2099767076 1728041205 4562846939 3428766136 2752907652 9894495125 2984367441 2169450459 4650349790 6617609566 5200567544 5251248195 5871855102 7618126077 0475330618 0741099274 1131236388 5361021298 7161122496 4362464779 9320389584 6547169147 1857621195 6141770705 4102686290 0620895186 0454223841 4929865736 4812736698 4766151173 1111357118 8769563425 5494258564 8111686979 1656089368 0551484192 6417811589 3098019287 8367119914 3168285815 6614435576 4647896910 4495515085 3166251750 7399986772 0180900677 7706937796 2143123169 5732614148 9690538784 5610191357 0609047560 7667972743 1612802286 5734312126 7285749944 9933563930 7492970157 4302157200 0205736259 2775659569 4499864837 6509061250 3157753157 4031089287 1731233246 2363674854 7927642602 8581253054 4513757387 5381087400 0680833565 9435300303 0217602385 9715246479 2905126907 1275055528 0062060175 8748280036 3868082663 4542720378 6110240049 5973858514 8265333211 2768322231 9890395315 2687773436 0394834931 3303108830 6562486713 7140125977 6088417386 7088186928 4648163830 4033390138 4607951754 9180295602 3631800084 7152317200 3720384631 6580617355 3242089956 4441570186 1442238690 4533399186 4151782246 6465153565 8301212621 6191796054 8835241221 8508117748 6918733761 6895675778 6574407337 6975182491 4994815410 0327310196 5000766191 0023443381 1900965334 5102840507 2653461116 1894860705 6921852246 0373759853 4144723887 2054110556 4470493107 8206904566 4761807471 8027634507 8540144036 5092981731 0115553617 8166925429 7335096504 1837244042 6555924399 1982667479 3011557240 4591833532 4454916342 7500729501
6529800864 4110183554 0968996769 7076921163 0954360736 7567503510 5587144892 5911773639 4856501249 4759881511 3783540614 9850215979 4026871649
189318253 4538007300 5509636221 9213830868 3953857839
34 5669620343 9123168615 7461556188 8473980781
251500 8194272391 9138812780 6187112631
2219493651 2292310439
1631939280 4747560403
12632712 2913591601
60997 4711356921
46733 7762551899
6347 7872214611
5581 3065501901
3231 1439700193
433 0548063241
196 6718426833
174 2486749669
78 4323812287
42 1266597949
21 5596876447
17 2445939317
13 9410082477
4 4321414047
9475947517

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

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 F2 < N < F3 and N ≡ 1 (mod F), we can let N = c2·F2 + c1·F + 1.

Brillhart, Lehmer and Selfridge's Theorem shows that N is prime if and only if c12-4·c2 is not a square.

Here, c12-4·c2 is ≡ 48 (mod 64) and therefore cannot be a square and N is prime.