Logo

The Data Daily

Difficult Probability Problem: Distribution of Digits in Rogue Systems

Difficult Probability Problem: Distribution of Digits in Rogue Systems

I recently posted a table summarizing probabilistic properties of digits in various number representation systems, see here.  The topic is already rather difficult for well-behaved systems (those listed in my table) but some systems are rogue, and do not have these nice statistical properties. Here we focus on one of these less known systems, though at this point it is not sure if it is well-behaved or not: It is part of this problem to figure this out.

The problem is about the representation of a real number x greater or equal to 1, in some base b, using the following infinite product rather than the familiar (well-behaved) representation such as in the decimal or binary system: 

The a(k)'s are called the (binary) digits of x in that system. The case b = 2 has been used to compute logarithms using the Feynman's algorithm, see here. Not all x's can not be represented that way. A necessary condition for x to have such a representation in base b is

The last inequality becomes an equality when all the a(k)'s are equal to 1. Also not all b's work. If b > 2, not all x's can be represented this way. The algorithm to compute, for a given x, the digits a(k)'s, is as follows:

2. Compute a(k) for k > 0 with the formula3. Compute w(k) using  

Here are some difficult problems:

Note that the algorithm to generate the digits becomes very unstable after 30 iterations or so. You may need to use high precision computing or some other technique to analyze long-term behavior, see here.

For illustration purposes, here are the first few digits a(k) for x = Pi / 2 and b = 2, yielding 8 correct decimals when expressed in traditional base 10: 

For the well-behaved number representation systems, there is always an algorithm of the form

to compute the digits, see here.  The function g is associated with the equilibrium distribution (see same reference.)

Is it also the case here? I suspect that the answer is no. Are there bivariate functions g and h that would work, for instance something like x(k) = g(x(k-1), x(k-2)) ? This could lead to a stochastic integral equation, perhaps easily solvable, as for the well-behaved systems.

Finally, there are other number representation systems that are not well-behaved, for instance Egyptian fractions.

For related articles from the same author, click here or visit www.VincentGranville.com. Follow me on Twitter at @GranvilleDSC or on LinkedIn.

Images Powered by Shutterstock