n / k の和
A006218(http://oeis.org/A006218/list)
と比較し、答え合わせしてみる。
def s(n)
(1..n).inject(0){|s, k| s + n / k}
end
def A006218(n)
(0..n).map{|i| s(i)}
end
ary = A006218(59)
# OEIS A006218のデータ
ary0 =
[0,1,3,5,8,10,14,16,20,23,27,29,35,37,41,45,50,52,
58,60,66,70,74,76,84,87,91,95,101,103,111,113,119,
123,127,131,140,142,146,150,158,160,168,170,176,
182,186,188,198,201,207,211,217,219,227,231,239,
243,247,249]
# 一致の確認
p ary == ary0
0 件のコメント:
コメントを投稿
注: コメントを投稿できるのは、このブログのメンバーだけです。