DouKaku?

1以上200未満の整数nのうち、 exp(pi * sqrt(n))がほとんど整数であるようなnを求めるコードを書いてください。 なお、expは底がeである指数関数 - Wikipedia、 piは円周率、sqrtは平方根です。また「ほとんど整数である」とは 整数からプラスマイナス0.0001の範囲にあることとします。

Pythonで34行のスクリプトを書いて得られた出力の例が下のようになります。

37 199148647.999978
58 24591257752.000000
67 147197952743.999999
163 262537412640768744.000000 

この問題は光成さんに教えて頂いた e^{π*sqrt{163}}≒26253741640768744 が元になっています。ご協力ありがとうございました。

ext:long-float-digits の存在を知ったので使ってみたくなりまして。

1
2
3
4
5
6

(setf (ext:long-float-digits) 192)
(loop
  for x from 1 below 200
  as y = (exp (* pi (sqrt (float x 1L0))))
  if (< (abs (- y (round y))) 1e-4)
  do (format t "~D ~F~%" x y))