Primality Certificate for (8662^3037-1)/8661

Andy Steward11,955 digits29 October 2005
Originally by A.A.D.Steward 2005

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

From Factorisation
86622 · 61 · 71
Φ28663
Φ33 · 25012969
Φ45 · 53 · 283133
Φ67 · 13 · 824413
Φ11629520134278949 · 3777655147787100914203783
Φ12181 · 31102416793753
Φ2223 · 2927 · 29173 · 16013760930641 · 75597366690341651
Φ2347 · 139 · 967 · 379087 · 3597569 · 1095043295519 · 3924057665321 · 349439876584511 · 32796220431269349522729553111
Φ3367 · 14983 · 7589543400319 · p60
Φ441786754179937 · p67
Φ46599 · p84
Φ6620593 · p75
Φ6942643 · 947509 · c163
Φ92829 · 1381 · 6073 · 241791089 · p156
Φ1328713 · 1122397 · 5525170069 · c138
Φ138277 · 4003 · 6763 · 29947 · 190027 · c154
Φ253174571 · 35344607 · c854
Φ276c347
Φ50623 · 1013 · c862
Φ7593037 · 9109 · c1726
Φ1012c1733
Φ151845025399 · 578819473 · c1717
Φ30361004917 · 2644357 · p3453

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 :

478 1317472161 1136066193 9671678080 8588434171 4153429246 8082054479 2818273744 5267255547 0319339641 3082758748 8789783519 5673516730 0283428324 6193709160 0556558975 0880276334 4868553823 3483860502 8241938121 7843171485 9315241287 7345506169 0164794094 0672904927 1089722633 7997293173 5329812845 1671003756 6374287333 6316259947 6850813406 2771353789 9459114575 5488602480 9136320014 5778981651 0532173299 8942428349 9816185773 8661430362 2664740566 2705447844 0724032631 1902978129 4826063046 5541542553 7016794178 8714649573 4636538341 1383975233 2356764706 4074582290 2792782291 7272822378 7452750930 6740260695 0345595612 0627160297 7717201662 6754004949 8990821929 4382089459 8572443673 3283891397 1274968777 3669568723 6847937606 2110814054 8394846564 7163889135 6527971353 1300016802 7879050468 3596247732 3372544611 1786287216 8833871685 0037256150 9735910135 7160106030 5299409764 5604369153 4713695256 9420480520 3031446691 1408328513 4814023147 3854765698 4103805785 2841157503 6670053345 0640292925 6383505155 7131248955 4227456205 7070530786 1758796182 4019321294 7661643056 6580605286 6190899451 7391082218 8949732996 7442059022 2919368055 8548286528 0107872127 4064061580 8640700013 9437948169 0880433092 1801124586 9786003457 3642458654 0326447325 6677266133 9704061166 1667470833 3795484975 8016972539 9125120781 6717524048 2808225296 1742750197 3283920515 6882244481 9904410743 6161183543 2974827191 5225079634 3100384195 5455373201 6735879398 1326592499 8781102931 0033905059 0350994396 4102404708 4271966200 3175308760 0472959318 8799160326 1422232552 9183339864 5187656339 2409645032 9188961472 7800527236 0682324162 2055634970 9227058484 7879831165 7103136804 5016576709 7020511880 8137514161 8138717261 4876725822 7291074231 2761488113 2309152174 8846657170 0471028391 5412480511 0285673013 0479171127 7580435961 5712715740 3060675131 2216577824 2612608311 9680475527 0997343202 1432171026 1586985836 9087906091 7004605841 1755327192 4515526028 7542769653 7435870998 9900599024 6944180630 4440817399 3087458603 1245604999 4036452427 6325281970 4938832446 8916936410 3519298312 8879312126 2008375899 0909662834 6753876911 7212348261 8611582148 5856233135 8185919245 7624789004 4570370825 8639314612 0855595411 4189850294 6223187188 2191551166 9851494733 6261348616 7369297393 2783681306 8800820394 8203007284 3217721135 3234733128 2426350499 5469825684 3498994052 7001076953 4965387576 7960604834 5477815320 5532970565 5652404325 2501355386 4452092792 5238622701 4984873499 3076608915 8956602213 6700692875 7476187380 4612240027 4981040311 7247700629 3330573306 1555809911 2617283371 2697170875 5251437774 6336173173 4759065014 6957315568 1580749765 4141940717 2474763457 7525210698 6429243153 5589745461 9498067842 8653604821 1673663252 1364874520 5139865023 2700047660 2551367936 9682986891 8740774014 9810123853 7097161528 7232840034 1981360021 5117272948 2601159483 5973893032 5407281167 3921563534 4643205199 6165764758 5862665671 1292809327 1039708605 4372930737 8595144487 9443096837 5577122070 3918372706 2321607044 2355111788 4960533468 2887158774 3619203686 2607992235 1710322133 4927706369 7724750126 8733563831 8074712775 1784369455 6231892440 3624350142 2251557718 1662105287 0499722540 8710577954 4082895192 3960849732 9683591937 9673628610 6198898014 3676107073 3569159045 8201272646 9584545649 2260422639 8669204478 2687547786 5599724634 4363859074 5963331061 9027495704 7955059455 9826218081 6874069899 8577204731 9407480858 4744443375 5221851870 4623620715 0201576064 4232051937 9523253247 8856689234 3725833429 6670569780 9842434322 8834093354 8249243723 9505131270 5149453034 7876404378 3692309512 1657991132 5425259353 0872889383 7525857002 0494048666 8162991252 5553774930 7078387221 9321992928 0010427376 7142933615 2316745300 0910128220 6318914629
107055 2835286822 9064064942 9122864522 5780085689 6826166559 7550252737 1015688719 2709105891 3477964110 9249612295 6575547780 2635067961 2737884705 8301709288 2788449621
7081 4629824864 3643841338 0437043791 4518517120 4405702450 6511069459 2939404096 6228312817
27459 5934713008 7143772512 2626377043 6290264599 1214432646 4384482266 8365713503
3164453 9153157992 5185169051 1881517339 3955786174 6702645485 9975607053
7420347021 6406681417 2687380514 6591461057 7130378610 1799502041
327962204 3126934952 2729553111
37776 5514778710 0914203783
7559736 6690341651
62952 0134278949
34943 9876584511

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

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 ≡ 53 (mod 64) and therefore cannot be a square and N is prime.