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]