张宇/王式安 考研数学 概率论争议题 【Python验证版】

张宇/王式安 概率论争议题 【Python版】

之前写过一篇关于这道题的Java验证版，但是当时的Java代码比较混乱，在这里补充一个Python代码。同样地，由于技术水平限制，这篇代码也存在很大的局限性，欢迎交流指正~

原文链接： 张宇/王式安 考研数学 概率论争议题 【Java验证版】

一、题目

本题为1987年实考题
设两箱内装有同种零件，第一箱装50件，其中有10 件一等品，
第二箱装30件，其中有18件一等品，先从两箱中任挑一箱，
再从此箱中前后不放回任取两个零件，求：
		（1）先取出的零件是一等品的概率p．
		（2）在先取出的是一等品的条件下，后取的仍是一等品的条件概率q．

二、争议

如果你看到这篇文章，估计已经多少听说过这个争议，在此简单介绍即可。
张宇和王式安第一问答案一样，无争议，为0.4。
第二问，张宇认为，应该用条件概率公式，答案约为0.4855.
但王式安先生认为，应该用全概率公式，答案为0.3849.

三、个人结论

本题思想简单，题目中也无明显歧义，简单重复验证即可。
个人觉得，应该是王式安先生错了，错误的原因在于忽略了“第一次取的是一等品”这一事件已经发生。
（应考虑贝叶斯，执果索因）

四、代码验证

1.代码思想

1.随机*抽取times次
    *此处的随机，我们约定为两个箱子中任取一个箱子进行实验并记录为样本
2.统计成功次数并计算概率

2.实现代码

'''
设两箱内装有同种零件，第一箱装50件，其中有10 件一等品，
第二箱装30件，其中有18件一等品，先从两箱中任挑一箱，
再从此箱中前后不放回任取两个零件，求：
		（1）先取出的零件是一等品的概率p．
		（2）在先取出的是一等品的条件下，后取的仍是一等品的条件概率q．
'''
#初始化两箱零件
#按照题意，实际上设前XX件为正品即可，但为了模拟真实场景，采用random生成
#0 为次品 1 为一等品
import random

def defQuality(box,list_random):
    for i in list_random:
        box[i]=1


box1=[0 for i in range(50)]
box2=[0 for i in range(30)]
list_random_1=random.sample(range(0,50,1),10)
list_random_2=random.sample(range(0,30,1),18)
defQuality(box1,list_random_1)
defQuality(box2,list_random_2)


#初始化完成
#求在先取出的是一等品的条件下，后取的仍是一等品的条件概率q．
list_box=[box1,box2]

def takeOut(box):#取出操作
    list_random_takeOut=random.sample(range(0,len(box),1),2)
    #主要是为了保证取到不相等的两个随机数,暂不重点考虑性能
    #如果您有更好的办法，请您在评论区指点一二，感激不尽
    firT=list_random_takeOut[0]
    secT=list_random_takeOut[1]
    if box[firT]==1:
        if(box[secT]==0):
            return 1,0
        if(box[secT]==1):
            return 1,1
    else:
        return 0,0

def test(times):#按次数测试
    isQualityCount,fir_is_quality=0,0
    for i in range(0,times):
        temp1,temp2=takeOut(list_box[random.randint(0,1)])
        fir_is_quality+=temp1
        isQualityCount+=temp2
    return isQualityCount/fir_is_quality

print("\nPlease input the times you want to test:")
inputTimes=int(input())
print("A total of %d tests results a success rate of %f."%(inputTimes,test(inputTimes)))

3.测试结果

A total of 1000 tests results a success rate of 0.515385.

A total of 10000 tests results a success rate of 0.495892.

A total of 100000 tests results a success rate of 0.485057.


显然，测试结果更倾向于张宇的答案，随着房本越大，概率点越靠近0.485这个答案

水平有限，如有错误，请您批评指正。