擅长:python、mysql、java
<p>因为您使用的是scipy,我想我应该提到scipy已经实现了统计分布。还要注意,当n这么大时,二项分布很好地近似于正态分布(或泊松,如果p很小)。在</p>
<pre><code>n = 450000
p = .5
k = np.array([17., 225000, 226000])
b = scipy.stats.binom(n, p)
print b.pmf(k)
# array([ 0.00000000e+00, 1.18941527e-03, 1.39679862e-05])
n = scipy.stats.norm(n*p, np.sqrt(n*p*(1-p)))
print n.pdf(k)
# array([ 0.00000000e+00, 1.18941608e-03, 1.39680605e-05])
print b.pmf(k) - n.pdf(k)
# array([ 0.00000000e+00, -8.10313274e-10, -7.43085142e-11])
</code></pre>