[Z]K.I.S.S.Random Genrator “保持简单”随机数发生器
2013-03-17 15:08
218 查看
(假)随机数发生器是Monte Carlo算法的基础,自然大家格外关注。
近几年,日本某君发明的Mersenne Twister[Makoto Matsumoto]可谓其中翘楚。
不过,它并非尽善尽美。一则它不适合用于数据加密,二则不能通过一些统计检验。
三嘛就是源程序比较长,不是一眼就能看明白的。
反过来,上世纪90年代后期就被Marsaglia[Marsaglia]发明出来的Keep-It-Simple-And-Stupid
算法却挺有趣而简短,而且能打通DIEHARD的各种检验[DIEHARD]。
下面是源代码:
unsigned int x = 123456789,
y = 362436000,
z = 521288629,
c = 7654321; /* Seed variables */
unsigned int KISS()
{
unsigned long long t, A = 698769069ULL;
x = 69069*x+12345;
y ^= (y<<13); y ^= (y>>17); y ^= (y<<5);
t = (A*z + c);
c = (t >> 32);
z = t;
return x+y+z;
}
居然就只是把“线性同余”,“移位轮转”和“带记忆乘法”这3种基本的随机数发生法一起用,
便获得很好的效果。的确够巧的,从这一点上说,比Mersenne Twister要有趣多了。
而且通过它,我们一次就学会3种随机数发生法。抽出其中任何1种,都可以单独使用,仿佛一把瑞士军刀。
[Makoto Matsumoto]Mersenne Twister Home Page http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html [Marsaglia]George_Marsaglia Wiki Page
en.wikipedia.org/wiki/George_Marsaglia
[DIEHARD]The Marsaglia Random Number CDROM http://stat.fsu.edu/pub/diehard/
近几年,日本某君发明的Mersenne Twister[Makoto Matsumoto]可谓其中翘楚。
不过,它并非尽善尽美。一则它不适合用于数据加密,二则不能通过一些统计检验。
三嘛就是源程序比较长,不是一眼就能看明白的。
反过来,上世纪90年代后期就被Marsaglia[Marsaglia]发明出来的Keep-It-Simple-And-Stupid
算法却挺有趣而简短,而且能打通DIEHARD的各种检验[DIEHARD]。
下面是源代码:
unsigned int x = 123456789,
y = 362436000,
z = 521288629,
c = 7654321; /* Seed variables */
unsigned int KISS()
{
unsigned long long t, A = 698769069ULL;
x = 69069*x+12345;
y ^= (y<<13); y ^= (y>>17); y ^= (y<<5);
t = (A*z + c);
c = (t >> 32);
z = t;
return x+y+z;
}
居然就只是把“线性同余”,“移位轮转”和“带记忆乘法”这3种基本的随机数发生法一起用,
便获得很好的效果。的确够巧的,从这一点上说,比Mersenne Twister要有趣多了。
而且通过它,我们一次就学会3种随机数发生法。抽出其中任何1种,都可以单独使用,仿佛一把瑞士军刀。
[Makoto Matsumoto]Mersenne Twister Home Page http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html [Marsaglia]George_Marsaglia Wiki Page
en.wikipedia.org/wiki/George_Marsaglia
[DIEHARD]The Marsaglia Random Number CDROM http://stat.fsu.edu/pub/diehard/
相关文章推荐
- K.I.S.S.Random Genrator/ “保持简单”随机数发生器
- JavaScript当中网页刷新后数值(text、radio、checkbox、select、textarea)保持不变(输入域数据保持)的简单方法
- 粘贴代码更简单的方法(保持VS.NET中的格式)
- 简单概率dp-hdu-4487-Maximum Random Walk
- 一个简单的注册表单验证保持原来输入的内容实例
- 保持简单苹果发展早期的设计理念
- 保持简单----纪念丹尼斯•里奇(Dennis Ritchie
- 保持简单----纪念丹尼斯里奇(Dennis Ritchie) (zz.is2120.BG57IV3)
- 使用scatter与random.normal函数绘制散点图简单实例
- 分享一个简单的页面刷新保持菜单不变的好东东
- 一个简单的注册表单验证保持原来输入的内容实例
- 保持简单苹果发展早期的设计理念
- 保持简单----纪念丹尼斯·里奇(Dennis Ritchie)
- 如何避免被颠覆掉的命运?保持你产品的简单
- 产品保持简单的重要性
- 保持简单苹果发展早期的设计理念
- HttpClient简单请求和session的保持
- 保持专注,保持简单
- Java中Random的简单使用说明
- 让你的伪造保持简单