3

我使用下面的 Stream 和惰性 val 编写了惰性筛算法的以下实现:

def primes(): Stream[Int] = {
   lazy val ps = 2 #:: sieve(3)
   def sieve(p: Int): Stream[Int] = {
       p #:: sieve(
            Stream.from(p + 2, 2).
             find(i=> ps.takeWhile(j => j * j <= i).
                     forall(i % _ > 0)).get)
  }
  ps
}

以及使用(可变)ListBuffer的以下实现:

import scala.collection.mutable.ListBuffer
def primes(): Stream[Int] = {
    def sieve(p: Int, ps: ListBuffer[Int]): Stream[Int] = {
        p #:: { val nextprime =
            Stream.from(p + 2, 2).
            find(i=> ps.takeWhile(j => j * j <= i).
                 forall(i % _ > 0)).get
            sieve(nextprime, ps += nextprime)
         }
    }       
    sieve(3, ListBuffer(3))}

当我做 primes().takeWhile(_ < 1000000).size 时,第一个实现比第二个快 3 倍。对此有何解释?

我编辑了第二个版本:它应该是 sieve(3, ListBuffer(3)) 而不是 sieve(3, ListBuffer()) 。

4

2 回答 2

6

好吧,我的猜测是这一行:

find(i=> ps.takeWhile(j => j * j <= i).forall(i % _ > 0)).get

On ListBuffertakeWhile创建一个临时集合(它越来越大)。同时,Stream由于其非严格性,避免这样做。一旦forall失败,它就会停止计算takeWhile.

于 2010-12-26T17:05:12.903 回答
2

没有真正回答这个问题,但因为我花了一些时间对各种组合进行基准测试......

如果在内部循环中使用并避免使用Iterator,您可以获得更好的性能,以最小化内存分配。ArrayBuffertakeWhile

def primes2(): Stream[Int] = {
  def sieve(p: Int, ps: ArrayBuffer[Int]): Stream[Int] = {
    def hasNoDivisor(prime_? :Int, j: Int = 0): Boolean = {
      val n = ps(j)
      if (n*n > prime_?) true
      else if (prime_? % n == 0) false else hasNoDivisor(prime_?, j+1)
    }
    p #:: { 
      val nextprime = Iterator.from(ps.last + 2, 2).find(hasNoDivisor(_)).get
      sieve(nextprime, ps += nextprime)
    }
  }     
  sieve(3, ArrayBuffer(3))
}

这是一个带有Iterator代替的版本Stream,它更快,primes3().toStream如果需要,您可以随时使用来获取 Stream。

def primes3() = List(2,3).iterator ++ new Iterator[Int] {
  val ps = ArrayBuffer[Int](3)
  def hasNoDivisor(prime_? :Int, j: Int = 0): Boolean = {
    val n = ps(j)
    if (n*n > prime_?) true
    else if (prime_? % n == 0) false else hasNoDivisor(prime_?, j+1)
  }
  def hasNext = true
  def next() = {
    val nextprime = Iterator.from(ps.last + 2, 2).find(hasNoDivisor(_)).get
    ps += nextprime
    nextprime
  }
}

结果:

primes : warming...
primes : running...
primes : elapsed: 3.711
res39: Int = 283145
primes2: warming...
primes2: running...
primes2: elapsed: 1.039
res40: Int = 283145
primes3: warming...
primes3: running...
primes3: elapsed: 0.530
res41: Int = 283146

我还尝试用几个循环替换from,find和,这样更快,但更难理解。hasNoDivisorwhile

于 2010-12-27T06:21:44.783 回答