Wednesday, November 17, 2021

[ysfdlwnw] factoring generalized Fibonacci, lag 24

we factor terms of a generalized Fibonacci sequence a[n] = a[n-23] + a[n-24] , initialized with only one 1, in Pari/GP.

(previously, with lag 19.)

allocatemem(16*10^6)
addprimes(...)
addprimes(...)
len=24; a=List([]); for(i=2,len,listput(a,0)); listput(a,1); for(i=1,8410,x=a[len]+a[len-1]; listpop(a); listinsert(a,x,1); gettime; y=factorint(x); t=gettime; print("out ",i," time ",t," ",x," ",y))

the two addprimes commands are given below.  they are the 2nd and 3rd largest prime factors, respectively, of the most difficult terms to factor among the first 8410 terms.  with both addprimes, the computation takes about 3 hours, down from the original several months of computing.

without allocatemem, pollard rho phase runs out of memory on term 6524 (GP/Pari Version 2.9.4).

here is a function that prints out any given term by index.  (it could be made more efficient, exponentiating a matrix.)

jumpto(n) = my(a,i,x,len); a=List([]); len=24; for(i=1,len-1,listput(a,0)); listput(a,1); for(i=1, +oo, x=a[len]+a[len-1]; listpop(a); listinsert(a,x,1); if(i==n,break)); x

for reference, here are the first 531 terms of this sequence:

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 3 3 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 4 6 4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 5 10 10 5 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 6 15 20 15 6 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 7 21 35 35 21 7 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 8 28 56 70 56 28 8 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 9 36 84 126 126 84 36 9 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 10 45 120 210 252 210 120 45 10 1 0 0 0 0 0 0 0 0 0 0 0 0 1 11 55 165 330 462 462 330 165 55 11 1 0 0 0 0 0 0 0 0 0 0 0 1 12 66 220 495 792 924 792 495 220 66 12 1 0 0 0 0 0 0 0 0 0 0 1 13 78 286 715 1287 1716 1716 1287 715 286 78 13 1 0 0 0 0 0 0 0 0 0 1 14 91 364 1001 2002 3003 3432 3003 2002 1001 364 91 14 1 0 0 0 0 0 0 0 0 1 15 105 455 1365 3003 5005 6435 6435 5005 3003 1365 455 105 15 1 0 0 0 0 0 0 0 1 16 120 560 1820 4368 8008 11440 12870 11440 8008 4368 1820 560 120 16 1 0 0 0 0 0 0 1 17 136 680 2380 6188 12376 19448 24310 24310 19448 12376 6188 2380 680 136 17 1 0 0 0 0 0 1 18 153 816 3060 8568 18564 31824 43758 48620 43758 31824 18564 8568 3060 816 153 18 1 0 0 0 0 1 19 171 969 3876 11628 27132 50388 75582 92378 92378 75582 50388 27132 11628 3876 969 171 19 1 0 0 0 1 20 190 1140 4845 15504 38760 77520 125970 167960 184756 167960 125970 77520 38760 15504 4845 1140 190 20 1 0 0 1 21 210 1330 5985 20349 54264 116280 203490 293930 352716 352716 293930 203490 116280 54264 20349 5985 1330 210 21 1 0 1 22 231 1540 7315 26334 74613 170544 319770 497420 646646 705432 646646 497420 319770 170544 74613 26334 7315 1540 231 22 1 1 23

the asymptotic growth rate of this lag-24 sequence is 1.029939697 (FindRoot[0==1+x-x^24, {x,1}, WorkingPrecision->10]).  log(2)/log(x) = 23.496 ; log(10)/log(x) = 78.1 .

for our arbitrarily chosen cutoff for factorizations taking too long (20-30 seconds, depending on computer), the first term for which induced an entry to addprimes was term 5185, 216 bits, 65 digits, which factors into 5 * (32-digit prime) * (33-digit prime).  the 32-digit prime is 35883310089669970781089136839889, the first entry in the addprimes commands below.

no two terms induced duplicated entries to our addprimes lists.

