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