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 多线程应用介绍
Dec 19 Python
python装饰器初探(推荐)
Jul 21 Python
详解python 发送邮件实例代码
Dec 22 Python
python django使用haystack:全文检索的框架(实例讲解)
Sep 27 Python
python数据封装json格式数据
Mar 04 Python
对Python中for复合语句的使用示例讲解
Nov 01 Python
Python实现KNN(K-近邻)算法的示例代码
Mar 05 Python
python pytest进阶之fixture详解
Jun 27 Python
Python对接六大主流数据库(只需三步)
Jul 31 Python
django-crontab 定时执行任务方法的实现
Sep 06 Python
Python代码生成视频的缩略图的实例讲解
Dec 22 Python
对tensorflow 中tile函数的使用详解
Feb 07 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调用三种数据库的方法(3)
2006/10/09 PHP
php cookie名使用点号(句号)会被转换
2014/10/23 PHP
php实现将上传word文件转为html的方法
2015/06/03 PHP
PHP中Trait及其应用详解
2017/02/14 PHP
phpStudy 2016 使用教程详解(支持PHP7)
2017/10/18 PHP
PHP正则判断一个变量是否为正整数的方法
2019/02/27 PHP
理解Javascript_08_函数对象
2010/10/15 Javascript
jquery的index方法实现tab效果
2011/02/16 Javascript
Json序列化和反序列化方法解析
2013/12/19 Javascript
JavaScript包装对象使用详解
2015/07/09 Javascript
Javascript中replace()小结
2015/09/30 Javascript
分享几种比较简单实用的JavaScript tabel切换
2015/12/31 Javascript
原生JS简单实现ajax的方法示例
2016/11/29 Javascript
探究JavaScript中的五种事件处理程序方式
2016/12/07 Javascript
利用JS实现简单的日期选择插件
2017/01/23 Javascript
jQuery实现一个简单的轮播图
2017/02/19 Javascript
vue 使用post/get 下载导出文件操作
2020/08/07 Javascript
echarts实现晶体球面投影的实例教程
2020/10/10 Javascript
python简单的函数定义和用法实例
2015/05/07 Python
在Django中使用Sitemap的方法讲解
2015/07/22 Python
详解Python中映射类型的内建函数和工厂函数
2015/08/19 Python
理解Python垃圾回收机制
2016/02/12 Python
详解python使用Nginx和uWSGI来运行Python应用
2018/01/09 Python
修复CentOS7升级Python到3.6版本后yum不能正确使用的解决方法
2018/01/26 Python
Pytorch中Tensor与各种图像格式的相互转化详解
2019/12/26 Python
Python3爬虫发送请求的知识点实例
2020/07/30 Python
python 实时调取摄像头的示例代码
2020/11/25 Python
CSS中的字体大小设置属性总结
2016/05/24 HTML / CSS
巴黎一票通:The Paris Pass
2018/02/10 全球购物
2014年文员工作总结
2014/11/18 职场文书
匿名信格式范文
2015/05/27 职场文书
董事长助理工作总结2015
2015/07/23 职场文书
2016年小学端午节活动总结
2016/04/01 职场文书
HTML+CSS制作心跳特效的实现
2021/05/26 HTML / CSS
深入解析Apache Hudi内核文件标记机制
2022/03/31 Servers
利用Python脚本写端口扫描器socket,python-nmap
2022/07/23 Python