Wednesday, December 28

The one-billionth Kaprekar number

When I published the first 11 million Kaprekar numbers on September 1, I asked what the one-billionth Kaprekar number was. Today I fashioned my answer. This was a troublesome compute that included a six-week run calculating nothing: I had used in my program the old Mathematica NthSubset directive but had neglected to first load the Combinatorica package in which it resides. Ouch!

Sunday, December 25

Digit pick-up

Consider the multiplication table for the number 125. For each product, we will pick up any digits thereof we have not already collected and add them to a growing list. When will we have collected all ten digits?

 1  125  +1 +2 +5 => {1,2,5}
 2  250        +0 => {0,1,2,5}
 3  375     +3 +7 => {0,1,2,3,5,7}
 4  500
 5  625        +6 => {0,1,2,3,5,6,7}
 6  750
 7  875        +8 => {0,1,2,3,5,6,7,8}
 8 1000
 9 1125
10 1250
11 1375
12 1500
13 1625
14 1750
15 1875
16 2000
17 2125
18 2250
19 2375
20 2500
21 2625
22 2750
23 2875
24 3000
25 3125
26 3250
27 3375
28 3500
29 3625
30 3750
31 3875
32 4000        +4 => {0,1,2,3,4,5,6,7,8}
33 4125
34 4250
35 4375
36 4500
37 4625
38 4750
39 4875
40 5000
41 5125
42 5250
43 5375
44 5500
45 5625
46 5750
47 5875
48 6000
49 6125
50 6250
51 6375
52 6500
53 6625
54 6750
55 6875
56 7000
57 7125
58 7250
59 7375
60 7500
61 7625
62 7750
63 7875
64 8000
65 8125
66 8250
67 8375
68 8500
69 8625
70 8750
71 8875
72 9000        +9 => {0,1,2,3,4,5,6,7,8,9}

So, it takes 72 multiplications before we've bagged all ten digits. In base ten, it turns out that 72 is the largest number of multiplications required for any given starting number.

What is the largest number of multiplications required in base b? Tomas Rokicki has worked them out for thousands of bases: It's the second number on each line in his list. I felt that these maxima (k) deserved to be graphed. Noting that in composite bases k is divisible by (b-1), I decided to do another graph of k/(b-1). Then, to explicitly list the k/(b-1) straight-line points, I did the families, with a question at the end.

Sunday, November 13

Meloe yellow



Oops. Just over four minutes after taking a picture of an oil beetle at my usual spot above the river, I accidently stepped on it. At the time I thought that the exuded bright-yellow mass was liquid, but an examination of the photo revealed it to be eggs. Also visible in the image is the cantharidin-containing ochre (generally described as 'yellow') fluid exuding from the joints (detail above). Samuel Maunder's 1848 description of Meloe is here and a more modern treatment, here. I noticed my first oil beetle in the local, private cemetery three years ago. A large grassy section of this cemetery sports an extensive covering of ground-bee dwellings and the beetles have taken full advantage, being this fall every bit as bountiful as were the bees in the spring.

Wednesday, November 9

Scientific American, weakly

