Jacobi符号计算的一种实现

来源:互联网 发布:苹果cms采集插件 编辑:程序博客网 时间:2024/06/07 15:55

        Jacobi符号是Legendre符号的推广,在素性检验以及因子分解等诸多方面极为有用,所以先写一个备着,写其他算法时方便直接调用。


(defun jacobi (n m)
  (let ((e 0)
    (a1 0)(a n)
    (s 0)(b m)
    (n1 0))
    (cond ((= a 0) 0)
      ((= a 1) 1)
      (t
       (if (/= 0 (mod a 2)) (progn (setf a1 a) (setf e 0))
         (loop while (= 0 (mod a 2)) do
         (setf a (/ a 2)) (setf a1 a) (incf e)))
       (cond ((= 0 (mod e 2)) (setf s 1))
         ((or (= 1 (mod b 8)) (= 7 (mod b 8)))
          (setf s 1))
         ((or (= 3 (mod b 8)) (= 5 (mod b 8)))
          (setf s -1)))
       (if (and (= 3 (mod b 4)) (= 3 (mod a1 4))) (setf s (- 0 s)))
       (setf n1 (mod b a1))
       (if (= 1 a1)
           s
         (* s (jacobi n1 a1)))))))

        

原创粉丝点击