Factorial number system

In combinatorics, the factorial number system, also called factoradic, is a mixed radix numeral system adapted to numbering permutations. It is also called factorial base, although factorials do not function as base, but as place value of digits. By converting a number less than n! to factorial representation, one obtains a sequence of n digits that can be converted to a permutation of n in a straightforward way, either using them as Lehmer code or as inversion table^[1] representation; in the former case the resulting map from integers to permutations of n lists them in lexicographical order. General mixed radix systems were studied by Georg Cantor.^[2] The term "factorial number system" is used by Knuth,^[3] while the French equivalent "numération factorielle" was first used in 1888.^[4] The term "factoradic", which is a portmanteau of factorial and mixed radix, appears to be of more recent date.^[5]

Definition

The factorial number system is a mixed radix numeral system: the i-th digit from the right has base i, which means that the digit must be strictly less than i, and that (taking into account the bases of the less significant digits) its value to be multiplied by (i − 1)! (its place value).

Radix	8	7	6	5	4	3	2	1
Place value	7!	6!	5!	4!	3!	2!	1!	0!
Place value in decimal	5040	720	120	24	6	2	1	1
Highest digit allowed	7	6	5	4	3	2	1	0

From this it follows that the rightmost digit is always 0, the second can be 0 or 1, the third 0, 1 or 2, and so on (sequence A124252 in the OEIS). The factorial number system is sometimes defined with the 0! place omitted because it is always zero (sequence A007623 in the OEIS).

In this article, a factorial number representation will be flagged by a subscript "!", so for instance 341010_! stands for 3₆4₅1₄0₃1₂0₁, whose value is

• = 3×5! + 4×4! + 1×3! + 0×2! + 1×1! + 0×0! 
• = ((((3×5 + 4)×4 + 1)×3 + 0)×2 + 1)×1 + 0
• =  463₁₀.

(Note that the place value is one less than the radix position, which is why these equations begin with 5!.)

General properties of mixed radix number systems also apply to the factorial number system. For instance, one can convert a number into factorial representation producing digits from right to left, by repeatedly dividing the number by the place values (1, 2, 3, ...), taking the remainder as digits, and continuing with the integer quotient, until this quotient becomes 0.

For example, 463₁₀ can be transformed into a factorial representation by these successive divisions:

• 463 ÷ 1 = 463, remainder 0
• 463 ÷ 2 = 231, remainder 1
• 231 ÷ 3 = 77, remainder 0
• 77 ÷ 4 = 19, remainder 1
• 19 ÷ 5 = 3, remainder 4
• 3 ÷ 6 = 0, remainder 3

The process terminates when the quotient reaches zero. Reading the remainders backward gives 341010_!.

In principle, this system may be extended to represent fractional numbers, though rather than the natural extension of place values (−1)!, (−2)!, etc., which are undefined, the symmetric choice of radix values n = 0, 1, 2, 3, 4, etc. after the point may be used instead. Again, the 0 and 1 places may be omitted as these are always zero. The corresponding place values are therefore 1/1, 1/1, 1/2, 1/6, 1/24, ..., 1/n!, etc.

Examples

Here are the first twenty-four numbers, counting from zero.

The table on the left shows permutations, and inversion vectors^[6] (which are reflected factorial numbers) below them. Another column shows the inversion sets. The digit sums of the inversion vectors (or factorial numbers) and the cardinalities of the inversion sets are equal (and have the same parity as the permutation). The digit sums form the sequence  A034968.

Permutohedron graph
showing permutations and their inversion vectors
(compare version with factorial numbers)
The arrows indicate the bitwise less or equal $\leq$ relation.

decimal factorial 0 0 1 10 2 100 3 110 4 200 5 210 6 1000 7 1010 8 1100 9 1110 10 1200 11 1210 12 2000 13 2010 14 2100 15 2110 16 2200 17 2210 18 3000 19 3010 20 3100 21 3110 22 3200 23 3210

For another example, the greatest number that could be represented with six digits would be 543210_! which equals 719 in decimal:

5×5! + 4×4! + 3x3! + 2×2! + 1×1! + 0×0!.