BoingBoing alerted me, last Friday, to the free accessing (this month only) of Scientific American weeklies from 28 August 1845 to 25 December 1909. A weeks-between lookup gives 3357 potential issues, but I think eight of those did not see publication. Of the remaining 3349 (that's four more than claimed by Scientific American), I managed to find 3246 complete, readable issues at the site, some hidden behind missing or misdirecting links.

I have been an aficionado (and eventual collector) of the magazine since discovering this issue near the end of 1967. Having a substantial, freely accessible pdf-archive of the old weeklies was, for me, a little like finding money on the street. I spent all day Saturday downloading, as a result of which I woke up Sunday morning with a debilitating lower-back issue (it became extremely painful to maintain an upright position) from which I did not recover until yesterday! I spent all of last night (until about 4:30 this morning, more carefully minding my posture) completing the heist.

Wednesday, October 26

Felix Aloysius Sullivan

Some three years ago I mentioned in my old blog, Felix, a relative of my wife. Today, I received a surprise email from a David Wilhelm about his relative, Kaiser Wilhelm. Wanting to share with David information about Felix, I did a search looking for my blog mention... and found instead this wonderful, previously unknown short biography!

Wednesday, October 12

Zillions

In response to the perhaps-little-known fact that the last prime in alphabetical order is two vigintillion, two undecillion, two trillion, two thousand, two hundred, ninety-three [Donald E. Knuth, Allan A. Miller: A Programming and Problem-Solving Seminar, page 12 (page 20 of the ftp-available pdf), June 1981 (the solution dating to October 1980)], I recently provided MathFun (28 September 2011) with:
____________________________________________________________
Bill Gosper:
"Yeah, what about 2*10^4679+3, (two zillion and three)-?"

Phil Carmody:
"It seems to precede 2*10^4772+2003 (two zillion two thousand and three)"

It's instructive here to look at the original (without using 'zillion') answer: 2*10^63 + 2*10^36 + 2*10^12 + 2*10^3 + 2*10^2 + 93 (two vigintillion, two undecillion, two trillion, two thousand, two hundred, ninety, three). The 93 is forced because it is the only two-digit number that will make the entire quantity prime, but 23 (twenty, three) would have been a more desirable choice to force the written number to the back of the dictionary. Just 2 (two) would have been even better but a number ending in 2 cannot be prime.

So, let's prepend the 'two zillion' (2*10^n) and force the presumably-optimal 23 ending.

 2*10^n   (two zillion)
+2*10^63  (two vigintillion)
+2*10^36  (two undecillion)
+2*10^12  (two trillion)
+2000     (two thousand)
+200      (two hundred)
+20       (twenty)
+3        (three)

Of course, n >=66 because we have perfectly good, widely-dictionaried names for numbers less than this. A quick run has n = {77, 113, 116, 158, 342, 464, 468, 565, 2171, 2274, 2340, 3347, 5724, ...}. That might have been the end of it, but I noticed that there is nothing to prevent our indefinite 'zillion' from being repeated. And every repetition forces the written number further back in the dictionary. So one begins to look at the divisors of n to see which combination results in the most number of greater-than-10^66 (zillion) amounts. For example, 10^5724 = 10^1431 * 10^1431 * 10^1431 * 10^1431, so, taking a zillion to be 10^1431, our large prime starts out with "two zillion zillion zillion zillion, two vigintillion, ..." Our n = 2340 solution allows for 30 written zillions (of 10^78 each).
____________________________________________________________

The term after n = 5724 is n = 39960. The resulting 39961-digit probable prime is now on Henri & Renaud Lifchitz's Probable Prime Records site.

Thursday, September 1

Wednesday, August 31

Spoonful

David 'Honeyboy' Edwards died on Monday, so yesterday I listened to a bit of what I had of his in iTunes (mostly the Delta Bluesman album). One of the songs was Just a Spoonful, which brought me to list all of the 'Spoonful' covers that I had: about a dozen. Cream's short six-and-a-half-minute version from Fresh Cream played next and I was reminded that they had a long version on their Wheels of Fire double-album which, back in the day, I once owned. But the album was missing from my iTunes library, so I downloaded it — and am listening to it today. Catherine actually brought into our relationship the Fresh Cream album (along with The Incredible String Band's Wee Tam & The Big Huge), acquired back then from Kay Owen, her British cousin. Much appreciated music, even after all these years.

Tuesday, August 16

The largest base-24 left-truncatable prime

10594160686143126162708955915379656211582267119948391137176997290182218433

As previously aspired to, I have now replicated Martin Fuller's 2008 determination of all base-24 left-truncatable primes. I don't know if Martin thought to store the largest of these as part of his program, but I don't see it anywhere on the Net and more specifically not in the a-file for A103443, where it belongs. It's the 53-digit number {17, 22, 19, 14, 19, 15, 10, 3, 10, 7, 1, 13, 18, 13, 9, 22, 15, 22, 23, 15, 14, 3, 13, 3, 20, 19, 17, 10, 6, 1, 20, 17, 9, 2, 18, 15, 12, 10, 3, 21, 11, 8, 16, 15, 4, 4, 4, 15, 11, 11, 7, 10, 17}.

Bookends

