隣り合う素数の和
オンライン整数列大辞典の
A164653(http://oeis.org/A164653/list)
と比較し、答え合わせしてみる。
require 'prime'
def A164653(n)
ary = [1, 3]
ary += Prime.each(n).each_cons(2).map{|c| c[0] + c[1]}
end
ary = A164653(508 / 2 + 3)
# OEIS A164653のデータ
ary0 =
[1,3,5,8,12,18,24,30,36,42,52,60,68,78,84,90,100,
112,120,128,138,144,152,162,172,186,198,204,210,
216,222,240,258,268,276,288,300,308,320,330,340,
352,360,372,384,390,396,410,434,450,456,462,472,
480,492,508]
# 一致の確認
p ary == ary0
0 件のコメント:
コメントを投稿
注: コメントを投稿できるのは、このブログのメンバーだけです。