尽管在计算机中并没有一个真正的随机数发生器,但是可以做到使产生的数字的重复率很低,以至于它们看起来是随机的。实现这一功能的程序叫做伪随机数发生器。

有关如何产生随机数的理论有许多,这里不讨论这些理论及相关的数学知识。因为讨论这一主题需要整整一本书的篇幅。这里要说的是,不管你用什么办法实现随机数发生器,你都必须给它提供一个被称为“种子(seed)”的初始值,而且这个值最好是随机的,或者至少是伪随机的。“种子”的值通常是用快速计数寄存器或移位寄存器来生成的。

本文中,笔者将介绍c语言所提供的随机数发生器的用法。现在的c编译程序都提供了一个基于一种ANSI标准的伪随机数发生器函数,用来生成随机数。Microsoft和Borland都是通过rand()和srand()函数来支持这种标准的,它们的工作过程如下:

  1. 首先,给srand()提供一个“种子”,它是一个unsignde int类型,其取值范围是从0到65,535 ;
  2. 然后,调用rand(),它会根据提供给srand()的“种子”值返回一个随机数(在0到32,767之间);
  3. 根据需要多次调用rand(),从而不断地得到新的随机数;
  4. 无论什么时候,你都可以给srand()提供一个新的“种子”,从而进一步“随机化"rand()的输出结果。

这个过程看起来很简单,问题是如果你每次调用srand()时都提供相同的“种子”值,那么你将会得到相同的“随机”数序列。例如,在以17为“种子”值调用srand()之后,在你首次调用rand()时,你将得到随机数94;在你第二次和第三次调用rand()时,你将分别得到26,602和30,017。这些数看上去是相当随机的(尽管这只是一个很小的数据点集合),但是,在你再次以17为“种子”值调用srand()之后,在对rand()的前三次调用中,所得到的返回值仍然是94、26,602和30,017,并且此后得到的返回值仍然是在对rand()的第一批调用中所得到的其余的返回值。因此,只有再次给srand()提供一个随机的“种子”值,才能再次得到一个随机数。

下面的例子用一种简单而有效的办法来产生一个相当随机的“种子”值——当天的时间值。

