Does anyone have an elegant/slick proof ?
Like in the answer given to the link you quote by Kiril, the fastest proof is given by fixed point iteration. The function you iterate, is:
$$g_c(x)=c^x$$
The fixed points of $g_c(x)$ are given in terms of the Lambert $W$ function as:
$$x_0=\frac{W(-\ln(c))}{-\ln(c)}$$
The fixed point condition (for attraction hence convergence) is therefore:
$$|g'(x_0)|\le 1\Rightarrow$$
$$|-W(-\ln(c))|\le 1 (*)$$
Now consider the function:
$$m(x)=x\cdot\exp(x)$$
$$(*)\Rightarrow W(-\ln(c))\in [-1,1]\Rightarrow$$
$$m(W(-\ln(c)))\in m([-1,1])\Rightarrow$$
$$m(W(-\ln(c)))\in [-e^{-1},e]\Rightarrow$$
$$-\ln(c)\in[-e^{-1},e]\Rightarrow$$
$$c\in[e^{-e},e^{1/e}]$$
Fixed point iteration, therefore, takes care of all cases where $c\in(e^{-e},e^{1/e})$ (open interval). The two end points have to be checked separately, because the fixed point condition is inconclusive there, being possibly $\pm1$.
The upper end point ($e^{1/e}$) is easy using standard Calculus (Consider the sequence $a_n=g_{e^{1/e}}^{(n)}(1)$ and prove that it converges to $e$, etc.). The lower end point ($e^{-e}$), is a bit tricky, but can also be done using Calculus, and you are done.