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