import math from sympy.ntheory import isprime def primefactor(x): res = 0 for i in range(2, math.ceil(math.sqrt(x))): if not isprime(i) or x % i != 0: continue res = i return res print(primefactor(600851475143))