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中各种方法的运作原理
Jun 15 Python
Python制作钉钉加密/解密工具
Dec 07 Python
基于python select.select模块通信的实例讲解
Sep 21 Python
python2 与python3的print区别小结
Jan 16 Python
Python实现简易版的Web服务器(推荐)
Jan 29 Python
PyQt5实现下载进度条效果
Apr 19 Python
Python实现读取SQLServer数据并插入到MongoDB数据库的方法示例
Jun 09 Python
解决django 新增加用户信息出现错误的问题
Jul 28 Python
Python帮你识破双11的套路
Nov 11 Python
Django继承自带user表并重写的例子
Nov 18 Python
python 实现端口扫描工具
Dec 18 Python
详解Python爬虫爬取博客园问题列表所有的问题
Jan 18 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
漫威DC御用漫画家去世 他的表情包曾走红网络
2020/04/09 欧美动漫
现磨咖啡骗局!现磨咖啡=新鲜咖啡?现磨咖啡背后的猫腻你不懂!
2019/03/28 冲泡冲煮
用PHP调用数据库的存贮过程
2006/10/09 PHP
PHP中如何防止外部恶意提交调用ajax接口
2016/04/11 PHP
PHP实现ASCII码与字符串相互转换的方法
2017/04/29 PHP
laravel框架邮箱认证实现方法详解
2019/11/22 PHP
基于jquery的lazy loader插件实现图片的延迟加载[简单使用]
2011/05/07 Javascript
前端必备神器 Snap.svg 弹动效果
2014/11/10 Javascript
jQuery使用drag效果实现自由拖拽div
2015/06/11 Javascript
理解Javascript的动态语言特性
2015/06/17 Javascript
js 动态给元素添加、移除事件的实现方法
2016/07/19 Javascript
JS去掉字符串前后空格或去掉所有空格的用法
2017/03/25 Javascript
Angular 4.x 动态创建表单实例
2017/04/25 Javascript
js评分组件使用详解
2017/06/06 Javascript
React Native react-navigation 导航使用详解
2017/12/01 Javascript
Nuxt.js之自动路由原理的实现方法
2018/11/21 Javascript
Vue 2.0双向绑定原理的实现方法
2019/10/23 Javascript
jquery使用echarts实现有向图可视化功能示例
2019/11/25 jQuery
如何利用node转发请求详解
2020/09/17 Javascript
[01:01:29]2018DOTA2亚洲邀请赛 4.4 淘汰赛 VP vs Liquid 第一场
2018/04/05 DOTA
python简单实现基数排序算法
2015/05/16 Python
scrapy爬虫完整实例
2018/01/25 Python
在pycharm 中添加运行参数的操作方法
2019/01/19 Python
Python 类方法和实例方法(@classmethod),静态方法(@staticmethod)原理与用法分析
2019/09/20 Python
浅谈css3中的前缀
2016/07/20 HTML / CSS
Stuart Weitzman欧盟:美国奢华鞋履品牌
2017/05/24 全球购物
工地资料员岗位职责
2013/12/31 职场文书
机关保密承诺书
2014/06/03 职场文书
党的群众路线教育实践活动个人对照检查剖析材料
2014/09/23 职场文书
公司费用报销管理制度
2015/08/04 职场文书
2016暑期社会实践新闻稿
2015/11/25 职场文书
2016教师节感恩话语
2015/12/09 职场文书
小学语文新课改心得体会
2016/01/22 职场文书
2019财务毕业实习报告
2019/06/27 职场文书
NGINX 权限控制文件预览和下载的实现原理
2022/01/18 Servers
win10系统xps文件怎么打开?win10打开xps文件的两种操作方法
2022/07/23 数码科技