生成不重复随机数的一段代码

小蚂蚁树

学习了这里的很多内容。很多高手。。。
灰袍法师的等概率遍历法应该是直接符合概率的定义的，剩余多少值域，剩余多少要取的值，去匹配rnd的值。不过，从实现来说，可能会受到rnd精度的影响？

如果VBA里面rnd的最大有效数字15位，就是≥10^-15，要在1～10的20次方的值域里面取10个（初始概率是10^-19），那么前面的值段就无法满足了。只能等到后来（第10^16个数以后？），被除数越来越小，商越来越大，才可能大于rnd的值。随机性受到精度的影响。

我用excel里面的RAND()公式，几万个值的最小值也只是10^-5级别的。VBA可能也差不多？那么精度就差得更多了。

另外，洗牌法，其中随机位置的范围应该是递减的。开头递增，末尾不变。

另外，检验是否随机，我觉得一个简单的必要条件是（不是充分条件），生成的数的平均值是否接近理论平均值（不管是允许重复还是不允许重复），就是值域首尾两个数的算术平均值。

小蚂蚁树

学习了灰袍法师的另外一个随机步长的方法。思路很巧妙。         
不过，这种方法好像有一个缺陷，就是，在取的数比较多，即“空档”比较少的情况下，会倾向于出前面的值。影响了随机性。

小蚂蚁树

比如在1到1000的范围内取900个随机数，理论上的平均值应该是500.5，而用这种方法得到的平均值一直是467左右。而1到500每次都是出现的。

我觉得是这段代码的原因：

1. For i = 1 To split_place - 1
   
2.     my_share = Round(arr(i) / sum * total_share, 0)
   
3.     total_share = total_share - my_share
   
4.     sum = sum - arr(i)
   
5.     arr(i) = arr(i - 1) + 1 + my_share
   
6.     loop_count = loop_count + 1
   
7. Next i
   

复制代码

这句：Round(arr(i) / sum * total_share, 0)，其中total_share是100，sum（假如不放大的话）≈450，而arr(i)（假如不放大的话）是在0到1之间。那么，只有当sum减小到<200的时候，才可能使得(arr(i) / sum * total_share)>0.5，即 my_share=1，否则 my_share 一直都是0。
sum从450减少到200，因为arr(i)平均值是0.5左右，所以 i 要到500以后，才可能使得 sum 变成200。  450-500*0.5=200

我添加了测试代码：

1. Dim SumEnough As Boolean '本句添加
   
2. For i = 1 To split_place - 1
   
3.     my_share = Round(arr(i) / sum * total_share, 0)
   
4.   If my_share > 0 And SumEnough = False Then '本句添加
   
5.     MsgBox i & "," & arr(i) & "," & sum & "," & total_share & "," & my_share
   
6.     SumEnough = True
   
7.   End If
   
8.     total_share = total_share - my_share
   
9.     sum = sum - arr(i)
   
10.     arr(i) = arr(i - 1) + 1 + my_share
    
11.     loop_count = loop_count + 1
    
12. Next i
    

复制代码

然后运行发现，每次都是 i 到500以上，才第一次出现 my_share >0。也就是说，此前500多个数是必然出现的。因此导致了随机性受到破坏。
当然，要按比例对应到“空档”，“空档”是整数，而比例是精度有限、“密度”不定的小数，两者之间有矛盾，可能无法避免这样的问题。

上善若水11

请教：如果生成的数据要求符合平均值±标准差，如何调整代码？

plu2

向各位大神学习！

bdggch

需要好好看，慢慢消化。

hsyn01

谢谢你的分享，值得学习学习。

EXCEL20022002

灰袍法师 发表于 2011-1-24 16:23
洗牌算法（记住：不叫跳蚤算法）的弱点是占用 max-min大小的内存，
max-min＝1亿的时候当然会内存不足

更好的附件在27楼和29楼,,请问同学, 29的附近下载不了,方便重新上传下吗?

EXCEL20022002

小蚂蚁树 发表于 2016-6-22 21:43
比如在1到1000的范围内取900个随机数，理论上的平均值应该是500.5，而用这种方法得到的平均值一直是467左右 ...

同学你好,  29楼的法师的代码 怎么下载不了,谢谢如果你有文件,谢谢上传下

EXCEL20022002

同学你好,  29楼的法师的代码 怎么下载不了,谢谢如果你有文件,谢谢上传下