Regular number

Formally, a regular number is an integer of the form ${\displaystyle 2^{i}\cdot 3^{j}\cdot 5^{k}}$, for nonnegative integers ${\displaystyle i}$, ${\displaystyle j}$, and ${\displaystyle k}$. Such a number is a divisor of ${\displaystyle 60^{\max(\lceil i\,/2\rceil ,j,k)}}$. The regular numbers are also called 5-smooth, indicating that their greatest prime factor is at most 5.[2] More generally, a k-smooth number is a number whose greatest prime factor is at most k.[3]

The first few regular numbers are[2]

Several other sequences at the On-Line Encyclopedia of Integer Sequences have definitions involving 5-smooth numbers.[4]

Although the regular numbers appear dense within the range from 1 to 60, they are quite sparse among the larger integers. A regular number ${\displaystyle n=2^{i}\cdot 3^{j}\cdot 5^{k}}$ is less than or equal to some threshold ${\displaystyle N}$ if and only if the point ${\displaystyle (i,j,k)}$ belongs to the tetrahedron bounded by the coordinate planes and the plane

$${\displaystyle i\ln 2+j\ln 3+k\ln 5\leq \ln N,}$$

as can be seen by taking logarithms of both sides of the inequality ${\displaystyle 2^{i}\cdot 3^{j}\cdot 5^{k}\leq N}$. Therefore, the number of regular numbers that are at most ${\displaystyle N}$ can be estimated as the volume of this tetrahedron, which is

$${\displaystyle {\frac {\log _{2}N\,\log _{3}N\,\log _{5}N}{6}}.}$$

Even more precisely, using big O notation, the number of regular numbers up to ${\displaystyle N}$ is

$${\displaystyle {\frac {\left(\ln(N{\sqrt {30}})\right)^{3}}{6\ln 2\ln 3\ln 5}}+O(\log N),}$$

and it has been conjectured that the error term of this approximation is actually ${\displaystyle O(\log \log N)}$.[2] A similar formula for the number of 3-smooth numbers up to ${\displaystyle N}$ is given by Srinivasa Ramanujan in his first letter to G. H. Hardy.[5]

AO 6456, a table of reciprocals of regular numbers from Seleucid Uruk, copied from an unknown earlier source

In the Babylonian sexagesimal notation, the reciprocal of a regular number has a finite representation. If ${\displaystyle n}$ divides ${\displaystyle 60^{k}}$, then the sexagesimal representation of ${\displaystyle 1/n}$ is just that for ${\displaystyle 60^{k}/n}$, shifted by some number of places. This allows for easy division by these numbers: to divide by ${\displaystyle n}$, multiply by ${\displaystyle 1/n}$, then shift.[6]

For instance, consider division by the regular number 54 = 2¹3³. 54 is a divisor of 60³, and 60³/54 = 4000, so dividing by 54 in sexagesimal can be accomplished by multiplying by 4000 and shifting three places. In sexagesimal 4000 = 1×3600 + 6×60 + 40×1, or (as listed by Joyce) 1:6:40. Thus, 1/54, in sexagesimal, is 1/60 + 6/60² + 40/60³, also denoted 1:6:40 as Babylonian notational conventions did not specify the power of the starting digit. Conversely 1/4000 = 54/60³, so division by 1:6:40 = 4000 can be accomplished by instead multiplying by 54 and shifting three sexagesimal places.

The Babylonians used tables of reciprocals of regular numbers, some of which still survive.[7] These tables existed relatively unchanged throughout Babylonian times.[6] One tablet from Seleucid times, by someone named Inaqibıt-Anu, contains the reciprocals of 136 of the 231 six-place regular numbers whose first place is 1 or 2, listed in order. It also includes reciprocals of some numbers of more than six places, such as 3²³ (2 1 4 8 3 0 7 in sexagesimal), whose reciprocal has 17 sexagesimal digits. Noting the difficulty of both calculating these numbers and sorting them, Donald Knuth in 1972 hailed Inaqibıt-Anu as "the first man in history to solve a computational problem that takes longer than one second of time on a modern electronic computer!" (Two tables are also known giving approximations of reciprocals of non-regular numbers, one of which gives reciprocals for all the numbers from 56 to 80.)[8][9]

Although the primary reason for preferring regular numbers to other numbers involves the finiteness of their reciprocals, some Babylonian calculations other than reciprocals also involved regular numbers. For instance, tables of regular squares have been found[6] and the broken tablet Plimpton 322 has been interpreted by Neugebauer as listing Pythagorean triples ${\displaystyle (p^{2}-q^{2},\,2pq,\,p^{2}+q^{2})}$ generated by ${\displaystyle p}$ and ${\displaystyle q}$ both regular and less than 60.[10] Fowler and Robson discuss the calculation of square roots, such as how the Babylonians found an approximation to the square root of 2, perhaps using regular number approximations of fractions such as 17/12.[9]

In music theory, the just intonation of the diatonic scale involves regular numbers: the pitches in a single octave of this scale have frequencies proportional to the numbers in the sequence 24, 27, 30, 32, 36, 40, 45, 48 of nearly consecutive regular numbers.[11] Thus, for an instrument with this tuning, all pitches are regular-number harmonics of a single fundamental frequency. This scale is called a 5-limit tuning, meaning that the interval between any two pitches can be described as a product 2ⁱ3^j5^k of powers of the prime numbers up to 5, or equivalently as a ratio of regular numbers.[12]

