Over at JD2718, Jonathan asks the following question:

[A conference] presenter posed a problem that required finding three numbers that multiplied to make 72. The list included 1, 8, 9 and 2, 2, 18 and 3, 4, 6 and several other triples (groups of three numbers).

So I wondered, could we look at 72 and see how many factor triples it has? In general, how many factor triples does a number have?

Answering this question involves a lot of nice mathematics!

There are two different things this question could mean: it could be asking for the number of *ordered triples* (*a*, *b*, *c*) of positive integers such that *a* × *b* × *c* = *n*; or it could be asking for the number of *multisets* {*a*, *b*, *c*} (so that order doesn’t matter). Equivalently, we could impose a particular order, like asking that . For example, the number 6 has two unordered factor triples (1×1×6 and 1×2×3) but nine ordered factor triples (1×1×6, 1×6×1, 6×1×1, 1×2×3, 1×3×2, 2×1×3, 2×3×1, 3×1×2, and 3×2×1). As is often the case in combinatorics, the unordered factorization problem is much trickier, so let’s start with the ordered case.

### The ordered case

Following Jonathan, first imagine that *n* is just a power of a single prime, say *n* = *p*^{k} for some prime *p* and nonnegative integer *k*. Then all of the factors *a*, *b*, and *c* must also be powers of *p*. In fact, if we set *a* = *p*^{x}, *b* = *p*^{y}, and *c* = *p*^{z}, then our original equation *a* × *b* × *c* = *n* gets transformed into the equation *x* + *y* + *z* = *k*, to be solved for nonnegative integers *x*, *y*, *z*. This is the classic setting for *stars and bars* (the good kind, not the treason-in-defense-of-slavery kind); as the Wikipedia article explains, the number of solutions is exactly the binomial coefficient . For example, when *n* = *p* is a prime number (so *k* = 1), there are solutions (1×1×*p*, 1×*p*×1, and *p*×1×1), while if *n* = *p*^{2} is the square of a prime then there are solutions (1×1×*p*^{2}, 1×*p*^{2}×1, and *p*^{2}×1×1, 1×*p*×*p*, *p*×*p*×1, and *p*×1×*p*).

Now what if *n* is divisible by more than one prime? Well, then the situations for the different primes are completely independent: for each prime, we have to choose how many of its factors go in the first coordinate, how many go in the second coordinate, and hand many go in the third coordinate. So, for example, since 6 = 2×3, we have the following table of possibilities:

(1, 1, 2) | (1, 2, 1) | (2, 1, 1) | |
---|---|---|---|

(1, 1, 3) | 1×1×6 | 1×2×3 | 2×1×3 |

(1, 3, 1) | 1×3×2 | 1×6×1 | 2×3×1 |

(3, 1, 1) | 3×1×2 | 3×2×1 | 6×1×1 |

So that recovers exactly the nine solutions we had before. In general, if is the prime factorization of *n*, then the number of ordered factor triples is .

### The unordered case

To a first approximation, one should expect there to be six times as many ordered triples as unordered triples, since a typical unordered triple like {1, 2, 3} can be ordered in six different ways. However, there are always triples like {1, 1, 6} (with two elements equal) that only have three different orderings, and even triples like {2, 2, 2} (with all three elements the same) that only have one ordering. So the trick is going to be to figure out how many of each kind there are.

#### A single prime

Let’s start with the case of a single prime again. As before, we’re studying the nonnegative integer solutions of the equation the equation *x* + *y* + *z* = *k*. The unordered solutions of this equation, or equivalently the solutions such that , are the integer partitions of *k* having at most three parts; however, rather than try to count these directly (which would lead to sequence A001399 in the OEIS), it will be more profitable to follow the ideas behind Burnside’s lemma and count ordered solutions based on how much duplication they have.

We know there are ordered solutions in total. If *k* happens to be divisible by 3, then there is one solution with *x* = *y* = *z*, and otherwise there are none. Now let’s think about solutions with *y* = *z*: in this case, the equation becomes *x* + 2*y* = *k*. Counting solutions of this two-variable equation is easy: *y* can be any nonnegative integer up to (but not larger than) *k*/2, and for each such *y* there is exactly one value of *x* that works. So the number of solutions in this case is , where that bracket notation represents the floor function (rounding down to the nearest integer)<!– and we add one because the smallest allowable *y* is 0–>. But we have to be a bit careful here: if *k* is divisible by 3, one of the solutions we’ve just counted is the solution with *x* = *y* = *z*.

