Looking for an explanation or insight on a form of the inverse of a restricted gamma function

I posted this yesterday asking about how to find an inverse of a restricted gamma function. To put it concisely, I was looking at the gamma function in the positive reals after restricting it to be greater than its minimum to make it 1:1 where an inverse does exist (there's a graph in the last post so you can see what I mean better).

After doing more research into finding an inverse, I found this post/answer where the high precision answer is pretty much exactly what I want. The problem is that I don't understand it, how it's derived, where it came from, or anything. Does anybody have any insight on the derivation of this or where it came from or anything? Any help will be appreciated.

There is a magificent approximation. If you look at this question of mine which was about the solution of $$n!=a^n\, 10^k$$ @robjohn proposed (with very detailed explanations) $$\bbox[5px,border:2px solid #C0A000]{n\sim ea\exp\left(\operatorname{W}\left(\frac k{ea}\log(10)-\frac1{2ea}\log(2\pi a)\right)\right)-\frac12}\tag{4}$$ where $$W(.)$$ is Lambert function.
So, for your case, make $$a=1$$ and $$k=\frac{\log (n!)}{\log (10)}$$ to get $$n \sim e \exp\left(W\left(\frac{2 \log (n!)-\log (2 \pi )}{2 e} \right)\right)-\frac 12$$
Just for the fun of it, try with $$n=123$$ the result would be $$122.9999299$$ !
If we know that the result should be an integer, then a good inverse of $$\Gamma(n+1)$$ is given by $$n\sim e\exp\left(\operatorname{W}\left(\frac1{e}\log\left(\frac{\Gamma(n)}{\sqrt{2\pi}}\right)\right)\right)+\frac12$$