C++学习笔记-随机数

随机函数Rand SRand

原理

计算机常用线性同余法产生随机数，LCG( Linear Congruential Generator)用递推的方法求随机序列

$$X(n+1) = ( a * x(n) + c ) \% M$$

其中：$0<=a < M, 0<=c< M, 0<= x(0)< M.$ 为了计算迅速，M一般取2的幂，求余可用移位实现。
为了尽可能让序列周期长，最大为M，a和c一般互斥，与M也互斥。

周期为M的充分必要条件：

$$\begin{equation} \begin{split} T(k)&=(x*a^k）+(a^(k-1)*c)+(a^(k-2)*c)+...+a*c+c \\ T(k) \% M &= x (k=M) \\ T(k)\%M &\neq x (0<k<M)\\ \end{split} \end{equation}$$

只要证明x=0~M-1中任意一个上式成立即可。很容易推断只要有一个x周期为M，那么其他所有x周期也为M，序列周期为M。
由x=0和x=1容易推出充分条件：$a^M \% M = 1$. a和M必然互斥(反正法).
必要条件：TODO

C++函数

C++标准库里有两个随机函数，一个设置随机种子 x(0)，一个产生伪随机数。

int Rand(void)

产生随机函数，使用的seed是上一次调用SRand时指定的seed

void srand (unsigned int seed)

Initialize random number generator
The pseudo-random number generator is initialized using the argument passed as seed.

For every different seed value used in a call to srand, the pseudo-random number generator can be expected to generate a different succession of results in the subsequent calls to rand.

Two different initializations with the same seed will generate the same succession of results in subsequent calls to rand.

If seed is set to 1, the generator is reinitialized to its initial value and produces the same values as before any call to rand or srand.

In order to generate random-like numbers, srand is usually initialized to some distinctive runtime value, like the value returned by function time (declared in header ). This is distinctive enough for most trivial randomization needs.

Ref：http://www.cplusplus.com/reference/cstdlib/srand/

使用的时候先指定随机种子srand(time(NULL))，再调用rand().rand()里会将这一次得到的随机数作为下一次的种子。
程序每次启动，计算机自己初始化的种子相同，得到的随机序列也相同。主要目的为了方便调试。因此想得到真正的随机序列，自己初始化一下种子。