惊人的答案:平均要取多少个(0,1)中的随机数才能让和超过1
icon2 Brain Storm | icon4 2010-08-08 13:16| icon349 Comments | 本文内容遵从CC版权协议 转载请注明出自matrix67.com

    数学常数最令人着迷的就是,它们常常出现在一些看似与之毫不相干的场合中。 随便取一个 0 到 1 之间的数,再加上另一个 0 到 1 之间的随机数,然后再加上一个 0 到 1 之间的随机数⋯⋯直到和超过 1 为止。一个有趣的问题:平均需要加多少次,才能让和超过 1 呢?答案是 e 次。

  

 
    为了证明这一点,让我们先来看一个更简单的问题:任取两个 0 到 1 之间的实数,它们的和小于 1 的概率有多大?容易想到,满足 x+y<1 的点 (x, y) 占据了正方形 (0, 1)×(0, 1) 的一半面积,因此这两个实数之和小于 1 的概率就是 1/2 。类似地,三个数之和小于 1 的概率则是 1/6 ,它是平面 x+y+z=1 在单位立方体中截得的一个三棱锥。这个 1/6 可以利用截面与底面的相似比关系,通过简单的积分求得:

       ∫(0..1) (x^2)*1/2 dx = 1/6

  

 
    可以想到,四个 0 到 1 之间的随机数之和小于 1 的概率就等于四维立方体一角的“体积”,它的“底面”是一个体积为 1/6 的三维体,在第四维上对其进行积分便可得到其“体积”

       ∫(0..1) (x^3)*1/6 dx = 1/24

    依此类推, n 个随机数之和不超过 1 的概率就是 1/n! ,反过来 n 个数之和大于 1 的概率就是 1 - 1/n! ,因此加到第 n 个数才刚好超过 1 的概率就是

       (1 - 1/n!) - (1 - 1/(n-1)!) = (n-1)/n!

    因此,要想让和超过 1 ,需要累加的期望次数为

       ∑(n=2..∞) n * (n-1)/n! = ∑(n=1..∞) n/n! = e

 
来源:http://www.mostlymaths.net/2010/08/and-e-appears-from-nowhere.html

