跑上百万次代码,终于把三门问题验证了
在帅地玩编程公众号后台回复 剑指offer，送你一本剑指offer书籍
来源：小浩算法
作者：程序员浩哥
三门问题（Monty Hall problem）亦称为蒙提霍尔问题、蒙特霍问题或蒙提霍尔悖论，出自美国的电视游戏节目Let's Make a Deal。今天为大家进行完整分析。
话不多说，直接看题目。
三门问题：参赛者的面前有三扇关闭着的门，其中一扇的后面是天使，选中后天使会达成你的一个愿望，而另外两扇门后面则是恶魔，选中就会死亡。
当你选定了一扇门，但未去开启它的时候，上帝会开启剩下两扇门中的一扇，露出其中一只恶魔。（上帝是全能的，必会打开恶魔门）随后上帝会问你要不要更换选择，选另一扇仍然关着的门。
按照常理，参赛者在做出最开始的决定时，对三扇门后面的事情一无所知，因此他选择正确的概率是1/3，这个应该大家都可以想到。
接下来，主持人排除掉了一个错误答案（有恶魔的门），于是剩下的两扇门必然是一扇是天使，一扇是恶魔，那么此时无论选择哪一扇门，胜率都是1/2，依然合乎直觉。
所以你作为参赛者，你会认为换不换都无必要，获胜概率均为1/2。但是，真的是这样吗？
正确的答案是，如果你选择了换，碰见天使的概率会高达2/3，而不不换的话，碰见天使的概率只有1/3。怎么来的？
我们用一个很通俗的方法，能让你一听就懂。首先刚开始选择的一扇门的概率为1/3，而另外两扇门的总概率为2/3。
现在上帝打开了其中一扇为恶魔的门，我们知道这个门后面不会再有天使，所以相当于这部分概率被第三个门持有。
剩下的那扇门的概率（2/3）相当于刚开始选择的门（1/3）的二倍。所以我们得换。
如果还没有听懂。我们可以假设有一百扇门，里边有99只都是恶魔。现在你随机选择一扇门，选择到天使的概率是1/100。
这时，上帝打开其中的98扇，里边都是恶魔。这时候就相当于99/100的概率都集中在了另一扇门里。自然，我们需要选择换。
为了验证结果，我用代码跑了一百万次。什么？用贝叶斯分析分析！太俗（请隔壁去找李永乐老师），咱们还是直接上代码。
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}
跑了一百万次，结果当然不让我们失望！
推荐阅读