2016-02-19 81 views
1

我想做一个没有利用明显的数学攻击的筛子。我想强暴它。我的算法是基于这样一个概念,即筛选器对很多不是素数的检查进行了很多检查,只是返回操作的结果来检查这些数据,而不是找出素数。我认为一些卡迈克尔数字证明它对于非常大的东西是无效的。我可能是错的。我继续检查范围内的数字,并遵循从Wikipedia给出的基本算法。Eratosthenes变体的筛选

def primes(n) 
    nums = (2..n) 
    not_prime = [] 
    primes = [] 
    nums.to_a.each_with_index do |it, idx| 
     primes << it unless not_prime.include?(it) 
     p primes.last 
     p nums.step(primes.last).to_a 
     nums.step(primes.last).each_with_index do |num, idx| 
     next if idx == 0 
     not_prime << num 
     end 
    end 

    primes 
end 

当我的范围不会将一行:

nums.step(primes.last).each_with_index 

比第一个(2)等数字,我得到关闭的情况-X错误(配合名单,我相信上)。例如,找到所有非素数两个倍数,但是对于三的倍数,范围上的步骤返回2,5,811等,它们都是1。

我想弄清楚使用Range对象或转换为Array的解决方案,但我喜欢我的(错误的)解决方案的简洁性。有人认为他们可以帮我解决这个问题吗?

编辑:

我修好了!该解决方案创建了一个全新的范围来迭代,而不是采用原始范围。见下文。向JörgW Mittag大声呼喊,鼓励创建一个新的范围,而不是试图摆弄我正在尝试做的原始不变对象。有时候圆孔中的方形钉听起来好多了。

def primes(n) 
    nums = (2..n) 
    not_prime = [] 
    primes = [] 
    nums.to_a.each_with_index do |it, idx| 
     next if not_prime.include?(it) 
     primes << it 
     ((primes.last)..n).step(primes.last).each_with_index do |num, idx| 
     next if idx == 0 || num == primes.last 
     not_prime << num 
     end 
    end 

    primes 
end 
+0

另一个解决方案是只创建自己的自定义数据结构的房子我的Range对象并执行操作,但我喜欢用stdlib的项目和原语的人工挑战...... – noname

+1

您遍历范围'2..n'在步骤'3'中,除了'2,5,8,...'之外还有什么可能*? –

+0

我明白这是错的。我该如何解决这个错误? – noname

回答

0
def primes(n) 
    nums = (2..n) 
    not_prime = [] 
    primes = [] 
    nums.to_a.each_with_index do |it, idx| 
     next if not_prime.include?(it) 
     primes << it 
     ((primes.last)..n).step(primes.last).each_with_index do |num, idx| 
     next if idx == 0 || num == primes.last 
     not_prime << num 
     end 
    end 

    primes 
end