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