Project1

标题: 随机数发生器Mersenne Twister 照抄版 [打印本页]
作者: viktor    时间: 2012-11-8 14:24
标题: 随机数发生器Mersenne Twister 照抄版
本帖最后由 viktor 于 2012-11-8 14:29 编辑

嗯……其实ruby1.7以上的srand/rand函数本身就是用的mersenne twister产生的强随机数
其实ruby1.9自带了Random类，但是XP用的ruby1.81没有

总之 反正照着标准C程序mt19937ar.c抄了一份ruby的，可以制定随机数种子 和初始向量的随机数产生器。
用法
mt = MT19937.new(seed)    # seed是随机数种子
iv = [123, 456, 789, 1234, 5678] ￥ iv是初始向量，想填几个数字都行，自由发挥
mt.init_array iv # 可以在用种子初始化的基础上，再定义初始向量，增强随机性

随机数：
mt.randi

另外底下的常数除了birthday是原作者的生日（作者自重）以外其他都不要改了。后果貌似很严重……

  2. #============
  3. # rpg.blue
  4. # Mersenne Twister copycat version.
  5. # original: mt19937ar.c
  6. # author: viktor
  7. #============
  8. class MT19937
  9.   N = 624
  10.   M = 397
  11.   A = 0x9908b0df
  12.   MASK_S32 = 0x7fffffff
  13.   MASK_U32 = 0xffffffff
  14.   PRIMITIVE = 1812433253
  15.   BIRTHDAY = 19650218
  16.   MAGIC = [166_4525, 156_608_3941, 0x9d2c5680, 0xefc60000]
  17.   
  18.   attr_reader :seed
  19.   attr_reader :iv

  21.   def init_mt(s)
  22.     @mt[0] = s & MASK_U32
  23.     N.times{|i|@mt[i+1]=(PRIMITIVE * (@mt[i] ^ (@mt[i] >> 30))+i+1) & MASK_U32}
  24.     @life = N
  25.   end
  26.   
  27.   def init_array(iv)
  28.     @iv = iv   
  29.     i, j = 1, 0
  30.     [iv.length, N].max.times {
  31.       @mt[i] = (@mt[i] ^ (( @mt[i-1] ^ (@mt[i-1] >> 30 )) * MAGIC[0]))
  32.       @mt[i] =(@mt[i]+iv[j]+j) & MASK_U32
  33.       i+=1; j=(j+1) % iv.length
  34.       if i>=N
  35.         @mt[0] = @mt[N-1] # 折叠数组
  36.         i=1
  37.       end      
  38.     }
  39.     (N-1).times {
  40.       @mt[i] = (@mt[i] ^ (( @mt[i-1] ^ (@mt[i-1] >> 30 )) * MAGIC[1]))
  41.       @mt[i] =(@mt[i]-i) & MASK_U32
  42.       i+=1; if i>=N
  43.         @mt[0] = @mt[N-1]
  44.         i=1
  45.       end
  46.     }
  47.     @mt[0] = 0x80000000
  48.   end
  49.   
  50.   def randi
  51.     if @life >= N
  52.       # start a new round, update table mt
  53.       N.times {|k|
  54.         tmp = (@mt[k] & 0x80000000) | (@mt[(k+1)%N] & MASK_S32)
  55.         @mt[k] = @mt[(k+M)%N] ^ (tmp>>1) ^ (tmp[0]*A)
  56.       }
  57.       @life=0
  58.     end
  59.     ret = @mt[@life]
  60.     ret ^= (ret >> 11);
  61.     ret ^= (ret << 7) & MAGIC[2];
  62.     ret ^= (ret << 15) & MAGIC[3];
  63.     ret ^= (ret >> 18);
  64.     @life+=1
  65.     ret
  66.   end
  67.   
  68.   def initialize(seed=nil)
  69.     @mt=[]
  70.     srand(seed)
  71.   end
  72.   
  73.   def srand(seed=nil)
  74.     last_seed = @seed
  75.     seed = rand(1<<32) if seed==nil
  76.     @seed = seed
  77.     init_mt(@seed)
  78.     return last_seed
  79.   end
  80. end


  83. # test
  84. mt = MT19937.new(MT19937::BIRTHDAY)
  85. iv = [0x123, 0x234, 0x345, 0x456]
  86. mt.init_array iv
  87. out = []
  88. 1000.times {out << mt.randi}
  89. p out
复制代码
作者: 1243852    时间: 2012-11-8 15:34
其实，我确实不懂脚本。
其实，我确实看不懂你说的是什么，但是，应该很赞。
总之，反正我看不懂。



欢迎光临 Project1 (https://rpg.blue/) Powered by Discuz! X3.1