Question

我正在研究Euler Project问题，这是问题五：

最大素数问题3 13195的主要因子是5、7、13和29。

600851475143的最大素数是多少？

我得到了工作代码：

def factor(x, f=2):


    while x >= f*f:
        while x % f == 0:
            x = int(x/f)
        factor = f
        f += 1

    print(f'x = {x},\nlast factor = {factor}') # print for debug only 
    return max(x, factor)

factor（19 * 19 * 19 * 19 * 19 * 19 * 19 * 19 * 19 * 1999989899）

x = 33170854034208712，最后一个因子= 182128674

33170854034208712

有人知道为什么这不能产生正确答案吗？

Answer 1

我只是从头开始思考这个问题。下面是我的代码，看起来很像你的

def largest_prime_divisor(x, f=2):
    while f**2 <= x:
        if x % f == 0:
            x //= f
        else:
            f+=1
    return x

print(largest_prime_divisor(25)) # 5
print(largest_prime_divisor(13195))  # 29

n = 19*19*19*19*19*19*19*19*19*1999989899
print(largest_prime_divisor(n))  # 577 this is prime

Answer 2

您应该使用整数除法运算符//而不是浮点除法运算符/

除了您的代码是正确的

            while x % f == 0:
                x = int(x//f)
            factor = f

输出

x = 577,
last factor = 541

除法运算符差异

Python 3.7.1 (default, Nov  6 2018, 18:46:03) 
[Clang 10.0.0 (clang-1000.11.45.5)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> 
>>> 
>>> print(645372136089564734321//19)
33966954531029722859
>>> print(645372136089564734321/19)
3.396695453102972e+19

最大数量的最大素数失败

2 个答案: