160517

Ruby

Goldbach's conjecture(2)

ペアの数を数えてみた。、
オンライン整数列大辞典の
A061358（http://oeis.org/A061358/list）
と比較し、答え合わせしてみる。

require'prime'

def A061358(n)
  (0..n).map{|i| Prime.each(i / 2).count{|j| (i - j).prime?}}
end
ary = A061358(105)

# OEIS A061358のデータ
ary0 =
[0,0,0,0,1,1,1,1,1,1,2,0,1,1,2,1,2,0,2,1,2,1,3,0,
 3,1,3,0,2,0,3,1,2,1,4,0,4,0,2,1,3,0,4,1,3,1,4,0,5,
 1,4,0,3,0,5,1,3,0,4,0,6,1,3,1,5,0,6,0,2,1,5,0,6,1,
 5,1,5,0,7,0,4,1,5,0,8,1,5,0,4,0,9,1,4,0,5,0,7,0,3,
 1,6,0,8,1,5,1]
# 一致の確認
p ary == ary0