Fixed para.lsp primep primep*

sasagawa888 · May 22, 2024 · 574cdbd · 574cdbd
1 parent 41a8012
commit 574cdbd
Showing 1 changed file with 11 additions and 10 deletions.
diff --git a/tests/para.lsp b/tests/para.lsp
@@ -3,15 +3,15 @@
 ;(mp-create 5)
 ;(primep* 100000000000031)
 (defun primep* (n)
-    (if (= (mod n 2) 0)
-        nil
-        (let* ((limit (isqrt n))
-               (span (div limit 5)))
-            (mp-part (coprimep n 3 span)
-                     (coprimep n (near-odd span) (* 2 span ))
-                     (coprimep n (near-odd (* 2 span)) (* 3 span))
-                     (coprimep n (near-odd (* 3 span)) (* 4 span))
-                     (coprimep n (near-odd (* 4 span)) limit)))))
+    (cond ((= n 2) t)
+          ((= (mod n 2) 0) nil)
+          (t (let* ((limit (isqrt n))
+                    (span (div limit 5)))
+                 (mp-part (coprimep n 3 span)
+                      (coprimep n (near-odd span) (* 2 span ))
+                      (coprimep n (near-odd (* 2 span)) (* 3 span))
+                      (coprimep n (near-odd (* 3 span)) (* 4 span))
+                      (coprimep n (near-odd (* 4 span)) limit))))))
 
 (defun near-odd (n)
     (if (= (mod n 2) 0)
@@ -25,6 +25,7 @@
           (t (coprimep n (+ s 2) e))))
 
 (defun primep (n)
-    (cond ((= (mod n 2) 0) nil)
+    (cond ((= n 2) t)
+          ((= (mod n 2) 0) nil)
           (t (coprimep n 3 (isqrt n)))))