151231(2)

Ruby


p(n^3 | 和因子は立方数)

151109(3)分でp(n | 和因子は立方数) を求めたが、

そのn^3 番目について求めてみた。
オンライン整数列大辞典の
A259792（http://oeis.org/A259792/list）
と比較し、答え合わせしてみる。

def A259792(n)
  return [1] if n == 0
  ary = (1..n).map{|i| i * i * i}
  m = ary[-1]
  ps = Array.new(m + 1){0}
  ps[0] = 1
  a_ary = [1]
  ary.each{|num|
    (num..m).each{|i|
      ps[i] += ps[i - num]
    }
    a_ary << ps[num]
  }
  a_ary
end
ary = A259792(26)

# OEIS A259792のデータ
ary0 =
[1,1,2,5,17,62,258,1050,4365,18012,73945,301073,
 1214876,4852899,19187598,75070201,290659230,
 1113785613,4224773811,15866483556,59011553910,
 217410395916,793635925091,2871246090593,
 10297627606547,36620869115355,129166280330900]
# 一致の確認
p ary == ary0