def problem007():
"""prime number theorem says # primes < n = n / ln n + B, B = 1.0836 --> nPrimes ~ n / ln n
Solving in reverse: e^(-1 * ProductLog(-1, -(1/a))) --> a ~ 116684
"""
n = 10001
r = range(116684)
print r[-1]
l = tools.filter_primes(r)
print l[n]
def problem010c():
l = tools.filter_primes(range(2,limit))
print sum(l)
