Archive for 2015年9月|Monthly archive page

projecteuler25

 newLISP の length を数値に適用すると、整数部分の桁数が返ります。
 ということで、projecteuler問題25は、フィボナッチ数が 1000桁になるのは何番目かという問題。
 1000桁ですから大整数の出番。そして、前述のlengthを使えば、

(silent)
(setq fibo  '(1L 1L))
(while (< (length (fibo -1)) 1000)
  (push (+ (fibo -2) (fibo -1)) fibo -1))
(println (length fibo))

 これを実行すれば

> 
4782

 と答えが得られます。
 silent は余計な表示を抑えるための関数。なくても答えは得られます。
 リストの長さが求める答えになっているので、スクリプトの最後にも length を使っています。
 フィボナッチ数列は変数 fibo に入っていますから、中身を覗くと

> (0 12 fibo)
(1L 1L 2L 3L 5L 8L 13L 21L 34L 55L 89L 144L)
> (fibo -1)
1070066266382758936764980584457396885083683896632151665013235203375314520604694040621889147582489792657804694888177591957484336466672569959512996030461262748092482186144069433051234774442750273781753087579391666192149259186759553966422837148943113074699503439547001985432609723067290192870526447243726117715821825548491120525013201478612965931381792235559657452039506137551467837543229119602129934048260706175397706847068202895486902666185435124521900369480641357447470911707619766945691070098024393439617474103736912503231365532164773697023167755051595173518460579954919410967778373229665796581646513903488154256310184224190259846088000110186255550245493937113651657039447629584714548523425950428582425306083544435428212611008992863795048006894330309773217834864543113205765659868456288616808718693835297350643986297640660000723562917905207051164077614812491885830945940566688339109350944456576357666151619317753792891661581327159616877487983821820492520348473874384736771934512787029218636250627816L
> (length (fibo -1))
1000
> 

 こんな感じ。
 末尾に L が付いているのは、大整数の証。
 以上、如何でしょうか?