-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathproblem7.py
50 lines (42 loc) · 1.13 KB
/
problem7.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
import sys
import math
from datetime import datetime
if len(sys.argv) < 3:
print("Error: missing argument(s)\nusage: 'python",sys.argv[0]," $index, $step'")
sys.exit(0)
try:
index = int(sys.argv[1])-1
step = int(sys.argv[2])
except ValueError:
print("Arguments must be decimal integer values, instead received '%s' and '%s'" % (sys.argv[1],sys.argv[2]))
sys.exit(0)
primes = [2,3,5,7]
def generatePrimesSmallerThan(limit):
print("Generating all primes smaller than %8d " % limit, end='')
time = datetime.today()
numbers = list(range(limit-step,limit))
for prime in primes:
i=0
while i < len(numbers):
if numbers[i]%prime==0:
numbers.pop(i)
else:
i+=1
while 0<len(numbers):
primes.append(numbers.pop(0))
if primes[-1]==1 :
primes.pop(-1)
continue
toRemove = primes[-1]
while(toRemove<limit):
try:
numbers.remove(toRemove)
except:
pass
toRemove += primes[-1]
print("(%.3f s elapsed)" % (datetime.today()-time).total_seconds())
currentLimit=step
while len(primes)<index+1:
generatePrimesSmallerThan(currentLimit)
currentLimit+=step
print("%d is prime number %d." % (primes[index], index+1))