Deje $I_i$ ser el indicador de la variable aleatoria que denota si la i-ésima punto en $S$ es pareto óptimo. A continuación, el número de pareto-óptima puntos en $S$ $M = \sum_{i=1}^{n} I_i$ $E[M] = \sum_{i=1}^n E[I_i]$ por la linealidad de las expectativas. Desde $I_i$ es un binario con valores de variable aleatoria, también tenemos $E[I_i] = P(I_i = 1)$. Poner esto juntos, llegamos $E[M] = \sum_{i=1}^n P(I_i = 1)$.
Para evaluar $P(I_i=1)$, tenemos que averiguar la probabilidad de que no hay ningún punto por encima y a la derecha de la i-ésima punto. Esto puede ser hecho por primera acondicionado en la i-ésima del punto de ubicación dentro de la plaza de la unidad y, a continuación, exigiendo que todos los restantes $(n-1)$ puntos no se encuentran en el pequeño rectángulo de arriba y a su derecha. Si la i-ésima punto es $(x_i,y_i)$, queremos que ninguno de los restantes $(n-1)$ puntos en el rectángulo de lados a$1-x_i$$1-y_i$. La probabilidad de que un punto al azar, se encuentra en este rectángulo es simplemente su área y de la probabilidad de que ninguno de los $(n-1)$ puntos se encuentran en este rectángulo es, por tanto,$(1-(1-x_i)(1-y_i))^{n-1}$.
\begin{eqnarray}
P(I_i=1 \mid \text{ith pt} = (x_i,y_i)) &=& (1-(1-x_i)(1-y_i))^{n-1} \\
P(I_i=1) &=& \int_0^1 \int_0^1 (1-(1-x_i)(1-y_i))^{n-1} dx_i dy_i \\
&=& \int_0^1 \int_0^1 (1-uv)^{n-1} du dv
\end{eqnarray}
Esta integral es bastante fácil de calcular. Mediante la celebración de $v$ constante, podemos evaluar primero el interior de la integral como
$\int_0^1 (1-uv)^{n-1} du = \frac{1}{nv} (1-(1-v)^{n})$
Haciendo la sustitución de $v \rightarrow 1-v$, obtenemos
\begin{eqnarray}
P(I_i = 1) &=& \frac{1}{n} \int_0^1 \frac{1-v^n}{1-v} dv \\
&=& \frac{1}{n} \int_0^1 (1+v+\dots+v^{n-1}) dv \\
&=& \frac{1}{n} \left(1+\frac{1}{2}+\dots+\frac{1}{n}\right) = \frac{H_n}{n}
\end{eqnarray}
Finalmente, calculamos el $E[M] = \sum_{i=1}^n P(I_i=1) = H_n$ donde $H_n$ es el n-ésimo número Armónico.