Python使用random和tertools模块解一些经典概率问题


Posted in Python onJanuary 28, 2015

random 模块中的常用函数

random()

返回一个位于区间 [0,1] 内的实数;

uniform(a, b)

返回一个位于区间 [a,b] 内的实数;

randint(a, b)

返回一个位于区间 [a,b] 内的整数;

choice(sequence)

返回一个位于 sequence 中的元素,其中,sequence 为一个有序序列,如 list、string 或者 tuple 等类型;

randrange([start], stop[, step])

等效于 choice(range([start], stop[, step]));

shuffle(sequence [, random])

无返回值,用于打乱 sequence 中元素的排列顺序;

sample(sequence, n)

返回一个由 n 个 sequence 中的元素组成的分片,其中,sequence 也可以是 set 类型。

利用 itertools 得到排列、组合

permutations(sequence, k))

从序列 sequence 中得到包含 k 个元素的所有排列。
combinations(sequence, k))

从序列 sequence 中得到包含 k 个元素的所有组合。

羊车门问题

有一个抽奖节目,台上有三扇关闭的门,一扇门后面停着汽车,其余门后都是山羊,只有主持人知道每扇门后面是什么。参赛者可以选择一扇门,在开启它之前,主持人会开启另外一扇门,露出门后的山羊,然后允许参赛者更换自己的选择。问题是:参赛者更换选择后能否增加赢得汽车的机会?

有很多时候,我们并不知道自己的理论分析正确与否,但如果知道概率论中的 大数定律,又碰巧懂一点编程,无疑可以利用计算机重复模拟事件以求解问题。该问题的 Python 3.x 解答程序如下:

from random import *
def once(doors = 3):  # 一次事件的模拟

 car = randrange(doors) # 一扇门后面停着汽车

 man = randrange(doors) # 参赛者预先选择一扇门

 return car == man # 参赛者是否最初就选择到车
h = 0 # 坚持选择赢得汽车的次数                    

c = 0 # 改变选择赢得汽车的次数 

times = int(1e6) # 重复实验的次数
for i in range(times): 

 if once(): h += 1

 else:  c += 1
print("维持选择:",h/times*100,"%\n改变选择:",c/times*100,"%")

运行结果:

维持选择: 33.268 %
改变选择: 66.732 %

扑克牌问题

概率论给我们带来了很多匪夷所思的反常结果,条件概率尤其如此。譬如:

四个人打扑克,其中一个人说,我手上有一个 A。请问他手上有不止一个 A 的概率是多少?
四个人打扑克,其中一个人说,我手上有一个黑桃 A。请问他手上有不止一个 A 的概率又是多少?

from random import *
cards = [i for i in range(52)]

counter = [0, 0, 0, 0]
def once(): # 0 表示黑桃 A

 global cards

 ace = set(sample(cards, 13)) & {0,1,2,3}

 return len(ace), 0 in ace
for i in range(int(1e6)):

 a, s = once() # a 表示 A 的个数, s 表示是否有黑桃 A

 if a:

  counter[1] += 1

  if s: counter[3] += 1

 if a > 1:

  counter[0] += 1

  if s: counter[2] += 1
print('情况一:', counter[0]/counter[1], '\n情况二:', counter[2]/counter[3])

运行结果:

情况一: 0.3694922900321386
情况二: 0.5613778028656186

有趣的事情出来了:如果这个人宣布了手中 A 的花色,他手中持有多个 A 的概率竟然会大大增加。可这又该如何理解呢?

一个家庭中有两个小孩,已知其中一个是女孩,求另一个小孩也是女孩的概率

网络上每一次有人发帖提出与条件概率有关的悖论时,总会引来无数人的围观和争论,哪怕这些问题的实质都是相同的。本题目无疑是争论的最多的问题之一。

说起来网上的分析都像模像样,一些原本都迷糊的人被人讲的晕头转向,一会觉得这个对,一会又觉得那个对。现在我不给你分析那些道理,就用计算机来模拟问题,让你直接得到结论,而毋须明白个中缘由。

from random import * # 0 表示女孩,1 表示男孩
family = (lambda n :[{randrange(2),randrange(2)} for i in range(n)])(int(1e6))
both = family.count({0}) # 都是女孩的家庭数

exist = len(family) - family.count({1}) # 有女孩的家庭数
print(both/exist)

运行结果:
0.33332221770186543

没有那些深奥的分析过程,寥寥数行代码就得到了问题的答案,想必这也是计算机引入数学计算与证明的好处。

生日悖论

每个人都有生日,偶尔会遇到与自己同一天过生日的人,但在生活中这种缘分似乎并不常有。我们猜猜看:在 50 个人当中出现这种缘分的概率有多大,是 10%、20% 还是 50%?

