漫画：跑上百万次代码验证三门问题

三门问题（Monty Hall problem）亦称为蒙提霍尔问题、蒙特霍问题或蒙提霍尔悖论，出自美国的电视游戏节目Let's Make a Deal。今天为大家进行完整分析。

话不多说，直接看题目。

01三门问题

三门问题：参赛者的面前有三扇关闭着的门，其中一扇的后面是天使，选中后天使会达成你的一个愿望，而另外两扇门后面则是恶魔，选中就会死亡。

当你选定了一扇门，但未去开启它的时候，上帝会开启剩下两扇门中的一扇，露出其中一只恶魔。（上帝是全能的，必会打开恶魔门）随后上帝会问你要不要更换选择，选另一扇仍然关着的门。

02普通人的直觉

按照常理，参赛者在做出最开始的决定时，对三扇门后面的事情一无所知，因此他选择正确的概率是1/3，这个应该大家都可以想到。

接下来，主持人排除掉了一个错误答案（有恶魔的门），于是剩下的两扇门必然是一扇是天使，一扇是恶魔，那么此时无论选择哪一扇门，胜率都是1/2，依然合乎直觉。

所以你作为参赛者，你会认为换不换都无必要，获胜概率均为1/2。但是，真的是这样吗？

03题目分析

正确的答案是，如果你选择了换，碰见天使的概率会高达2/3，而不不换的话，碰见天使的概率只有1/3。怎么来的？

我们用一个很通俗的方法，能让你一听就懂。首先刚开始选择的一扇门的概率为1/3，而另外两扇门的总概率为2/3。

现在上帝打开了其中一扇为恶魔的门，我们知道这个门后面不会再有天使，所以相当于这部分概率被第三个门持有。

剩下的那扇门的概率（2/3）相当于刚开始选择的门（1/3）的二倍。所以我们得换。

如果还没有听懂。我们可以假设有一百扇门，里边有99只都是恶魔。现在你随机选择一扇门，选择到天使的概率是1/100。

这时，上帝打开其中的98扇，里边都是恶魔。这时候就相当于99/100的概率都集中在了另一扇门里。自然，我们需要选择换。

04贝叶斯证明 代码证明

为了验证结果，我用代码跑了一百万次。什么？用贝叶斯分析分析！太俗（请隔壁去找李永乐老师），咱们还是直接上代码。

 1func main() {
 2    //换门遇见天使的次数和不换门遇见天使的次数
 3    changeAngelCount, unchangeAngelCount := 0, 0
 4    for i := 0; i < 1000000; i++ {
 5        //门的总数
 6        doors := []int{0, 1, 2}
 7        //天使门和选中的门
 8        angelDoor, selectedDoor := rand.Intn(3), rand.Intn(3)
 9        //上帝移除一扇恶魔门
10        for j := 0; j < len(doors); j++ {
11            if doors[j] != selectedDoor && doors[j] != angelDoor {
12                doors = append(doors[:j], doors[j+1:]...)
13                break
14            }
15        }
16        //统计
17        if selectedDoor == angelDoor {
18            unchangeAngelCount++
19        } else {
20            changeAngelCount++
21        }
22    }
23    fmt.Println("不换门遇见天使次数:", unchangeAngelCount, "比例：", (float32(unchangeAngelCount) / 1000000))
24    fmt.Println("换门遇见天使次数:", changeAngelCount, "比例：", (float32(changeAngelCount) / 1000000))
25}

跑了一百万次，结果当然不让我们失望！

所以，今天的问题你听明白了吗？评论区留下你的想法吧！

"要进学习群的小伙伴记得加我，和小浩做朋友"

左侧为私人名片，右侧为公号名片

加我之前怎么得关注一下吧

长按扫一下

这个二维码是我用PS+Python做的

酷不酷炫

关注领取一套算法教程+《图解算法》

发现了吗？近日小浩已经将所有的广告都关掉了！也就是说，所有的文章，对我一点收入都不会有。我唯一想要的，就是大家可以顺手帮我转发，足够我高兴一整天。这也是支持我一直写下去的动力，感谢你。