Yesterday, nobody contributed anything.  :-(
The official answer is: ln(x * (abs(floor(sqrt(-x))) + x)).