日志
快速求因子
|||
(defun factor (n)
(defun f (z) (+ (* z z) 1))
(setq x 1 y x k 0 i 1 d 1)
(while (= d 1)
(setq k (+ 1 k))
(setq x (f x))
(setq d (gcd (- x y) n))
(cond
((and (< d n) (> d 1)) (princ d)(princ))
((= d n) (princ "此数是质数!"))
(t (if (= k i) (progn (setq y x) (setq i (lsh i 1)))))
)
)
)
(defun f (z) (+ (* z z) 1))
(setq x 1 y x k 0 i 1 d 1)
(while (= d 1)
(setq k (+ 1 k))
(setq x (f x))
(setq d (gcd (- x y) n))
(cond
((and (< d n) (> d 1)) (princ d)(princ))
((= d n) (princ "此数是质数!"))
(t (if (= k i) (progn (setq y x) (setq i (lsh i 1)))))
)
)
)
