C++ 随机数rand()和srand()

       标准函数库中函数rand()可以生成0~RAND_MAX之间的一个随机数,其中RAND_MAX 是stdlib.h中定义的一个整数,它与系统有关。 例如在我的机器上,RAND_MAX=32767。可以使用 printf("%ld",RAND_MAX);查看。

调用rand()可以生成一个随机数,但我们往往需要的是一个有上下界的随机数。

例如 要产生一个[a,b]之间的随机整数x

x = rand()%(b-a+1)+a
产生区间[a,b]上的随机实数
z = ((double)rand()/RAND_MAX)*(b-a) + a

rand()简介:   

   生成伪随机数的函数包含在标准库 <cstdlib>中,所以首先需在头文件中包含进来。

    注意:在C语言中包含在标准库<stdlib.h>中。
    函数一:int   rand(void); 
    从srand   (seed)中指定的seed开始,返回一个[seed,   RAND_MAX(0x7fff))间的随机整数。 

    函数二:void   srand(unsigned   seed); 
    参数seed是rand()的种子,用来初始化rand()的起始值。

    用法:它需要提供一个种子,这个种子会对应一个随机数,如果使用相同的种子后面的rand()函数会出现一样的随机数。如: srand(1); 直接使用1来初始化种子。不过为了防止随机数每次重复常常使用系统时间来初始化,即使用 time函数来获得系统时间,它的返回值为从 00:00:00 GMT, January 1, 1970 到现在所持续的秒数,然后将time_t型数据转化为(unsigned)型再传给srand函数,即: srand((unsigned) time(&t)); 还有一个经常用法,不需要定义time_t型t变量,即: srand((unsigned) time(NULL)); 直接传入一个空指针,因为你的程序中往往并不需要经过参数获得的t数据。srand((int)getpid()); 使用程序的ID(getpid())来作为初始化种子,在同一个程序中这个种子是固定的。

    可以认为rand()在每次被调用的时候,它会查看: 
    1) 如果用户在此之前调用过srand(seed),给seed指定了一个值,那么它会自动调用 srand(seed)一次来初始化它的起始值。 
    2) 如果用户在此之前没有调用过srand(seed),它会自动调用srand(1)一次。 

   根据上面的第一点我们可以得出: 
   1) 如果希望rand()在每次程序运行时产生的值都不一样,必须给srand(seed)中的seed一个变值,这个变值必须在每次程序运行时都不一样(比如到目前为止流逝的时间)。 
   2) 否则,如果给seed指定的是一个定值,那么每次程序运行时rand()产生的值都会一样,虽然这个值会是[seed,   RAND_MAX(0x7fff))之间的一个随机取得的值。 
   3) 如果在调用rand()之前没有调用过srand(seed),效果将和调用了srand(1)再调用rand()一样(1也是一个定值)。 
举几个例子,假设我们要取得0~6之间的随机整数(不含6本身):

例一,不指定seed: 

for(int   i=0;i <10;i++){
    ran_num=rand()   %   6; 
    cout < <ran_num < < "   "; 
} 
每次运行都将输出:5   5   4   4   5   4   0   0   4   2 

例二,指定seed为定值1: 

srand(1); 
for(int   i=0;i <10;i++){
    ran_num=rand()   %   6; 
    cout < <ran_num < < "   "; 
} 
每次运行都将输出:5   5   4   4   5   4   0   0   4   2 
跟例子一的结果完全一样。 

例三,指定seed为定值6: 
srand(6); 
for(int   i=0;i <10;i++){
    ran_num=rand()   %   6; 
    cout < <ran_num < < "   "; 
} 
每次运行都将输出:4   1   5   1   4   3   4   4   2   2 
随机值也是在[0,6)之间,随得的值跟srand(1)不同,但是每次运行的结果都相同。 

例四,指定seed为当前系统流逝了的时间(单位为秒):time_t   time(0): 
#include   <ctime> 
//… 
srand((unsigned)time(0)); 
for(int   i=0;i <10;i++){
    ran_num=rand()   %   6; 
    cout < <ran_num < < "   "; 
} 

第一次运行时输出:0   1   5   4   5   0   2   3   4   2 
第二次:3   2   3   0   3   5   5   2   2   3 
总之,每次运行结果将不一样,因为每次启动程序的时刻都不同(间隔须大于1秒?见下)。 

