151231(3)

Ruby


p(n^4 | 和因子は四乗数)

151109(4)分でp(n | 和因子は四乗数) を求めたが、

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

def A259793(n)
  return [1] if n == 0
  ary = (1..n).map{|i| 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 = A259793(23)

# OEIS A259793のデータ
ary0 =
[1,1,2,7,36,253,1886,14800,118238,955639,7750456,
 62777522,506272363,4056634991,32252971687,
 254209569990,1985108901344,15352968310930,
 117579612410477,891596419221856,6694250497509934,
 49768995849050468,366423320400440927,
 2671969175372760210]
# 一致の確認
p ary == ary0