5-limit musical scales other than the familiar diatonic scale of Western music have also been used, both in traditional musics of other cultures and in modern experimental music: Honingh & Bod (2005) list 31 different 5-limit scales, drawn from a larger database of musical scales. Each of these 31 scales shares with diatonic just intonation the property that all intervals are ratios of regular numbers.[12] Euler's tonnetz provides a convenient graphical representation of the pitches in any 5-limit tuning, by factoring out the octave relationships (powers of two) so that the remaining values form a planar grid.[12] Some music theorists have stated more generally that regular numbers are fundamental to tonal music itself, and that pitch ratios based on primes larger than 5 cannot be consonant.[13] However the equal temperament of modern pianos is not a 5-limit tuning,[14] and some modern composers have experimented with tunings based on primes larger than five.[15]

In connection with the application of regular numbers to music theory, it is of interest to find pairs of regular numbers that differ by one. There are exactly ten such pairs ${\displaystyle (x,x+1)}$ and each such pair defines a superparticular ratio ${\displaystyle {\tfrac {x+1}{x}}}$ that is meaningful as a musical interval. These intervals are 2/1 (the octave), 3/2 (the perfect fifth), 4/3 (the perfect fourth), 5/4 (the just major third), 6/5 (the just minor third), 9/8 (the just major tone), 10/9 (the just minor tone), 16/15 (the just diatonic semitone), 25/24 (the just chromatic semitone), and 81/80 (the syntonic comma).[16]

In the Renaissance theory of universal harmony, musical ratios were used in other applications, including the architecture of buildings. In connection with the analysis of these shared musical and architectural ratios, for instance in the architecture of Palladio, the regular numbers have also been called the harmonic whole numbers.[17]

Algorithms for calculating the regular numbers in ascending order were popularized by Edsger Dijkstra. Dijkstra (1976, 1981) attributes to Hamming the problem of building the infinite ascending sequence of all 5-smooth numbers; this problem is now known as Hamming's problem, and the numbers so generated are also called the Hamming numbers. Dijkstra's ideas to compute these numbers are the following:

• The sequence of Hamming numbers begins with the number 1.
• The remaining values in the sequence are of the form ${\displaystyle 2h}$, ${\displaystyle 3h}$, and ${\displaystyle 5h}$, where ${\displaystyle h}$ is any Hamming number.
• Therefore, the sequence ${\displaystyle H}$ may be generated by outputting the value 1, and then merging the sequences ${\displaystyle 2H}$, ${\displaystyle 3H}$, and ${\displaystyle 5H}$.

This algorithm is often used to demonstrate the power of a lazy functional programming language, because (implicitly) concurrent efficient implementations, using a constant number of arithmetic operations per generated value, are easily constructed as described above. Similarly efficient strict functional or imperative sequential implementations are also possible whereas explicitly concurrent generative solutions might be non-trivial.[18]

In the Python programming language, lazy functional code for generating regular numbers is used as one of the built-in tests for correctness of the language's implementation.[19]

A related problem, discussed by Knuth (1972), is to list all ${\displaystyle k}$-digit sexagesimal numbers in ascending order (see #Babylonian mathematics above). In algorithmic terms, this is equivalent to generating (in order) the subsequence of the infinite sequence of regular numbers, ranging from ${\displaystyle 60^{k}}$ to ${\displaystyle 60^{k+1}}$.[8] See Gingerich (1965) for an early description of computer code that generates these numbers out of order and then sorts them;[20] Knuth describes an ad hoc algorithm, which he attributes to Bruins (1970), for generating the six-digit numbers more quickly but that does not generalize in a straightforward way to larger values of ${\displaystyle k}$.[8] Eppstein (2007) describes an algorithm for computing tables of this type in linear time for arbitrary values of ${\displaystyle k}$.[21]

Heninger, Rains & Sloane (2006) show that, when ${\displaystyle n}$ is a regular number and is divisible by 8, the generating function of an ${\displaystyle n}$-dimensional extremal even unimodular lattice is an ${\displaystyle n}$th power of a polynomial.[22]

As with other classes of smooth numbers, regular numbers are important as problem sizes in computer programs for performing the fast Fourier transform, a technique for analyzing the dominant frequencies of signals in time-varying data. For instance, the method of Temperton (1992) requires that the transform length be a regular number.[23]

Book VIII of Plato's Republic involves an allegory of marriage centered on the highly regular number 60⁴ = 12,960,000 and its divisors (see Plato's number). Later scholars have invoked both Babylonian mathematics and music theory in an attempt to explain this passage.[24]

Certain species of bamboo release large numbers of seeds in synchrony (a process called masting) at intervals that have been estimated as regular numbers of years, with different intervals for different species, including examples with intervals of 10, 15, 16, 30, 32, 48, 60, and 120 years.[25] It has been hypothesized that the biological mechanism for timing and synchronizing this process lends itself to smooth numbers, and in particular in this case to 5-smooth numbers. Although the estimated masting intervals for some other species of bamboo are not regular numbers of years, this may be explainable as measurement error.[25]

• Table of reciprocals of regular numbers up to 3600 from the web site of Professor David E. Joyce, Clark University.
• RosettaCode Generation of Hamming_numbers in ~ 50 programming languages