关于time_t   time(0):
time_t被定义为长整型,它返回从1970年1月1日零时零分零秒到目前为止所经过的时间,单位为秒。比如假设输出:
cout < <time(0); 
值约为1169174701,约等于37(年)乘365(天)乘24(小时)乘3600(秒)(月日没算)。
另外,关于ran_num   =   rand()   %   6,
将rand()的返回值与6求模是必须的,这样才能确保目的随机数落在[0,6)之间,否则rand()的返回值本身可能是很巨大的。 
一个通用的公式是: 
要取得[a,b)之间的随机整数,使用(rand()   %   (b-a))+   a   (结果值将含a不含b)。 
在a为0的情况下,简写为rand()   %   b。 

最后,关于伪随机浮点数
用rand()   /   double(RAND_MAX)可以取得0~1之间的浮点数(注意,不同于整型时候的公式,是除以,不是求模),举例: 
double   ran_numf=0.0; 
srand((unsigned)time(0)); 
for(int   i=0;i <10;i++){
    ran_numf   =   rand()   /   (double)(RAND_MAX); 
    cout < <ran_numf < < "   "; 
} 
运行结果为:0.716636,0.457725,…等10个0~1之间的浮点数,每次结果都不同。 

如果想取更大范围的随机浮点数,比如1~10,可以将 
rand()   /(double)(RAND_MAX)   改为   rand()   /(double)(RAND_MAX/10) 
运行结果为:7.19362,6.45775,…等10个1~10之间的浮点数,每次结果都不同。 
至于100,1000的情况,如此类推。

计算机产生随机序列rand()(线性同余方法)

众所周知,计算机产生随机序列实际上是一个有很长周期的序列(通过线性同余算法),这就是伪随机数。随机数产生时需要一个seed(种子),才能产生随机数。

在Pascal 和 Basic 中,都可以使用 randomize;来初始化种子。而在C中,需要用srand(unsigned int)函数给予一个种子。 例如srand(time(NULL));

一般来说在OI中,依靠随机数的算法较为稳定,它不受随机序列的影响(难道有人看过默认的随机序列,专门设计出一组卡随机的测试数据????)。所以srand调用与否没有太大关系。况且time()函数定义在 time.h 中,某些评测系统禁止使用time.h。我一般都是srand(rand());。

线性同余方法(LCG)是个产生伪随机数的方法。
它是根据递归公式:
其中A,B,M是产生器设定的常数。
LCG的周期最大为M,但大部分情况都会少于M。要令LCG达到最大周期,应符合以下条件:
   1. B,M互质;
   2. M的所有质因子的积能整除A − 1;
   3. 若M是4的倍数,A − 1也是;
   4. A,B,N0都比M小;
   5. A,B是正整数。

使用线性同余算法产生伪随机数 转自http://www.blog.edu.cn/user1/20989/archives/2005/989163.html C语言中有个random(n)函数,可以产生0——n-1之间的伪随机数,rand()函数可以产生0——32767之间的伪随机数,一般需要配合使用srand(long)给出随机种子,或者使用randomize()函数来根据系统时间指定随机种子。

     常用的产生伪随机数的算法是线性同余法,下面是我写的代码,相信通过代码,大家可以很容易地明白伪随机数原理:

     算法:an+1=b*an+c mod m;a0=d。d为种子。
int d;/*种子*/
long My_Rand(long b,long c,long m)
{
   return  d=(b*d+c)%m;/*产生随机数,并记录,用来产生下一个伪随机数*/
}

long gcd1(long m,long n)
{/*求最大公约数,算法中的b和m一般取互素*/
   long r;
   while((r=m%n)!=0)
   {
     m=n;
     n=r;
   }
   return n;
}

long gcd(long m)
{/*计算与m互素的b*/
   long i=2;
   while(gcd1(i,m)!=1)
      i++;
   return i;
}

void main()
{
   long b,c,m;
   int count;
   int i;
   m=32767;
   b=gcd(m);
   c=5001;
   d=3;
   count=800;
   for(i=0;i<count;i++)/*产生800个伪随机数*/
   {
      if(i%10==0)
           printf("
");
      printf("%8ld",My_Rand(b,c,m));
   }
}
原文地址:https://www.cnblogs.com/tham/p/6827415.html