Marsaglia XORshift随机数算法

引理:
二进制矩阵相乘中加法为异或。
给定非空向量beta和n阶矩阵T,序列beta*T,beta*T^2,beta*T^3,...的秩为2^n-1的充要条件是矩阵T为非奇异矩阵.


L是左移位操作,y=y^(y<<a)表示为y=y*(E+L^a),y=y^(y>>b)表示为y=y*(E+R^b),令T=(E+L^a)(E+R^b),n=32或64,找不到这样的非奇异矩阵。
但是令T=(E+L^a)*(E+R^b)*(E+L^c)能找到很多。对于n=32,T=(E+L^a)*(E+R^b)*(E+R^c)有81个可选的(a,b,c)三元组。对于n=64,T=(E+L^a)*(E+R^b)*(E+L^c)有275个可选三元组。还有变种的xor128和mwc周期可以更长。
该算法相当于一个线性反馈移位寄存器,是最快的非密码安全随机数生成器(肯定比线性同余法快得多),大约只要10个时钟周期。java的IdentityHashCode就是用这个算法,不过如果要做加密或蒙特卡洛就不行了。


参考资料:
C++ 快速随机数生成器:http://blog.jobbole.com/94034/
wiki: https://en.wikipedia.org/wiki/Xorshift
XORShift RNGS: http://www.jstatsoft.org/v08/i14/paper
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值