Schizophrenic thinking

Some years ago I wrote up something about schizophrenic numbers. I read about these numbers in one of Clifford Pickover’s books, but he didn’t really get into what made them work the way they do, so I figured it out on my own. I posted it to the web. Now I can’t find a copy in my files or on my computer, and even the Internet Archive fails me — it didn’t archive it.

So I reconstructed it, more or less. I think.

The schizophrenic numbers are defined as follows. To start with we define

f(x) = \begin{cases} 0& \mbox{if } m=0\\10f(m-1)+m& \mbox{if } m>0\end{cases}\\=\sum\limits_{k=0}^{m} 10^{m-k}k

The first several f(m) are

f(0)=0 \\ f(1)=1\\f(2)=12\\f(3)=123\\f(4)=1234\\f(5)=12345\\f(6)=123456\\f(7)=1234567\\f(8)=12345678\\f(9)=123456789\\f(10)=1234567900\\f(11)=12345679011

and so on. Now the mth schizophrenic number is

S(m)=\sqrt{f(2m+1)}

Are any of the f(2m+1) perfect squares (for m>0)? No, they aren’t, as will become evident below. Therefore S(0) is the only rational in the set. The rest are irrational, so should be non-repeating decimals.

But let’s look at S(20)=\sqrt{12345679012345679012345679012345679012341}. Its value is

S(20)=111111111111111111111.111111111111111111...

Wat.

It can’t be a rational number, unless it’s an integer, and it’s not an integer. What’s going on? Let’s get some more significant figures.

S(20)=111111111111111111111.11111111111111111109005...

Oh, okay, that’s better. It’s not a repeating decimal after all…

S(20)=111111111111111111111.1111111111111111110900555555\\5555555555555555555555555555555...

Wat.

S(20)=111111111111111111111.1111111111111111110900555555\\555555555555555555555555555555535605416...

Um, okay?

S(20)=111111111111111111111.1111111111111111110900555555\\555555555555555555555555555555535605416666666666666666666...

Wat.

Seems not to be able to make up its mind whether to be rational or not. Hence schizophrenic. Here’s about 500 decimal places, broken up into “rational” and “irrational” pieces:

S(20)=111111111111111111111.111111111111111111 \\ 0900\\5555555555555555555555555555555555555\\3560541\\666666666666666666666666666666666\\2886115347\\22222222222222222222222222222\\1326704128428819\\444444444444444444444444\\2068634941610546874\\999999999999999999999\\32467614881946461588541\\6666666666666666\\465564871268691522639973958\\333333333333\\2714065492129693624196136474609374\\999999\\804414573486517125197501992119683159722159\\214377070903714104597488983578152126715487\\497795899993976792215472407658894849931527\\56767644661366482...

(I calculated that by hand. HA HA HA no actually I used a Python script, with the mpmath library.) The repeating pieces get shorter and shorter, the non-repeating pieces get longer and longer, and then (apparently) the repetition stops (aside from random repeats) and it finally starts behaving like we expect irrational numbers to do.

I mean, sure, we all know there exist an infinite number of irrationals that have arbitrarily large (but finite) numbers of repeating digits in them. But who’d expect behavior like this from the square root of an integer? What’s going on?

Well, observe that if x = 0.0123456790123456790123... then

10^9x = 12345679.0123456790123456790123...=12345679+x \\ \Rightarrow x=12345679/999999999=1/81

So f(m) is a little smaller than 10^{m+1}/81. In particular you can see if you think about it, and you can prove inductively using the recurrence formula for f(m), that

f(m)=10^{m+1}/81-(10+9m)/81

which means

S(m)=\sqrt{10^{2m+2}/81-(19+18m)/81} \\ =\frac{10^{m+1}}{9}\sqrt{1-\frac{19+18m}{10^{2m+2}}}

Now you can use the Taylor series expansion

\sqrt{1-y} = 1-y/2-y^2/4-y^3/8-5y^4/128....

Here y is (19+18m)/10^{2m+2}, which for large m is a bunch of zeros after the decimal point, followed by a few digits. In the expansion the terms after the first are powers of y multiplied by an integer and divided by powers of 2, which for higher powers are lots of zeros after the decimal point, followed by some digits. The infinite sum gets divided by 9 — which turns the first term into 0.11111... and each of the following terms into zeros after the decimal point, followed by some digits, followed by some digit repeating to infinity — and then multiplied by 10^{m+1}, which just shifts the decimal point. Look at the first few terms for S(5):

111111.111111... \\ -0.0000060555555...\\-0.000000000000000165013888888...\\-0.0000000000000000000000000089932569444444...\\-0.0000000000000000000000000056207855902777777...

For the first three terms, there are enough zeros in the subsequent term to preserve some of the repeating digits in the sum. The fifth term, though, scribbles over the repeating 4s in the fourth term, and so we only see three sets of repeating digits in S(5):

S(5)=111111.111105055555555390541666657673409...

But for larger m, each term has more leading zeros due to the power of 10^{2m+2} in the denominator, and more terms contribute repeating digits before the scribble-over point is reached. (And you get more and more 1s before and after the decimal point; in other words, you get closer and closer to a power of 10 times 1/9, so no, none of the corresponding f(m) are perfect squares.) I referred somewhat facetiously to “rational” and “irrational” pieces above, but it’s sort of true; the repeating digits are where you can see that there is a finite number of rational terms in the Taylor series contributing to the sum.

Advertisements

6 thoughts on “Schizophrenic thinking

  1. In the definition you have the function equal to a summation but in the summation the integers are the reverse of the function as defined.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s