随机函数rand()算法

今天学习了一下随机函数rand的算法。

这个算法叫做线性同余算法(linear congruential generator (LCG))。

不同的编译器取的常数不同，可以参考wiki：http://en.wikipedia.org/wiki/Linear_congruential_generator

以下程序可以输出和系统一样的随机数。

#include <stdio.h>
#include <stdlib.h>


#define _A 214013LL
#define _B 2531011LL


int a;


void mysrand(int x)
{
    a = x;
}


int myrand()
{
    return ((a = a * _A + _B)>>16) & 0x7fff;
}


int main()
{
    mysrand(0);
    srand(0);
    while (1)
    {
        printf ("%d, %d
", rand(), myrand());
        getchar();
    }
    return 0;
}