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目录操作之python遍历文件夹后将结果存储为xml
Jan 27 Python
跟老齐学Python之集合的关系
Sep 24 Python
Python判断文件和文件夹是否存在的方法
May 21 Python
python线程池(threadpool)模块使用笔记详解
Nov 17 Python
Python实现的txt文件去重功能示例
Jul 07 Python
Python使用matplotlib 模块scatter方法画散点图示例
Sep 27 Python
下载官网python并安装的步骤详解
Oct 12 Python
python元组和字典的内建函数实例详解
Oct 22 Python
Python算法的时间复杂度和空间复杂度(实例解析)
Nov 19 Python
Django异步任务线程池实现原理
Dec 17 Python
python 解决tqdm模块不能单行显示的问题
Feb 19 Python
在CentOS7下安装Python3教程解析
Jul 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
PHP语法速查表
2007/01/02 PHP
在Laravel框架里实现发送邮件实例(邮箱验证)
2016/05/20 PHP
php pdo操作数据库示例
2017/03/10 PHP
解决FireFox下[使用event很麻烦]的问题
2006/11/26 Javascript
Firefox getBoxObjectFor getBoundingClientRect联系
2008/10/26 Javascript
JavaScript 精粹读书笔记(1,2)
2010/02/07 Javascript
Jquery Ajax的Get方式时需要注意URL地方
2011/04/07 Javascript
js自动闭合html标签(自动补全html标记)
2012/10/04 Javascript
javascript的字符串按引用复制和传递,按值来比较介绍与应用
2012/12/28 Javascript
javascript数组去重3种方法的性能测试与比较
2013/03/26 Javascript
解析dom中的children对象数组元素firstChild,lastChild的使用
2013/07/10 Javascript
如何书写高质量jQuery代码(使用jquery性能问题)
2014/06/30 Javascript
JS使用ajax方法获取指定url的head信息中指定字段值的方法
2015/03/24 Javascript
JavaScript实现横向滑出的多级菜单效果
2015/10/09 Javascript
js html5 css俄罗斯方块游戏再现
2016/10/17 Javascript
js编写简单的计时器功能
2017/07/15 Javascript
微信小程序动态增加按钮组件
2018/09/14 Javascript
JS中使用react-tooltip插件实现鼠标悬浮显示框
2019/05/15 Javascript
layui扩展上传组件模拟进度条的方法
2019/09/23 Javascript
layer ui 导入文件之前传入数据的实例
2019/09/23 Javascript
Vue computed 计算属性代码实例
2020/04/22 Javascript
Python程序员鲜为人知但你应该知道的17个问题
2014/06/04 Python
Python遍历目录的4种方法实例介绍
2015/04/13 Python
微信小程序跳一跳游戏 python脚本跳一跳刷高分技巧
2018/01/04 Python
python实现雪花飘落效果实例讲解
2019/06/18 Python
python自动结束mysql慢查询会话的实例代码
2019/10/27 Python
查看端口并杀进程python脚本代码
2019/12/17 Python
基于python使用tibco ems代码实例
2019/12/20 Python
Python压缩模块zipfile实现原理及用法解析
2020/08/14 Python
Myprotein加拿大官网:欧洲第一的运动营养品牌
2018/01/06 全球购物
拉斯维加斯酒店、演出、旅游、俱乐部及更多:Vegas.com
2019/02/28 全球购物
新员工培训个人的自我评价
2013/10/09 职场文书
政治思想表现评语
2014/05/04 职场文书
学习十八届四中全会依法治国心得体会
2014/11/03 职场文书
捐款感谢信
2015/01/20 职场文书
2015年财务经理工作总结
2015/05/13 职场文书