So far we have the following cases: if *k* is a multiple of 3, there is 1 solution (*x*, *x*, *x*) and solutions (*x*, *y*, *y*) with , while if *k* is not a multiple of 3 then there are solutions (*x*, *y*, *y*) with . And of course, for each ordered solution (*x*, *y*, *y*), there is also a solution (*y*, *x*, *y*) and a solution (*y*, *y*, *x*). So if *k* is divisible by 3, we’ve accounted for ordered solutions, while if *k* is not divisible by 3 then we’ve accounted for solutions. Finally, all the other ordered solutions have all entries different, so they come in groups of six rearrangements, leading to unordered solutions if *k* is divisible by 3 or to unordered solutions if *k* is not divisible by 3. (It isn’t obvious that these numbers are integers, is it? But they have to be, as a consequence of our proof. One could prove it directly by carefully considering all possibilities for the remainder upon dividing *k* by 6.) If all we wanted was the total number of unordered solutions, we could add the three numbers from the previous paragraph up: for example, when *k* is divisible by 3, we get unordered solutions. However, as we’ll see next, for purposes of combining different primes together, it’s really important that we know how many of each symmetry type there are.

#### Products of many primes

Having dealt with powers of a single prime, let’s move to the general case. Suppose that . First, let’s assume that *n* is not a perfect cube (so it cannot be written as *a* × *a* × *a* = *n* for integer *a*). As we saw above, there are a total of *ordered* factor triples for *n*. Almost all of these come in groups of six; the exceptions are the triples (*a*, *a*, *b*), with only three rearrangements. To get such a triple, you need that simultaneously every single prime must divide the first two factors the same number of times. As we saw when studying single primes, there are exactly ways to distribute the powers of *p*_{1} so that the first two factors get the same number of powers. Similarly, there are ways to distribute the powers of *p*_{2}, and so on. Taking all the different primes into consideration, this gives factor triples of the form (*a*, *a*, *b*). Since *n* is not a perfect cube, in each of these factorizations we have , and so there exactly three permutations (*a*, *a*, *b*), (*a*, *b*, *a*), (*b*, *a*, *a*) of each such triple. That leaves other ordered factor triples. We argued at the beginning of this paragraph that all of these triples come in groups of six rearrangements, so the total number of *unordered* factor triples is what we get by dividing this number by six and adding it to what we already had, namely

which we can simplify somewhat to

Finally, in the case that *n* *is* a perfect cube, we have to adjust the count slightly: among those ordered triples whose first two entries are equal, there is now one in which all three entries are equal. Making the necessary adjustments, we get in this case that there are

factor triples.

For example, consider the factor triples of 1000 = 2^{3} × 5^{3}. There are 10 ways to distribute the powers of 2 (three rearrangements of 1 × 1 × 8, six rearrangements of 1 × 2 × 4, and 2 × 2 × 2) and similarly 10 ways to distribute the powers of 5, so a total of 100 *ordered* factor triples. Of these, there are four in which the first two entries are equal: 1 × 1 × 1000, 2 × 2 × 250, 5 × 5 × 40, and 10 × 10 × 10. The first three of these have three permutations each, so this accounts for 1 + 3 × 3 = 10 of the 100 ordered factor triples. The remaining 90 triples come in 15 groups of six; their unordered representations are 1 × 2 × 500, 1 × 4 × 250, 1 × 5 × 200, 1 × 8 × 125, 1 × 10 × 100, 1 × 20 × 50, 1 × 25 × 40, 2 × 4 × 125, 2 × 5 × 100, 2 × 10 × 50, 2 × 25 × 40, 4 × 5 × 50, 4 × 10 × 25, 5 × 8 × 25, and 5 × 10 × 20. Thus, the total number of unordered factorization triples of 1000 is

as expected. For the original example 72 = 2^{3} × 3^{2} (so not a cube), we get

unordered factor triples.

#### TL;DR: OEIS

This sequence is number A034836 in the OEIS. My impression is that the conjectured formula there is equivalent to what I’ve worked out above, but I haven’t checked that carefully.