随机数的实现原理

凌风探梅

于 2016-04-25 09:47:42 发布

阅读量607

点赞数

分类专栏： C C++

C++ 同时被 2 个专栏收录

40 篇文章

订阅专栏

37 篇文章

订阅专栏

http://en.wikipedia.org/wiki/Random_number_generation

The GLIBC random number generator 讲了GLIBC中random()函数的实现原理，并给出了一段代码：

 
        #include <stdio.h>
       
        #define MAX 1000
       
        #define seed 1
       
        main() {
       
        int 
         r[MAX];
       
        int 
         i;
       
        r[0] = seed;
       
        for 
         (i=1; i<31; i++) {
       
        r[i] = (16807LL * r[i-1]) % 2147483647;
       
        if 
         (r[i] < 0) {
       
        r[i] += 2147483647;
       
        }
       
        }
       
        for 
         (i=31; i<34; i++) {
       
        r[i] = r[i-31];
       
        }
       
        for 
         (i=34; i<344; i++) {
       
        r[i] = r[i-31] + r[i-3];
       
        }
       
        for 
         (i=344; i<MAX; i++) {
       
        r[i] = r[i-31] + r[i-3];
       
        printf
        (
        "%d\n"
        , ((unsigned 
        int
        )r[i]) >> 1);
       
        }
       
        }