Clearly the next factorial number representation after 543210_! is 1000000_! which designates 6! = 720₁₀, the place value for the radix-7 digit. So the former number, and its summed out expression above, is equal to:

6! − 1.

The factorial number system provides a unique representation for each natural number, with the given restriction on the "digits" used. No number can be represented in more than one way because the sum of consecutive factorials multiplied by their index is always the next factorial minus one:

$\sum _{i=0}^{n}{i\cdot i!}={(n+1)!}-1.$

This can be easily proved with mathematical induction, or simply by noticing that ${\displaystyle \forall i,i.i!=(i+1-1).i!=(i+1)!-i!}$ : subsequent terms cancel each other, leaving the first and last term (see Telescoping series)

However, when using Arabic numerals to write the digits (and not including the subscripts as in the above examples), their simple concatenation becomes ambiguous for numbers having a "digit" greater than 9. The smallest such example is the number 10 × 10! = 36288000₁₀, which may be written A0000000000_!, but not 100000000000_! which denotes 11! = 39916800₁₀. Thus using letters A–Z to denote digits 10, 11, 12, ..., 35 as in other base-N make the largest representable number 36! − 1 = 371993326789901217467999448150835199999999₁₀. For arbitrarily greater numbers one has to choose a base for representing individual digits, say decimal, and provide a separating mark between them (for instance by subscripting each digit by its base, also given in decimal, like 2₄0₃1₂0₁). In fact the factorial number system itself is not truly a numeral system in the sense of providing a representation for all natural numbers using only a finite alphabet of symbols.

Permutations

There is a natural mapping between the integers 0, ..., n! − 1 (or equivalently the numbers with n digits in factorial representation) and permutations of n elements in lexicographical order, when the integers are expressed in factoradic form. This mapping has been termed the Lehmer code (or inversion table). For example, with n = 3, such a mapping is

decimal	factorial	permutation
010	000!	(0,1,2)
110	010!	(0,2,1)
210	100!	(1,0,2)
310	110!	(1,2,0)
410	200!	(2,0,1)
510	210!	(2,1,0)

The leftmost factoradic digit 0, 1, or 2 is chosen as the first permutation digit from the ordered list (0,1,2) and is removed from the list. Think of this new list as zero indexed and each successive digit dictates which of the remaining elements is to be chosen. If the second factoradic digit is "0" then the first element of the list is selected for the second permutation digit and is then removed from the list. Similarly if the second factoradic digit is "1", the second is selected and then removed. The final factoradic digit is always "0", and since the list now contains only one element it is selected as the last permutation digit.

The process may become clearer with a longer example. For example, here is how the digits in the factoradic 4041000_! (equal to 2982₁₀) pick out the digits in (4,0,6,2,1,3,5), the 2982nd permutation of the numbers 0 through 6.

                                 4041000! → (4,0,6,2,1,3,5)
factoradic:  4          0                        4        1          0          0        0!
             |          |                        |        |          |          |        |
    (0,1,2,3,4,5,6) -> (0,1,2,3,5,6) -> (1,2,3,5,6) -> (1,2,3,5) -> (1,3,5) -> (3,5) -> (5)
             |          |                        |        |          |          |        |
permutation:(4,         0,                       6,       2,         1,         3,       5)

A natural index for the group direct product of two permutation groups is the concatenation of two factoradic numbers, with two subscript "!"s.

           concatenated
 decimal   factoradics            permutation pair
    010     000!000!           ((0,1,2),(0,1,2))
    110     000!010!           ((0,1,2),(0,2,1))
               ...
    510     000!210!           ((0,1,2),(2,1,0))
    610     010!000!           ((0,2,1),(0,1,2))
    710     010!010!           ((0,2,1),(0,2,1))
               ...
   2210     110!200!           ((1,2,0),(2,0,1))
               ...
   3410     210!200!           ((2,1,0),(2,0,1))
   3510     210!210!           ((2,1,0),(2,1,0))

Fractional values

Unlike single radix systems whose place values are baseⁿ for both positive and negative integral n, the factorial number base cannot be extended to negative place values as these would be (−1)!, (−2)! and so on, and these values are undefined. (see factorial)

One possible extension is therefore to use 1/0!, 1/1!, 1/2!, 1/3!, ..., 1/n! etc. instead, possibly omitting the 1/0! and 1/1! places which are always zero.

