rand函数的使用及实现原理

最近很想知道c语言中random是怎么实现的，查了一些资料，总结了一下：

参考：http://topic.csdn.net/u/20091124/14/2598eca6-b66a-47d6-9603-61f669fe0334.html

http://jonny131.iteye.com/blog/188829以及维基百科

======================================================================

首先有两个函数：

rand（产生随机数）
表头文件: #include<stdlib.h>
定义函数 :int rand(void)

srand（设置随机数种子的值）
表头文件：#include<stdlib.h>
定义函数：void srand (unsigned int seed);

rand的内部实现是用线性同余法做的，它不是真的随机数，而是一个伪随机数，是根据一个数（可称为“种子”）为基准以某个递推公式推算出来的一系数，当这系列数很大的时候，就符合正态公布，从而相当于产生了随机数，这不是真正的随机数，当计算机正常开机后，这个种子的值是定了的，除非你破坏了系统，为了改变这个种子的值，C提供了 srand()函数，它的原形是void srand( int a) 功能是初始化随机产生器既rand()函数的初始值，即使把种子的值改成a; 只不过是因为其周期特别长，所以有一定的范围里可看成是随机的，rand()会返回一随机数值，范围在0至RAND_MAX 间。在调用此函数产生随机数前，必须先利用srand()设好随机数种子，如果未设随机数种子，rand()在调用时会自动设随机数种子为1。rand ()产生的是假随机数字，每次执行时是相同的。若要不同,以不同的值来初始化它.初始化的函数就是srand()。

线性同余法求伪随机数的公式：
X(n+1) = (a *X(n) + c) mod m
公式中m > 0 , 0<=a<m, 0<c<m,  0<=X(0)<m。 方程求解过程详细参考：维基百科_线性同余方程
公式主要用到递归方法得到模除后的伪随机数,从X(0)到X(n),每一项都只依赖于前面一项，X序列将产生m个不同的数，然后重复。例如：

int rand(unsigned int seed)  
{  
   return (unsigned int)(seed * 1103515245 + 12345) % 32768;  
}  

rand()函数返回0至RAND_MAX之间的随机整数值，RAND_MAX是32767(int），即双字节（16位数）。若用unsigned int 双字节是65535，四字节是4294967295的整数范围。
0~RAND_MAX每个数字被选中的机率是相同的。

范例：

/* 产生介于1 到10 间的随机数值，此范例未设随机数种子，完整的随机数产生请参考
srand（）*/
#include<stdlib.h>
main()
{
	int i,j;
	for(i=0;i<10;i++)
	{
		j=1+(int)(10.0*rand()/(RAND_MAX+1.0));
		printf("%d ",j);
	}
}

  执行：
9 4 8 8 10 2 4 8 3 6
9 4 8 8 10 2 4 8 3 6 //再次执行仍然产生相同的随机数

srand()用来设置rand()产生随机数时的随机数种子。参数seed必须是个整数，通常可以利用geypid()或time(0)的返回值来当做seed。如果每次seed都设相同值，rand()所产生的随机数值每次就会一样。

范例

/* 产生介于1 到10 间的随机数值，此范例与执行结果可与rand（）参照*/
#include<time.h>
#include<stdlib.h>
main()
{
	int i,j;
	srand((int)time(0));
	for(i=0;i<10;i++)
	{
		j=1+(int)(10.0*rand()/(RAND_MAX+1.0));
		printf(" %d ",j);
	}
}

  执行：与rand范例比较
5 8 8 8 10 2 10 8 9 9
2 9 7 4 10 3 2 10 8 7

当然也可以用"int x = rand() % 10;"来生成 0 到 10 之间的随机数，但这种方法不如上边的“j=(int)(n*rand()/(RAND_MAX+1.0))”较好一些。

我想可能因为如果没有设置种子的值的话前边的方法产生的值都会是0。

int main(void)
{
  int i;
  time_t t;
  srand((unsigned) time(&t));
  printf("Ten random numbers from 0 to 99\n\n");
  for(i=0; i<10; i++)
  printf("%d\n", rand() % 100);
  return 0;
}

除以上所说的之外，补充一点就是srand这个函数一定要放在循环外面或者是循环调用的外面，否则的话得到的是相同的数字。

MSDN中的例子。

// crt_rand.c
// This program seeds the random-number generator
// with the time, then displays 10 random integers.
//
#include <stdlib.h>
#include <stdio.h>
#include <time.h>
int main( void )
{
  int i;
   
  // Seed the random-number generator with current time so that
  // the numbers will be different every time we run.
  //
  srand( (unsigned)time( NULL ) );//randomize();
  // Display 10 numbers.
  for( i = 0; i < 10;i++ )
  printf( " %6d", rand() );
  printf("\n");
  // Usually, you will want to generate a number in a specific range,
  // such as 0 to 100, like this:
  int RANGE_MIN = 0;
  int RANGE_MAX = 100;
  for (i = 0; i < 10; i++ )
  {
  int rand100 = (((double) rand() /
  (double) RAND_MAX) * RANGE_MAX + RANGE_MIN);
  printf( " %6d", rand100);
  }
 }

我们可以产生可以预见的随机序列，那我们如何才能产生不可预见的随机序列呢？利用srand((unsign)(time(NULL))是一种方法，因为每一次运行程序的时间是不同的，对了，你知道time() 函数的功能是返回从1970/01/01到现在的秒数的吧，可能这个起始时间不正确，你查一下对不对吧，C还提供了另一个更方便的函数， randomize()原形是void randomize()，功能是用来始初rand() 的种子的初始值，而且该值是不确定的，它相当于srand((unsign)(time(NULL)) 不过应注意的是randomize()的功能要通过time来实现所以在调用它时头文件要包含time.h罢了( 可以在TC2.0中运行通过，在VC6.0不能运行)。

测试：

对于rand函数，我用大数据测试了一下，看看是否等概率。

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <string.h>
int main( void )
{
  int i;
  int test[100000];
  memset(test, 0, sizeof(test));
  srand( (unsigned)time( NULL ) );//randomize();
  // Display 10 numbers.
  for( i = 0; i < 100000000;i++ )
    test[(int)rand()]++;
  for(i = 0; i < 10; ++i)
  {
    printf("%d\n", test[i]);
  }
  return 0;
}

   执行：

3052
2968
3095
3033
3107
3008
3034
3071
3068
3167

平均每个数出现的概率约为3000/100000000 = 0.00003, 近似等于1/32767 = 0.0000305...看来还是挺符合的。