python-numpy--求概率

import numpy as np         #导入numpy
from numpy import random   #导入random

#分析，一共有12个星座，一个宿舍8个人，问一个宿舍有三个人的星座相同的概率

a=random.randint(12,size=(100,8))     #创建一个100行8列的二维数组，数据由0到11随机组成
print(a)  #输出二维数组
aa=0

for i in range(100): #循环100次
    v=np.argmax(np.bincount(a[i]))  #判断二维数组中第行出现次数最多的数
   # print(v)
    z=list(a[i]).count(v)      #将第行转换为列表ii，并求出出现最多的次数
   # print(z)
    if z>=3:     #判断，如果出现次数大于或者等于三，表示宿舍里有三个人或三人以上星座相同
        aa+=1     #加一

print(f'出现的概率是{aa/100}')    #一共测试100次