How can I efficiently approximate the sum of Bernoulli random variables for any number of summands in partial sum?