An alternative approach, which doesn't use the explicit formula in terms of eigenvalues. Instead, find the characteristic polynomial. See that, for any $\lambda$,
$$p(\lambda) = \det (B - \lambda I) = -\lambda^3 + \lambda^2 \mathrm{Tr} \, B - \lambda x + \det B$$
for some number $x$. Here we merely have $\lambda = -1$.
The number $x$ can be computed in a variety of ways. Using exterior algebra is one method: construct the $3 \times 3$ matrix $B_2$, which acts on $3 \times 1$ column vectors corresponding to elements of $\bigwedge^2 \mathbb R^3$. Then $x = \mathrm{Tr} \, B_2$.
The relationship between $B_2$ and $B$ is explicitly
$$B_2 (a \wedge b) = B(a) \wedge B(b)$$
for any vectors $a, b$.
Now, use a common inversion identity:
$$B^{-1}(a) = \star B_2^T(\star a)/\det B$$
where $\star$ is the usual Hodge dual. This means we can write $B_2$ as
$$B_2(a \wedge b) = \star (B^T)^{-1}(\star [a \wedge b]) \det B$$
You can verify now (e.g. by breaking into a basis) that $\mathrm{Tr} \, B_2 = \det B \, \mathrm{Tr} \, (B^T)^{-1}$. By the arguments given in other answers, this is merely $\det B \, \mathrm{Tr} \, B$, and as a result, we have
$$p(\lambda) = \det(B - \lambda I) = -\lambda^3 + \lambda^2 \mathrm{Tr} \, B - \lambda [\det B \, \mathrm{Tr} \, B] + \det B$$
For $\lambda =-1$, and since $\det B = 1$, the result follows.