hypergeometric.py 311 B

123456789101112
  1. #超几何分布
  2. from fractions import Fraction
  3. from math import comb
  4. N=11 #总数
  5. M=2 #次品数
  6. n=3 #抽取数
  7. m=max(0,n-N+M)
  8. r=min(n,M)
  9. print("EX=",Fraction(n*M,N))
  10. print("DX=",Fraction(n*M*(N-M)*(N-n),N*N*(N-1)))
  11. for k in range(m,r+1):
  12. print("P(X=",k,")=",Fraction(comb(M,k)*comb(N-M,n-k),comb(N,n)))