150726

Ruby


Prime numbers of the form 33…331

掲題について興味深いサイトを見つけた。

（http://stdkmd.com/nrr/3/33331.htm）

なお、オンライン整数列大辞典にも
A051200（http://oeis.org/A051200/list）
があるが、初項が異なる。

最初の8項なら、数分で求まる。

require 'prime'

m = 18

(0..m).each{|i|

  n = (10 ** i - 7) / 3

  p [i, n] if n.prime?

}

出力結果

[2, 31]

[3, 331]

[4, 3331]

[5, 33331]

[6, 333331]

[7, 3333331]

[8, 33333331]

[18, 333333333333333331]