大数上的最大素因子失败

Question

我在研究Euler项目的问题，这是第五个问题：

最大素因子 问题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)

系数（19*19*19*19*19*19*19*19*19*19*19999899）

x=33170854034208712， 最后一个系数=182128674

33170854034208712

有人知道为什么这不能给出正确的答案吗？你知道吗

Answer 1

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

除此之外，你的代码是正确的

            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

Answer 2

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

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