151227(2)

Ruby

離散力学系における軌道の計算(2)

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

def f(lambda, x)
  lambda * x * (1 - x)
end

def A079271(n)
  x = 1 / 3.to_r
  ([x] + (1..n).map{x = f(4, x)}).map{|i| i.numerator}
end
ary = A079271(7)

# OEIS A079271のデータ
ary0 =
[1,8,32,6272,7250432,1038154236987392,
 3383826162019367796397224108032,
 674838593766753484487654913831820720085359667709963001167872]
# 一致の確認
p ary == ary0