Without use of the Lambert W function (you will be hard put to find implementations outside of Mathematica, Maple, or Matlab) there is a very simple technique you can use to find the real solutions to
[tex]x - \log x = c,\,\,c>=1[/tex]
Fixed point iteration can be used to find both solutions. Note that the above is equivalent to [itex]x=\log x + c[/itex]. Writing this as a fixed point iteration scheme,
[tex]x_{n+1} = \log x_n + c[/tex]
Simply start with some initial guess [itex]x_0[/itex] and iteratively generate [itex]x_1, x_2,[/itex] etc. using the above. If the sequence [itex]{x_0, x_1, x_2, \cdots}[/itex] does converge, it will converge at one of the solution points. An initial guess of 1 will always converge. It will take a while if c is fairly close to 1. Fixed point iteration is very simple but it can be very slow to converge.
In fact, this will only find the larger of the two solutions. So how to find the other solution? Since exp(x) is a bijection, we can go from [itex]x-\log x = c[/itex] to [itex]\exp(x)/x = \exp(c)[/itex], or [itex]x = \exp(x-c)[/itex]. This suggests another fixed point iteration scheme,
[tex]x_{n+1} = \exp(x_n-c)[/tex]
If this sequence converges it will converge to the smaller of the two solutions, and once again, an initial value of 1 will always result in a convergent sequence. Once again the convergence is rather slow when c is just slightly larger than 1.