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