python - 将蟒蛇长成四肢

Question

我正在尝试在 Python 中编写蒙哥马利乘法的实现，我需要一个等效于 GMP 的mpz_getlimbn()Python longs 但我似乎一辈子都找不到。

任何帮助将不胜感激。

编辑

我已经实现了以下内容，但我得到了 GMP 中没有出现的肢体的索引超出范围错误。

def unpack(x, b):
    if gmpy2:
        return [long(x) for x in gmpy2.unpack(gmpy2.mpz(x), b)]

    b = 2 ** b
    r = []
    while x:
        x, temp = divmod(x, b)
        r.append(temp)
    return r

Answer 1

我修改了你的unpack()，它似乎对我有用。如果您仍然收到错误，请发布完整的错误。

>>> import gmpy2
>>> 
>>> def unpack(x, b):
...     try:
...         return [x for x in gmpy2.unpack(gmpy2.mpz(x), b)]
...     except NameError:
...         b = 2 ** b
...         r = []
...         while x:
...             x, temp = divmod(x, b)
...             r.append(temp)
...         return r
... 
>>> unpack(123456**7, 15)
[mpz(0), mpz(0), mpz(4096), mpz(25855), mpz(24508), mpz(31925), mpz(15111), mpz(10775)]
>>> del(gmpy2)
>>> unpack(123456**7, 15)
[0, 0, 4096, 25855, 24508, 31925, 15111, 10775]

使用时gmpy2，我将结果mpz显示gmpy2为已使用。

Python 的长整数类型使用存储 15 位或 30 位的肢体。sys.int_info将为您的系统提供详细信息。

顺便说一句，我坚持gmpy2，很高兴看到有人使用unpack().