rand函数的实现原理
rand函数的实现原理
rand函数产生的是伪随机数,也就是说它不是一个真实的随机数。
那么伪随机数是怎么实现的呢?原理大概如下:
如果约定:a1=f(seed),an+1=f(an)那你可以行到一个序列:a1,a2,a3...an,那么要制作一个伪随机函数rand,只需要让它每调用一次就返回序列的下一个元素就行。
就相当于第1次调用rand返回a1,第2次返回a2,…,第n次返回an,这样每次调rand都能拿到一个不同的数,只要整个序列的规律不明显,整个函数看起来就是随机的。
现在计算机上的rand函数都是用这样的原理实现的,这里的seed被称为“随机数种子”。
但这里有一个问题,如果seed不变,那我们每次调用rand函数获取的序列都是相同的。这就会造成有的程序跑一遍退出后,再重新跑一遍,两次的输出结果是相同的。所以我们还需要一个接口去设置seed值,这个接口就是srand函数。
linux下的rand是用类似下面的代码实现的:
static unsigned long next = 1; /* RAND_MAX assumed to be 32767 */ int myrand(void) { next = next * 1103515245 + 12345; return((unsigned)(next/65536) % 32768); } void mysrand(unsigned seed) { next = seed; }
myrand、mysrand分别对应rand和srand,但实际的rand实现会复杂一些。
使用rand函数的一个问题
有些人使用rand函数时,因为初始化时没有调用srand会造成程序每次运行输出结果都相同:http://ask.csdn.net/questions/174871
如果你没有调用srand设置随机数种子,seed的默认值会是0,而seed为0时所决定的序列是固定的,而第一次调用rand()就是返回这个固定序列里的第1个元素,那它的值也是固定的,自然你的程序每次输出都一样了。
所以正确的写法应该是程序初始化时用srand设置不同的随机数种子(只需要设置一次),例如srand(time(NULL)),但要注意,time(NULL)的值是隔1秒才改变一次的,必要情况下可以考虑使用精度更高的时间函数,如gettimeofday。
下面这段程序,只要你不是在同一秒内执行两次,每次输出结果都是不一样的:
#include <stdio.h> #include <stdlib.h> #include <time.h> int main() { srand(time(NULL)); // 设置随机数种子 for (int i = 0; i < 10; i++) { printf("%u\n", rand()); } getchar(); return 0; }