2012-09-18 42 views
0

我使用ECL和newlisp運行面函數。關於newlisp中的階乘

ECL:

>(defun fac (n) (if (= n 1) 1 (* n (fac (- n 1))))) 
>(fac 20) 
22432902008176640000 
>(fac 30) 
2265252859812191058636308480000000 
>(fac 40) 
815915283247897734345611269596115894272000000000 
... 

newlisp

>(define (fac n) (if (= n 1) 1 (* n (fac (- n 1))))) 
>(fac 20) 
22432902008176640000 
>(fac 30) 
-8764578968847253504 

爲什麼newlisp與ECL返回不同的結果呢?

回答

1

newLISP被設計成小(〜260K),並且你會不會在默認情況下得到的事情之一是真正的大整數的支持,以及必要的自動轉換普通大小和大整數之間。如果你經常需要整數,如1000!在你的工作中,你可能選擇了錯誤的語言。但是,如果你想使用newLISP其他原因,和真的需要計算100!然後安裝GMP庫,並編寫代碼如下:

(load "gmp.lsp") 
(define (factorial-gmp num) 
(if (= num 0) 
    "1" 
    (GMP:* (string num) (factorial-gmp (- num 1))))) 
2

newLISP中最大的64位整數是9,223,372,036,854,775,807。

2

嘗試浮點:

(define (fac n) (if (= n 1) 1 (mul n (fac (sub n 1))))) 

(fac 20) => 2.432902008e+18 
+0

(FAC 100)=> 9.332621544e + 157,這是偉大的!但是,newlisp不能計算(fac 1000),但ECL確實! –

+0

順便說一句,爲什麼'+ - * /'不像CL那樣具有多態性? –

+0

該手冊告訴您如何定義這些運算符,以便它們可以使用不精確的FP編號以及精確的整數... :) – cormullion

2

由於10.4.8版本newLISP具有內置支持(bigint)爲無限精度整數。這使得GNU GMP模塊gmp.lsp不再適用。但由於BIGINT的工作方式,你會看到比典型的遞歸一個完全不同的階乘函數:

(define (fac n) (apply * (map bigint (sequence 1 n)) 2)) 
(fac 1000) 

> 
(lambda (n) (apply * (map bigint (sequence 1 n)) 2)) 
402387260077093773543702433923003985719374864210714632543799910429938512398629020592044208486969404800479988610197196058631666872994808558901323829669944590997424504087073759918823627727188732519779505950995276120874975462497043601418278094646496291056393887437886487337119181045825783647849977012476632889835955735432513185323958463075557409114262417474349347553428646576611667797396668820291207379143853719588249808126867838374559731746136085379534524221586593201928090878297308431392844403281231558611036976801357304216168747609675871348312025478589320767169132448426236131412508780208000261683151027341827977704784635868170164365024153691398281264810213092761244896359928705114964975419909342221566832572080821333186116811553615836546984046708975602900950537616475847728421889679646244945160765353408198901385442487984959953319101723355556602139450399736280750137837615307127761926849034352625200015888535147331611702103968175921510907788019393178114194545257223865541461062892187960223838971476088506276862967146674697562911234082439208160153780889893964518263243671616762179168909779911903754031274622289988005195444414282012187361745992642956581746628302955570299024324153181617210465832036786906117260158783520751516284225540265170483304226143974286933061690897968482590125458327168226458066526769958652682272807075781391858178889652208164348344825993266043367660176999612831860788386150279465955131156552036093988180612138558600301435694527224206344631797460594682573103790084024432438465657245014402821885252470935190620929023136493273497565513958720559654228749774011413346962715422845862377387538230483865688976461927383814900140767310446640259899490222221765904339901886018566526485061799702356193897017860040811889729918311021171229845901641921068884387121855646124960798722908519296819372388642614839657382291123125024186649353143970137428531926649875337218940694281434118520158014123344828015051399694290153483077644569099073152433278288269864602789864321139083506217095002597389863554277196742822248757586765752344220207573630569498825087968928162753848863396909959826280956121450994871701244516461260379029309120889086942028510640182154399457156805941872748998094254742173582401063677404595741785160829230135358081840096996372524230560855903700624271243416909004153690105933983835777939410970027753472000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000L
1

從版本10.4.8 newLISP有內置的支持(BIGINT)無限精度的整數。

BIGINT遞歸階乘:

> (define (fac n) (if (= n 1) 1 (* (bigint n) (fac (- n 1))))) 
> (fac 20) 
2432902008176640000L 
> (fac 30) 
265252859812191058636308480000000L 
> (fac 40) 
815915283247897734345611269596115894272000000000L