結果

問題 No.1140 EXPotentiaLLL!
コンテスト
ユーザー Common Lisp
提出日時 2024-11-03 02:29:39
言語 Common Lisp
(sbcl 2.6.3)
コンパイル:
sbclc _filename_
実行:
sbcl --script Main.fasl
結果
AC  
実行時間 1,642 ms / 2,000 ms
コード長 2,157 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 473 ms
コンパイル使用メモリ 45,728 KB
実行使用メモリ 71,312 KB
最終ジャッジ日時 2026-05-06 02:00:38
合計ジャッジ時間 16,811 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 12
権限があれば一括ダウンロードができます
コンパイルメッセージ
; compiling file "/home/judge/data/code/Main.lisp" (written 06 MAY 2026 02:00:21 AM):

; wrote /home/judge/data/code/Main.fasl
; compilation finished in 0:00:00.060

ソースコード

diff #
raw source code 

(declaim (optimize (speed 3) (safety 0)))
(defun sieve-of-atkin (lim)
  (declare (type fixnum lim))
  (let* ((slim (sqrt lim))
         (p (make-array (+ lim 60) :element-type 'fixnum :initial-element 0)))
    
    (declare (type short-float slim))
    (loop for v from 1 to (floor slim 2) do
          (let ((x (1+ (* 4 v v)))
                (y 8))
            (declare (type fixnum x))
            (declare (type fixnum y))
            (loop while (<= x lim) do
                  (when (not (= (mod x 12) 9))
                        (setf (aref p x) (logxor (aref p x) 1)))
                  (setf x (+ x y)
                        y (+ y 8)))))
    (loop for v from 1 to (floor (/ slim (sqrt 3))) by 2 do
          (let ((x (+ (* 3 v v) 4))
                (y 12))
            (declare (type fixnum x))
            (declare (type fixnum y))
            (loop while (<= x lim) do
                  (when (= (mod x 12) 7)
                        (setf (aref p x) (logxor (aref p x) 1)))
                  (setf x (+ x y)
                        y (+ y 8)))))
    (loop for v from 2 to (floor (/ slim (sqrt 2))) do
          (let ((x (1- (* 2 v (1+ v))))
                (y (- (* 4 v) 8)))
            (loop while (and (<= x lim) (<= 0 y)) do
                  (when (= (mod x 12) 11)
                        (setf (aref p x) (logxor (aref p x) 1)))
                  (setf x (+ x y)
                        y (- y 8)))))
    (loop for n from 5 to (floor slim) do
          (when (/= (aref p n) 0)
                (loop for z from (* n n) below lim by (* n n) do
                      (setf (aref p z) 0))))
    (setf (aref p 2) 1
          (aref p 3) 1)
    p))

(defun main (&rest argv)
  (declare (ignorable argv))
  (let* ((q (read))
         (p (sieve-of-atkin 5000001)))
    (declare (type fixnum q))
    (dotimes (_ q)
      (let* ((a (read))
             (s (read)))
        (declare (type fixnum a))
        (declare (type fixnum s))
        (format t "~d~%" (if (= (aref p s) 1)
                             (if (zerop (mod a s))
                                 0
                                 1)
                             -1))))))

(main)
0