49 条回复

  • 楼层: 沙发 | | question1663 说:

    BRILLIANT!

  • 楼层: 板凳 | | voices 说:

    佩服!

  • 楼层: 地毯 | | Just :: Solo 说:

    惊人的答案:平均要取多少个(0,1)中的随机数才能让和超过1...

    数学常数最令人着迷的就是,它们常常出现在一些看似与之毫不相干的场合中。 随便取一个 0 到 1 之间的数,再加上另一个 0 到 1 之间的随机数,然后再加上一个 0 到 1 之间的随机数⋯⋯........

  • 楼层: 地板 | | wecing 说:

    我靠,这个答案太悬疑了……

  • 楼层: 地下室 | | Tweets that mention Matrix67: My Blog » Blog Archive » 惊人的答案:平均要取多少个(0,1)中的随机数才能让和超过1 -- Topsy.com 说:

    [...] This post was mentioned on Twitter by lileding, 阳阳猪的 GR Share. 阳阳猪的 GR Share said: 惊人的答案:平均要取多少个(0,1)中的随机数才能让和超过1 http://is.gd/e8vSJ [...]

  • 楼层: 地基 | | Phil 说:

    您真的用 macbook 了?

  • 楼层: 地壳 | | 掌柜的马甲 说:

    e又是这它啊...

  • 楼层: 地幔 | | Kane 说:

    突然发现源地址 http://www.mostlymaths.net 竟然被墙...

  • 楼层: 地核 | | Eagle_Fantasy 说:

    妙哉 无处不在的e

  • 楼层: 10楼 | | 白左 说:

    看见了mac疑似物

  • 楼层: 11楼 | | 白左 说:

    啊呀,所以说这是为什么呢,为什么宇宙中有一个数是大量规律的不变常数呢?
    一个可能的答案是,为了减轻模拟我们这个宇宙的超级计算机的负荷,所有的积分常数和偏移常数都是一样的

    超意识A:喂,给这个基准常数取一个数吧,整数?
    超意识B:iyada,整数很容易让原始人生疑的,取个无理数吧。
    (超意识B在单向加密软件中输入了他的名字,很快,一个长度为无限的无理数MD5i码就生成好了,前三位是2.71。超意识B随手将它复制进变量池最后一个空缺的位置,回车)

    下一秒,宇宙诞生

  • 楼层: 12楼 | | 严酷的魔王 说:

    不错不错~

  • 楼层: 12a楼 | | shenpeng 说:

    刚刚用这个模型去求e的值,两次400000次重复试验试验,结果分别是2.71525和2.7154275,收敛好慢啊……

  • 楼层: 14楼 | | 青年 说:

    数学果然很美~

  • 楼层: 15楼 | | DarkRaven 说:

    另有一个经典题目:
    一串物品,每一个都有一个喜好度,一个一个给你看,你只能选择一次.
    最佳方案是先看前n/e个,选出最大的,然后之后看到一个比这个最大值大的就选.
    额,怎么证明呢?

  • 楼层: 16楼 | | DarkRaven 说:

    @白左:
    无论在什么宇宙,只要数学系统相同,结果就是e,而数学是人类的创造!

  • 楼层: 17楼 | | sexla 说:

    laiguangguangle,实在看不懂

  • 楼层: 18楼 | | biohu 说:

    很神奇...

  • 楼层: 19楼 | | 剑影若兰 说:

    严格来说应该是满足[0,1]均匀分布的随机数吧。。。

  • 楼层: 20楼 | | cbkid 说:

    博主也用mac了?linux没人了吗?

  • 楼层: 21楼 | | ncmooc 说:

    这个答案令人深思。
    惊现mac疑似物。

  • 楼层: 22楼 | | blade 说:

    猜的是3。。。还挺近的。。。一开始想随机取一个期望是0.5,所以大于两个的话基本就行。。。没想到这么复杂:D

  • 楼层: 23楼 | | snaily 说:

    15L的问题也很有趣,求证!

  • 楼层: 24楼 | | snaily 说:

    15L的问题有证明不?

  • 楼层: 25楼 | | 花荣 说:

    再引申一下,这些随机数加在一起,和是多少?
    1.359?

  • 楼层: 26楼 | | 呼吸 说:

    看似毫不相干 我也这么想 幸会幸会

  • 楼层: 27楼 | | Ming 说:

    一定要添加一条:
    随机数是均匀分布在[0,1]上,否则结论不成立

    我很喜欢你对于'n维立方体一角'体积的证明,很简单明了,得到1/n!的结果比我的笨方法好多了.
    但是最后一步有点多余了,因为E(n|S(n)>1)=E(n|S(n)<1)=E(n|S(n)=1).

  • 楼层: 28楼 | | wuzhengkai 说:

    归纳的很精妙!

  • 楼层: 29楼 | | onlytest 说:

    这个是什么软件...

  • 楼层: 30楼 | | Boleyn Su 说:

    首先,将题目转化为平均要取多少个1到n-1中的随机数才能让和超过n。显然当n趋于无穷大时,这两个问题是等效的。
    然后,解:
    设f(x)表示平均取多少次才能大于x.
    显然 f(x)=0 x0来说,设之前最后一次取得了i,那么i在1到n-1之间,且取到任意一个数的概率为1/(n-1).
    所以 f(x)=sum{(1+f(x-i))/(n-1);1<=i<=n-1}
    =1+sum{f(x-i);1<=i<=n-1}/(n-1)
    =1+sum{f(t);0<=t<=x-1}/(n-1)
    令 F(x)=sum{f(i);0<=ioo) f(n)=(1+1/(n-1))^(n-1)*(1+1/(n-1))
    =e*1
    =e

  • 楼层: 31楼 | | Boleyn Su 说:

    无语 XXXX为被和谐内容 这个HTML隔绝方式 也太奇特了吧
    首先,将题目转化为平均要取多少个1到n-1中的随机数才能让和超过n。显然当n趋于无穷大时,这两个问题是等效的。
    然后,解:
    设f(x)表示平均取多少次才能大于x.
    显然 f(x)=0 x<0
    f(x)=1 x=0
    对于x>0来说,设之前最后一次取得了i,那么i在1到n-1之间,且取到任意一个数的概率为1/(n-1).
    所以 f(x)=sum{(1+f(x-i))/(n-1);1<=i<=n-1}
    =1+sum{f(x-i);1<=i<=n-1}/(n-1)
    =1+sum{f(t);0<=t<=x-1}/(n-1)
    令 F(x)=sum{f(i);0<=i<=x}
    则 f(x)=1+F(x-1)/(n-1) ......[1]
    所以 f(x-1)=1+F(x-2)/(n-1) ......[2]
    [1]式-[2]式得 f(x)-f(x-1)=(F(x-1)-F(x-2))/(n-1)
    所以 (f(x)-f(x-1))*(n-1)=f(x-1)
    所以 f(x)=f(x-1)*n/(n-1)
    =(n/(n-1))^x*f(0)
    =(n/(n-1))^x
    所以 f(n)=(n/(n-1))^n
    lim(n->+oo) f(n)=(1+1/(n-1))^(n-1)*(1+1/(n-1))
    =e*1
    =e

  • 楼层: 32楼 | | math 说:

    加到第 n 个数才刚好超过 1 的概率就是

    (1 - 1/n!) - (1 - 1/(n-1)!) = (n-1)/n!
    这步怎么感觉不对啊,减是啥意思

  • 楼层: 33楼 | | sugarzh 说:

    关于15L的问题,魔方俱乐部里面有人讨论过这个
    http://bbs.mf8.com.cn/viewthread.php?tid=22929

    结论在主贴,证明在14L

  • 楼层: 34楼 | | sugarzh 说:

    关于15L的问题,中国魔方俱乐部有人讨论过
    http://bbs.mf8.com.cn/viewthread.php?tid=22929

    结论在14L有证明

  • 楼层: 35楼 | | sugarzh 说:

    - - 不小心风怒了。。。第一次卡了,刷新也刷不出来

  • 楼层: 36楼 | | Nice Guy 说:

    从0到1取一个数小于1的概率也是1吗?!

  • 楼层: 37楼 | | Nice Guy 说:

    n维空间与n-1维空间的运算这样不对吧

  • 楼层: 38楼 | | aLex 说:

    15L本来是选秘书的题吧,大学概率书上有,很经典那,选老婆是不是也适用呀哈哈

  • 楼层: 39楼 | | Matrix67: My Blog » Blog Archive » 趣题:平均要取到第几个随机数才会让序列第一次下降 说:

    [...] 题目来源: Mind Your Decisions 大家有兴趣的话可以看看我之前发过的一个类似的题目,这两个问题似乎是有一腿。 Posted in Brain Storm Tags: 证明, 趣题, [...]

  • 楼层: 40楼 | | paramecium 说:

    昨天刚学了*中心极限定理*手痒想试试其威力的工科男解法:

    对每一个随机数字xi服从均匀分布U(0,1),均值E(xi)=1/2,方差D(xi)=1/12,令X=∑(i=1..n)xi,由中心极限定理,Zn = (X-∑E(xi)) / √∑D(xi)服从标准正态分布N(0,1),P{X>1} = P{Zn>(1-n/2)/√(n/12)} = Φ((1-n/2)/√(n/12))。令(1-n/2)/√(n/12) (7+√33)/2=6.372,即平均*最少*要取6.372个随机数,X=Σxi>1,置信度超过99.87%(=99.74%+0.13%)。

  • 楼层: 41楼 | | paramecium 说:

    发现少了字,再发一遍。

    对每一个随机数字xi服从均匀分布U(0,1),均值E(xi)=1/2,方差D(xi)=1/12,令X=∑(i=1..n)xi,由中心极限定理,Zn = (X-∑E(xi)) / √∑D(xi)服从标准正态分布N(0,1),P{X>1} = P{Zn>(1-n/2)/√(n/12)} = Φ((1-n/2)/√(n/12))。令(1-n/2)/√(n/12)[小于号]-3σ=-3,解得n[大于号](7+√33)/2=6.372,即平均*最少*要取6.372个随机数,X=Σxi>1,置信度超过99.87%(=99.74%+0.13%)。

    “尖括号”里面的字被吃掉了,强大的html过滤...

  • 楼层: Answer to Life, the Universe, and Everything | | paramecium 说:

    还有,31L的问题与原问题并不等价,最多只能说两个问题的答案刚好相等。M牛的问题是(0,1)范围内的实数集(不可数集),而{1/n,2/n..1}是(0,1)范围内的有理数集的一个子集(可数集),不可数集和可数集明显是不能一一对应的,31L的解法把所有的无理数都忽略掉了(比如√2/2)。

  • 楼层: 43楼 | | yh 说:

    15L那个问题若干年前的mm群光棍节模拟赛里有。。

  • 楼层: 44楼 | | Math Oh! Computer » Blog Archive » 随机取数 说:

    [...] 这样就不对了。 最后在matrix67上面看到了这个问题的方法,而且答案也是十分的神奇:。 [...]

  • 楼层: 45楼 | | ZWJ 说:

    如果n个数地和大于一常数K,那n的期望和K是什么函数关系呢?好像是分段的。用体积法可以证明K<1时,n=e^K。1<K<2时,n=e^K-(K-1)*e^(K-1)。后面的维数太高,几何直观出不来了。跪求高人啊!

  • 楼层: 46楼 | | yang_bigarm 说:

    这个解法太高明了,我以前做过这个问题,先求均匀分布之和的表达式f[x](这一步就巨复杂了),然后再算f[x]从0到1积分(又是一堆巨麻烦的积分),可以得到相同的答案。

  • 楼层: 47楼 | | mmaa01 说:

    我曾经考虑过一个问题:
    一个数4,可以分解为4个相同的数之和吗?
    一个数5,可以分解为5个相同的数之和吗?
    一个数6,可以分解为6个相同的数之和吗?
    结果发现都是1.
    同样又想
    一个数4,可以分解为4个相同的数之和吗?
    一个数5,可以分解为5个相同的数之和吗?
    一个数6,可以分解为6个相同的数之和吗?
    答案是可以。但是每个分解都不一样。分别是4的4次根、5的5次根、6的6次根。大小不一样。
    自然联想到哪个数的分解是最大的呢,结果是3的分解最大。
    那么一个实数n可不可以考虑它是n(个)相等的数之乘积呢。就是n的1/n次方。
    发现在实数域里面,通过求导可以得到,最大的那个分解的数是e。
    用上面的语言描述就是e可以分解为e(个)相同的数之积。这个数是所有实数中按这个规律分解得到的最大的一个数1.4446678610097661336583391085964。(用xp自带的计算器算了一下,里面没有e,自己倒腾了几下不知道对不对)

  • 楼层: 48楼 | | orbea jersey 说:

    这个理念不错 我很喜欢!

  • 楼层: 49楼 | | AirLee 说:

    另一个问题:
    定义运算Ranint(1,x),x∈N,x>1,是从1~x中随机取一个整数,设Xn+1=Ranint(1,Xn),X0=m时,若Xn=1,则期望E(n)=Hm≈ln(m),Hm为调和级数.
    拓展:1、运算改为Ran(0,x),x∈R:从(0,x)中随机取一个实数,仿上迭代过程,当Xn<1时结束,E(n)=?
    2、运算改为Ran(0,kx),k>0.k>e时,Xn有增加的趋势,k<e,Xn有减少的趋势,猜想E(Xn/X0)≈(k/e)^n,n趋向于正无穷.

    两个拓展我都没搞出来,我觉得跟这个问题有渊源,求大牛解答.

您也随便说几句吧:

您可以在 Gravatar 设置您的头像。