# include <stdlib. h>
# include <stdio. h>
# include <sys/types. h>
# include <sys/timeb. h>
void main (void){
int i ;
unsigned int seedVal;
struct_timeb timeBuf ;
_ftime (&timeBuf) ;
seedVal = ( ( ( ( (unsigned int)timeBuf, time & 0xFFFF) +
(unsigned int)timeBuf, millitm) ^
(unsigned int)timeBuf, millitm) ;
srand ((unsigned int)seedVal) ;
for(i=O;i<lO;++i)
printf (" %6d\n" ,rand ( ) ) ;
}

上例先是调用_ftime()来检索当前时间,并把它的值存入结构成员timeBuf.time中,当前时间的值从1970年1月1日开始以秒计算。在调用了_ftime()之后,在结构timeBuf的成员millitm中还存入了在当前那一秒已经度过的毫秒数,但在DOS中这个数字实际上是以百分之一秒来计算的。然后,把毫秒数和秒数相加,再和毫秒数进行一次异或运算。你可以对这两个结构成员施加更多的逻辑运算,以控制seedVal的取值范围,并进一步加强它的随机性,但上例所用的逻辑运算已经足够了。

注意,在前面的例子中,rand()的输出并没有被限制在一个指定的范围内,假设你想建立一个彩票选号器,其取值范围是从1到44。你可以简单地忽略掉rand()所输出的在该范围之外的值,但这将花费许多时间去得到所需的全部(例如6个)彩票号码。假设你已经建立了一个令你满意的随机数发生器,它所产生的随机数据范围是从0到32,767(就象前文中提到过的那样),而你想把输出限制在1到44之间,下面的例子就说明了如何来完成这项工作:

int i ,k ,range ;
int rain, max ;
double j ;
min=1; /* 1 is the minimum number allowed */
max=44; /* 44 is the maximum number allowed */
range=max-min; /* r is the range allowed; 1 to 44 */
i=rand(); /* use the above example in this slot */
/* Normalize the rand() output (scale to 0 to 1) */
/* RAND_MAX is defined in stdlib, h */
j= ((double)i/(double)RAND_MAX) ;
/* Scale the output to 1 to 44 */
i= (int)(j * (double)range) ;
i+ =min;

上例把输出的随机数限制在1到44之间,其工作原理如下:

  1. 得到一个在O到RAND_MAX(32,767)之间的随机数,把它除以RAND_MAX,从而产生一个在0到1之间的校正值;
  2. 把校正值乘以所需要的范围值(在本例中为43,即44减去1),从而产生一个在O到43之间的值;
  3. 把该值和所要求的最小值相加,从而使该值最终落在正确的取值范围——1到44之内。

你可以用不同的min和max值来验证这个例子,你会发现它总是会正确地产生在新的rain和max值之间的随机数。

C语言生产随机数的方法的更多相关文章

  1. 用c语言产生随机数的方法

    用c语言产生随机数的方法 在C语言中,rand()函数可以用来产生随机数,但是这不是真正意义上的随机数,是一个伪随机数,是根据一个数,我们可以称它为种子,为基准以某个递推公式推算出来的一系数,当这系列 ...

  2. c语言产生随机数的方法

    在C语言中,rand()函数可以用来产生随机数,但是这不是真真意义上的随机数,是一个伪随机数,是根据一个数,我们可以称它为种子,为基准以某个递推公式推算出来的一系数,当这系列数很大的时候,就符合正态公 ...

  3. 【转】linux shell实现随机数多种方法(date,random,uuid)

    在日常生活中,随机数实际上经常遇到,想丢骰子,抓阄,还有抽签.呵呵,非常简单就可以实现.那么在做程序设计,真的要通过自己程序设计出随机数那还真的不简单了.现在很多都是操作系统内核会提供相应的api,这 ...

  4. C语言中随机数相关问题

    用C语言产生随机数重要用到rand函数.srand函数.及宏RAND_MAX(32767),它们均在stdlib.h中进行了声明. int rand(void);//生成一个随机数 voidsrand ...

  5. Python语言上机题实现方法(持续更新...)

    Python语言上机题实现方法(持续更新...) 1.[字符串循环左移]给定一个字符串S,要求把S的前k个字符移动到S的尾部,如把字符串"abcdef"前面的2个字符'a'.'b' ...

  6. 深入理解Java 8 Lambda(语言篇——lambda,方法引用,目标类型和默认方法)

    作者:Lucida 微博:@peng_gong 豆瓣:@figure9 原文链接:http://zh.lucida.me/blog/java-8-lambdas-insideout-language- ...

  7. 现有语言不支持XXX方法

    史上最强大的IDE也会有bug的时候哈,今天遇到这个问题特别郁闷,百度了下,果然也有人遇到过这个问题 解决方法: 1.调用的时候参数和接口声明的参数不一致(检查修改) 2.继承接口中残留一个废弃的方法 ...

  8. 【erlang】erlang几种生成随机数的方法

    erlang有三个生产随机数的办法 random:uniform(). 这个函数是erlang库random模块提供的.一般都采用这个. 1> random:uniform(). 0.44358 ...

  9. java 产生随机数的方法

    有三种方法: Math.random():这个方法返回一个[0.0, 1.0)的一个随机double型数.它实际是调用Random类的nextDouble()方法.只不过Math类使用的是一个静态随机 ...

随机推荐

  1. hdu 畅通工程再续

    http://acm.hdu.edu.cn/showproblem.php?pid=1875 #include <cstdio> #include <cstring> #inc ...

  2. selenium太有爱,我已离不开!!!

    自动化测试,超有用. PROXY,PLUGIN,PROFILE,WINDOWS HANDLE个个搞定!!! from selenium import webdriver from selenium.c ...

  3. Unity 3D中的菜单项

    1.File菜单:主要是包含项目和场景的创建.保存以及输出等功能.2.Edit(编辑)菜单:只要包括对场景进行一系列的编辑以及环境设置操作等命令.3.Assets(资源)菜单:掌握资源在Unity中的 ...

  4. PSAM读卡芯片TDA8007BHL开发

    WWT:Work Waiting Time ATR:Answer To Reset,复位应答 etu =F/Df 1.     PSAM概述和应用 PSAM(PurchaseSecure Access ...

  5. logstash tag使用说明

    zjtest7-frontend:/usr/local/logstash-2.3.4/config# cat stdin04.conf input { stdin { } } filter { # d ...

  6. NOI2013 UOJ122 向量内积

    神题...... 还是大神讲得比较清晰~orz http://dffxtz.logdown.com/posts/197950-noi2013-vector-inner-product 启发题:poj3 ...

  7. Delphi 调试 通过BreakPoint

    1.打个断点, 如下图 2. 在断点上,邮件,如下图 3. 弹出一个窗体 ,如下图 在 condition 中写条件就可以了.  这样就可以按你假设的条件来进行了,方便.

  8. easyui获取日期datebox中的值

    <input type="text" class="easyui-datebox" id="CTIME" style="wi ...

  9. hdu 4983 Goffi and GCD(欧拉函数)

    Problem Description Goffi is doing his math homework and he finds an equality on his text book: gcd( ...

  10. python高级编程:缓存

    # -*- coding: utf-8 -*-__author__ = 'Administrator'#缓存"""对于运行代价很高的函数和方法结果,可以进行缓存,只要:1 ...