习题1.27

先写代码

(defn square [x] (* x x))
(defn expmod[base exp m]
  (cond (= exp 0) 1
        (even? exp) (mod (square (expmod base (/ exp 2) m)) m)
        :else (mod (* base (expmod base (- exp 1) m)) m)
        )
  )



(defn fermat-test[n]
  (defn try-it [a]
    (cond (= a n) (println "test end" n);
          (= (expmod a n n) (mod a n)) (try-it (+ 1 a))
          :else (println "no") 
          )
    )
  (try-it (+ 1 0))
)

执行结果如下

输出结果显示，对于每一个小于561的数字，他的561次方 和他与561的模相同，也就是通过了测试。

 

对于一个非素数，执行结果显示为no，表明没有逃过费曼检查测试。