151227(5)

Ruby

Stern's diatomic series(2)

第10 ** i 項を求めてみた。

def A002487(n)
  a, b = 0, 1
  ary = [a, b]
  i = 1
  while i < n
    a, b = b, b + a - 2 * (a % b)
    ary << b
    i += 1
  end
  ary[0..n]
end

N = 8
ary = A002487(10 ** N)
(0..N).each{|i| p ary[10 ** i]}

出力結果
1
3
7
11
43
127
191
1249
897