For the first part, look at the convexity of the function at the extremum.
For the second part of the question : if we consider $$F = c\Big(x\log(zy) + (1-x)\log(z(1-y))\Big)= c\Big(x\log(y) + (1-x)\log((1-y))+\log(z)\Big)$$ I suppose that we can forget $c$ and $z$ as long as we search for extremum conditions (I assume $c\neq 0$ and $z\neq 0$) . So, the critical points do not depend on $c$ or $z$.
Now, if we compute the partials, we have $$F'_x=\log (y z)-\log ((1-y) z)=\log (y )-\log (1-y)$$ $$F'_y=\frac{x}{y}-\frac{1-x}{1-y}$$ Setting the first equal to $0$ leads to $y=\frac{1}{2}$ and replacing in the second it cancels for $x==\frac{1}{2}$ too.
This just to confirm that you are totally right. It is very good you notice points like this one.
Your intuition proves to be correct, but you would need to give something a bit more formal in the solution to an exam problem. Your proposed constraint isn't strong enough since there is nothing in the inequality that forbids any one (or even two) of the edge measurements to exceed 8 meters.
The given constraint requires that the solution(s) for the maximum cost lie inside the cube $ \ [ 0 \ , \ 8 ] \ \times \ [ 0 \ , \ 8 ] \ \times\ [ 0 \ , \ 8 ] \ $ . This constraint is essential since the one for the volume of the metal box, $ \ xyz \ = \ 8 \ $ , alone would also any one or two of the edge measurements to "run to infinity". The "hyperbolic" surface is open, so the cost function, which is a linear combination of pairs of length variables, could be made arbitrarily large.
the open surface xyz = 8 with the location of the minimal cost indicated
$ \ \ $
With the "maximum lengths" constraints, the volume constraint surface intersects the planes $ \ x \ = \ 8 \ \ , \ \ y \ = \ 8 \ \ , \ \ z \ = \ 8 \ $ at three hyperbolic curves $ \ yz \ = \ 1 \ \ , \ \ xz \ = \ 1 \ \ , \ \ xy \ = \ 1 \ $ , respectively. This now requires the variables to lie in the intervals $ \ \frac{1}{8} \ \le \ x \ \le \ 8 \ $ and similarly for $ \ y \ $ and $ \ z \ $ in order to work on all three curves simultaneously.
the constraint planes added, with the intersection hyperbolas visible and the location of the maximal cost indicated
$ \ \ $
The Lagrange-multiplier method will only locate the one point where the cost function $ \ C \ = \ 8 \ xy \ + \ 4 \ xz \ + \ 2 \ yz \ $ can be tangent to $ \ xyz \ = \ 8 \ $ , which occurs for $ \ C \ = \ 48 \ $ , as you found. There are no other places where the "cost surface" will be tangent to the open volume constraint surface for any value of $ \ C \ $ .
Instead, we are left to check the curves on the three individual planes. We find on the plane $ \ z \ = \ 8 \ $ , for instance, that the cost function becomes $ \ C \ = \ 8 \ xy \ + \ 4 \ x \ \cdot \ 8 \ + \ 2 \ y \ \cdot \ 8 \ $ . Together with the hyperbolic equation $ \ xy \ = \ 1 \ $ , we have $ \ C \ = \ 8 \ \cdot 1 \ + \ 32 \ x \ + \ 16 \ \frac{1}{x} \ $ ; but this reduced function has a positive second derivative at its "critical point", so this is not a maximum value for the cost function. We obtain similar results on the other two planes, $ \ y \ = \ 8 \ $ and $ \ z \ = \ 8 \ $ .
This leaves us to check the three hyperbolic curves at each of their two endpoints:
z = 8 :
$$ \ C \ = \ 8 \ (1) \ + \ 32 \ \cdot \ \frac{1}{8} \ + \ 16 \ \cdot \ 8 \ \ \ , \ \ \ C \ = \ 8 \ (1) \ + \ 32 \ \cdot \ 8 \ + \ 16 \ \cdot \ \frac{1}{8} \ \ ;$$
y = 8 :
$$ \ C \ = \ 8 \ \cdot \ \frac{1}{8} \ \cdot \ 8 \ \ + \ 4 \ ( 1 ) \ + \ 2 \ \cdot \ 8 \ \cdot \ 8 \ \ \ , \ \ \ C \ = \ 8 \ \cdot \ 8 \ \cdot \ 8 \ + \ 4 \ (1) \ + \ 2 \ \cdot \ 8 \ \cdot \ \frac{1}{8}\ \ ;$$
x = 8 :
$$ \ C \ = \ 8 \ \cdot \ 8 \ \cdot \ \frac{1}{8} \ + \ 4 \ \cdot \ 8 \ \cdot \ 8 \ + \ 2 \ (1) \ \ \ , \ \ \ C \ = \ 8 \ \cdot \ 8 \ \cdot \ 8 \ + \ 4 \ \cdot \ 8 \ \cdot \ \frac{1}{8} \ + \ 2 \ (1) \ \ . $$
It is clear from this that in order to maximize the cost, we want the "product of two lengths equal to 1" factors to be multiplied by the lowest coefficients, $ \ 4 \ $ and $ \ 2 \ $ . Consequently, we need to have the product $ \ xy \ = \ 8 \ \cdot \ 8 \ $ , that is, as we recognized intuitively, the gold faces of the box should be squares with side $ \ 8 \ $ meters. This is either the fourth or the sixth of the equations above, giving a maximal cost of $ \ 518 \ $ dollars. (In fact, there are only three distinct cost values at the six indicated endpoints. In the diagram above, we see that this occurs because pairs of "cost hyperbolas" meets at just three locations.)
Best Answer
In fact you got quite a good "guess". When it comes to finding maxima and minima of a function $f:\mathbb{R}^2 \to \mathbb{R}$, you begin by finding the critical points which, as you said, are obtained by solving the system: $$\begin{cases} f'_x(x_0, y_0) = 0 \\ f'_y(x_0, y_0) = 0\end{cases}$$ The solution(s), $(x_0, y_0)$, is either going to be maxima, minima or saddlepoint. You determine which one it is by studying the Hessian matrix evaluated at that point.
I'm sure you will get to it when you start doing multivariable calculus. Great guess, however!