I have N bernoulli dependent variables, X1, X2, ..., XN, and Xi ~ B(1, pi), pi is already konwn. And Y=X1+X2+...+XN, I need to get the distribution of Y. If Xi, Xj is independent, then I can use the simulation, generate X1,...,XN using the bernoulli distribution, and get the Y; repeat 10000...