What part is not clear? The symmetry automatically pulls you into other two situations cyclically.
Let us take Lagrange Multiplier (as others have also done).
I take the unified Lagrangian combining object and constraint functions of
Volume and Area together.
I also choose it such that in $ V - A \cdot \lambda \tag{1}$ $\lambda$ would be
physically a linear dimension for a side of a rectangular parallelepiped,except
for a constant factor.
$ x y z - ( x y + y z + z x) \lambda \tag{2}$
Partial differentiation with respect to x gives $ y z - ( y+z) \lambda =0,\, \dfrac1y + \dfrac1z= \dfrac{1}{\lambda} \tag{3}$
Remember that when number of independent variables are more than 2, partial differentiation should be done with respect to each variable.
So similarly by cyclic symmetry, $ z x -( z + x)\lambda =0 , \,\dfrac1z + \dfrac1x= \dfrac{1}{\lambda} \tag{4} $
and
$ xy -( x+y)\lambda =0 , \,\dfrac1x + \dfrac1y= \dfrac{1}{\lambda} \tag{5} $
Summming up the three and halving,
$ \dfrac1x + \dfrac1y + \dfrac1z= \dfrac{3}{2\lambda} \tag{6} $
Subtracting from this the second part of $ (3), (4), (5)$ we get
$ \dfrac1x = \dfrac{1}{2 \lambda} \tag{7}, x = 2 \lambda $
that gives you
$ x = y = z = 2 \lambda = a , $ say.
So finally $ V = a^3$ and $ A = 6 a^2.$