首页 > 解决方案 > 简单地将 1 加到一个大数上不起作用?(Python 3.9)

问题描述

注意:我在 Python 方面没有那么丰富的经验,因此我的代码可能没有它应该/应该的那么好。

我正在尝试创建一个工具来帮助计算某种形式的数字的代数因子(参见https://en.wikipedia.org/wiki/Aurifeuillean_factorization)。这主要是作为测试/学习经验,但是在尝试计算参数“c”时遇到了问题,该参数定义为 2^(2k+1)+1。添加步骤对我不起作用。我只是将返回值作为 2^129,而不是我希望得到的 2^129+1。这是 Python 本身的问题,还是我在这方面犯了某种错误。

代码:

import math


def make_aurifeuille_factors(base, exponent):
    if base == 2 and exponent % 4 == 2:
        k = (exponent - 2) / 4
        c = int(1 + 2 ** (2*k + 1))
        d = int(2 ** (k + 1))
        L = c + d
        M = c - d

        return int(k), int(c), int(d), int(L), int(M)


def gcd(a, b):
    return int(math.gcd(a, b))


print(make_aurifeuille_factors(2, 258))

标签: pythonprime-factoring

解决方案


k = (exponent - 2) / 4k生成 a ,这float意味着您可能会在计算中引入数值错误。使用整数除法从一开始就留在int世界中:

def make_aurifeuille_factors(base, exponent):
    if base == 2 and exponent % 4 == 2:
        k = (exponent - 2) // 4
        c = 1 + 2 ** (2*k + 1)
        d = 2 ** (k + 1)
        L = c + d
        M = c - d

        return k, c, d, L, M 

推荐阅读