La fórmula principal que se necesita aquí es
$(ST)^{-1} = T^{-1}S^{-1}, \tag 0$
que se cumple para dos invertibles cualesquiera $N \times N$ matrices $S$ y $T$ En lo que sigue aplicaremos repetidamente este conocido resultado.
Dado que
$(X + XB)^{-1} = BX^{-1} + B, \tag 1$
observamos vía (0) que el lado izquierdo puede escribirse
$(X + XB)^{-1} = (X(I + B))^{-1} = (I + B)^{-1}X^{-1}; \tag 2$
por lo que (1) se convierte en
$(I + B)^{-1}X^{-1} = BX^{-1} + B, \tag 3$
de donde
$(I + B)^{-1}X^{-1} - BX^{-1} = B, \tag 4$
o
$((I + B)^{-1} - B)X^{-1} = B; \tag 5$
volviendo a tomar las inversas utilizando (0) se obtiene
$X((I + B)^{-1} - B)^{-1} = B^{-1}, \tag 6$
y después de un poco más de álgebra ahora podemos aislar $X$ :
$X = B^{-1}((I + B)^{-1} - B) = B^{-1}(I + B)^{-1} - B^{-1}B$ $= B^{-1}(I + B)^{-1} - I = ((I + B)B)^{-1} - I = (B + B^2)^{-1} - I. \tag 7$