梅森旋转算法

梅森旋转算法（Mersenne twister）是一个伪随机数发生算法。由松本真和西村拓士在1997年开发，基于有限二进制字段上的矩阵线性递归。可以快速产生高质量的伪随机数，修正了古典随机数发生算法的很多缺陷。

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31

def _int32(x):     return int(0xFFFFFFFF & x)   class MT19937:     def __init__(self, seed):         self.mt = [0] * 624         self.mt[0] = seed         self.mti = 0         for i in range(1, 624):             self.mt[i] = _int32(1812433253 * (self.mt[i - 1] ^ self.mt[i - 1] >> 30) + i)         def extract_number(self):         if self.mti == 0:             self.twist()         y = self.mt[self.mti]         y = y ^ y >> 11         y = y ^ y << 7 & 2636928640         y = y ^ y << 15 & 4022730752         y = y ^ y >> 18         self.mti = (self.mti + 1) % 624         return _int32(y)         def twist(self):         for i in range(0, 624):             y = _int32((self.mt[i] & 0x80000000) + (self.mt[(i + 1) % 624] & 0x7fffffff))             self.mt[i] = (y >> 1) ^ self.mt[(i + 397) % 624]               if y % 2 != 0:                 self.mt[i] = self.mt[i] ^ 0x9908b0df