フェルマーテストで素数判定ははやくなる。
(use srfi-19) (use srfi-27) (define (even? n) (= (remainder n 2) 0)) (define (square x) (* x x)) (define (fermat-test n) (define (try-it a) (= (expmod a n n) a)) (try-it (+ 1 (random-integer (- n 1))))) (define (fast-prime? n times) (cond ((= times 0) #t) ((fermat-test n) (fast-prime? n (- times 1))) (else #f))) (define (expmod base exp m) (cond ((= exp 0) 1) ((even? exp) (remainder (square (expmod base (/ exp 2) m)) m)) (else (remainder (* base (expmod base (- exp 1) m)) m)))) (define (timed-prime-test n) (newline) (display n) (start-prime-test n (current-time))) (define (start-prime-test n called-time) (if (fast-prime? n 100) (report-prime (time-difference (current-time) called-time)) #f)) (define (report-prime elapsed-time) (display " *** ") (display elapsed-time) #t )