Primality Certificate for (2781^3121-1)/2780

Andy Steward10,746 digits13 August 2001
Originally by A.A.D.Steward 2001

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

From Factorisation
27813 · 3 · 3 · 103
Φ22 · 13 · 107
Φ37 · 19 · 58171
Φ42 · 37 · 104513
Φ55 · 11 · 251 · 911 · 991 · 4801
Φ6307 · 25183
Φ82 · 113 · 137 · 337 · 5732513
Φ1033151 · 1803645511
Φ125821 · 54541 · 188401
Φ13313 · 547 · 5779998263261 · 216325983824119616975243
Φ153511 · 1018640381032915313965111
Φ162 · 17 · 241 · 433 · 7937 · 14970337 · 8486646209
Φ202801 · 8580581 · 12919681 · 11521957261
Φ2473 · 49010039303328366986479177
Φ2613 · 53 · 3797 · 4758833 · 24130003 · 712095219277331467943
Φ3031 · 61 · 941791 · 49620271 · 40500356971
Φ3979 · 9283 · 27269115823 · 59547611676441409 · 72908487184696252129 · p30
Φ4041 · 881 · 3041 · p48
Φ4897 · 1971889 · p47
Φ52157 · 521 · 22694929509574921 · p62
Φ606781 · 148021 · p47
Φ6553301041 · c158
Φ781093 · 9444441241 · 10561697930083 · 3210465558223274327308980337 · p30
Φ80401 · p108
Φ104c166
Φ1203456560641 · 309924109123211521 · c84
Φ130131 · p164
Φ1562423929 · 167955373 · c151
Φ1951171 · 64867141 · c320
Φ208c331
Φ2402161 · 352081 · 538149265921 · 882330087118170651841 · c179
Φ26073061 · 60186881 · 38425632281401 · c305
Φ3121406090891009377 · c316
Φ39083071 · c326
Φ520395850001 · 4565980121 · 16685029569521 · c630
Φ6241249 · c659
Φ7802341 · 57949321 · 1847879281 · c641
Φ10402081 · 71761 · 17330561 · 160080961 · 14820591761 · c1289
Φ156065521 · 1111045190815561 · c1303
Φ31203121 · 61511443384561 · 24443067188429761 · p2612

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 :

29 8506683592 1389906102 8777819224 4934278707 8084621398 8359930371 2990335890 6867784655 9867520028 9139530980 1023495434 3455574212 6985557403 5898742077 8327269619 3075795040 6420382247 6399881824 9820913089 7881102797 4166680715 9628676777 5351079284 1786853358 2471909950 1618753165 7965750837 5291409102 8676771815 1680671992 0014188985 3214081620 5716394662 4972368408 9206734905 7034857658 0882048678 7370302674 9884240016 5972236392 3435224570 2513791303 4757901434 1892726978 4204553912 6427668181 2558493795 3462922816 9586170517 5609243960 4825509350 4216920880 1368656518 9610467307 7637034263 5734875379 0950672527 7078278170 3015990423 1287331910 0290608679 3264316842 3571240207 6197322057 8666344898 5156782367 9645957719 8558855684 7027967614 6207452180 7751729354 1387895540 7083179755 5839021294 0540849655 6770219229 2270185133 5329044587 9030855174 2365818214 2488253269 2971245375 1965158455 5597234227 1706694880 6136975551 7908028933 7846513153 9411806070 1260156948 8490112729 3113951427 9494463140 5738713680 6969181947 1683611737 2195325741 2802586049 7380770028 9856648038 1639496681 5428271004 8987724513 1376565321 7605782044 4084245466 9885653495 1855921160 1470486587 3190933388 9101749414 5138065413 4773260259 3618379834 1994810540 1312952010 6678193112 3188485826 4075789397 0365155025 4941351881 7311270511 0166001557 6097551584 6818113176 1167215967 1292234707 6247957288 9358938053 1766993128 9988934557 3773079726 0200819842 5827431310 9509196764 7090154256 5275633766 7813401069 3736369193 9723858639 5115180532 8984218293 4042970220 0829028780 8948505561 7472785417 0794666512 9893107257 0406540315 7998972673 8676012556 9288078551 6446468798 7607272670 1980046581 0242983986 5621364878 8847337395 2812192198 0302080554 7185085215 6720827054 6514391034 4978071458 9263553913 2528305639 3859803809 2488957266 4583548629 2430879716 5756518031 5543108589 8472415008 6949059364 8897218735 2136623539 3380076132 7212693618 9075608458 1696827718 4095648148 1217759280 8454301027 7718336674 4469864276 1619600362 6567884568 8313074782 6327865287 8163483664 1776092272 5490995162 4829674105 2845676245 2724147219 6262826459 8371112707 1065447894 9525041607 2036700191 5103179679 4975633091 1601640784 1707415205 6502782717 0506448653 2726054906 6079632572 9330831058 2105116371 3724261443 5002578825 8884874936 9830143832 2949684940 7537681565 5485801594 5905585285 8299774293 1802182424 8349521332 3319385204 1679438822 3183232829 6130664157 1635457709 8169501785 3231786094 1081594338 6398276335 8112587435 0546230908 7068170587 6474171235 6822091291 5236086795 3466186407 7765399179 6862853000 3787650828 9288058227 6218354436 2960092094 2620661174 4513132973 9447260277 0069712848 2585795654 5428098574 5720775690 5884237135 2416237292 6874579117 4784500368 1970032208 4525754772 8859129245 7973814348 5262477745 7000614970 1495784481
1601 5197789387 0895531122 7689683089 3883750962 0038076700 3065945986 0895070006 2396972988 5084736905 2185825924 5990784512 2367581276 2505969665 3521477332 0140116331 5023554891
40858956 5070271230 8499761567 4005573124 2594139020 7363164513 8829236015 0191801455 3035173333 9249831395 5603272881
24 6693361442 7841131492 1199998831 6908381557 1641076323 8906639253
11653050 7152288579 1559669169 5139278283 1717275881
6692087 6916146556 4217956980 8915015918 0809477777
1275260 1819427016 0321396644 2887720088 1846453441
5272767320 9367224356 0220409691
1308825299 0099478504 2590365607
32104655 5822327432 7308980337
490100 3930332836 6986479177
10186 4038103291 5313965111
2163 2598382411 9616975243
8 8233008711 8170651841
7 1209521927 7331467943
7290848718 4696252129
30992410 9123211521
5954761 1676441409
2444306 7188429761
2269492 9509574921
140609 0891009377
111104 5190815561
6151 1443384561
3842 5632281401
1668 5029569521
1056 1697930083
577 9998263261
53 8149265921
4 0500356971
2 7269115823
1 4820591761
1 1521957261
9444441241
8486646209
4565980121
3456560641
1847879281
1803645511
395850001

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

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