Python使用random和tertools模块解一些经典概率问题

combinations(sequence, k))
从序列 sequence 中得到包含 k 个元素的所有组合。

def once(doors = 3):  # 一次事件的模拟
 car = randrange(doors) # 一扇门后面停着汽车
 man = randrange(doors) # 参赛者预先选择一扇门
 return car == man # 参赛者是否最初就选择到车

h = 0 # 坚持选择赢得汽车的次数                   
c = 0 # 改变选择赢得汽车的次数
times = int(1e6) # 重复实验的次数

for i in range(times):
 if once(): h += 1
 else:  c += 1

print("维持选择：",h/times*100,"%\n改变选择：",c/times*100,"%")

cards = [i for i in range(52)]
counter = [0, 0, 0, 0]

def once(): # 0 表示黑桃 A
 global cards
 ace = set(sample(cards, 13)) & {0,1,2,3}
 return len(ace), 0 in ace

for i in range(int(1e6)):
 a, s = once() # a 表示 A 的个数， s 表示是否有黑桃 A
 if a:
  counter[1] += 1
  if s: counter[3] += 1
 if a > 1:
  counter[0] += 1
  if s: counter[2] += 1

print('情况一：', counter[0]/counter[1], '\n情况二：', counter[2]/counter[3])

family = (lambda n :[{randrange(2),randrange(2)} for i in range(n)])(int(1e6))

both = family.count({0}) # 都是女孩的家庭数
exist = len(family) - family.count({1}) # 有女孩的家庭数

print(both/exist)

counter, times = 0, int(1e6)
for i in range(times):
 if len({randrange(365) for i in range(50)}) != 50: # 存在同一天生日的人
  counter += 1

print('在 50 个人中有相同生日的概率为：',counter/times)