2013-08-04 28 views
1

這是我製作的主要發電機。我應該能夠產生最多20000的素數。它應該生成提供給它的素數作爲參數。但是,它只能達到11並停止D :.任何人都可以解釋這裏有什麼問題?爲什麼我的素數發生器發生故障?

def find_primes(limit): 
    prime_holder = [2, 3, 5 ,7] 
    divided_pass = 0 
    for i in range(11, 20000): 
     for j in range(0, len(prime_holder)): 
      if i%prime_holder[j] != 0: 
       divided_pass += 1 
     if divided_pass == len(prime_holder): 
      prime_holder.append(i) 
      divided_pass = 0 
     if len(prime_holder)-1 == limit: 
      break 
    return prime_holder 

my_primes = find_primes(50) 
for x in my_primes: 
    print x; 
raw_input() 
+1

如果在該範圍內產生的所有質數,你應該使用[埃拉托色尼的篩(HTTP:/ /en.wikipedia.org/wiki/Sieve_of_Eratosthenes)。 –

+0

我製作了基於SoE的腳本。但是這個劇本是我第一個創作的劇本。這不起作用,我不明白爲什麼。 – Stormboy

回答

2

您需要設置divided_pass0爲每一個新i

def find_primes(limit): 
    prime_holder = [2, 3, 5 ,7] 
    for i in range(11, 20000): 
     divided_pass = 0 
     for j in range(0, len(prime_holder)): 
      if i%prime_holder[j] != 0: 
       divided_pass += 1 
     if divided_pass == len(prime_holder): 
      prime_holder.append(i) 
     if len(prime_holder)-1 == limit: 
      break 
    return prime_holder 

>>> print find_primes(50) 
[2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233] 

P.S - 列表是迭代在Python,所以你不需要做for j in range(len(...))和可以做for j in prime_holder: if i % j != 0。在給定範圍內尋找素數有更好的實現(你應該研究這些)。

您也可以使用for-else結構,其中else零件在沒有遇到break時執行。現在你的代碼被減少到了。

def find_primes(limit): 
    prime_holder = [2, 3, 5 ,7] 
    for i in range(11, 20000): 
     for j in prime_holder: 
      if i%j == 0: 
       break 
     else: 
      prime_holder.append(i) 
     if len(prime_holder)-1 == limit: 
      break 
    return prime_holder 
+2

「for-else」的完美用法。 – nneonneo

0

你應該清楚每個週期divided_pa​​ss的價值:

def find_primes(limit): 
    prime_holder = [2, 3, 5 ,7] 
    divided_pass = 0 
    for i in range(11, 20000): 
     for j in range(0, len(prime_holder)): 
      if i % prime_holder[j] != 0: 
       divided_pass += 1 
     if divided_pass == len(prime_holder): 
      prime_holder.append(i) 
     divided_pass = 0 
     if len(prime_holder)-1 == limit: 
      break 
    return prime_holderdivided_pass