近几个月,网络上最火的就是锦鲤女孩,是国庆假期某支付平台官微搞了个带有锦鲤图的转发抽奖活动,一名网名叫@信小呆的呆萌小姐姐抽中了这个巨额大奖,小姐姐也因此又被广大锦鲤爱好者们当成了“中国活体锦鲤”,疯狂转发。

学习Python中有不明白推荐加入交流裙
                号:735934841
                群里有志同道合的小伙伴,互帮互助,
                群里有免费的视频学习教程和PDF!

自此之后,朋友圈,以及各大媒体论坛盛行各种“锦鲤”。作为程序员的我们,自然应该想一想,如何完成“锦鲤”的抽取?

实际这就是一个抽样问题。可以抽象为如下问题:

  • 随机的选取容量为N的数组中的k个元素,要求是不能重复选取,并且不能删除数组中的元素,只能够进行交换。其中 k≤N 。

看到这个问题,你想到了什么?

我先想到的就是抽签算法。当由 k 个人抽 K 张签,无论先后顺序,每个人抽中的概率都是1/K 。同理, k 个人抽 N 张签,无论先后顺序,每个人抽中的概率都是1/N 。

可以简单说明一下:

1、当k=1时,由于是从 N 张签中抽取,所以抽中的概率是1/N,成立。

2、当k=2时,在剩下的 N-1 个中随机选:1/(N-1),由于第1次没有选中它, 而是在另外N-1个中选:(N-1)/N,因此概率为:(N-1)/N * 1/(N-1) = 1/N。

3、当k=3时,概率为 (N-1)/N * (N-2)/(N-1) * 1/(N-2) = 1/N。

4、重复上述流程,到k=N。

既然如此,我们可以对N个数,进行k次抽签操作,算法代码如下:

import random
def SelectRandomK(L, N, k):for i in range(0, k):# 产生i到N-1间的随机数j = random.randint(i, N - 1)L[i], L[j] = L[j], L[i]

这个算法实现的缺点就是依赖了数据总数N。如果不知道N有没有办法呢?

那就是蓄水池算法。蓄水池算法是大数据抽样常用的算法,在不知道数据总数目的情况下可以完成随机抽样。

先从最简单的蓄水池抽样算法说起,即蓄水池中数据的数目为1。先把第一个数据以概率1/1放入蓄水池,第二个数据以1/2的概率替换蓄水池中数据,第三个数据以1/3的概率替换蓄水池中数据,第k个数据以1/k的概率替换蓄水池中数据,如此重复,直至遍历完所有的数据。

这样完成后,每个数据被抽中的概率是相等的,即使不知道数据总数目。下面,就用数学归纳法完成证明:

只需要证明当遍历至第k行时,前面k行中的任意一行被抽取的概率均为1/k。

证明:(1)当i=1时,第一行被抽取的概率为1/1,成立。

(2)假设当i=k时成立,则前面k行中的任意一行被抽取的概率为1/k。现证明i=k+1时成立。

当i=k+1时,第k+1行替换原有样本的概率为1/(k+1),所以第k+1行被抽取的概率是1/(k+1)。

前面k行任意一行被抽取的概率为 1/k*k/(k+1)=1/(k+1),

即当i=k+1时成立。证毕。

Python代码实现如下:

def SelectRandom1(L):# 计数器num = 2for item in L[1:]:if random.random() < 1.0/num:L[0], L[num - 1] = L[num - 1], L[0]num += 1

那么,如何以等概率选择k个数呢?跟单个数类似,实现方法如下:

先把前k个数据放入蓄水池,对第k+1个数据,我们以 k/(k+1)概率决定是否要把它放入蓄水池,放入时随机的选取一个作为替换项。对第n个数据,我们以 k/n概率决定是否要把它放入蓄水池,放入时随机的选取一个作为替换项。这样一直做下去,直至遍历完所有的样本空间。可以证明,对于任意的样本空间N,对每个数据的选取概率都为k/N。

也可以通过数学归纳法完成证明:

只需要证明当遍历至第n(n>=k)行时,前面n行中的任意一行被抽取的概率均为k/n。

证明:(1)当i=k时,前面k行被抽取的概率为k/k=1,成立。

(2)假设当i=n时成立,则前面n行中的任意一行被抽取的概率为k/n。现证明i=n+1时成立。

当i=n+1时,第n+1行替换原有样本的概率为k/(n+1),所以第n+1行被抽取的概率是k/(n+1)。

前面n行任意一行被抽取的概率为

k/n*(k/(n+1)*(k-1)/k+(n+1-k)/(n+1))=k/n*(n/(n+1))=k/(n+1)

即当i=n+1时成立。证毕。

Python代码实现如下:

def SelectRandomK(L, k):# 计数器num = k + 1for item in L[k:]:if random.random() < float(k/num):# 产生0到k-1之间的随机数j = random.randint(0, k - 1)L[num - 1], L[j] = L[j], L[num - 1]num += 1

关于随机抽样算法,你深入理解了吗? 你明白朋友圈“锦鲤”该如何抽取了吗?

