As I said in a comment, there is a solution involving the generalized Lambert function. If you look here, C. B. Corcino and R. B. Corcino proved
that in the asymptotic expression of the generalized Bell numbers
one needs to solve the equation $$x\,e^x + r x=n$$ which is exactly your case. The remaining of the paper addresses the more general case of $$e^{-cx}=a \frac{ \prod_{i=1}^n (x-t_i)} {\prod_{i=1}^m (x-s_i) }$$
Loooking at the posted equation, from a purely numerical point of view, Newton method will converge quite fast using $t_0=1$ observed by inspection. The successive iterates would be
$$\left(
\begin{array}{cc}
0 & 1.000000000 \\
1 & 1.063499184 \\
2 & 1.060100663 \\
3 & 1.060090320
\end{array}
\right)$$ which is the solution for ten significant figures.
We also could use $[1,n]$ Padé approximants built around $t=1$ and get explicit approximations of the solution. For example
$$t_{(1)}=1+\frac{2 (3-e) (2 e-1)}{2+e+5 e^2}\approx 1.0599978$$
$$t_{(2)}=1+\frac{3 (3-e) \left(2+e+5 e^2\right)}{2 \left(-3+9 e+15 e^2+8 e^3\right)}\approx 1.0600921$$
$$t_{(3)}=1+\frac{8 (3-e) \left(-3+9 e+15 e^2+8 e^3\right)}{24-21 e+291 e^2+265 e^3+65 e^4}\approx 1.0600903$$