这个东西,对于数学上的分析我看明白了90%.但并未做深入研究.这个代码是抄来的,核心部分是抄来的.感觉不错.
生成固定的随机数序列,只要种子一样的话.这个随机数生成器比较透明,而且较容易改动.很不错.
每次实现代码的时候,都会发现很多问题.而且,一旦写上了代码,时光就飞逝了...很不错.
好吧,帖出来吧.又是来网吧,又是边下电影边写Blog,呵呵.很不错.
/*10-35-11-02-08-20.11.c -- 第十章第三十五题*/
#include
#define FALSE (0)
#define TRUE (1)
#define A (48271L)
#define M (2147486647L)
#define Q (M / A)
#define R (M % A)
unsigned long int Seed ;
int main (void) ;
int eatEnter (void) ;
void initialize (unsigned long int initial_value) ;
double random (void) ;
int main (void)
{
unsigned long int initial_value ;
char ch ;
printf ("Please input a initial value (>= 0, enter 'q' to quit.):") ;
while (scanf ("%lu", &initial_value) != 0)
{
printf ("Initial value has been initialized is : %lu/n", initial_value) ;
eatEnter () ;
initialize (initial_value) ;
puts ("Enter 'q' to input a new initial value and enter any key else to generate a random value again.") ;
while ((ch = getchar ()) != 'q')
{
if (TRUE == eatEnter ())
break ;
printf ("It is %lf./n", random ()) ;
}
printf ("Please input a initial value (>= 0, enter 'q' to quit.):") ;
}
return 0 ;
}
int eatEnter (void)
{
if ('q' == getchar ())
return TRUE ;
else
return FALSE ;
}
void initialize (unsigned long int initial_value)
{
Seed = initial_value ;
}
double random (void)
{
unsigned long int temp_seed ;
temp_seed = A * (Seed % Q) - R * (Seed / Q) ;
if (temp_seed >= 0)
Seed = temp_seed ;
else
Seed = temp_seed + M ;
return (double) Seed / M ;
}
