Approximate binomial by normal (#877)

Approximate binomial by normal (#877)

Use μ=np\mu = np, σ=np(1−p)\sigma = \sqrt{np(1-p)} with continuity correction.