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