def compute():
i = 647
limit = 1000000
run = 0
i_run = 0
primes = set(prime_sieve(limit))
while run < 4:
factors = list_divisors(i)
factors = [f for f in factors if (f % 2 == 1 or f == 2) and f in primes]
if len(factors) == 4:
if run == 0:
i_run = i
run += 1
else:
run = 0
i += 1
return i_run
def compute():
primes = set(prime_sieve(100000))
rad = sorted([(i, prod([y for y in list_divisors(i) if y in primes])) for i in range(1, 10**5 + 1)],
key=lambda x: x[1])
return rad[9999][0]