朋友圈“锦鲤”盛行,如何抽取“锦鲤”?Python放大招了~相关推荐

  1. python朋友圈自动点赞_基于AirTest+Python的ios自动化测试demo(微信朋友圈无限点赞)...

    AirTest相比Appuim有个好处就是可以对GUI图片进行捕捉和最新版本支持WebView(目前Appuim不支持iOS12的WebView进行Xpath抓取) AirTest环境搭建可参考以下链 ...

  2. python微信朋友圈刷图_10分钟用Python做个微信朋友圈抽奖九宫格

    最近在朋友圈看到个好玩的抽奖九宫格: 随便点开一个: 哈哈,笑出猪叫!这个玩法利用了朋友圈发长图会显示中间局部的设定,搞笑之余也为朋友圈广告营销贡献了新思路. 制作图的过程也不复杂,长图对应部分设计好 ...

  3. 朋友圈如何秀一把!用Python一键生成炫酷九宫格图片!

    作为一个男同胞来说,为了给女朋友拍一张美美的照片,着实需要花费很大的时间和精力,不仅仅需要从众多的图片中精心挑选,而且还需要有着超强的图片精修能力,才能得到一张张达到女友要求的图片,真心不容易啊- 朋 ...

  4. python发朋友圈_10分钟教你用Python发一个高逼格的朋友圈

    程序猿声 你与千万程序猿在一起 01 前言 Hello~各位小伙伴们大家好.现在大家是越来越离不开手机,离不开微信了.每天打开手机的第一或者第二件事就是赶紧打开朋友圈看看有什么好玩的东西.偶尔忍不住了 ...

  5. python训练营 朋友圈_教你如何用python来弄一个高逼格的朋友圈

    小编今天要给大家介绍一个Python库: PIL(Python Image Library) 下面我们用一个实际的例子 看看50行python代码可以做什么神奇的事情 这是小编发的一个朋友圈 切图前是 ...

  6. python朋友圈表白_备战520|Python花式表白的几种姿势

    大家好,520它又要来了 所以今天的主题是粉色的 为了各位秃头程序员不再头疼 本文给大家介绍几种用Python表白的姿势 绝不是画个爱心曲线那么简单- . ? 属于TA的词云 用Python将你们的  ...

  7. 朋友圈装死,微博蹦迪|Python教你如何掌握女神情绪变化

    很多人都是在朋友圈装死,微博上蹦迪. 微信朋友圈已经不是一个可以随意发表心情的地方了,微博才是! 所以你不要傻傻盯着女神的朋友圈发呆啦! 本文教你如何用Python自动通知女神微博情绪变化,从今天开始 ...

  8. 微信朋友圈的实现(图片查看,手势放大、缩小)

    效果图: 简单的就这样,实现效果的方法有多种,可以使用ListView嵌套GridView,也可以直接使用RecycleView,这里只说RecycleView的实现方法. 难一点的地方就是每个ite ...

  9. python是个什么东西刷爆朋友圈_五星好评!这门Python神作刷爆朋友圈

    原标题:五星好评!这门Python神作刷爆朋友圈 2019年,有人说是Python元年,背靠AI和大数据,连续半年拿下TOP1编程语言,这些都给了Python的底气,所以无论是薪资还是招聘需求,都越来 ...

最新文章

  1. redis Could not connect to Redis at 127.0.0.1:6379: Connection refused 问题解决
  2. Unity Standard Assets 简介之 Cameras
  3. mysql sql诊断建议_MySQL诊断调优常用SQL语句
  4. java实现图形界面输入半径求圆面积_测试开发工程师系列之Android自动化测试Appium(Python)连载(7)安卓图形界面...
  5. JavaFX 2.0和Scala,例如牛奶和饼干
  6. 解决win2003安装exchangeServer后关机慢的方法
  7. STM32 学习总结2 ----利用中断来控制按键点灯、捕获功能练习
  8. 应用Composite模式
  9. 视频转gif怎样操作?如何快速在线生成gif动图?
  10. mysql数据库基础知识--一分钟让你数据库入门(sql基础语法篇)
  11. Java 内存模型(Java Memory Model)
  12. 在JavaScript中实现继承的几种方式
  13. 软件工程——软件维护
  14. flex布局 justify-content 新属性值 space-evenly
  15. Unity LookAt 2D锁定
  16. BuildPack:无需编写 Dockerfile,新一代的企业镜像打包工具
  17. 实例:用C#.NET手把手教你做微信公众号开发(8)--普通消息处理之链接(普通消息终结篇)
  18. python实现绘制ROC曲线寻找指标最佳诊断临界值
  19. 连接mysql工具_非常实用的mysql可视化连接工具,开发测试必备软件,快收藏吧...
  20. linux怎么查看内容并显示行号,linux中查看文件时显示行号

热门文章

  1. c51单片机矩阵键盘1602计算器_【Proteus】单片机配合矩阵键盘LCD1602制作简易计算器...
  2. 将Android替换为Linux,android:在AndroidO下将selinux设置为permissive
  3. qq离开状态经典语句
  4. Docker Desktop卸载后Hyper-V服务驻留与第三方虚拟机产生冲突
  5. JAVA语言中的this关键字用法
  6. 编译mavros报错提示:Your workspace is not extending any other result space, but it is set to use a `linked`
  7. laravel map方法踩坑记录
  8. win10共享打印机(连接共享打印机的两种方法,IP地址变化的连接方法)
  9. nginx重定向导致端口消失无法正常访问
  10. Geomagic Touch(USB版本)在ROS下的配置和使用