久々にプログラムです.
ラマヌジャンのτ関数の計数をある程度計算してみようと思います.
とはいっても今回は難しいことはしないです.
内容としては,
- どこまでの値がほしいかを入力します.
- 24乗の係数をパスカルの三角形から計算します.せいぜい24までなのでこの方法でも十分速いと思います.
- 各に対しとすると,動的計画法チックに解くことができます.
- 指数値が最初の入力値を超えた場合は計算を終了し,次のnについて計算します.
- nが最初の値を超えた場合は終了し,結果を出力します.
という流れ.
一応途中で打ち切りはしていますが,若干効率悪めなので1000程度では高速ですが,10000とかになると若干厳しいものがあります.
プログラム自体はここ
github.com
1000まで出力した結果はここ
github.com
です.