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实现socket客户端和服务端简单示例
Feb 24 Python
Python实现批量把SVG格式转成png、pdf格式的代码分享
Aug 21 Python
简单谈谈python中的多进程
Nov 06 Python
Python中模块string.py详解
Mar 12 Python
详解用python实现简单的遗传算法
Jan 02 Python
Python爬虫获取图片并下载保存至本地的实例
Jun 01 Python
python读写配置文件操作示例
Jul 03 Python
python递归法实现简易连连看小游戏
Mar 25 Python
使用python快速实现不同机器间文件夹共享方式
Dec 22 Python
Python中zip()函数的解释和可视化(实例详解)
Feb 16 Python
Python中BeautifulSoup通过查找Id获取元素信息
Dec 07 Python
Python使用tkinter实现小时钟效果
Feb 22 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的范围解析操作符(::)的含义分析说明
2011/07/03 PHP
什么情况下可以不写PHP的闭合标签“?>”
2014/08/28 PHP
PHP后端银联支付及退款实例代码
2017/06/23 PHP
php更新cookie内容的详细方法
2019/09/30 PHP
PHP快速导出百万级数据到CSV或者EXCEL文件
2020/11/27 PHP
HR vs ForZe BO3 第一场 2.13
2021/03/10 DOTA
javascript实现左右控制无缝滚动
2014/12/31 Javascript
JavaScript使用Prototype实现面向对象的方法
2015/04/14 Javascript
angularjs自定义ng-model标签的属性
2016/01/21 Javascript
BootStrap中关于Select下拉框选择触发事件及扩展
2016/11/22 Javascript
thinkjs之页面跳转同步异步操作
2017/02/05 Javascript
简单的JS控制button颜色随点击更改的实现方法
2017/04/17 Javascript
React学习笔记之条件渲染(一)
2017/07/02 Javascript
vue-cli配置文件——config篇
2018/01/04 Javascript
Node.js中的cluster模块深入解读
2018/06/11 Javascript
vue2.0的虚拟DOM渲染思路分析
2018/08/09 Javascript
React 源码中的依赖注入方法
2018/11/07 Javascript
ckeditor一键排版功能实现方法分析
2020/02/06 Javascript
Vue + ts实现轮播插件的示例
2020/11/10 Javascript
介绍Python中的fabs()方法的使用
2015/05/14 Python
Python 中 -m 的典型用法、原理解析与发展演变
2019/11/11 Python
Python3 使用selenium插件爬取苏宁商家联系电话
2019/12/23 Python
python使用QQ邮箱实现自动发送邮件
2020/06/22 Python
Python代码注释规范代码实例解析
2020/08/14 Python
详解python模块pychartdir安装及导入问题
2020/10/22 Python
python 自定义异常和主动抛出异常(raise)的操作
2020/12/11 Python
Shopty西班牙:缝纫机在线销售
2018/01/26 全球购物
华为慧通笔试题
2016/04/22 面试题
通用自荐信范文
2014/03/14 职场文书
慈善晚会策划方案
2014/05/14 职场文书
宣传口号大全
2014/06/16 职场文书
亲子阅读的活动方案
2014/08/15 职场文书
领导干部群众路线对照检查材料
2014/11/05 职场文书
2015年劳动部工作总结
2015/05/23 职场文书
使用Python的开发框架Brownie部署以太坊智能合约
2021/05/28 Python
nginx服务器的下载安装与使用详解
2021/08/02 Servers