Back on May 7, I created a 9 times 48969-digits-each, type-2 Belgian-numbers html-file to elegantly display solutions to Eric Angelini's second type of Self Belgian Numbers. On May 17, I published a list of thirteen type-2 Belgian primes plus six more probable primes, the largest composed of 27867 digits. By May 24, I had determined for what number-length (after length 1) all nine type-2 Belgian-number templates again produced solutions: length 1899283.

All that needed to be done now was to find all solutions to length 1899283 so that these might be as elegantly displayed as my initial length-48969 set. The calculation took about a week for each of the nine templates and, on July 31, I html-colourized the result. Only problem was, the numbers would not display properly in any of the browsers that I tried. The far-right digits did not align, even though all nine numbers have the same number of digits and they are rendered in a monospace font.

Thankfully, Firefox 6 has just been released and it did display the numbers correctly! The length-1899283 solutions turned out to be terms 3594728-3594736 of A107070 and because these nine solutions along with the first nine solutions form "bookends" to all the other solutions within the nine templates, that's what I ended up naming the over-80-MB file.

Friday, July 8

Kaiser Wilhelm

Just now, I fixed the birth year and the 1921 stated age in the Wikipedia entry for baseballer Irvin K. Wilhelm — yet again! This time I added references.

Wednesday, June 29

Wednesday, June 15

A076623 decomposed

I've been busy with a number of ongoing projects. I'm quite proud of my "number of left-truncatable primes by digit length" which I published yesterday and extended today. I still want to get actuals for bases 73, 32, 28, 79, 55, 83, 34, 89, 51, 65, 38, and (finally) 24, but it'll take a couple of months.

Tuesday, May 3

Aliquot sequence 46758

I have gotten back to pursuing the fate of aliquot sequence 46758. On April 7, I found the 53-digit factor of term #3237, and early this morning, the 44-digit factor of term #3238. Terms 3239 (this used to be my Beer Store employee number) to 3249 proved relatively easy to factorize, so now I am working on the 113-digit composite of term #3250.

Sunday, May 1

iMac interruptus

There was a power interruption this morning. This is the first time that my wife's iMac has rebooted since I started it running back in January. I make that to be 100 days.

Saturday, April 9

ABC

Three 3-digit numbers, a, b, and c, have no digits in common and a^b ± c are both prime.

Wednesday, March 30

444085 Ruth-Aaron pairs

I contacted Donovan Johnson and he was kind enough to share his list of Ruth-Aaron numbers less than 10^12. So I added the additional terms to my compilation: 444085 factored pairs @ 66.5 MB.

Tuesday, March 29

150000 Ruth-Aaron pairs

I am done with my Ruth-Aaron pairs compilation for now. I finished around 3:30 a.m. this morning: 150000 factored pairs @ 22.5 MB. I discovered that the number of terms less than 10^11 was not 106770 as indicated here, but 106670 as indicated here (I had not been aware of this sequence before today). I also now know that the Herman Baer Ruth-Aaron number 174024968568 sits at index 149492. Mr. Baer estimates here that its index is ~2050, but reduces that guess to ~1123 just a few weeks later.

Sunday, March 20

Friday, March 18

28, 11, 12, 13, ... (reprise)

Back on January 4, I asked: What might be this sequence's final term? After asking Lars Blomberg this question two days ago, this morning he replied: 3008889. Additionally, he provided me with all 1124577 terms. My own calculation — after 73 days of computation — had only achieved 481185 terms. Why do I bother?

Tuesday, March 15

Ides of March

I had given today as the date when I might have finished my Ruth-Aaron pairs compilation. Well, I have another two weeks to go: Running the third Mathematica did take away from the iMac's available clock cycles. I've been keeping busy working with Eric Angelini's colorless green ideas. More specifically, I'm looking for minimal loops of a given length.

Tuesday, March 8

Thursday, March 3

Snow big deal (reprise)

I had submitted the photo of the squirrel that I shared here last month to the Toronto Star's "camera club" (topic: snow) and was delighted, yesterday, to see it online. Here it is in the context of the other nine winners. I looked for the promo in today's print edition but it wasn't there. :(

Tuesday, March 1

250 komets, mostly done

