150803

Ruby


4p - 1 型のヘーグナー数の性質

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

include Math

@heegner = [7, 11, 19, 43, 67, 163]

# ほとんど整数であることの確認

p @heegner.map{|i| (E ** (PI * sqrt(i)) - 744) ** (1.0 / 3)}

def A199743()

  @heegner.map{|i| ((E ** (PI * sqrt(i)) - 744) ** (1.0 / 3)).round}

end

ary = A199743()

# OEIS A199743のデータ

ary0 = [15,32,96,960,5280,640320]

# 一致の確認

p ary == ary0

出力結果

[14.929943081975248, 31.998093332227423, 95.99999195891692, 959.9999999999179, 5

279.999999999993, 640319.9999999987]

true