150816(2)

Ruby


円周率

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

require 'bigdecimal/math'
include BigMath

def A000796(n)
  m = n * 2
  BigDecimal(PI(m).to_s).to_s("F").scan(/\d/).map(&:to_i)[0..n + 1]
end
ary = A000796(103)

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