While Lars Blomberg went about extending the columns (up and down), I colored in the cells of our spreadsheet (on the basis of the ratio of a cell's value to that of its predecessor):


•     f(n)/f(n-1) < 1   pale orange
•     f(n)/f(n-1) = 1   yellow
• 1 < f(n)/f(n-1) < 2   cyan
•     f(n)/f(n-1) = 2   light green
• 2 < f(n)/f(n-1) < 3   blue

Hmm. I can't get the 'blue' color here to match the one that I used in the spreadsheet. The tops of the 250 komets remain colorless, pending calculations that show their predecessors to be larger than they are, in which case they'll end up 'pale orange'.

Sunday, February 27

A missed opportunity

I rarely go anywhere without my camera and even at home it is always close by. Last night I lent the Nikon to my daughter (for today) and this morning I woke up to winter wonderland: delicate heaps of snow on thin branches — already only a memory just a few hours later. I'll recall my maplepan captured 2 March 2003 (saved photo #52 of my previous Coolpix).

Saturday, February 26

Tuesday, February 15

Monday, February 14

Komets and planits, oh my

Lars Blomberg and I have been busy populating entries in a Google docs spreadsheet that extrapolates "komet" k2 (A038834 <-> A038807) to k6, k11, k13, ... k781: 250 kometary path fragments (and 5 planits) in all. I have had a (much smaller) version of this table as long ago as 1997 but it is certainly nice to see it updated, extended, and out in the open. Almost all of the really large numbers in the chart were calculated by Lars.

Monday, February 7

60000 Ruth-Aaron pairs

I have added another 10000 Ruth-Aaron pairs to my list: 60000 factored pairs @ 7.6 MB. I was tempted to run the third Mathematica on this effort — cutting my expected completion date by a couple of weeks — but, on Saturday night, I started Mathematica 8 on a different venture: duplicating the extension of A038807 sent me by Lars Blomberg last summer.

Saturday, February 5

Reality check

By this morning, Charles Greathouse and Donovan Johnson had calculated, respectively, the sixth and seventh terms of the sequence that I mentioned in my last entry: A185584. At that point, my own effort had only reached 2 billion. Ordinarily, I might have let things run — or started a new search beginning with the larger value — but I thought better of it.

A number of folk thought it was a good idea to prepend zero to the sequence — which, in my mind, lacked logic — and I objected*. Nevertheless, it was added to the submitted sequence — but later dropped, after I pointed out that the closely related A064510 lacked the zero.

*The crux of my objection was that the primary sense of the word "sum" is an aggregate of two or more things. Aggregating one thing is already a stretch but I'll buy it. Aggregating no things? { } != {0} but (apparently) Apply[Plus, { }] == 0.

Thursday, February 3

An unexpected bonus

Claudio Meller today asked Sequence Fanatics about an extension to the sequence 1, 130, 1860, 148480, 3039520, ... So I fired up Mathematica 8 on my wife's iMac and started a search for the next term. I had indicated here on Saturday my expectation that doing so would degrade the performance of Mathematica 6 and Mathematica 7, concurrently running their respective Ruth-Aaron numbers compilations, but when I consulted the computer's Activity Monitor, I found — to my surprise, and delight — all three Mathematica kernels running at ~100%. The effect — and I wonder if it is just the appearance of a processing-power benefit — may be due to hyper-threading.

Wednesday, February 2

Snow big deal

Here in Toronto, they certainly oversold the snowmageddon that was to befall us today. Nevertheless, our school boards — concerned that it had been over 12 years since the last one — decided to declare a snow day. The kids were happy, at least. My photo was taken 50 days ago.

Monday, January 31

50000 Ruth-Aaron pairs

This 5000 Ruth-Aaron-pairs addition to my compilation (50000 factored pairs @ 6.3 MB) arrived much faster than my previous updates, because two days of the calculation were done on the new, much faster iMac. From now on, therefore, I'll update only when I have collected another 10000 pairs. My current aim for the compilation is ~150000 pairs (up to ~175 billion) @ ~20 MB.

Sunday, January 30

A small premium

I have run that program now for 24 hours. Whereas Mathematica was checking 1.5 billion numbers per day in the greater-than-100-billion search-range, this search has managed to look at 1.8 billion numbers in its first day (it's at 27 billion) — because factoring smaller numbers takes less time. In spite of that small premium, it'll still take until mid-March to complete the run to 100 billion. This is, nevertheless, a significant improvement over my previous estimate of August-or-so, had I not had access to the new iMac.

Saturday, January 29

Double take

One of my pre-retirement, beer-store acquaintances emailed me yesterday: Are you using Mathematica 8 yet? The copies of Mathematica 4 and 5 on my old machine are Power-PC only, so I can't run them on my wife's Intel. Mathematica 6 is what is collecting the Ruth-Aaron numbers greater than 100 billion. This evening, I installed Mathematica 7 on her iMac and had it running the same program — but starting at 25.2 billion (which is where I was at). Each Mathematica runs at close to 100% CPU because the iMac's "i5" chip is dual-core. Running Mathematica 8 on top of that would only serve to degrade the performance. ;)

Thursday, January 27

How fast?

I've now had my wife's new iMac collecting, in Mathematica, Ruth-Aaron numbers (greater than 100 billion) for one full week. In that period of time, it has looked at some 10.5 billion candidates — roughly one million every minute. This is more than four times the rate that I have been getting on my old Mac — in a range where the numbers are five times the size. So, I could reach 100 billion (from my current position of 24.5 billion) by mid-March (at the latest) if I were to transfer the job to my wife's iMac today. Instead, I will let her computer continue its current compilation for another two weeks.

45000 Ruth-Aaron pairs

I've added another 5000 Ruth-Aaron pairs to my list: 45000 factored pairs @ 5.7 MB.

Tuesday, January 25

4000555

I have now collected more than 275000 terms of A179066 and have encountered the first few values thereof exceeding 1124577 (the assumed length of this sequence) — which, therefore, lack a symmetrical counterpoint in the x=y mirror. Term #261379 is the very first of these with a value of 4000555.

Thursday, January 20

Transition

It's been a long time coming but I finally got my wife to replace her "Bondi Blue" with a new "mid-2010" iMac. She insisted on the 21.5" screen but I maxed out every other spec. This will allow me to run Mathematica in the background and give me a handle on just how much faster this machine is compared to my old workhorse. I have started it collecting Ruth-Aaron numbers greater than 10^11.

Tuesday, January 18

My secret journal

When I started my "blahgodo" blog on 6 April 2008, I deliberately disabled its search and comment capabilities. It sat on my computer with a single link in my "chesswanks" portal, hidden in my avatar's glasses. I notified a handful of friends/acquaintances about the blog and eventually ended up with a handful of followers (who, by and large, were not the people that I had notified). The first Google mention of the word "blahgodo" did not appear until 14 August 2008.

By contrast, Glad Hobo is both searchable and commentable, but I have yet to tell a single soul about it. And although it is tempting, I will not, for I am curious just how long it will take to get that out-of-the-blue first comment.

Monday, January 17

A large Ruth-Aaron pair

A large (3109-digit) Ruth-Aaron pair from Jens Kruse Andersen in 2006:

1905497580199509189284170340160617729454507519284447380894032125401299216605231959404746406491456076087537452500615450864475294509574851694872558965147637458142847851169925504384597156839186342774703794218835079914787224201023704381910357892050633144667126701025695274917267244904046631547412451217764892682665820564602733117589573504948536227821345923192129480709161020261811449211910222044348747482379430725999743632221747352369895912918537806024983087876413015656051286630362514544336250655742275490390123627033642602816846082773618047618057241712864339666641628923304519793358493637984057056708431504856716666486852783064558823303214311031327843985329990250941661759058076234800283993403507076244375580959674432176165987722177157610949025014540262241440456225977481329826275439637821113978467713347592277272778761877369486620895052353323089466275186762154730760340318912354139790806897575696138703430821296308462103569017153398443267481786206210886767724818966146946844173334119114448743319742017799652250972014097975370944402009702825892737428030948767100797625686574016828306896510177585513807643247111028436387799038506440379364429459573363671897648038471055229959043865825572200500964733710963640202420287896153160330276281618227487298617540255476280074559350892254762502210817894443167720935846502569560506110667691599491208731919024952421363870761711803501239185538574069391983342719955717058882890896368138193494147004587331151843148675367549228967086029644320447859065120962952089596374447116233085434132873847609496753472985903393433789231694813099478361000748611580864826593011061847641838713556713471254725932258368864283683692200692371209106534770572307218871042916413593183485395365506379734347460063367332536222904080159915268050488463591666677883273886773515619453824283477025504454432123743394046125013460971239497800890371700898189508814245802605172076033911743962026904069806100895805280160105496873001315827715390516280823634234028772631187703911152987130309050986582352435463047787954537597743446653459394894961821259826264830962239473426783366979045957590674800827386475515258839150265832778437479906407032823425469516491021728001379715431899719366154512563502547285212591255838316026666455949752922243730064842786597562244544978884895727223273891803087198867866259079665555303285891331560572045662554268535966709593656589236064020830654852199844137321000857229185204802690120413827152016690707913406032458904557410391839906976201600478728787648878095918260077763451151802737157974388818693751187019698506269372899366791461115456256745354823759379125968739181772530377198570431567908347972211264615201563021468158000344014411840455964750289454658586942925561481991489057183205484114267948734129461009557038116373961305283946779805357275258052988081142338786052205011092957765100556971046168725405699196504876638805413437488325425290215402290526524610163462014214603181590049079455254375903521028239539201511806075458657292398588631028289834824303308800555229520881653153370337142941704610130531032478966096949380681751035597648875605646925429060966480912553804692548277672433233670259
and the adjacent
1905497580199509189284170340160617729454507519284447380894032125401299216605231959404746406491456076087537452500615450864475294509574851694872558965147637458142847851169925504384597156839186342774703794218835079914787224201023704381910357892050633144667126701025695274917267244904046631547412451217764892682665820564602733117589573504948536227821345923192129480709161020261811449211910222044348747482379430725999743632221747352369895912918537806024983087876413015656051286630362514544336250655742275490390123627033642602816846082773618047618057241712864339666641628923304519793358493637984057056708431504856716666486852783064558823303214311031327843985329990250941661759058076234800283993403507076244375580959674432176165987722177157610949025014540262241440456225977481329826275439637821113978467713347592277272778761877369486620895052353323089466275186762154730760340318912354139790806897575696138703430821296308462103569017153398443267481786206210886767724818966146946844173334119114448743319742017799652250972014097975370944402009702825892737428030948767100797625686574016828306896510177585513807643247111028436387799038506440379364429459573363671897648038471055229959043865825572200500964733710963640202420287896153160330276281618227487298617540255476280074559350892254762502210817894443167720935846502569560506110667691599491208731919024952421363870761711803501239185538574069391983342719955717058882890896368138193494147004587331151843148675367549228967086029644320447859065120962952089596374447116233085434132873847609496753472985903393433789231694813099478361000748611580864826593011061847641838713556713471254725932258368864283683692200692371209106534770572307218871042916413593183485395365506379734347460063367332536222904080159915268050488463591666677883273886773515619453824283477025504454432123743394046125013460971239497800890371700898189508814245802605172076033911743962026904069806100895805280160105496873001315827715390516280823634234028772631187703911152987130309050986582352435463047787954537597743446653459394894961821259826264830962239473426783366979045957590674800827386475515258839150265832778437479906407032823425469516491021728001379715431899719366154512563502547285212591255838316026666455949752922243730064842786597562244544978884895727223273891803087198867866259079665555303285891331560572045662554268535966709593656589236064020830654852199844137321000857229185204802690120413827152016690707913406032458904557410391839906976201600478728787648878095918260077763451151802737157974388818693751187019698506269372899366791461115456256745354823759379125968739181772530377198570431567908347972211264615201563021468158000344014411840455964750289454658586942925561481991489057183205484114267948734129461009557038116373961305283946779805357275258052988081142338786052205011092957765100556971046168725405699196504876638805413437488325425290215402290526524610163462014214603181590049079455254375903521028239539201511806075458657292398588631028289834824303308800555229520881653153370337142941704610130531032478966096949380681751035597648875605646925429060966480912553804692548277672433233670260

Saturday, January 15

Friday, January 14

Automobile number

An article describing Thomas A. Edison's Latest Invention appeared in Scientific American (the weekly journal of practical information) one hundred years ago today. Right beside it is a piece called The Modern Pleasure Electric Vehicle. (William Hudson writes: "The electric pleasure car is coming into its own. For certain kinds of service it is ideal. Its simple and responsive control will always be its most remarkable feature. The gasoline car has perhaps obscured the development which the electric has been undergoing in recent years.") Several other car-related articles follow in this special automobile number: Some Remarkable Mechanical Road Guides, Why You Can Buy a Good Car for Little Money, The Commercial Motor Truck vs. the Horse, A Few Shop Jobs on an Old Car, and Repainting the Old Car at Home. There follows a nice two-page spread showing twenty-five motor cars, complete with vital statistics (except for some damaged entries at the bottom of my pages). There's an article titled The Small, Inexpensive Garage. Then one about Automobile Cylinder Lubricating Oils and, on the page opposite, an advertisement for that new Edison storage battery (pictured) made by the Edison Storage Battery Company.

Thursday, January 13

The big picture

I had Mathematica create a very large (4834 x 4766 pixels) graph of A179066. Most browsers will shrink the image to fit window size — which is how it looks best — but then one can expand it to actual size to better appreciate detail. The highest-up points in this part of the graph — there are higher-ups (>200000) that are not shown — are these: {97860,121125}, {98065,121114}, {98067,121134}, {98076,121143}, {98470,121123}, {98506,121132}, {98560,121141}, {98605,121177}, {98607,121152}, {98650,121213}, {98670,121215}, {98704,121222}, {98706,121224}, {98740,121231}, {98760,121233}. Their x=y mirror images are on the far right.

I had not really intended for Glad Hobo to be only recreational mathematics (as it has been so far) but rather a less-constrained, different-emphasis successor to my "blahgodo" blog. Tomorrow I will publish something different.

Tuesday, January 11

685/252

= 11 + 22 - 33/44*55 + 66/77 - 88/99 - 110 + 121

This is a pretty good approximation to the number e.

96319/30660

= 3*7/15 + 1/292/1 + 1/1 + 2/1/3 + 1/14

The fraction in the title approximates π, and the 3, 7, 15, 1, 292, 1, 1, 1, 2, 1, 3, 1, 14 are the first thirteen terms of π's simple continued fraction (= 80143857/25510582).

Monday, January 10

Fractals everywhere

Perhaps it should not come as a surprise that sequence A179066 has fractal properties. I have put here two images: First, of values less than 7500, then of values less than 75000. The document has shrunk the pictures for easier comparison: Click on each one to look at the graphs full-size.

Saturday, January 8

Progress and prognosis

Last night I reached term #35000 in my evaluation of A179066. I had been trying a few program improvements on yet a fourth copy of Mathematica to see if I could speed things up a little. Most of these didn't really pan out — the overhead of the new code canceling any hoped-for benefits — but one obvious implementation did shave some time off the calculations, so I aborted my run and started the new routine from where I left off. Now I am about to pass term #70000 and — putting one of my two Ruth-Aaron-pairs tabulations on hold in order to devote more CPU muscle to this — tomorrow morning I should be passing term #100000.

My expectation is that — any day now — someone else will have calculated all 1124577 terms and that will be the end of my effort. In particular, Maximilian Hasler — who co-authored A179066 with Eric Angelini — thinks that "it should not take that long".

Thursday, January 6

Computing power

How difficult might it be to calculate term #1124577 of A179066? For me, very. I've been running my program now for a couple of days and have only managed some 25000 terms. I can probably see my way to 100000 terms but, short of a significantly more optimized version of my search routine (and even with that, the search beyond would prove taxing), there it will end. Remember also that I am running two other copies of Mathematica in the background (my ongoing Ruth-Aaron pairs tabulation) and these will eat up available computer-clock cycles, reducing the Mathematica kernel's CPU usage to some 60% of the available 2 MHz. Have I mentioned yet that I could really use a modern replacement for my seven-and-a-half-year-old machine?

Wednesday, January 5

Proof by contradiction

Yesterday I asked for the "final term" of a sequence. It may not be immediately obvious that there is for this sequence such a thing, but consider this: A179066 is defined in such a way that term #1124578 must have a digital root of 1 (the digital root of 1124578) and yet be composed of only the digits 0, 3, 6, and 9 (those digits that are not used in 1124578). But numbers made up of only 0s, 3s, 6s, and 9s will be evenly divisible by 3 or, alternatively, leave a remainder of 0, 3, or 6 when divided by 9 (making their digital roots 3, 6, or 9), thus contradicting our already established fact that the number that we are looking for has a digital root of 1. As defined, term #1124578 cannot exist and term #1124577, therefore, is likely the final term of this sequence.

Tuesday, January 4

28, 11, 12, 13, ...

Eric Angelini came up with this: Each successive term is the smallest positive integer not yet used, where both the term and its index have the same digital root (remainder when divided by 9) but do not have in common any of each other's base-ten digits. What might be this sequence's final term?

Monday, January 3

35000 Ruth-Aaron pairs

I've just added another 5000 Ruth-Aaron pairs to my ongoing tabulation thereof (bringing the total now to 35000 pairs). The project will likely consume the better part of this year. I realize that this would go much, much faster using a modern computer (my dual 2 GHz PowerPC G5 is now well over seven years old) but, in the absence of anyone else taking this on, at least I've got a decent start on it.

Sunday, January 2

1+2+3/45/6/7*89

Using the digits 1, 2, 3, 4, 5, 6, 7, 8, 9 in order, replace each comma with one of +, -, *, /, or a concatenation of adjacent numbers. Parentheses or other punctuation are not allowed and, as usual, multiplications and divisions are to be evaluated before additions and subtractions. My title today evaluates to a quantity fairly close to π. Can you come up with an expression that is even closer?

Saturday, January 1

Countdown

Just how likely was my pre-midnight countdown? Consider the ordering of 10, 9, 8, 7, 6, 5, 4, 3, 2, 1 and allow the commas to be replaced by one of +, -, *, /, or a concatenation of adjacent numbers (with the added proviso of no concatenation after the 10 to force the generally accepted start of the countdown). We exclude parentheses, so the mathematical statements must be evaluated in the standard order of operations (i.e., multiplications and divisions evaluated before additions and subtractions). Of the 1562500 (4*5^8) possible numerical outcomes, 486501 are unique. Of these, only 51249 are integers (31373 positive), the largest being 10*987654321. The smallest positive integers that do not appear are 1992, 2102, 2133, etc., so every year from 1993 to 2101 will have at least one solution: The largest number of solutions in this interval for a given year is 13 and the average, 616/109. Only two of the solutions start with 10/ (1997 = 10/9*87*6*5 - 43*21 and 2021 = 10/9*8*765/4 + 321). While 1992 is missing from our integer solutions, it may be approximated by the nearby 10 + 9 + 8 + 76*543/21 (as well, if we allow concatenation after the 10, it's 1098 + 765 + 4*32 + 1). You might like to look for one of the four solutions that evaluate to 377/120, the closest that we can get to the irrational number π.

An interesting bonus in last night's countdown exemplar is the symmetry of the operations: +, *, *, /, *, *, +. How common might that be? Out of the 616 solutions in the 1993-2101 interval, I count 12 that qualify:

1994 = 10*9 - 8*7 + 654*3 - 2*1
2001 = 10 - 9*87 + 65*43 - 21
2011 = 10 + 9*8*7/6*5*4 + 321
2020 = 10 - 9 + 8*7 + 654*3 + 2 - 1
2023 = 10*987/6 + 54/3*21
2024 = 10 + 9 + 8*7*6*5 + 4 + 321
2040 = 10 + 9 + 8*7 + 654*3 + 2 + 1
2072 = 10 + 9*8 + 7 + 654*3 + 21
2088 = 10*9*87/6/5/4*32*1
2093 = 10*98 + 7*6*5 + 43*21
2098 = 10 - 98*7 + 65*43 - 21 or 10 - 9*8 + 7 - 6 + 5*432 - 1