from math import sqrt
a=[2]
y=3
def more_primes(x):
for n in xrange(y, x, 2):
is_prime = 1
for b in a:
if n%b==0:
is_prime = 0
break
if b*b > n:
break
if is_prime:
a.append(n)
print n,
print ""
if x%2==0:
return (x+1)
else:
return (x+2)
y=more_primes(40) ; print(y)
y=more_primes(42) ; print(y)
print(a)
y=more_primes(12) ; print(y)
y=more_primes(42) ; print(y)
ZnJvbSBtYXRoIGltcG9ydCBzcXJ0CgphPVsyXQp5PTMKCmRlZiBtb3JlX3ByaW1lcyh4KToKICAgIGZvciBuIGluIHhyYW5nZSh5LCB4LCAyKToKICAgICAgICAgICAgaXNfcHJpbWUgPSAxCiAgICAgICAgICAgIGZvciBiIGluIGE6CiAgICAgICAgICAgICAgICBpZiBuJWI9PTA6CiAgICAgICAgICAgICAgICAgICAgaXNfcHJpbWUgPSAwCiAgICAgICAgICAgICAgICAgICAgYnJlYWsKICAgICAgICAgICAgICAgIGlmIGIqYiA+IG46CiAgICAgICAgICAgICAgICAgICAgYnJlYWsKICAgICAgICAgICAgaWYgaXNfcHJpbWU6CiAgICAgICAgICAgICAgICBhLmFwcGVuZChuKQogICAgICAgICAgICAgICAgcHJpbnQgbiwKICAgIHByaW50ICIiCiAgICBpZiB4JTI9PTA6CiAgICAgICAgcmV0dXJuICh4KzEpCiAgICBlbHNlOgogICAgICAgIHJldHVybiAoeCsyKQogICAgICAgIAp5PW1vcmVfcHJpbWVzKDQwKSAgICAgOyBwcmludCh5KSAKeT1tb3JlX3ByaW1lcyg0MikgICAgIDsgcHJpbnQoeSkKcHJpbnQoYSkKCnk9bW9yZV9wcmltZXMoMTIpICAgICA7IHByaW50KHkpCnk9bW9yZV9wcmltZXMoNDIpICAgICA7IHByaW50KHkp
3 5 7 11 13 17 19 23 29 31 37
41
41
43
[2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41]
13
13 17 19 23 29 31 37 41
43