Quant Memo
Foundational

Expected Maximum of Several Dice

Working out the expected value of the largest roll among several dice using the complement trick — via 1 minus the CDF — a technique that comes up constantly for 'expected max' and 'expected best-of-n' questions.

Problem. Roll nn fair six-sided dice. What's the expected value of the largest number showing? For n=1n=1 it's obviously 3.53.5 (the standard expected value of one die), but the answer for n=2n=2, n=3n=3, and beyond isn't obvious by symmetry the way single-die expectation is — and trying to sum "probability of each value times that value" directly gets messy fast. There's a cleaner route.

The complement trick: E[X]=kP(Xk)E[X] = \sum_k P(X \geq k)

For a random variable XX taking positive integer values from 11 to 66, there's a useful identity:

E[X]=k=16P(Xk).E[X] = \sum_{k=1}^{6} P(X \geq k) .

In plain English: instead of weighting each outcome by its value and summing, you can equivalently sum, for each threshold kk from 1 to 6, the probability that XX reaches at least that threshold. This is easier here because "the max of nn dice is at least kk" has an easy complement: "the max is at least kk" fails only if every die shows less than kk, i.e., every die is in {1,,k1}\{1, \dots, k-1\}, which happens with probability (k16)n\left(\frac{k-1}{6}\right)^n. So

P(maxk)=1(k16)n.P(\max \geq k) = 1 - \left(\frac{k-1}{6}\right)^n .

Worked example 1: two dice

For n=2n=2: P(max1)=1P(\max \geq 1) = 1 (trivially, since max is always 1\geq 1), P(max2)=1(1/6)2=35/36P(\max \geq 2) = 1 - (1/6)^2 = 35/36, P(max3)=1(2/6)2=32/36P(\max \geq 3) = 1 - (2/6)^2 = 32/36, P(max4)=1(3/6)2=27/36P(\max \geq 4) = 1-(3/6)^2 = 27/36, P(max5)=1(4/6)2=20/36P(\max \geq 5) = 1-(4/6)^2 = 20/36, P(max6)=1(5/6)2=11/36P(\max \geq 6) = 1-(5/6)^2 = 11/36. Summing:

E[max]=36+35+32+27+20+1136=161364.472.E[\max] = \frac{36+35+32+27+20+11}{36} = \frac{161}{36} \approx 4.472 .

Compare to the single-die expectation of 3.53.5: taking the max of two dice pulls the expectation up noticeably, which matches intuition — the larger of two rolls tends to land above the midpoint.

Worked example 2: how it scales with more dice

The same formula extended to n=3n=3 dice gives E[max]4.958E[\max] \approx 4.958, and for n=4n=4, E[max]5.27E[\max] \approx 5.27. As nn grows, the expected maximum climbs toward 6 but never reaches it, since there's always some (shrinking) chance every die falls short. This "expected max grows and saturates near the top of the range" pattern is the same shape you'd expect from order statistics in general — the maximum of more samples from a bounded distribution creeps toward the upper bound but the increments shrink (going from n=1n=1 to n=2n=2 adds about 0.970.97; from n=3n=3 to n=4n=4 adds only about 0.310.31).

Distribution · binomial
mean 3.00123456outcomes (k) →
mean 3.00std dev 1.22peak at k = 3

Use the explorer above to build intuition for how a bounded outcome's distribution compresses toward one tail as you effectively "select the best of several trials" — the same compression that drives E[max]E[\max] toward 6 as nn grows, even though no single die's distribution changes.

For the expected maximum (or any order statistic) of nn independent bounded random variables, use E[X]=kP(Xk)E[X] = \sum_k P(X \geq k) and compute each P(maxk)P(\max \geq k) via the complement — "all dice below kk" — which is usually far easier than summing kP(max=k)k \cdot P(\max = k) directly. For nn six-sided dice, P(maxk)=1(k16)nP(\max \geq k) = 1 - \left(\frac{k-1}{6}\right)^n.

This complement trick isn't specific to dice — it's the standard tool for "expected best of nn" questions across finance too, such as the expected best price you'll get quoting to nn competing counterparties, or the expected maximum drawdown-free run length in a simulated path.

Related concepts

Practice in interviews

Further reading

  • Ross, A First Course in Probability, ch. 7
ShareTwitterLinkedIn