在python或它的包装器中是否有任何drand48()等效?

时间:2011-09-02 17:17:40

标签: python random

我正在将一些c代码转换为python。我想知道在python中是否有任何drand48()等效(完全相同!)或它的包装器?

非常感谢

1 个答案:

答案 0 :(得分:7)

您正在寻找random.random

>>> import random
>>> help(random)
>>> random.random()
0.8423026866867628

它不使用与drand48相同的生成器。 drand48系列使用48位线性同余生成器,而Python随机模块使用优越的Mersenne Twister算法。

如果您希望drand48完全相同的输出,则可以在Python中实现它。

# Uncomment the next line if using Python 2.x...
# from __future__ import division
class Rand48(object):
    def __init__(self, seed):
        self.n = seed
    def seed(self, seed):
        self.n = seed
    def srand(self, seed):
        self.n = (seed << 16) + 0x330e
    def next(self):
        self.n = (25214903917 * self.n + 11) & (2**48 - 1)
        return self.n
    def drand(self):
        return self.next() / 2**48
    def lrand(self):
        return self.next() >> 17
    def mrand(self):
        n = self.next() >> 16
        if n & (1 << 31):
            n -= 1 << 32
        return n   

但是,输出将远远低于Python的标准随机数生成器。 SVID 3在1989年宣布rand48系列函数已过时。