With this method, all rational numbers have a terminating expansion, whose length in 'digits' is less than or equal to the denominator of the rational number represented. This may be proven by considering that there exists a factorial for any integer and therefore the denominator divides into its own factorial even if it does not divide into any smaller factorial.

By necessity, therefore, the factoradic expansion of the reciprocal of a prime has a length of exactly that prime (less one if the 1/1! place is omitted). It can also be proven that the last 'digit' or term of the representation of a rational with prime denominator is equal to the difference between the numerator and the prime denominator.

There is also a non-terminating equivalent for every rational number akin to the fact that in decimal 0.24999... = 0.25 = 1/4 and 0.999... = 1, etc., which can be created by reducing the final term by 1 and then filling in the remaining infinite number of terms with the highest value possible for the radix of that position.

In the following selection of examples, spaces are used to separate the place values, otherwise represented in decimal. The rational numbers on the left are also in decimal:

• $1/2=0.0\ 1_{!}$
• $1/3=0.0\ 0\ 2_{!}$
• $1/4=0.0\ 0\ 1\ 2_{!}$
• $1/5=0.0\ 0\ 1\ 0\ 4_{!}$
• $1/6=0.0\ 0\ 1_{!}$
• $1/11=0.0\ 0\ 0\ 2\ 0\ 5\ 3\ 1\ 4\ 0\ A_{!}$
• $2/11=0.0\ 0\ 1\ 0\ 1\ 4\ 6\ 2\ 8\ 1\ 9_{!}$
• $9/11=0.0\ 1\ 1\ 3\ 3\ 1\ 0\ 5\ 0\ 8\ 2_{!}$
• $1/15=0.0\ 0\ 0\ 1\ 3_{!}$
• $237/360=0.0\ 1\ 0\ 3\ 4_{!}$

There are also a small number of constants that have patterned representations with this method:

• $e=1\ 0.0\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1\ 1..._{!}$
• $e^{-1}=0.0\ 0\ 2\ 0\ 4\ 0\ 6\ 0\ 8\ 0\ A\ 0\ C\ 0\ E..._{!}$
• $\sin(1)=0.0\ 1\ 2\ 0\ 0\ 5\ 6\ 0\ 0\ 9\ A\ 0\ 0\ D\ E..._{!}$
• $\cos(1)=0.0\ 1\ 0\ 0\ 4\ 5\ 0\ 0\ 8\ 9\ 0\ 0\ C\ D\ 0..._{!}$
• $\sinh(1)=1.0\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0..._{!}$
• $\cosh(1)=1.0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1\ 0\ 1..._{!}$

See also

• Primorial number system
• Combinatorial number system (also called combinadics)
• Steinhaus–Johnson–Trotter algorithm, an algorithm that generates Gray codes for the factorial number system

References

1. Knuth, D. E. (1973), "Volume 3: Sorting and Searching", The Art of Computer Programming, Addison-Wesley, p. 12, ISBN 0-201-89685-0
2. Cantor, G. (1869), Zeitschrift für Mathematik und Physik, 14.
3. Knuth, D. E. (1997), "Volume 2: Seminumerical Algorithms", The Art of Computer Programming (3rd ed.), Addison-Wesley, p. 192, ISBN 0-201-89684-2.
4. Laisant, Charles-Ange (1888), "Sur la numération factorielle, application aux permutations", Bulletin de la Société Mathématique de France (in French), 16: 176–183.
5. The term "factoradic" is apparently introduced in McCaffrey, James (2003), Using Permutations in .NET for Improved Systems Security, Microsoft Developer Network.
6. Weisstein, Eric W. "Inversion Vector". MathWorld.

• Mantaci, Roberto; Rakotondrajao, Fanja (2001), "A permutation representation that knows what "Eulerian" means" (PDF), Discrete Mathematics and Theoretical Computer Science, 4: 101–108.
• Arndt, Jörg (2010). Matters Computational: Ideas, Algorithms, Source Code. pp. 232–238.

External links

• A Lehmer code calculator Note that their permutation digits start from 1, so mentally reduce all permutation digits by one to get results equivalent to those on this page.
• Factorial number system