2013-08-04 62 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