the final term done, the 8410th, was 354 bits, 107 digits, and was the most time-consuming factorization: 257 hours.  11 * (49-digit prime) * (58-digit prime) .  the 49-digit prime is 1506291006111119117205211837954859047440030121999 , the final entry in the first addprimes command below (second largest prime factors).

goal was factorization of a random-ish collection of numbers of many different sizes.  in retrospect, there were better ways to sample a geometrically growing random set of integers (future post tykznwfp).  previously, factoring random numbers 128-256 bits.

the small numbers scattered among the addprimes lists below, some as small as 18 digits, represent instances of the elliptic curve method being unlucky.  future work: identify where ECM was lucky.


\\ second largest prime factors, 884 primes
addprimes([35883310089669970781089136839889, 153505251588962425237, 65638786155986809, 6335151847969010091037, 83760219098371044014977, 108580968225874351874999416991, 25148004015165244945393, 346874642039053903, 141905839647428613187, 155950072352408546336434929029, 31695074109394258271, 7075949561986955161111, 125483384255560494440256649, 409560837590847550933746593, 12473249520758578507143390773231, 2810264938578751397147, 17385895598475846559133959, 19006523816425424693, 16206094373800117798177637, 2501912890817549778663555882541, 71551955299832764348843, 191065982914609136798378183, 33042768649358639137672369, 1240706693728224290647, 81155427312671365848286210846009, 43163174962449998794753514311, 110035050983875374718606602113, 26307854174268046560883297205454731, 2420787406638495969223, 439423190482007024750843389, 162357542812059374124689948099, 4203003726922737668816771, 34877148612585284194278833, 191338612293976252796780621209, 66842928216881960549977905644131, 5161861367102510064212163797, 417892158488830677096314973979, 159300167253469173179, 107898937199391510017500103, 3016347922027349373419373167, 123666967923718613633, 17655251116100366309, 1203491122519025797, 4613211177207116570971, 19723720405744196003, 194878241936273643855517523, 40653974502953972903137, 94026581851621796717, 40988261561385405099932107, 18871835948325629576963893, 552622890722408992857229519, 1146734147003352508341824410891, 464218964329931177186826001, 19325463480002724973102829, 185131052392684573047972238889219, 704219834788168094084743, 360105654228634026854133871, 35144238726511396515401, 58353547941244794347, 238243554125224366348028382559, 39876464302417338577036103647, 1017476210558008476833, 203402490652059275985826939, 866989487436993264259049535061, 6608967975692051760082730839, 66973401804976882714049357, 650105045402896843952273, 1974649814742640599115901633328797, 5338046368915062952759217, 13531402087340150879216419315519, 743764443065750615064397, 2280015094453148244653, 4081421335632058138897492132181, 2489197589701671349, 6702390421178620686281, 1140370359494428247938622546349683, 8855293109537402154509, 38446847658423228196343, 364889503687309647501600549649547, 873741357802646848032187961911, 9752949346243296231811, 41308466177712936023944784254027, 85464556443785093, 361335558416328913013438957, 160690400832726349514142713, 23432601626315867, 41529121437393908664342499, 22261274218525066162222862907218057, 205391996226575615407, 1527655959295362914173082609, 661558506163603653397801, 5081184794696793315681503, 1944170353778077410960584901653, 60043778859191504003362410173, 2628170515316281810423480175599, 223678765464483589658779, 3639576972170554849508095601, 19902449794680151551443786796634757, 1540753250181253612070147716339, 239684517499807765054038961, 3374878201013512482766667, 3525855290745536644687859203, 284235995965269093828800099, 149579326541296101390489450329111, 15641810332203294019951, 173667712823578344209242231928131, 1924630846143174370239659452957, 2383004447765599568322454254397, 311944908026286851, 844023876970433615543, 2588266488728386930818224878283109853, 1549375183978647196943742959627, 95961836785874982136067453456471, 50922135911096107408117049, 47791666744379988463510230754057, 42710349504521258305820052211, 27114912709565088403948453, 3188221413877264553939412099839, 131326053847804569152758575330791, 385757482876481726501, 90102216412371430926861941, 429146086666430796946308293, 530547325948505522193805171, 153856507570538281234181, 168857542989049482459103, 14683486239752177614829, 5161827257408654722511977057, 332462535739503406258419260832239, 33624436710443037576569095781899, 914795309714143746365951, 997494115250597025961961, 58344607021256458975011394207, 803998489011048627969735607, 172636793427747629005570277, 8846506736153187752655899, 442835399372669800254280711057607, 55631414814783489646969480388248496701, 1442444510961613852817990077, 3212165308500447730666515431489147, 5503137884444636306224433, 54166242077332034608129026095203, 4842893494044004359133915303, 7431165711739942869601741, 6420460859384758037260815890683, 136671282868298358848891025353079243983, 871462660121042952962743582651, 91354844304282074600448336292103033713, 347771791214109926109463, 54802097354616633434697539514209, 140022872608738886604592886641, 108937285167552944843455824352581739, 2384061519099881653738156849, 6879330955267671109103, 100598712949879180066157, 22265304911646173453571259174008372419813, 607814192622859935017515917229, 165581786573241387429781, 1488325233857100929288154811, 186423992743441023996049, 307730963774079817894031, 1735758232010784023888411, 166379280512578165294916153, 523139845158004481492863, 131830440065868720078283, 2973433744285732447910913703821551, 67061346279547306138445909, 136996063484078854337155563828218284229, 1037176865141656123713303004631, 885967681963076310749, 47922441186914985233, 22930859609123598922980077, 390969596197593615342653, 240320452086552278431597, 15958267839658131535801399372777, 286316911705026911992574788261930859, 37477253347495984286299600451, 11966109851192479625908149551, 4809808791257642410489125373, 1887627648082311436098139, 6435957759614472319, 211511854005771607555009, 85977932470638309842653967, 80843867967554827734928759, 89796451811415988626490175581, 3897813258342946773825988012246668143, 11269415900608682254116893601301367, 30110300314646907503, 10034799545048863011286323693233, 4344318680657622521345221771, 5777991916510513627296759322921007, 4182001258276275642649039845959, 921561239704253584210429499893, 756661704960482446966063, 52235499951374127730627321, 128616786702360134938431153620200761623, 3373939617681324936107, 1896548581445534867480827841023042095761, 868643344872479348237211337969, 543734900819272346561884889, 22722627437188979271541, 157960573236092703995677067273231, 13892358884749929685668522889217754738157, 13909817747505572273, 47131452342179671180534373, 117088324288704235631009117, 562029965967075806303, 6323863428884693022162966773, 48343215099300749666783310287959, 6595574760921953501220751593457, 1446895022568931426117335679447, 106904929311577916095331, 173352915648297197175007, 6403471970134703903, 240652019453800653474494161, 470900224298406252365549487024553, 17550741045672163676579, 166644645131221497203, 1539937318644616316507, 83027490999511593555211, 79444446773167411889, 97055550232347657379916647, 11480392486727455456717, 166371821130542178114452665073, 132920420876154559635518160341848931909, 10463703112223400390475463989487, 1067823602987332237247057, 723825574673307664484750917, 109141706745011649809, 1939213815838953321329079256446547, 50463982824348025417253, 12678119149100743400451349, 3625720310314251085103, 6327081263150458177705754921507, 14623520896916829675866983, 2730113561836606187697853, 522923328881073510167741347722851983, 56236846768392642046665013, 135854441250155471919063121607, 37797194180329014416429, 2304814157787285735588973329941, 4842209641580282719415699293, 4101692701586040006189713, 770365612728908836626312017, 11283709456532905026184339, 1153504277887012780393387296153361, 12187214811100302067110139, 43621816193965972890791994447017, 232681262031802957515298251288079, 255208943676890762518766462156947, 17600700964464556692385844263, 5301091324548809371285999609, 173974779718007120346017543, 6577623264048732178605011, 40649840033814571635173171, 2847126448036660093177521898399, 522345345881346732683, 209696457419579308571, 913878722749048980866709311250833, 1300278765965918280768991, 1744672074219865452295801351, 131775522164313275071325503, 104176952207104466475940790248177, 6408680102224474015017863111827, 21564917223530556259828417057619, 437496265771852169, 592859578708579260516902754193, 10581844894410039059967931, 1702247874065147040451, 63009656240536348144171, 6608449951425047711810639850220690669, 5645862913448749839432496303897507, 3467028549826624619757938362288091, 8234068779856298543634371185701254587, 34158311100191370467947487, 471958895439992338782642023, 59371065121384561678489797181701313, 2782824786840696841379707, 222622602099593970306998786222921, 12038290883270552281, 9176005280455208981671194863354197, 100513363315083059017951759, 149401849743030716317, 438456215459144864836253797, 117034103830863217833151294129039, 1475963780941091432839672594177151, 26836683751018181551691406341923, 60288709415185194278677, 439429147756871154648757, 32236938873448200208769398347227, 6743823091967885673561983, 7694699168068156129, 22530815998857279987782304144837857, 8801759897034751833272341, 420588224681033838971, 14568637705657498887358640653741, 5454545251824726757159, 374871220985324771269, 8863639101682509009538494878395607, 426555152359925631771529, 685533198952864122281, 132336208695067236264178257023, 5747762695605677631187644690521344229, 17862970951494907430852187698494198663, 5428755736101910605035301041, 8063440676826363842254971792697883, 17711848635873043111974647633, 4334182393626145798721, 637770386366142517223905438680301073, 63044715504247762767761, 9321767037697845683983688237687, 297579775708088499863411, 95759296264696105949, 180878909381375037271900533220116138307, 10037583285565567853760367, 1029933589893039071354707526136823307, 3258819054529941139264123403, 5031474200869084431805030287618671, 158664753903407378173, 456771276204267675428742825445369, 1422246560811626160584774407, 3282518178996641820461367157, 9389357041019051757176652862035697813, 9373147824580014750547, 333007201330915136966752223716595971, 999812166309283060745461909, 930673071840027943976585901546151, 420161481101669645723417015875169367251, 834354759244296226739115151031, 1048815696915202751228323845167, 32359639890055384751357360964162181, 2292553716056596592512558577, 24732087155106229312765831, 87098630216478436988797, 2209681081550049846574261, 2862716212203663780941, 121160653368080357412924760853, 76145329689486974244730714396459, 138513157184580209750285443, 50637181918065549146481992472091, 942467262593951323682716399046708522317, 580875894185119153827248699, 5070172343870660873401, 2124591069063730197943, 14199417523908927881444178421, 8883260271140956142679473371, 650417804596071044263081597312153, 12307108246639242300867776293, 201832153502026487427656509, 1083163751023366157699923573, 973747486295626729978616572523, 317487378200027218633, 225090236867882298559, 142952508522288220589683, 95758119248509685870839, 14354029427276297314847, 330156450099310307299, 12143401308482790179, 5958558216922879485931323539803115083, 355716003500919416423466902107, 1290718568930333315298660085043441, 686530451677268829357880891, 408936469005017921029, 175445724709201112857392146514424148239, 317360002351425484530121960189, 2950413490938589016486861, 339746507922044224546960498583, 64167532843269930242571551, 199850586871860672889, 5479370602154804442595391447, 99577804230583938402278717, 3783417794327744831182119132503774401481, 13716906991084595101262445575417, 210423165484804063078512118483, 287070536075729935712849893, 33354392518318350723403854648143603509, 4455416911295415107296760509, 6220744783686863373582165802351, 77078322001219834207453, 6319593289955487901664395373, 33041929410150928981527384029964716674657, 7491173285235961568422360237769, 40445781138099519448309, 5947825634217541676603317, 28821651286846446851399, 24361789764729067653204987377, 813905048046717450176200484449, 1266288693403558091384040799, 46196879599825181676896554747897, 9675909456748904779033291, 3262947943349521498056870257141589437, 464894071046782109854568374267523, 787694135491846132487064899, 273348465798208158621011050639, 33761132480132651221, 71894209226574036464690514313149437, 104700357430550162423, 246267016315135628104811081411739309941, 7932016655177758283753, 30926614350445225564440867947, 213412037011161451578860857743503, 63712303701792202235626061, 2447571327002012697723534023, 834822597761616506778361, 2241468911765114607311511229271, 3241975471471717141243193123, 20055648608999842837709, 4292576066236367802776131, 268216279863408780232922633, 2179845071386053446916529807723, 18989063737458776761015808985667, 10121521691862486959, 1094676847874303054029, 302777780803178373541, 28800750111265828110469, 47970993386990426958633227, 710434704998695207684709, 12630594675776662675597295556427883, 202108642204054895996889509, 126515872780531077203, 200938574198573328778685263, 11868250454875698311401, 748671626217620383954421033014079, 123780594388332199960026850828177, 663998937632753439815992754213, 17688711174421497924612277, 114452888743843605591791159, 58560103926017224782367588479218509733, 14104039112867371127557, 4367498550454211613743581, 327463917139191581351, 249854599271727784716539, 9248157882323843825522814922284429843871, 34718323748436568929847, 3400929577999714979821963021, 767203472585678935362889457, 5169672101606468219233, 1478772771057941641802338873050629, 5232203996916638896847652534230516642773, 1495708640030689696950708484847, 705343634020427845793429, 3515266964412990282553333, 2955116076892047231377679723653053, 132337729139126579361645280991, 1140127820828077387380643, 19218820615630093835921639, 93975735375425968461619, 13724020251335342734703, 30048255131768371208475820521424223, 39414472450837838033126977, 155452384505225798812917586112683806060817, 2174219165323794083760137, 241039486068992879647, 73541441334664663463957859791, 3657014080460910853359676281731, 34058056914728753082328163206484119, 38453211012909205073865313517202205081, 14667392460080523302022515816909, 1215774185969629905977, 160115450119079578717519439, 2438656526051769909224378527907, 5687032416882714608216693912051, 7241398218724059139549, 744767445662936292477120881, 50993158779209544870888007051991, 18446968168906556161079309606509, 53205249890353828002727, 10894901283148888892377881414088559317, 9963330594642802268261, 107662442665345144021882672630363, 689038151847807881136131, 226748687866130166078888873075981412751, 615014450020930276066231, 300116751465561367356592718609, 3049236577429871243542506529, 21241446614677454783014754812963145822707, 34612428299667760748034535834051766971, 8945371929768138545281502361651489436729, 139210028620526096050263716917419985157291, 6503094612087367469808631, 1383545790376939260864804752113, 42005547272261738326160890569683, 1789209225788028310259914947107739599, 112715641529106393785574028811, 842384872732152736367175367, 1540923207388513909664245897, 2383091726971524402401162173310506595441, 28850362865111518696088955780737, 15359457065194110825302994791443, 35419823197764527015407160804871103, 45468642822334328343431, 607295578125194987330394705105134869, 2073568366715474536915096693, 10321085044981906904916154080812791349114461643, 1797699677598692443, 110340206313961643020838311, 31701410074045099415941, 13249191969614692985173977527827, 553785851101732760851332159119207029, 85665976695938015762484802297602991, 59304014827032364000127151075572386600847519, 2145286856013414965596089702407, 280457524057024747269223, 22245977394939118613993, 23694442753702084299039843052091693621, 122918202052205553178652531483717, 41824093802372966773057357155467149, 482235364613000399157261886032931, 9914439018165582596983179228690226793, 5146009553343023936351559657031, 556526323268077048908877961941750529, 487754558599575542381636864894089949, 55926805954186661041919857, 228863378704196159122290436957, 18484328114239734418212840626759, 117324538889764319441722300171, 1243782121301819789919616896108712253, 91030606596409632663597697025716981, 34187657545925254134677005160267, 41686962513790625138159233038748874066456383, 3258985795108900024912551634190417, 1267872172695978995825029253, 283351307719354207558481493530567, 44276063218875697089587843, 198398477735534078505129140099880556675506523, 158914403316435703196240904793381, 362636313132679879361562174594471121231, 879016894061064634891, 3858956474618975249989037731, 604405980852003261467579, 28379781758555540603021, 12805691260888561831249, 2226478802180475967142813, 1712258164121048693271073273183, 249485610783174500208956071606173121, 10239191571190919104657782461, 275387689232037472479601, 40666024844846045549712251, 1160758916984646331184622440373607133, 184215447715080621861121, 3644537933897601908374972087672897253, 92270856956565660424424025196436623, 6493130270946749249290972301, 1044616255730583331086136036505553821423, 95640255077895957269, 1341997514284766314388887045013174350621, 189300005241544458558181189, 23887504667865050910394624466341593166048433, 20322397076372690314253040958413277, 51585626763264606208594393, 7482392584187176019344076948223410769228503, 3119853068049990137984751940903, 30658874852753997359468718961, 4378356310712785710947203, 26546682305834264858839, 335896248798402738442571857134991, 16706511537336029339494843182770863, 23536353739460190937921, 58823183112031618446436269497, 176745274277466060775034437807300730321809, 48824302062757521545863801693, 95680915593387144745519039741546447, 43284071144949113834512925563, 1460470448368792843749077, 18197044237771816643117841476787990957882013, 697530060642379822035327494694181026023, 25616363804075521368386884292677, 3827599811784963639671, 171426499943854180340650196421421, 5700907983495487039687134533, 21994865419804958711026155826265281, 518177385546238771436219129947801, 2246896320935945629983509487036941, 1090607885919099265593008610661, 745248812539609488941, 362111833149082027509581, 471843820960541864524607, 285844085621094728799792333040019, 705501812286125536807947146003085950507090311, 1994611241565664930007, 52479913459488298075375499873, 170384895511874274865502599634754710437, 12009770237368173686946212203, 330960110137528821711211622659, 17935580932712022092662015227948563, 99654320441123231396059, 3721712480894499473161385902909, 132160002291672263083592925921113, 23883381259628552391599, 5302376876482059737428014337, 3536115623713258718227, 64978287251079088194843113932511128561466725973, 54758498872712708403217, 13578580319915743512190171, 7460917382998610900739808953701, 2047469506542764166951467, 8099827783370005667764301674511783, 57146068234154202665589505471100225604216421, 50917517315032731017, 242275401279836021227166992063, 275179932166104342662502569581297, 3282954612906649452231545023, 11092198042916866055144407, 29859092142133077216680623860601544232230537, 23853854325102529400915442783901, 440023165834835603530825751, 7001672425345324493959716092747, 3409325938371426514610100760929229412473, 253017641607238054023101329, 4729796119048578109867089765629683, 46756143082188976603031, 145291431293074766261, 8443418014194995113314349, 95615103208991966817898922347, 242686715066615174080422200941683772759, 2841862141747580495835106447, 403715854797958610907416711001721078411, 21365350977750811026070697, 5892196098277717080323, 10180023351079538025816250157831, 413614689897281866110689000195530492829687, 132496174723892163555792581, 498524146105634689991, 4521631426789830359774384821, 16914848621014053257401920597956329, 9641360917013256154987, 1928773857286678004170880387, 10662053909072513524848212228146816531, 1595087455974410144495171, 65776300234788772737358045160123443, 7232373042702871761759319206885094244503, 88893213548962948273132647879445381, 669278664781329239925826636402631, 4040444094299805826290516990841, 3444465702970145616457551387547, 1778769258081082004642116002536233, 348622083328871162021466089, 243268318015924573857547, 178392372518343998527814204569241, 2276034401389936251271, 185056110597579550887185833, 37367791973804673855530444441026969, 47561762328005391998687, 192311224492171856838750541, 33281890583034500455123, 3448690313313333590127968225071957, 497055834555138001708015269549774553, 92036706373679388693541156532437, 31269511572865748745377221579828867, 21431235855863949013037401, 88012986769655074055100801649271, 6672403474691001320559374197589719, 1198761917797965518899616696383, 116011435213066457261496693935975881253, 196771616692759324126843397, 10793685802847240093980949163181, 99887996776774987076758103903, 66670678615405677884425976200018351, 1429313583029781414220277657, 3861931023219513168049255570625227039, 28016640347729588482888553, 22453700047806142084526049473, 4335253345354436113392067331063, 5770813531690837224583386090048725221, 334932742798102448839197255618828975809, 4228195878835690750165779588168177478699, 113408930541705588109, 184015648212483069923531, 6994242278447475610655287039074678121, 30740606510391062938597470206999, 887421333089160232977287586013, 1928248625571209790661132159822799, 85343713010444188037110687, 50391478105170531547814063341, 14726460841684052503499, 99411093116216924986553, 45275547184932137873, 4682541688241014548637, 2216591925775169656421, 132102071960676330181, 943881373525768986803137, 18512775089416124012270092833787, 6059938718036085674767507802209, 151127955448485357405475906242841, 45969653725875793062704609, 477391130914957325683, 188213651289115112382127831877, 78114336155263841715619, 42463553769011608440763, 36106143139899513297522209, 68182872224174253040333281228974249020031, 28673652297516407366467, 446914759377569531551973, 38915537085953618973437, 18761563582187813647873, 16540139295362323138404691826712141134959, 133763001561943272983073439, 129913272182381414042211393765999096221, 62827928068709627329792667912227999, 2411298607046022759618727, 483773564744359950827572009115506601, 1619363184044663965603187911, 1864375622899642961069116492487, 193446815671835971415785991, 3337019093464932131859152229551393, 110854707743050267513793075250352961, 240232945778268907818414853841, 45174414887868935951780323541306407909, 183250064005228678926778542240482567, 340668288853214994103, 441115909446244662303209, 669686450073007411779127, 470085910720210199, 291006691714107447261669058193663708653530274159, 2257960311107943596177, 280709696477196049490034510376589, 1051298230603322717811623, 10394117338430379881580112183, 1910989083125104005004533653519121181, 2852252704886479447859618359, 2707293809864050028427407229392714114267507, 900886072950553047663789519762700121615107, 1262793052896800830393, 4590799086510476192537799446257893934010284852043, 44344904077234400977846507046642779, 43761129709959613608559559109455573812051, 47488611431494925174166789421, 26185140637728096100308103, 146642948098078226717938916623523, 34734756779404552630081716437, 365693405471529534563413, 918176666435578331544013, 4089416398734577672586282782787, 294740623099171751378957232637, 92980529964224371317717171862255314703, 25285508867908397408173, 331368882893583882514625835068713129515428098897, 402464418667689907202584203568651, 5450363361834747321617750342927999, 26530103926165834079102239965811, 542694418781450873617054849, 859968008523703555659444411968927, 47912662861480954625290726608586882577, 1726642344460471936312057212917, 77086415335294161998689, 4714512363849648934715517239, 113378155868574217740400405783, 11430512015022098279189259857437427, 36975436356925856640589, 3135805051489236957236477058736994307475717, 11669205397311561954056016463047652281599, 15481587445311809640779, 3933343539581316377008651321, 1685933908892909420376741598485176833511, 4621423274999920547548653788652899, 3833845252580373085118414348843, 42374943967615283852168070778684499, 71088058755428779681621789137657790645389277, 460806070492490393442576850973, 1171288899434457715309, 713556177743845631473, 65140521115272692352017689417387, 3661438546197358549687, 10246859229574848923613663353424491, 170265183833584770152995968398163839, 325846353879104261267558729, 411254209564978615894459, 29268468521387536791722527474698451, 3832057042868911600729040423, 5446061638611623381879, 616237238103693253361061002067563, 1039100872501128470865511024867, 669470205434048274882320881, 81023339233654089030761, 1564880648894348187028910791746980641, 8586370948563916871279062673694421, 95397859932945390342841723, 157420292457894276506851, 81846466611343796335241, 8925451425576790968375857522865480659, 87906306604729634731, 237773965120107416809, 1632993064891178855088229389948248678780399, 19443489010601526982943622972930019604260972271, 14009144301213225731175397721, 52019192895917833872209846396228234301885241, 226433244071099413468358332127213, 205119435484085499896369, 4617825207958447300285445154179376239, 123701537449403743263879664579, 52621668071905119582799, 5424495120882308208233884725067, 1284284337597829066447673423327, 296953480547178643445024185795112003, 591059507831414337806369, 253058113968831958525909395046206071563, 18807665488255615334294856829, 36739339654727301628030968702107, 1841208994616446056773, 27529569895070550460909667308225040351939, 83656913811004684871624833385601011, 12676048600137825786157, 33693655847906055951066350209, 78217354407132132720340890638005624816501, 10413791844382975050472705799, 489756852757918931396600106383, 51617843744565406048574047, 21014795823349263918049, 85768195214233239229809211, 1639050980447232698873396060571477310060468007, 58541388091004642091632500999398495437, 14122759423763556996418893061, 3830370097247654971150939452070709, 58263358895048138064236300898863, 1138136195974622931982751374973177, 27689403202380612862163217195946579, 2577907700164135714557114475205023, 13406563100657664173332909170509, 512215538227839479102694826854599617817, 58382528900413037074211, 95071703006261664423565748766873542357636597, 83665856147281885541919223, 320890778199943929231800530816639, 58493318988935277271241842215667, 93762703477799393010063981959, 1155611972401101700553, 27420324108264430088657767, 5599252503429489200364159471881, 575688050841903769182707680553529548743383079177, 115147422339996728173181, 17157257053530970825133, 231393386099920101772871507753, 31621847337480584576399362699, 17372424578650511804581317847, 1827967032641546651530156064637586783416227, 7453805052135432858678662352439, 59034202962431199923880242870529047, 1800090042559914778509923, 522807506260091624001449, 78164392953914761288865461, 56785567839991930609201188978666886524638821, 140502292495105886749634811397, 63955606112722155804045574677571, 29915051525468753000146731983791, 55136776122396999370101949218953, 51151920582744324329837498203693, 147570604311650625961635132359, 97908413726748317466956964701718035513, 442705757403837969832639936861, 9626785055476375982554046011060091193509176186201, 1523508277303890307680598295214198331796057098057429, 2228242700518039505481930114338907297763, 2295243319436406424136156216395082699, 21547340218502606918532608786669, 1020012110483594703106395821, 65689114086358895220347007687684269153131, 214572294771151294179372657857, 1234472517878661087305951070433707958699121, 41985198187351888134792316062487, 31529165404986041605281190913, 11369576702678545997280641642947935714586172006131, 1892837128501798979681179303421101, 102097740517890636222563404542620839920824363, 11761053241255319249168539979, 712557296843400168243454242367733323, 18427409688800345789, 2703988866742198762688823133, 362251876259637480414446519076330831479924634287, 69285730336986731641, 28227473468812082927261674754276123653607, 1455606615703131218389577, 104220107601997450923569, 105121769145431743652851, 1366952266934703697361, 63891103315135756401250319, 3081436746224903550148922544047286765639527, 424035699034895696813, 238249748208853116847627114283828076485769667, 2552398005399309046831, 1238426543009315467495119293, 10683567036703387214158691006010799561, 24124725763648233170537253509, 1069112251595276709800023561341578335445971001533, 357990206974732629654887, 1506291006111119117205211837954859047440030121999]);

\\ third largest prime factors, 44 primes
addprimes([59576718681627821, 98501123617245949, 138598798273932761, 5408979750324745901, 6462339564524872397, 20652875228488335017, 9163858315338862231, 356120421609113092153889, 6932943755493033953, 3354908076703056902939, 127948735601176844742109, 44161085023388023319, 225215623122648550052863, 20398034995569724763, 20611044148501285500533, 97586618754794392073, 26002217302461336127, 7785570944466431887, 8247006860941500257, 1655733203483381773445239, 2953069114940846260159, 374514648625536459379, 15059370546976585199536492065499, 643928788451607137, 33014245889922881, 8257029439741656395243, 7203813202136662187929, 1571908746467275931927, 94531644301368913421, 699184400849008886446013, 13525372779622820625647803208993, 93700146536472538850512807, 55253357281568457877, 464529691551452966431, 61962374346345162428939, 2973647431499584536625367, 1091823859608267013505531, 58349173047963275430283, 51694948793834014911443, 2021225824713875422650053219, 195067149129775878515527871831, 349793956030892222094153407, 4126129372002305119, 24511366861314949223]);

No comments :