結果

問題 No.181 A↑↑N mod M
ユーザー Common LispCommon Lisp
提出日時 2024-11-17 01:59:47
言語 Common Lisp
(sbcl 2.3.8)
結果
AC  
実行時間 10 ms / 5,000 ms
コード長 2,782 bytes
コンパイル時間 1,457 ms
コンパイル使用メモリ 30,848 KB
実行使用メモリ 22,144 KB
最終ジャッジ日時 2024-11-17 01:59:50
合計ジャッジ時間 1,548 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 8 ms
22,144 KB
testcase_01 AC 8 ms
22,144 KB
testcase_02 AC 8 ms
22,144 KB
testcase_03 AC 8 ms
22,144 KB
testcase_04 AC 8 ms
22,016 KB
testcase_05 AC 9 ms
22,144 KB
testcase_06 AC 8 ms
22,016 KB
testcase_07 AC 8 ms
22,016 KB
testcase_08 AC 7 ms
22,016 KB
testcase_09 AC 7 ms
22,144 KB
testcase_10 AC 8 ms
22,144 KB
testcase_11 AC 8 ms
22,144 KB
testcase_12 AC 8 ms
22,016 KB
testcase_13 AC 7 ms
22,144 KB
testcase_14 AC 8 ms
22,144 KB
testcase_15 AC 8 ms
22,016 KB
testcase_16 AC 8 ms
22,144 KB
testcase_17 AC 8 ms
22,144 KB
testcase_18 AC 7 ms
22,016 KB
testcase_19 AC 8 ms
22,016 KB
testcase_20 AC 8 ms
22,144 KB
testcase_21 AC 9 ms
22,144 KB
testcase_22 AC 8 ms
22,144 KB
testcase_23 AC 9 ms
22,144 KB
testcase_24 AC 9 ms
22,016 KB
testcase_25 AC 8 ms
22,144 KB
testcase_26 AC 9 ms
22,144 KB
testcase_27 AC 8 ms
22,144 KB
testcase_28 AC 8 ms
22,144 KB
testcase_29 AC 7 ms
22,016 KB
testcase_30 AC 7 ms
22,144 KB
testcase_31 AC 9 ms
22,144 KB
testcase_32 AC 8 ms
22,144 KB
testcase_33 AC 8 ms
22,016 KB
testcase_34 AC 10 ms
22,016 KB
testcase_35 AC 8 ms
22,144 KB
testcase_36 AC 8 ms
22,144 KB
testcase_37 AC 7 ms
22,016 KB
testcase_38 AC 8 ms
22,144 KB
testcase_39 AC 9 ms
22,016 KB
testcase_40 AC 9 ms
22,144 KB
testcase_41 AC 9 ms
22,144 KB
testcase_42 AC 8 ms
22,144 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
; compiling file "/home/judge/data/code/Main.lisp" (written 17 NOV 2024 01:59:47 AM):

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

ソースコード

diff # 

(defun prime-table (&optional (m 1000000000))
  (let* ((thre (ceiling (sqrt m)))
         (flg (make-array (1+ thre) :initial-element nil))
         (i 2)
         (res '()))
    (loop while (<= (* i i) thre) do
          (unless (aref flg i)
                  (loop for j from (* i i) to thre by i do
                        (setf (aref flg j) t)))
          (incf i))
    (loop for i from 2 to thre do
          (unless (aref flg i)
                  (push i res)))
    (nreverse res)))

(defvar *prime-table* (prime-table))

(defun totient (n_)
  (let ((n n_)
        (res n_))
    (dolist (p *prime-table*)
      (when (> (* p p) n)
            (return))
      (when (zerop (mod n p))
            (setq res (* (floor res p) (1- p)))
            (loop while (zerop (mod n p)) do (setf n (floor n p)))))
    (when (/= n 1)
          (setf res (* (floor res n) (1- n))))
    res))

(defun mpow (a_ p_ m_)
  (let ((a a_)
        (p p_)
        (m m_)
        (flg t)
        (res (mod 1 m_)))
    (loop unless (zerop p) do
          (when (oddp p)
                (setf res (* res a))
                (when (>= res m)
                      (setf flg nil
                            res (mod res m))))
          (when (= p 1)
                (return))
          (setf a (* a a))
          (when (>= a m)
                (setf flg nil
                      a (mod a m)))
          (setf p (ash p -1)))
    (cons flg res)))

(defun tetration (a b m)
  (labels ((rec (aa bb mm)
                (cond ((= 0 aa) (cons t (logxor (logand bb 1) 1)))
                      ((= 1 aa) (cons t 1))
                      ((= 1 mm) (cons nil 0))
                      ((= 0 bb) (cons t 1))
                      ((= 1 bb) (cons (< aa mm) (mod aa mm)))
                      (t (let* ((phi-m (totient mm))
                                (flg1pre (rec aa (1- bb) phi-m))
                                (flg1 (car flg1pre))
                                (pre (cdr flg1pre))
                                (flg t)
                                (res 0))
                           (if flg1
                               (let ((flgres (mpow (mod aa mm) pre mm)))
                                 (setf flg (car flgres)
                                       res (cdr flgres))
                                 (cons (and flg flg1) res))
                               (let ((flgres (mpow (mod aa mm) (+ pre phi-m) mm)))
                                 (setf flg (car flgres)
                                       res (cdr flgres))
                                 (cons (and flg flg1) res))))))))
    (mod (cdr (rec a b m)) m)))

(defun main (&rest argv)
  (declare (ignorable argv))
  (let* ((a (read))
         (n (read))
         (m (read)))
    (format t "~d~%" (tetration a n m))))

(main)
0