from random import *
counter, times = 0, int(1e6)

for i in range(times):

 if len({randrange(365) for i in range(50)}) != 50: # 存在同一天生日的人

  counter += 1
print('在 50 个人中有相同生日的概率为:',counter/times)

运行结果:
在 50 个人中有相同生日的概率为: 0.970109

在 50 个人中有相同生日的概率高达 97%,这个数字恐怕高出了绝大多数人的意料。我们没有算错,是我们的直觉错了,科学与生活又开了个玩笑。正因为计算结果与日常经验产生了如此明显的矛盾,该问题被称为「生日悖论」,它体现的是理性计算与感性认识的矛盾,并不引起逻辑矛盾,所以倒也算不上严格意义上的悖论。
Python 相关文章推荐
python使用marshal模块序列化实例
Sep 25 Python
Python函数式编程指南(三):迭代器详解
Jun 24 Python
解决python3中cv2读取中文路径的问题
Dec 05 Python
python中栈的原理及实现方法示例
Nov 27 Python
python 列表、字典和集合的添加和删除操作
Dec 16 Python
python使用rsa非对称加密过程解析
Dec 28 Python
tensorflow 实现自定义layer并添加到计算图中
Feb 04 Python
python 递归调用返回None的问题及解决方法
Mar 16 Python
Pytorch 解决自定义子Module .cuda() tensor失败的问题
Jun 23 Python
Django基于Models定制Admin后台实现过程解析
Nov 11 Python
python 实现图片批量压缩的示例
Dec 18 Python
Python提取PDF指定内容并生成新文件
Jun 09 Python
Python中的异常处理学习笔记
Jan 28 #Python
Python中集合类型(set)学习小结
Jan 28 #Python
Python中使用装饰器和元编程实现结构体类实例
Jan 28 #Python
Python实现Const详解
Jan 27 #Python
python映射列表实例分析
Jan 26 #Python
Python操作MySQL简单实现方法
Jan 26 #Python
Python中的is和id用法分析
Jan 26 #Python
You might like
php5.2时间相差8小时
2007/01/15 PHP
PHP批量删除jQuery操作
2017/07/23 PHP
[原创]PHP正则匹配中英文、数字及下划线的方法【用户名验证】
2017/08/01 PHP
javascript 模拟JQuery的Ready方法实现并出现的问题
2009/12/06 Javascript
jquery图片不完全按比例自动缩小的简单代码
2013/07/29 Javascript
jquery submit ie6下失效的原因分析及解决方法
2013/11/15 Javascript
jquery实现不同大小浏览器使用不同的css样式表的方法
2014/04/02 Javascript
原生JS绑定滑轮滚动事件兼容常见浏览器
2014/06/30 Javascript
JavaScript数组去重的3种方法和代码实例
2015/07/01 Javascript
创建自己的jquery表格插件
2015/11/25 Javascript
IE下JS保存图片的简单实例
2016/07/15 Javascript
Angularjs中使用layDate日期控件示例
2017/01/11 Javascript
Vue.js 2.0窥探之Virtual DOM到底是什么?
2017/02/10 Javascript
用户管理的设计_jquery的ajax实现二级联动效果
2017/07/13 jQuery
JavaScript实现时间表动态效果
2017/07/15 Javascript
vue 实现特定条件下绑定事件
2019/11/09 Javascript
详解三种方式在React中解决绑定this的作用域问题并传参
2020/08/18 Javascript
[02:29]DOTA2英雄基础教程 陈
2013/12/17 DOTA
[00:53]2015国际邀请赛 中国区预选赛一触即发
2015/05/14 DOTA
python实现基于信息增益的决策树归纳
2018/12/18 Python
Python button选取本地图片并显示的实例
2019/06/13 Python
原生python实现knn分类算法
2019/10/24 Python
Python Flask上下文管理机制实例解析
2020/03/16 Python
python中threading开启关闭线程操作
2020/05/02 Python
HTML5 文件域+FileReader 分段读取文件并上传到服务器
2017/10/23 HTML / CSS
Avène雅漾美国官方网站:敏感肌肤护理专家
2016/10/24 全球购物
Yahoo-PHP面试题4
2012/05/05 面试题
用友笔试题目
2016/10/25 面试题
投资意向书范本
2014/04/01 职场文书
倡导文明标语
2014/06/16 职场文书
分公司总经理岗位职责
2014/07/30 职场文书
2014年教师党员自我评议
2014/09/19 职场文书
领导干部群众路线教育实践活动剖析材料
2014/10/10 职场文书
2015年基层党支部工作总结
2015/05/21 职场文书
MySql存储过程之逻辑判断和条件控制
2021/05/26 MySQL
PHP中strval()函数实例用法
2021/06/07 PHP