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