伪随机数生成器
2015-06-18 10:39
369 查看
在所有的编程语言中,基本上都会有产生伪随机数的库函数。长期以来,我一直很好奇,这些函数,究竟是如何产生伪随机数的。最近通过stack overflow和维基百科,我大概了解了一下伪随机数的生成过程。这篇文章的目的,并不是帮助大家了解很多伪随机算法的原理、效率、历史,也不是为了深入某一个或者某几个伪随机数算法。这篇文章的意义,在于帮助大家消除对于伪随机数生成过程某种程度上的好奇心。对于伪随机生成器更广阔以及更深入的理解,请移步https://en.wikipedia.org/wiki/Random_number_generation。
事实上,可以生产伪随机数的算法有很多,比如说Mersenne twister,LCG(Linear
congruential generator),和XORShift,在这里我们只是简单介绍一下LCG和XORShift。
XORShift是一种优秀,短小并且简洁的随机数生成器,它通过了Diehard
tests,Diehard tests是一连串的统计测试,用来测试生成器产生的伪随机数的质量。它的一种实现形式如下:
上述公式中的三个常量是通过谨慎选择的,一种典型的选择方式如下:
其C++实现如下,代码中的BSD和MS可视为,根据具体操作系统平台,运行不同代码。
事实上,可以生产伪随机数的算法有很多,比如说Mersenne twister,LCG(Linear
congruential generator),和XORShift,在这里我们只是简单介绍一下LCG和XORShift。
XORShift是一种优秀,短小并且简洁的随机数生成器,它通过了Diehard
tests,Diehard tests是一连串的统计测试,用来测试生成器产生的伪随机数的质量。它的一种实现形式如下:
#include <stdint.h> /* These state variables must be initialized so that they are not all zero. */ uint32_t x, y, z, w; uint32_t xorshift128(void) { uint32_t t = x ^ (x << 11); x = y; y = z; z = w; return w = w ^ (w >> 19) ^ t ^ (t >> 8); }另外一种算法为LCG。LCG利用不连续的分段线性方程产生伪随机数。它是最有历史并且最广为人知的一种伪随机数生产算法,其背后理论相对简单,并且易于实现。其用来产生伪随机数的公式通常如下:
number = (previous_number * constant + other_constant) mod third_constant
上述公式中的三个常量是通过谨慎选择的,一种典型的选择方式如下:
number = (previous_number * 214013 + 2531011) mod 2^15
其C++实现如下,代码中的BSD和MS可视为,根据具体操作系统平台,运行不同代码。
#include <iostream> //-------------------------------------------------------------------------------------------------- using namespace std; //-------------------------------------------------------------------------------------------------- class mRND { public: void seed( unsigned int s ) { _seed = s; } protected: mRND() : _seed( 0 ), _a( 0 ), _c( 0 ), _m( 2147483648 ) {} int rnd() { return( _seed = ( _a * _seed + _c ) % _m ); } int _a, _c; unsigned int _m, _seed; }; //-------------------------------------------------------------------------------------------------- class MS_RND : public mRND { public: MS_RND() { _a = 214013; _c = 2531011; } int rnd() { return mRND::rnd() >> 16; } }; //-------------------------------------------------------------------------------------------------- class BSD_RND : public mRND { public: BSD_RND() { _a = 1103515245; _c = 12345; } int rnd() { return mRND::rnd(); } }; //-------------------------------------------------------------------------------------------------- int main( int argc, char* argv[] ) { BSD_RND bsd_rnd; MS_RND ms_rnd; cout << "MS RAND:" << endl << "========" << endl; for( int x = 0; x < 10; x++ ) cout << ms_rnd.rnd() << endl; cout << endl << "BSD RAND:" << endl << "=========" << endl; for( int x = 0; x < 10; x++ ) cout << bsd_rnd.rnd() << endl; cout << endl << endl; return 0; } //--------------------------------------------------------------------------------------------------
相关文章推荐
- Objective C面向对象
- 如何正确的使用SharedPreferences?
- Tomcat工作原理
- 返回一个二维整数数组中最大联通子数组的和
- android 关于listView 显示不全的问题
- ns3不使用Ipv4AddressHelper的替代方法
- iOS链接错误,Undefined symbols for architecture arm64:std::__1::basic_string<char, std::__1::char_traits<
- hdu 4463 有一条边必须加上 (2012杭州区域赛K题)
- C#获取程序的当前路径,启动和关闭一个Process
- eclipse-android-activity_main/fragment_main文件处理
- jquery在线预览PDF文件,打开PDF文件
- 每天一道算法题(14)——N个降序数组,找到最大的K个数
- 每天一道算法题(14)——N个降序数组,找到最大的K个数
- 读书笔记——读《构建之法:现代软件工程》第13~17章
- ThinkPHP里用U方法调用js文件实例
- php实现mysql数据库分表分段备份
- 配置tomcat服务
- SVN 中的 trunk branches tag 详解
- 【转载】2015 Objective-C 三大新特性 | 干货
- Ubuntu下apt-get安装软件时缺少依赖包解决方法