Probability distribution of product of integers I have a scoring system based on 5 factors with integer values from 1 to 5:
Score = A * B * C * D * E

So the Score can range from 1 to 3125. Each of the factors is uniformly distributed between 1 and 5. I want to find a way to partition the results into 5 ranges so that the probability of a score in each range is equal.
Because the final values for Score are only integers with prime factors 2, 3, and 5, the distribution of possible scores is weighted heavily toward low numbers. I did a simple graph in Excel and it looks like a Poisson or Gamma distribution, but that's just an uneducated guess.
Simply put, I want a way to say, "This score is between X and Y, so it is in the 2nd quintile of all scores."

The problem in this case is that since the numbers you're multiplying can be the same, there are many possible ways how to express a single product value. For example, there are 80 ways how to express 40 as the product of five numbers from 1 to 5. The chance that the boundaries of quintiles will fall between different product values is very small, and it is indeed this case. Just to test the idea, I made a small Haskell program to produce all possible values:

import Control.Applicative
import Data.List (sort)

scores = sort $ prod <$> xs <*> xs <*> xs <*> xs <*> xs where
    prod a b c d e = a * b * c * d * e
    xs = [1..5]

And confirmed the boundaries of the 1st and 2nd quintile contain the same value:
> scores !! 624
40
> scores !! 625
40

So if you put the boundary between 1st and 2nd quintile to be $\geq 40$, the 1st quintile will be too small. And if you put it to be $> 40$, it will be too large.
