有两个关于随机数几率的问题

stella

1）：是rand(301)得到300的几率大，还是rand(101)+rand(101)+rand(101)得到300的几率大？

2）：是rand(101)+rand(101)+rand(101)得到300的几率大，还是rand(101)*3得到300的几率大？

stella

用来作为随机装备中生成str_plus的一个算法，数值大概在3~400之间，rand(角色等级)*(rand(2)+1)+rand(角色等级)+rand(角色等级)，各位看看还有没有更好玩的算法？

金圭子

金圭子 发表于 2014-11-26 14:31
现在看到公式就头疼，说实话上次玩我之前推荐的那个游戏我突然发现我连傅里叶级数和拉格朗日中值法都忘光 ...

把高数放到游戏里面真是无力吐槽，不过Fourier级数和Lagrange定理还是很重要的

好像就是这个题目：



据说是里面最容易的一题了，不过我还是不知道他在说啥（唉）
其实这个游戏真的是蛮有意思的，配合了那画面和那背景音乐的这种台词，真的是非常让人哈哈哈哈哈哈哈哈……
http://www.66rpg.com/game/6738?s ... 1a374a433f596ad1440

（我会告诉你其实我是在刷积分么）

金圭子

RyanBern 发表于 2014-11-26 14:27
利用中心极限定理来产生服从正态分布的方法比较少见，因为中心极限定理描述的是极限行为，而n=5明显太小 ...

现在看到公式就头疼，说实话上次玩我之前推荐的那个游戏我突然发现我连傅里叶级数和拉格朗日中值法都忘光了，无比惊讶（好吧我当年小学数学好中学物理化学好而到了大学就只剩下计算机码代码了）。
另外n=5是小了点，不过好像图形还是有点正太的……嗯，换个大点的n来跑下……

n=10和n=20的比照

VB 代码复制

Private Sub Command7_Click()
  Me.Cls
  Dim x As Single, y As Single
  Dim i As Integer, j As Integer
  Dim a As Integer, b As Integer
  Dim arr(500) As Integer, arr2(500) As Integer
  Me.DrawWidth = 2
  For i = 1 To 30000
    x = 0
    For j = 1 To 10
      Randomize
      x = x + Rnd * 1000
    Next j
    a = Math.Round(x / 20)
 
    y = 0
    For j = 1 To 20
      Randomize
      y = y + Rnd * 1000
    Next j
    b = Math.Round(y / 40)
    arr(a) = arr(a) + 1
    arr2(b) = arr2(b) + 1
    Me.PSet (a * 10, 5000 - arr(a) * 10)
    Me.PSet (5000 + b * 10, 5000 - arr2(b) * 10)
  Next i
End Sub

Private Sub Command7_Click()


  Me.Cls


  Dim x As Single, y As Single


  Dim i As Integer, j As Integer


  Dim a As Integer, b As Integer


  Dim arr(500) As Integer, arr2(500) As Integer


  Me.DrawWidth = 2


  For i = 1 To 30000


    x = 0


    For j = 1 To 10


      Randomize


      x = x + Rnd * 1000


    Next j


    a = Math.Round(x / 20)


 


    y = 0


    For j = 1 To 20


      Randomize


      y = y + Rnd * 1000


    Next j


    b = Math.Round(y / 40)


    arr(a) = arr(a) + 1


    arr2(b) = arr2(b) + 1


    Me.PSet (a * 10, 5000 - arr(a) * 10)


    Me.PSet (5000 + b * 10, 5000 - arr2(b) * 10)


  Next i


End Sub

RyanBern

金圭子 发表于 2014-11-26 14:08
嗯，做了个模拟正太分布（好吧，翻了下百度发现“指数律”这个名字好像不是人家的学名，人家叫“正态分布” ...

利用中心极限定理来产生服从正态分布的方法比较少见，因为中心极限定理描述的是极限行为，而n=5明显太小，n取太大又不现实。
因此可以参考下面的方法（需要初等概率论的知识）。
这里假设计算机能生成(0,1)之间的随机数，即随机变量X~U(0, 1)，设F(x)为你要模拟的随机变量的分布函数（这里为正态分布的分布函数），那么，Y=F-1(X)就服从标准正态分布。其中F-1(x)为F(x)的广义逆（在正态分布的条件下，这个也是F的严格逆）
由于一维正态分布的分布函数F不容易算出，这里运用2维正态分布来处理一维的情况。

1.事实:如果X是(0, 1)均匀分布，那么Y=-ln(X)服从参数为1的指数分布

2.二维正态分布，相关系数为0，均值为(0, 0)方差为(1 ,1)的正态分布，其密度为exp(-(x^2+y^2)/2)/2pi，极坐标下为rexp(-r^2/2)/2pi，这说明r服从一个指数分布，θ服从(0. 2pi)上的均匀分布，所以

按照上述公式产生的X和Y都满足N(0, 1)，而且相互独立。
其中ξ1和ξ2服从(0, 1)均匀分布，且相互独立。

（科普一下，请勿认真）

金圭子

参考Command2写了个对比的：

VB 代码复制

