不合时宜的PASCAL问题

DoDo

请不要怪我在这个论坛中提出关于PASCAL的问题,但有人能够解答么?
这是别人考我的一道题:

1. 
   
2. 
   
3. 产生随机数的函数中
   
4. FUNCTION rand(VAR seed:real):real;
   
5. CONST
   
6.    a=93.0;
   
7.    m=8192.0;
   
8.    c=1.0;
   
9. BEGIN
   
10.    seed:=a*seed+c;
    
11.    seed:=round((seed/m-trunc(seed/m))*m);
    
12.    rand:=seed/m
    
13. END;
    
14. seed的初始量可以是小于8192.0的任一实数，那a和c是怎样确定的呢？
    

复制代码

无双

PASCAL也是一种设计语言
所以都可以讨论

a和c不是在变量定义部分定义好了吗
const后面

DoDo

他的意思是,那几个数的值是因何而定的

无双

c=1是为了不让seed=0  防止除0错误
a 和m是不是平方关系
总之要返回的值在那个范围内

DoDo

a和m不是刚好平方关系,m开平方约为91

我也认为a和m的值是无关紧要的,只要保证结果不溢出,也不因精度忽略产生问题,那么稍大稍小都无所谓.

但事情是这样的:我在班级的论坛主持编程版块,给大家讲基础程序设计,刚讲完QUICK BASIC,正在讲PASCAL,在一篇我写的教程之后,那个人跟了此帖.这段程序我记得是清华大学那本黄皮教材上的一个例子,但我那本书早被同学抢走了,便无从考证.而我当然不能按自己的猜测去回复同学,所以才到处发问,但都没有找到理想的答案.

无双

seed的初始量可以是小于8192.0的任一实数，那a和c是怎样确定的呢？

函数的输入输出是什么 c=1是为了避免出现 seed =0 情况

你可以把参数根据最大最小值可能出现的结果判断一下
然后再猜测

DoDo

唯一的发现是8192是32767的1/4

libinary

这个应该和写这个函数的人有关，取值好象没什么规定，
不过这个函数还是有点怪，一般随机数生成都是一个序列，这个函数只能生成一个，下一次调用还得把seed传回去。
简单的伪随机数生成一般用Xn = (Xn-1 * a + b) % c 其中a、b的值好象没什么规定(但是应该是有一定的技巧的，比如a，上面的函数就取了93，是一个质数)，c就是随机数的范围。
上面因为用了real，所以用seed:=round((seed/m-trunc(seed/m))*m);实现类似整除的功能。

libinary

这有一篇文章，中间有Borland的实现，可以参考一下。
http://www-900.ibm.com/developer ... playing/index.shtml

DoDo

Thank you, I'll show him the article.