def solve(): x = np.int64(600851475143) ps = eu.primes(int(np.sqrt(x))) return np.max([p for p in ps if x % p == 0])
def solve(): p = eu.primes(200000) return p[10001 - 1]