Dado un cubo de $n \times n \times n$ o $[0,\, n]^3$ queremos encontrar el avión $ax+by+cz=d$ que cruza
el mayor número de cubos unitarios dentro de $[0,\, n]^3$, y encontrar ese número.
Podemos hallar una sola unidad de cubo de la $[x_k,\, x_k+1] \times [y_j,\, y_j+1] \times [z_l,\, z_l+1]$, con $j,k,l \in [0, \, n-1]$.
Los cubos atravesado por el avión va a ser aquellos para los cuales
$$
\eqalign{
& ax_{\,k} + by_{\,j} + cz_{\,l} < d < a\left( {x_{\,k} + 1} \right) + b\left( {y_{\,j} + 1} \right) + c\left( {z_{\,l} + 1} \right)\quad \Rightarrow \cr
& \Rightarrow \quad d - \left( {a + b + c} \right) < ax_{\,k} + by_{\,j} + cz_{\,l} < d\quad \Rightarrow \cr
& \Rightarrow \quad {d \over {a + b + c}} - 1 < {{ax_{\,k} + by_{\,j} + cz_{\,l} } \over {a + b + c}} < {d \over {a + b + c}} \cr}
$$
Considere la posibilidad de $x_k$ como la realización de un uniforme discreta variable aleatoria $x$ sobre el apoyo a la $[0,\, n-1]$,
con una probabilidad de $1/n$, la media de $(n-1)/2$ y la varianza $(n^2-1)/12$ .
Mismo para $y, \, z$.
Su suma ponderada
$$
{{ax_{\,k} + by_{\,j} + cz_{\,l} } \over {a + b + c}}
$$
tendrá media, moda y mediana en $(n-1)/2$ y la varianza
$$
\sigma ^2 = {{a^2 + b^2 + c^2 } \over {\left( {a + b + c} \right)^2 }}\left( {{{n^2 - 1} \over {12}}} \right)
$$
Claramente, la menor es la varianza de la más grande es la porción de la pmf que satisface la desigualdad dada anteriormente, puesto que la desigualdad del medidor es constante en $1$.
Y la varianza es claramente mínimo para la igualdad de pesos.
Así llegamos a considerar la desigualdad
$$ \bbox[lightyellow] {
\left\{ \matriz{
x_{\,k} ,y_{\,j} ,z_{\,l} ,n,s \in \mathbb Z \hfill \cr
d \in \mathbb R \hfill \cr
0 \le x_{\,k} ,y_{\,j} ,z_{\,l} \le n - 1 \hfill \cr
d - 3 < x_{\,k} + y_{\,j} + z_{\,l} = s < d \hfill \cr} \right.
\etiqueta{1}}$$
Ahora, el número de puntos en el plano diagonal de un $m$-D cube
$$N_{\,b} (s,r,m) = \text{No}\text{. de soluciones a}\;\left\{ \begin{gathered}
0 \leqslant \text{integer }x_{\,j} \leqslant r \hfill \\
x_{\,1} + x_{\,2} + \cdots + x_{\,m} = s \hfill \\
\end{reunieron} \right.$$
está dada por
$$ \bbox[lightyellow] {
N_b (s,r,m)\quad \left| {\;0 \leqslant \text{números enteros }s,m,r} \right.\quad =
\sum\limits_{\left( {0\, \leqslant } \right)\,\,k\,\,\left( { \leqslant \,\frac{s}{i+1}\, \leqslant \,m} \right)}
{\left( { - 1} \right)^k \binom{m}{k}
\binom
{ s + m - 1 - k\left( {i + 1} \right) }
{ s - k\left( {i + 1} \right)}\ }
\etiqueta{2.un}}$$
como se explica en este post.
Por otra parte el número de puntos en o por debajo de la diagonal del plano son
$$ \bbox[lightyellow] {
\eqalign{
& M_b (s,r,m)\quad \left| {\;0 \le {\rm enteros }s,m,r} \right.\quad = \cr
& = \sum\limits_{0\, \le \,\,j\,\, \le \,s} {N_b (s,r,m)} = \cr
& = \sum\limits_{\left( {0\, \le } \right)\,\,k\,\,\left( { \le \,{s \over {i + 1}}\, \le \,m} \right)} {\left( { - 1} \right)^k \left( \matriz{
m \hfill \cr
k \hfill \cr} \right)\left( \matriz{
s + m - k\left( {i + 1} \right) \cr
s - k\left( {i + 1} \right) \cr} \right)} \cr}
\etiqueta{2.b}}$$
En este punto necesitamos la ayuda de una visualización gráfica para comprender el comportamiento de la desigualdad 1) wrt $N_b$
El dibujo representa los histogramas de $N_{\,b} (s,n-1,3)$ para $n=3$ e $n=4$.
$N_{\,b} (s,n-1,3)/n^3$ es el pmf de la suma de $s$ de los tres uniformes discretas variables aleatorias.
El dibujo muestra que la mayor parte de el histograma es interceptado cuando el medidor de ancho de $3$ de la desigualdad
es casi centrado alrededor de la media.
Que es en realidad, así que cuando n es impar, mientras que para incluso $n$ vamos a cambiar el calibre ligeramente a la izquierda (o a la derecha).
Por desgracia, la fórmula para $N_b$ sólo es válida para la integral parámetros (reescritura de la binomial a través de gamma produce una función discontinua).
Podemos eludir el de arriba y el uniforme de la desigualdad mediante la introducción de un fijo $1/2$ cambio de la media y, a continuación, la reescritura de la desigualdad como
$$
\eqalign{
& d - 3 < x_{\,k} + y_{\,j} + z_{\,l} = s < d\quad \Rightarrow \cr
& \Rightarrow \quad 3{{n - 1} \over 2} - 3/2 - 1/2 < s \le 3{{n - 1} \over 2} + 3/2 - 1/2\quad \Rightarrow \cr
& \Rightarrow \quad \left\lfloor {3{{n - 1} \over 2} - 3/2 - 1/2} \right\rfloor
< s \le \left\lfloor {3{{n - 1} \over 2} + 3/2 - 1/2} \right\rfloor \cr}
$$
y en general, para una dimensión $m$
$$ \bbox[lightyellow] {
\eqalign{
& d - m < x_{\,k} + y_{\,j} + z_{\,l} = s < d\quad \Rightarrow \cr
& \Rightarrow \quad m{{n - 1} \over 2} - m/2 - 1/2 < s \le m{{n - 1} \over 2} + m/2 - 1/2\quad \Rightarrow \cr
& \Rightarrow \quad \left\lfloor {{{mn - 1} \over 2}} \right\rfloor - m
< s \le \left\lfloor {{{mn - 1} \over 2}} \right\rfloor \cr}
\etiqueta{3}}$$
lo que conduce a
$$ \bbox[lightyellow] {
N(n,m) = M_b \left( {\left\lfloor {{{mn - 1} \over 2}} \right\rfloor ,\;n - 1,\;m} \right)
- M_b \left( {\left\lfloor {{{mn - 1} \over 2}} \right\rfloor - m,\;n - 1,\;m} \right)
\etiqueta{4}}$$
Los valores para los más pequeños de $m$ e $n$ dado por la fórmula se
que comprobar en directo la computación.