Correctly rounding expansion factors after scaling

The expansion factors after scaling are currently rounded by simply first taking the floor as the base value and the remainder as a probability of have an extra household. An array of random doubles is then generated and compared to these probabilities to decide whether to add this remaining household.

However, KM used the "Truncate-Replicate-Sample" method in his version. We should consider this maybe in the future.