This Markov Chain Monte Carlo simulation converges at large enough times to the uniform distribution over partitions of the number $N$. The limit shape of the distributions for large $N$ is $e^{-x/a}+e^{-y/a}=1$, where $a=\sqrt{6N}/\pi$.
N = redraws =
MC steps per redraw =
show limit shape show means
This work is licensed under a Creative Commons Attribution 3.0 United States License