蒙提霍尔悖论源码验证

蒙提霍尔悖论又称三门问题Monty Hall Problem,这个问题出自美国的电视游戏节目Let's Make a Deal。问题名字来自该节目的主持人蒙提·霍尔(Monty Hall)。在这个秀上有三扇门,其中有一扇门打开后可以获得一辆汽车,而其余的两扇门打开后则是山羊。

游戏规则

游戏的步骤如下:

  • 参加游戏的人选定其中的一扇门
  • 主持人选择另外两扇门中的一扇,且打开后只能是山羊,而不能是汽车
  • 给参加游戏的人一个是否交换剩下两扇门的机会
  • 游戏者做出决定,并打开门

这个悖论的关键在于:挑战者是应该选择交换还是选择坚持?这两种做法有没有区别?

从理论的角度分析

一般来说,大多数的挑战者在选择后会认为其实换不换无所谓,其概率是相等的。但是实际上,如果交换的话会有更大的机率获得汽车。现在我们从理论的角度去考虑这个问题:

首先,我们假设有三扇门,分别为门1,门2, 门3,且挑战者选择了门1。

很显然,汽车在每扇门后的概率为:

门1 门2 门3
对应门后有汽车的概率 1/3 1/3 1/3

按照游戏规则,接下来主持人需要为挑战者去除一个干扰项:

门1 门2 门3
对应门后有汽车的概率 1/3 1/3 1/3
去除一个干扰项 去除门2:1/2, 去除门3:1/2 只能去除门3:1 只能去除门2:1

这一步是由挑战者确定是否交换,对于不交换的情况,那只有门1是汽车的时候,挑战者才能拿到汽车的奖励,此时的概率为1/3,而对于交换的情况,只要不是门1的情况,都可以获得奖励,概率为2/3。如下表:

门1 门2 门3
对应门后有汽车的概率 1/3 1/3 1/3
去除一个干扰项 去除门2:1/2, 去除门3:1/2 只能去除门3:1 只能去除门2:1
不交换,获得汽车的概率 P(获得汽车) = 1/3 * 1/2 + 1/3 * 1/2 =1/3 P(获得汽车) = 0 P(获得汽车) = 0
交换后获得汽车的概率 P(获得汽车) = 0 P(获得汽车) = 1/3 * 1 = 1/3 P(获得汽车) = 1/3 * 1 = 1/3

由此可以看出,直觉中换不换概率相等的认识是错误的。

程序模拟源码

前面一节对蒙提霍尔悖论做了一些简单的分析,这一部分会用程序来模拟,并得出了相似的结果。源码如下:

# monty problem
# code is not optimal, but it can demonstrate the problem

import random
def open_the_door( doors, player_choice ):

  # choose which door should be open
  sample = [0,1,2]

  if doors[player_choice]:
    sample.remove(player_choice)
    open_door = random.choice(sample)
  else:
    for index, value in enumerate(doors):
      if index != player_choice and not value:
        open_door = index

  return open_door

def monty_problem( total_times, exchange = False ):

  bingo = 0

  for i in range(total_times):
    
    doors = [False, False, False]
    doors[random.randint(0,2)] = True
    player_choice = random.randint(0,2)
    open_door = open_the_door(doors, player_choice)

    if exchange:
      difference = set([0, 1, 2]) - set([player_choice, open_door])
      player_choice = list(difference)[0]

    if doors[player_choice]:
      bingo += 1

  print('P(exchange={}) = {}'.format(exchange ,bingo / total_times))

if __name__ == '__main__':
  monty_problem(100000, False)
  monty_problem(100000, True)

运行结果如下:

P(exchange=False) = 0.33178
P(exchange=True) = 0.66669

从这里可以看出,不交换和交换的概率分别接近在上一部分的推理结果。这一部分请自行验证。

总结

其实这个问题还可以从另外一个角度,凭直觉给出答案。比如说,现在不仅仅是只有3扇门,而是有100扇,甚至是1000扇门,那么在你选择了其中的一扇门后,打开其余不是汽车的98或998扇门,那么这个时候,你是换还是不换呢?

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 194,761评论 5 460
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 81,953评论 2 371
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 141,998评论 0 320
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,248评论 1 263
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 61,130评论 4 356
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,145评论 1 272
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,550评论 3 381
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,236评论 0 253
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,510评论 1 291
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,601评论 2 310
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,376评论 1 326
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,247评论 3 313
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,613评论 3 299
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,911评论 0 17
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,191评论 1 250
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,532评论 2 342
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,739评论 2 335

推荐阅读更多精彩内容