迁移C＃到Python - Random类

Question

我需要一些C＃代码移植到Python。原始代码使用Random类。迁移代码必须是周期准确的（即连续调用Next()必须产生在两个码相同的结果）。一些问题：

• Python中是否存在与C＃的Random相当的调用？
• Alliteratively，假设我可以修改这两个来源，有没有与C＃和Python的作品伪随机库？

Answer 1

我不知道任何可用于Python和C＃的库，它们都为这两个库生成相同的随机数。但是，您可能能够利用IronPython。默认Random实现IronPython和CPython中之间是不同的，但WichmannHill类没有。

您可以使用C＃来实例化的IronPython的WichmannHill类并获得尽可能CPython中相同的值相同的种子。或者，您可以通过翻译random.py中的Python代码，相对容易地在C＃中实现Wichmann-Hill算法。

另一种选择是采取Random的Mersenne Twister算法的CPython implementation，并将其转换为C＃以获得相同的结果。

Answer 2

我知道这是一个老问题，但我最终需要一个解决的办法。我最终在Python中实现了C＃的Random类。只要你不需要大于2147483647的随机数，它就会工作，我最终不需要那个功能，所以我没有实现它。

https://gist.github.com/BadStreff/541cf2e6953b3c666f83127a1d4f6a47

from ctypes import * 
# implemented from: 
# http://referencesource.microsoft.com/#mscorlib/system/random.cs,dec894a7e816e665 
class Random(object): 
    def __init__(self, seed): 
     self.seed = c_int(seed).value 
     self.MBIG = 2147483647 
     self.MMIN = -2147483648 
     self.MZ = 0 
     self.MSEED = 161803398 
     self.SeedArray = [0] * 56 

     if seed == self.MMIN: 
      subtraction = self.MBIG 
     else: 
      subtraction = abs(seed) 

     mj = c_int(self.MSEED - subtraction).value 
     self.SeedArray[55] = mj 
     mk = 1 
     for i in range(1, 55): 
      ii = (21 * i) % 55 
      self.SeedArray[ii] = mk 
      mk = mj - mk 
      if mk < 0: 
       mk += self.MBIG 
      mj = self.SeedArray[ii] 
     for k in range(1, 5): 
      for i in range(1, 56): 
       self.SeedArray[i] -= self.SeedArray[1 + (i + 30) % 55] 
       if self.SeedArray[i] < 0: 
        self.SeedArray[i] = c_int(self.SeedArray[i] + self.MBIG).value 
     self.inext = 0 
     self.inextp = 21 
     self.seed = 1 

    def InternalSample(self): 
     locINext = self.inext + 1 
     locINextp = self.inextp + 1 

     if locINext >= 56: 
      locINext = 1 
     if locINextp >= 56: 
      locINextp = 1 

     retVal = c_int(self.SeedArray[locINext] - self.SeedArray[locINextp]).value 
     if retVal == self.MBIG: 
      retVal -= 1 
     if retVal < 0: 
      retVal = c_int(retVal + self.MBIG).value 
     self.SeedArray[locINext] = retVal 
     self.inext = locINext 
     self.inextp = locINextp 
     return retVal 

    def Next(self, minValue=None, maxValue=None): 
     if minValue == None: 
      return self.InternalSample() 
     valRange = maxValue - minValue 
     if valRange <= self.MBIG: 
      return int(c_float(self.Sample() * valRange).value) + minValue 
     else: 
      return self.GetSampleForLargeRange() * valRange + minValue 

    def GetSampleRangeForLargeRange(self): 
     pass 

    def Sample(self): 
     s = self.InternalSample() 
     ret = c_double(s * c_double(1.0/self.MBIG).value).value 
     # print(f'sample: {s}\nret: {ret}') 
     return ret