Primality Certificate for (13803^3919-1)/13802

Andy Steward16,221 digits18 April 2010
Originally by A.A.D.Steward 2010
PrP: Bernardo Boncompagni 2010

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

From Factorisation
138033 · 43 · 107
Φ22 · 2 · 7 · 17 · 29
Φ3139 · 157 · 8731
Φ613 · 14654539
Φ6531307 · c2697
Φ1306412526618279224332911 · c2679
Φ19593919 · 3020779 · 21842851 · c5382
Φ3918144967 · 391801 · 822781 · 313714872871989739153 · p5362

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 :

22 8720841329 5243702777 5033523107 4284533349 2068985336 9260633122 3360805749 1738110943 9082214633 7931481260 2868158345 5400040986 7263429838 3932639021 3681710131 7602513582 2277357791 4737963071 4424560113 5276487356 1829602037 3498533576 8723715474 6094768755 6101154018 5145684916 1731755228 4595248095 5014171188 6931475123 7838892646 5437778440 3349312932 1964719914 9462217398 5440091044 3814231830 2621932010 8035558081 7604279969 0442059639 6924164411 0625635067 1590324188 9476441593 9703000965 6267890260 6707941625 0833045711 4749545162 9951146666 1763160690 1402086633 5506638446 6271396817 6086137672 7062767228 5184051018 3188309962 0900814001 2518231792 3839368368 5143356801 0525409918 7087795276 0529075509 2958667839 0260592877 6657984093 5143470012 9935082620 5546350764 9579134121 5898667688 5345103935 7574712516 4873196752 1464615821 3270866223 9983853793 4227592487 4816833230 8277300641 0157487940 2427454459 7998168704 7606500946 7863314070 5491905438 8880104427 4502320473 3643200683 8158670794 8963042076 5641924567 3810608150 2956028773 5685726216 4895173813 1969853655 9438290215 1078107598 4244945750 0629928315 5861560814 7703504754 4317738284 5229618716 0365807068 0160749274 8438071902 0060939142 5577653623 5011209841 0511519151 2429094553 9218244180 2392873555 8791096027 4474440250 2184583649 6174033472 3076457682 3226255355 7788834366 7409103656 5960001027 5313829938 7298888675 8987771001 4152692495 6537422741 1150788929 8149296626 8040519143 2841704431 3086316967 3046516230 5818651073 2035580093 6046167088 4377383192 7477879340 3407027813 3200878723 8204830239 1610583538 3272234547 3480636909 7702349514 5362444092 0659407669 1909669282 5548084300 6012142610 0684965484 0047189775 1233325931 5901305176 8952449372 9111016983 7278187024 8806914957 3517352382 7854923722 5576822883 2622192293 9691792088 9687882913 1544848469 2924988501 4604523962 8366071940 6848470864 3722150332 0115026183 9721211604 4890826982 4613189057 6062907025 6866027057 5307208527 3202513330 4679151386 1228598250 1918389913 8176084619 9475152671 3487131016 5696004389 9883711561 9939984578 1703559905 3444055836 2205879916 4785601184 5748864828 8151102004 5593030997 8467044212 6419023184 9566944290 5976197182 6574380276 1615100960 0983476383 4382086928 4772510662 7877278179 0816958157 2126835920 0730353164 2438842172 6373142248 8319910009 0298949018 3334356137 3693522907 9586976804 0229337524 8202891076 3187074142 1487695757 1702359962 7091473028 6905744826 3051299333 7738319686 6604946840 7744138088 6818948784 1932793542 0202591654 3556896891 6845908545 5370930906 8736925358 3117780297 9970772822 6543650460 3371231121 7722384473 4063105338 4408944674 4095375068 7145567008 7939414279 1955150428 0725420906 2762345131 4464550967 1555926724 8472072888 3400094887 2861015111 3006636088 7075900224 4980540663 2038449686 8242074841 2837016995 4677437924 7752844099 6135368210 4548814095 0804214854 9102885258 9914505897 1016508559 0504328437 7624113436 6377992162 5698686251 7855522311 1867318180 8785688021 6806011325 5896489355 9633017183 8792159180 2885953249 1650034708 9625008472 1420327395 8569568866 4711982444 9419812887 1942041918 8614743678 7408959382 5626982405 8131532189 0775326583 7914026019 3661944774 8373515822 1499981042 2713059469 0882057912 1799179924 4602195590 5991256050 5133919662 9942346639 0324899655 5243097102 7831526621 2892872925 4294648972 6960838087 3980074414 3902262907 7261161221 4475661307 6093504006 8885442413 1302427286 5069758392 4863684649 8426753975 7679318191 3474846239 3928993545 0081249601 4060026375 8007912974 8827950208 7946739934 0277620343 0981762764 2794508474 4268003785 5000002228 7226141529 9833172570 1979532619 0562907528 5757593074 4007622938 3455293212 6789655261 1822774352 5600106947 3776165636 2570112773 3232884606 6856539220 2765492886 3212871973 5025107309 1839030239 9731804585 2560406908 6628019046 8467246539 3784890533 3464735898 0473437301 4350116254 0469416312 5002488164 5369700705 6868263123 8642145243 4803631640 4309768132 8732441631 0422501173 8074941921 0461397927 3404341481 2677425913 0060068362 8662760574 4916871270 6352031275 6756484045 1126613082 3254751755 9735578493 8485896344 5354157477 8529291718 9359205521 7709047452 6570670884 4624681638 0485386330 5417625847 3880699291 6871149625 9308639743 4836210051 1976239379 7435420958 9772857854 8407875118 7859115819 8980597596 1593901441 2237268805 7320961333 0807770560 3747101739 7202652348 9903213948 4671010247 2651371417 8543417099 6022468463 2173993524 7553839249 2383428221 9633453102 8672839617 1324251867 9901196078 7596578837 3732528856 1122820414 1720440303 4359232662 6130461199 7265319240 5117901222 1181091581 2066825860 9772310048 4690920357 6846839106 2488994937 9944127159 4374465170 3767697782 9979765168 1824989078 9261570705 6128314639 5597601224 5713138203 1773224933 6948852196 7011110160 1760648216 8561282813 2875703704 1661591608 8791837386 7711211803 1590105198 1216770186 8662296145 3100471683 9765104845 8908441848 0626018679 0742868923 5869735538 2163760464 4862449788 5390148199 7867872234 2946526895 4990203347 0465417583 8684572645 0571528464 2710546883 5780267054 4284909874 3991187285 6130242546 9921492146 0247825023 6187227608 2656955219 1155232137 6366269502 2640129443 1693932568 1666536844 8213361847 3594540261 4739606255 4478065324 8123551932 3941695178 9744451613 6903170090 8527344852 0887229661 3443971066 2380526665 5932170663 5271161685 2924898088 0978002456 3387692273 4316403161 8756635266 0664254513 6479474875 2534838260 8978600980 7052326720 1006820066 3856811139 0662137052 6630536923 4485071800 1318214873 8437864253 3810585601 6212826279 2189779041 8952660224 0681294906 1217382074 6762377023 3306424471 8357310505 7982306446 7853640785 3074506410 8213893859 0875422954 1195370648 4831701379 1299449773 4951011324 9660740156 6762889301 6073131549 3484798792 9380733506 5938129371
4 1252661827 9224332911
3 1371487287 1989739153
21842851

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

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.

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

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

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

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