
Я хотел найти простые просты Eratosthenes подход. < /p>
Я написал реализацию этого в Python. Но это ужасно медленно. С, скажем, если я хочу найти все простые числа менее 2 миллионов. Это занимает> 20 минут. (Я остановил это в этот момент). Как я могу это ускорить? < /P>
def primes_sieve(limit):
limitn = limit+1
primes = range(2, limitn)
for i in primes:
factors = range(i, limitn, i)
for f in factors[1:]:
if f in primes:
primes.remove(f)
return primes
print primes_sieve(2000)
< /code>
Обновление: < /strong>
Я закончил профилировать в этом коде и обнаружил, что довольно много времени было потрачено на удаление элемента из списка. Вполне понятно, учитывая, что он должен пройти весь список (худший случай), чтобы найти элемент, а затем удалить его, а затем продумать список (может быть, какая-то копия продолжается?). В любом случае, я выбросил список для словаря. Моя новая реализация - < /p>
def primes_sieve1(limit):
limitn = limit+1
primes = dict()
for i in range(2, limitn): primes = True
for i in primes:
factors = range(i,limitn, i)
for f in factors[1:]:
primes[f] = False
return [i for i in primes if primes==True]
print primes_sieve1(2000000)
Подробнее здесь: https://stackoverflow.com/questions/393 ... mes-python