随机数生成算法

看到一些介绍随机数生成算法的文章，收集下来，有空深入了解下。

http://bufferpool.blogbus.com/logs/10182011.html

C语言中伪随机数生成算法实际上是采用了"线性同余法”。具体的计算如下：

Xi = (Xi-1 * A + C ) mod M 
其中A,C,M都是常数（一般会取质数）。当C=0时，叫做乘同余法。srand()函数置的seed实际上会作为X0被代入上式中，然后每次调用rand()函数都会用上一次产生的随机值来生成新的随机值。这样我们可以看出实际上用rand()函数生成的是一个递推的序列，一切值都来源于最初的seed。这就不难解释为什么当初始的seed取一样的时候，得到的序列都相同。也就解释了为什么我们需要调用系统时间这种变化的值来充当seed了。

 

比较有意思的是在CSDN上搜到了sjd163一篇文章讲到了一种通过交换来得到随机序列的方法。主要的依据是热力学第二定律的熵增加原理。主要的思想是将有序的序列通过一定的交换将有序变成无序。思想类似于扑克牌的洗牌。假设扑克牌是从小到大的顺序放好的，通过多次的洗牌之后扑克牌的顺序性就逐渐消失，变的随机。换句话说，系统的熵值变大。而且这种过程是不可逆的，只能从有序变成无序，不能从无序变成有序。

假设现在有一个数组，从小到大存放了一定数字。我们按照对插、倒排、交换、块分割等等方法对这个数组进行多轮的处理。原先包含在数组中的顺序关系就被我们逐渐打乱了。当然，为了使过程不可逆，我们可以在这些过程（比如交换）中引入简单伪随机值。通过n轮的处理之后，我们可以假设现在数组已经变成了随机序列。可以对结果进行概率分布的检验，如果检验不通过，重新进行处理，直到到达可以认为是随机序列的程度。

 

另外一篇是来自开源中国的，一个西安交大的学生，写了一个开源信号处理库（C++），实现了好几种随机算法

http://my.oschina.net/zmjerry/blog/3723