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