Private Sub Command6_Click()
  Me.Cls
  Dim x As Single, y As Single
  Dim i As Integer, a As Integer, b As Integer
  Dim arr(500) As Integer, arr2(500) As Integer
  Me.DrawWidth = 2
  For i = 1 To 30000
    Randomize
    x = Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000
    Randomize
    y = (Rnd * 5000 + Rnd * 5000 + Rnd * 5000 + Rnd * 5000 + Rnd * 5000) / 5
    a = Math.Round(x / 10)
    b = Math.Round(y / 10)
    arr(a) = arr(a) + 1
    arr2(b) = arr2(b) + 1
    Me.PSet (a * 10, 5000 - arr(a) * 10)
    Me.PSet (5000 + b * 10, 5000 - arr2(b) * 10)
  Next
End Sub

Private Sub Command6_Click()


  Me.Cls


  Dim x As Single, y As Single


  Dim i As Integer, a As Integer, b As Integer


  Dim arr(500) As Integer, arr2(500) As Integer


  Me.DrawWidth = 2


  For i = 1 To 30000


    Randomize


    x = Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000


    Randomize


    y = (Rnd * 5000 + Rnd * 5000 + Rnd * 5000 + Rnd * 5000 + Rnd * 5000) / 5


    a = Math.Round(x / 10)


    b = Math.Round(y / 10)


    arr(a) = arr(a) + 1


    arr2(b) = arr2(b) + 1


    Me.PSet (a * 10, 5000 - arr(a) * 10)


    Me.PSet (5000 + b * 10, 5000 - arr2(b) * 10)


  Next


End Sub

金圭子

嗯，做了个模拟正太分布（好吧，翻了下百度发现“指数律”这个名字好像不是人家的学名，人家叫“正态分布”，当然前面那个是我故意写错的囧），
跑了下果然还是很符合我预测的：

（我不会隐藏，就直接显示啦）

VB 代码复制

Private Sub Command4_Click()
  Me.Cls
  Dim x As Single
  Dim i As Integer, a As Integer
  Dim arr(500) As Integer
  Me.DrawWidth = 2
  For i = 1 To 30000
    Randomize
    x = Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000
    a = Math.Round(x / 10)
    arr(a) = arr(a) + 1
    Me.PSet (a * 10, 5000 - arr(a) * 10)
  Next
End Sub

Private Sub Command4_Click()


  Me.Cls


  Dim x As Single


  Dim i As Integer, a As Integer


  Dim arr(500) As Integer


  Me.DrawWidth = 2


  For i = 1 To 30000


    Randomize


    x = Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000 + Rnd * 1000


    a = Math.Round(x / 10)


    arr(a) = arr(a) + 1


    Me.PSet (a * 10, 5000 - arr(a) * 10)


  Next


End Sub

（好吧，上面公布了Command2，我这里就改成Command4吧，不过图我懒得换了，自己领悟吧）

再来个Command5：

VB 代码复制

Private Sub Command5_Click()
  Me.Cls
  Dim x As Single
  Dim i As Integer, a As Integer
  Dim arr(500) As Integer
  Me.DrawWidth = 2
  For i = 1 To 30000
    Randomize
    x = Rnd * 5000
    a = Math.Round(x / 10)
    arr(a) = arr(a) + 1
    Me.PSet (a * 10, 5000 - arr(a) * 10)
  Next
End Sub

Private Sub Command5_Click()


  Me.Cls


  Dim x As Single


  Dim i As Integer, a As Integer


  Dim arr(500) As Integer


  Me.DrawWidth = 2


  For i = 1 To 30000


    Randomize


    x = Rnd * 5000


    a = Math.Round(x / 10)


    arr(a) = arr(a) + 1


    Me.PSet (a * 10, 5000 - arr(a) * 10)


  Next


End Sub

金圭子

龙和许也 发表于 2014-11-26 12:41
啊！！！！！我不懂数学！我不懂数学！为啥在水区出现！

你玩过对面橙光的《本座倒要看看你能活几天》么：

http://www.66rpg.com/game/6738?s ... 1a374a433f596ad1440

我觉得是个很不错的游戏，宫斗和“高等”数学混合在一起的“笑果”很不错，第一次玩到皇后问那个问题的时候我简直笑疯了╮(╯▽╰)╭

金圭子

myownroc 发表于 2014-11-26 12:03
不过我还想说的是编程中的随机数是伪随机数，样本数量一大会有规律的……
于是单纯的rand(x)规律比较明显， ...

我擦，第一个也太明显了吧，典型的线性代数啊╮(╯o╰)╭
第二个倒是很符合指数律的。

嗯，拷贝到VB6实测了一下（我一直觉得vb6是个好东西，平时开个解密码玩配对算24点之类的直接用vb6写个程序跑下多快），
第一个点过去简直就是波浪一样的感觉（汗），
第二个倒是有点星云团的感觉（笑）。

很好奇Command2是啥……

先沉下去想想怎么写一个模拟指数律效果东东出来……

龙和许也

啊！！！！！我不懂数学！我不懂数学！为啥在水区出现！{:2_270:}