from functools import reduce

primes = []

i = 2

for i in range(2, 2000000):
    half = i / 2
    if all(p > half or i % p != 0 for p in primes):
        primes.append(i)

result = reduce(lambda x, y: x + y, primes)
print(primes)
print(result)

# answer 142913828922