Monday, March 23, 2026

A random 2000-digit emirp pair

Print[DateString[]];
c=0; While[c++; PrimeQ[IntegerReverse[r=RandomPrime[{10^1999,2*10^1999}]]]==False];
Print[{DateString[], c}]; r

Mon 23 Mar 2026 16:19:04

{Mon 23 Mar 2026 18:37:09, 305}

14855218387530938731264800899806518554251409448753318067281617721104517975019660570816502191121802141387619791604557638640603077054520997318307203718761757739257990957393586487037638214729490346907077085366489735302991920794042479840185960548800113792237378778189700214336498150535494706550605947647886952709343445154850786791334376804237058351294799668452160206131345032141325576891497636422462012076224432607954324056613347607288129441043897934814166919553140730183934337952784895386915184278773673374499327984879217657567333750360009291505056111578053084511773793397248168145190194638081779158080985373127557284411266286462358754778798450217253924795704562031795767586150594858361787012654211320033792723983297048174107728830987370130777266997038517867493302588134280307951436792574170079530221417396081625730302595728237532318481405855540975409890906976902320267307459218191770761431037983144020890302946754575512111129239578508080696921714584513835630889334839362187530001616420174589110874531666851407830633657753457085270606608334710585505853591828745779378101187982074326317494916222176946818501389870085531903573860851057307085813978075441446820544507426572075052156092931837527279166150509917972776167153657340031646179631733906184205918255895475754364175238057482645899534235102077183676796749323690141704894319808779172119291942162868677765137547565809194418652476017816016763486268007290212946536413065587272480288193061788482406655351702586154143963745149534333768316994754590468906962685037527510130840705726077739498669858651281009519862797180044407205030051742292526137434721833180012034655794438599828098375237119761007728650703770787960212849241788574933755087711157262236337487014829313684940715849983628697605605634086992019999606613109048737110152003061306410080888361354585657807374315692388968946536330992803388106574848702486355828678460394253967147711575515270818557921500146950778382811866201918398021331397625888200120115790700195112386497426010970976957364112704151652043

IntegerReverse[r]

34025615140721146375967907901062479468321159100709751102100288852679313312089381910266811828387705964100512975581807251557511774176935249306487682855368420784847560188330829903363564986988329651347370875658545316388808001460316030025101173784090131660699991029968043650650679682638994851704948631392841078473363226275111778055733947588714294821206978707730705682770016791173257389082899583449755643021008133812743473162529224715003050270444008179726891590018215685896689493777062750704803101572573058626960986409545749961386733343594154736934145168520715355660428488716039188208427278556031463564921209270086268436761061871067425681449190856574573156777686826124919291127197780891349840714109632394769767638177020153243599854628475083257146345757459855281950248160933713697164613004375635176167727971990505166197272573813929065125057027562470544502864414457087931858070375015806837530913558007898310581864967122261949471362347028978110187397754782819535850558501743380660607258075435775633603870415866613547801198547102461610003578126393843398803653831548541712969608080587593292111121557545764920309802044138973013416707719181295470376202320967960909890457904555850418481323573282759520303752618069371412203597007147529763415970308243188520339476871583079966277703107378903882770147184079238932729733002311245621078716385849505168576759713026540759742935271205489787745785326468266211448275572137358908085197718083649109154186184279339737711548035087511165050519290006305733376575671297848972399447337637787248151968359848725973343938103704135591966141843979834014492188270674331665042345970623442267021026422463679419867552314123054313160206125486699749215385073240867343319768705845154434390725968874674950605560749453505189463341200798187787373229731100884506958104897424049702919920353798466358077070964309492741283673078468539375909975293775716781730270381379902545077030604683675540619791678314120812119120561807506691057971540112771618276081335784490415245581560899800846213783903578381255841

Click on either number to check its primality.

Saturday, March 21, 2026

A random large emirp pair

One might think from my previous post that large base-ten emirps congregate near — or at least involve —  powers of ten. Well, record ones certainly do but that is surely an artifact of the convenience of searching for such integers in those locations, expressing them without having to show all of their digits, and even proving their primality.

I recently found that Mathematica has a RandomPrime function which can be configured to generate primes with a specific number of digits. By repeated application of it and checking each against the primality (most often lack-of-primality) of the integer created by reversing its decimal digits, I can now create random emirps.

While[PrimeQ[IntegerReverse[r = RandomPrime[{10^999, 10^1000}]]] == False]; r

3279288520829477235642879267287999159015788532533738855066046440495522398413879852102369701530100869580838076083149157548206285059151203914266414912431580164820025933969565023665242290427065104497454406152915882315883324964117545072539206755054402862798019415649139491569247916954175149446689650326543522910485350967102290630555006513883295579511863580034048864302120190250042483217978026251237207016901362308898161048199473494827584506196053490289923397647958362118557321048941198888848619470770742015547396263817662528411868482723797066437421483218035169802940903987225108903574471943089572864902145297776096342849868617007993544349744145780223402039236076990018520939944400626427281300176781542771736168164872520396475344293444882067084862153824121702601618059184428161120858343524815906494130051717507462817874659446095816309231120150861019479629358822063445459770400310760092233488204681637298530785169163900358943079863912831713129072815997830914582590098549219524593979172161769855800664064007

The reverse of this is:

7004604660085589671612719793954259129458900952854190387995182709213171382193689703498530093619615870358927361864028843322900670130040779545443602288539269749101680510211329036185906449564787182647057171500314946095184253438580211618244819508161062071214283512684807602884443924435746930252784618616371772451876710031827246260044499390258100996706329302043220875414479434453997007168689482436906777925412094682759803491744753098015227893090492089615308123841247346607973272848681148252667183626937455102470770749168488888911498401237558112638597467933299820943506916054857284943749918401618988032631096107027321526208797123842400520910212034688404300853681159755923883156005550360922017690535840192253456230569866449415714596197429651949319465149108972682044505576029352705457114694233885132885192516044547944015607240922425663205659693395200284610851342194146624193021519505826028457519413806708380859680010351079632012589783148932255940446406605588373352358875109519997827629782465327749280258829723

I did not think that I would be able to prove their primality, but factordb (click on either number to see its evaluation there) apparently has some "elves" who download the smallest probable primes in the database and run deterministic tests on them. I guess I got lucky.

Monday, March 16, 2026

Some recent large emirps

Let it be understood that all emirps come in pairs, say (p, q) where the number of (decimal) integer digits of p and of q are identical, but p < q. Since, in the following, we are dealing with record large integers, I will explicitly state the value of q, the larger of the pair, followed by a linked p in square brackets.

In 2007, Jens Kruse Andersen noted the 10007-digit 10^10006+941992101*10^4999+1 [p] as the then-largest-known emirp. Eighteen years later, Stephan Schöler managed to up this by four decimal digits with his 3867632931*10^10001+1 [p]. One month ago today, this was highlighted in an episode of Numberphile, precipitating (of course) a couple of new records:

Two days after the video, gamer Gelly Gelbertson found 10^10056+10^6692+10^5872+1 [p] (a 10057-digit term in OEIS A393530). Another two days and Vishwath Ganesan discovered 10^20000+518406362*10^9996+1 [p] (a 20001-digit emirp).

Update (March 19): Serge Batalov captured a near-repunit 10^22822-10^63-1 [p] emirp!

Email me if you spot any errors or have something to add.