150328

Ruby


Perrin Pseudoprime

1〜10＾6までで調べてみる。

require 'prime'

a, b, c = 2, 0, 3

i = 2

while i < 10 ** 6

  a, b, c = b + c, a, b

  i += 1

  if a % i == 0

    p i if !i.prime?

  end

end

出力結果

271441

904631