Yesterday, nobody contributed anything. :-( The official answer is: floor(ln(2 * x)) - 1 / (1 - 1 / x) ^ 2.