Doy una derivación geométrica, queremos contar integral de todos los puntos de la mentira en la superficie $|x|+|y|+|z|=P$. En realidad en el espacio 3D, en cada octante, la forma de la superficie es de forma triangular.
Por ejemplo, si $P=4$, entonces la forma en el primer octante sería
$$
\begin{array}[ccccccccc]
\ & & & & Q(0,0,4)& & & &\\
& & &D(0,1,3))& &D(1,0,3)& & &\\
& &D(0,2,2)& &S(1,1,2)& &D(2,0,2)& &\\
&D(0,3,1)& &S(1,2,1)& &S(2,1,1)& &D(3,0,1)&\\
Q(0,4,0)& &D(1,3,0)& &D(2,2,0)& &D(3,1,0)& &Q(4,0,0)\\
\end{array}
$$
donde $P$ denota los puntos que debe ser compartida por 4 octantes, $D$ denota los puntos que debe ser compartida por 2 octantes y $S$ denota los puntos sólo pertenece a este octantes.
Así que para el total de los octantes, el número de puntos con S
$$n_S=8\cdot\frac{(P-1)(P-2)}{2}=4P^2-12P+8$$
el número de puntos con D es
$$n_D=8\cdot\frac{3(P-1)}{2}=12P-12$$
el número de puntos con Q es
$$n_Q=8\cdot\frac{3}{4}=6$$
Por lo que el número total sería
$$n=n_S+n_D+n_Q=4P^2+2$$