Primality Certificate for (8781^2851-1)/8780

Andy Steward11,240 digits14 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 33.673517% factorization of N-1:

From Factorisation
87813 · 2927
Φ22 · 4391
Φ32083 · 37021
Φ55 · 11 · 31 · 61 · 61 · 101 · 9279401
Φ67 · 11013883
Φ101471 · 4041231972391
Φ15p32
Φ19761 · 2357 · p65
Φ255 · 401 · 702551 · 61996491771514944851 · p50
Φ30151 · 1831 · 2341 · 5821 · 9382897417979551681
Φ38191 · 9463 · 216829 · 3262867379 · p50
Φ50251 · 1387355351 · 1070617864321314851 · p50
Φ571007647 · 1977673 · 421029038330647 · 3137542451066807773411 · c94
Φ7529251 · 550858951 · 24911675216366871043411051 · c120
Φ95571 · 30781 · 44268449403709915861 · c258
Φ11416886735528053 · c129
Φ1502551 · 23007207283051 · c141
Φ1906782998376831 · p272
Φ285172026571 · 17649007710017311 · c544
Φ47526318211001 · 9587778662401 · c1397
Φ57043321 · 151051 · 8359621 · c552
Φ9505701 · 1325251 · 43038579752118751 · c1394
Φ14252851 · 62701 · c2832
Φ285010044760517681251 · p2824

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 :

2237 2986651447 9009499981 6505437298 6182888737 1345642621 2013978640 4937590501 3600284655 7140188274 0372621620 6629430757 2019900203 4652868095 6217259137 8266785521 3428531357 5269506317 4856478861 5000292886 0341622600 6195924110 6895361407 1313908621 3952663394 6448142443 3117391408 7845094280 3633443886 3534529314 2140332296 1369292381 9507287254 3612334169 2902904943 4118899764 9507002729 9333238075 3421694413 9831831514 7114910751 2913539906 2098661654 8933045661 7063707640 2685315291 6129864607 6139119458 2710057082 8691364609 5343343461 7583625679 9908281117 1008902203 4295505119 3575531647 1276996061 3753558699 4647979502 1133794988 4437198900 0848743719 5034269499 5500183759 6582114045 8595546588 2537671405 6039211903 3782612333 0602968785 9414655966 9648182334 1743043631 0310961482 9256398002 1335628853 8993653121 1081076094 3381970546 5180343367 7086749698 9872273401 5694620366 9617260065 7683015370 9689840186 6444859288 2495625208 8615837932 2186926937 6863833001 9721775419 9841450457 0023211606 8889584004 7973890361 7080200276 2821677319 0943936285 1314824896 7668264542 6007353131 8898003963 3829809029 5799094670 7602578367 6685092721 7315096578 4119729730 2376620323 8711167287 3383382761 7270127153 1460229695 7229369705 0462688170 4155120376 3706781918 4533605334 0577206552 9997867296 0399977788 4758137878 2593385659 1418661373 2606732103 5519043112 9887631197 6969760706 4550573682 0270470897 0329369926 3752467937 9259284089 8013003424 7443206672 7777666033 0809531936 3832899148 2942047949 9460776766 2009012161 8812242516 6695707426 0645351232 6173698140 5180888531 5416852545 4250372305 0900018799 1736752243 9004777080 4799803173 8255614484 8411613989 7173063254 0579726852 8650237678 0940499510 9404679003 0314743777 7468684757 5028799016 1611082468 9595043969 9565211116 2596153226 6191342269 1354636976 9738729673 8420845127 9143851657 2060347399 7959724632 6546942338 3588218910 5158385931 7618196558 0563998950 9461927301 9086304456 0052853949 1621829849 9634444788 8409097482 2846000847 1777924179 8862779579 4816195634 2757788271 8840850924 6505692100 3539538148 2017984348 9138891795 3547387459 9152131411 8395785328 9511161066 6807467559 7236390421 8247689246 1095168211 3645668562 2352236208 5661824073 4757992568 9575445374 7787805703 2804863159 9822648047 4500274674 2008599484 8618008945 8652639843 0646712462 7860977901 3638825922 2318948693 8194672722 4102041631 5443274191 0429930912 1328787327 5318665349 0561768549 2159922539 6532129782 1173229227 6163989744 8350597245 2235831563 0687292699 3222288387 9196943424 7142041241 6388573322 6020968141 3479723716 8707010613 6673863887 1534955608 4327157812 4504441630 2609225794 5739400701 3994824591 1858690549 2810752224 9065902992 2881758152 2914584805 6125401087 3244311009 4047438313 1997133718 7374765341 3822040801 2709485865 6106349114 7880329848 7502260657 4984461734 5242232242 1801635298 4566500377 8653671636 9928987828 8535017068 6004537132 1569522937 9430929258 0607701140 9213386050 9659782762 9202775744 2512157406 7150687573 7452500141 8231166692 5276017375 2768968541 3466315051
12 6997185867 5891538939 2168697879 8015662065 9676029515 9526022984 2122901325 4761311098 9432611569 3547397618 4261500550 8846758564 5909625940 5776561970 3490770236 6495540909 4721869838 6762465908 5687825940 7805954128 9961981396 5473090733 0362157780 5139835189 4494023973 6946969492 8792160751
53715 1522107756 5210782231 8399569629 8814924729 4210673883 9681720387
8505898485 6917359299 9021843247 0934473485 6581290501
7532916213 3649580310 3918564548 1500164205 3179249507
1992432523 2522170543 5180863830 1293775976 2371660351
35 3429141659 9354578561 6408168721
249116 7521636687 1043411051
31 3754245106 6807773411
6199649177 1514944851
4426844940 3709915861
938289741 7979551681
107061786 4321314851
4303857 9752118751
1764900 7710017311
1004476 0517681251
42102 9038330647
2300 7207283051
1688 6735528053
958 7778662401
678 2998376831
404 1231972391
2 6318211001
3262867379
1387355351
550858951
172026571
11013883
9279401
8359621
1977673
1325251
1007647
702551
216829
151051
62701
43321
37021
30781
29251
9463
5821
5701
4391
2927
2851
2551
2357
2341

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

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