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