結果

問題 No.721 Die tertia (ディエ・テルツィア)
コンテスト
ユーザー sansaqua
提出日時 2020-11-11 02:22:27
言語 Common Lisp
(sbcl 2.5.0)
結果
AC  
実行時間 11 ms / 2,000 ms
コード長 5,476 bytes
コンパイル時間 1,327 ms
コンパイル使用メモリ 43,452 KB
実行使用メモリ 22,656 KB
最終ジャッジ日時 2024-07-22 18:16:55
合計ジャッジ時間 1,635 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます
コンパイルメッセージ
; compiling file "/home/judge/data/code/Main.lisp" (written 22 JUL 2024 06:16:52 PM):

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

ソースコード

diff # 

(in-package :cl-user)
(eval-when (:compile-toplevel :load-toplevel :execute)
  (defparameter *opt*
    #+swank '(optimize (speed 3) (safety 2))
    #-swank '(optimize (speed 3) (safety 0) (debug 0)))
  #+swank (ql:quickload '(:cl-debug-print :fiveam :cp/util) :silent t)
  #+swank (use-package :cp/util :cl-user)
  #-swank (set-dispatch-macro-character
           #\# #\> (lambda (s c p) (declare (ignore c p)) `(values ,(read s nil nil t)))))
#+swank (set-dispatch-macro-character #\# #\> #'cl-debug-print:debug-print-reader)

(macrolet ((def (b)
             `(progn (deftype ,(intern (format nil "UINT~A" b)) () '(unsigned-byte ,b))
                     (deftype ,(intern (format nil "INT~A" b)) () '(signed-byte ,b))))
           (define-int-types (&rest bits) `(progn ,@(mapcar (lambda (b) `(def ,b)) bits))))
  (define-int-types 2 4 7 8 15 16 31 32 62 63 64))

(defconstant +mod+ 1000000007)

(defmacro dbg (&rest forms)
  #+swank (if (= (length forms) 1)
              `(format *error-output* "~A => ~A~%" ',(car forms) ,(car forms))
              `(format *error-output* "~A => ~A~%" ',forms `(,,@forms)))
  #-swank (declare (ignore forms)))

(declaim (inline println))
(defun println (obj &optional (stream *standard-output*))
  (let ((*read-default-float-format* 'double-float))
    (prog1 (princ obj stream) (terpri stream))))

;; BEGIN_INSERTED_CONTENTS
(defpackage :cp/date
  (:use :cl)
  (:export #:leap-year-p #:get-day-of-week #:date-to-jdn #:jdn-to-date))
(in-package :cp/date)

(declaim (inline leap-year-p))
(defun leap-year-p (year)
  (or (zerop (mod year 400))
      (and (zerop (mod year 4))
           (not (zerop (mod year 100))))))

;; Gauss's algorithm
;; https://en.wikipedia.org/wiki/Determination_of_the_day_of_the_week
(defun get-day-of-week (day month year)
  "Returns the day of the week (starting with Sunday) as an integer from 0 to 6.
Note that DAY and MONTH are 1-based."
  (declare ((integer 1 31) day)
           ((integer 1 12) month)
           (fixnum year))
  (let ((ms #.(coerce #(0 3 3 6 1 4 6 2 5 0 3 5)
                      '(simple-array (unsigned-byte 4) (*))))
        (leap-ms #.(coerce #(0 3 4 0 2 5 0 3 6 1 4 6)
                           '(simple-array (unsigned-byte 4) (*)))))
    (let ((m (aref (if (leap-year-p year) leap-ms ms)
                   (- month 1))))
      (declare ((integer 0 6) m))
      (mod (+ day
              m
              (* 5 (mod (- year 1) 4))
              (* 4 (mod (- year 1) 100))
              (* 6 (mod (- year 1) 400)))
           7))))

(declaim (inline date-to-jdn))
(defun date-to-jdn (day month year)
  "Converts a Gregorian calendar date to Julian day number."
  (declare ((integer 1 31) day)
           ((integer 1 12) month)
           (fixnum year))
  (+ (truncate (* 1461 (+ year 4800 (truncate (- month 14) 12))) 4)
     (truncate (* 367 (- month 2 (* 12 (truncate (- month 14) 12)))) 12)
     (- (truncate (* 3 (truncate (+ year 4900 (truncate (- month 14) 12)) 100)) 4))
     day
     -32075))

;; Reference: https://forum.arduino.cc/index.php?topic=557576.msg3802719#msg3802719
(declaim (inline jdn-to-date)
         (ftype (function * (values (integer 1 31) (integer 1 12) integer &optional))
                jdn-to-date))
(defun jdn-to-date (jdn)
  "Converts a Julian day number to Gregotian calendar date."
  (declare ((integer -32043) jdn)) ;; broken under this value
  (let* ((l (+ jdn 68569))
         (n (truncate (* 4 l) 146097))
         (l (- l (truncate (+ (* 146097 n) 3) 4)))
         (year (truncate (* 4000 (+ l 1)) 1461001))
         (l (+ (- l (truncate (* 1461 year) 4)) 31))
         (month (truncate (* 80 l) 2447))
         (day (- l (truncate (* 2447 month) 80)))
         (l (truncate month 11))
         (month (- (+ month 2) (* 12 l)))
         (year (+ year l (* 100 (- n 49)))))
    (values day month year)))

;; BEGIN_USE_PACKAGE
(eval-when (:compile-toplevel :load-toplevel :execute)
  (use-package :cp/date :cl-user))
(in-package :cl-user)

;;;
;;; Body
;;;

(defun main ()
  (let* ((line (read-line))
         (y (read-from-string line t nil :start 0 :end 4))
         (m (read-from-string line t nil :start 5 :end 7))
         (d (read-from-string line t nil :start 8 :end 10)))
    (multiple-value-bind (d2 m2 y2) (jdn-to-date (+ 2 (date-to-jdn d m y)))
      (format t "~4,'0D/~2,'0D/~2,'0D~%" y2 m2 d2))))

#-swank (main)

;;;
;;; Test and benchmark
;;;

#+swank
(progn
  (defparameter *lisp-file-pathname* (uiop:current-lisp-file-pathname))
  (setq *default-pathname-defaults* (uiop:pathname-directory-pathname *lisp-file-pathname*))
  (defparameter *dat-pathname* (uiop:merge-pathnames* "test.dat" *lisp-file-pathname*))
  (defparameter *problem-url* "https://yukicoder.me/problems/no/721"))

#+swank
(defun gen-dat ()
  (uiop:with-output-file (out *dat-pathname* :if-exists :supersede)
    (format out "")))

#+swank
(defun bench (&optional (out (make-broadcast-stream)))
  (time (run *dat-pathname* out)))

#-swank
(eval-when (:compile-toplevel)
  (when (or (> sb-c::*compiler-warning-count* 0)
            sb-c::*undefined-warnings*)
    (error "count: ~D, undefined warnings: ~A"
           sb-c::*compiler-warning-count*
           sb-c::*undefined-warnings*)))

;; To run: (5am:run! :sample)
#+swank
(5am:test :sample
  (5am:is
   (equal "2000/01/03
"
          (run "2000/01/01
" nil)))
  (5am:is
   (equal "2000/02/29
"
          (run "2000/02/27
" nil)))
  (5am:is
   (equal "2017/01/01
"
          (run "2016/12/30
" nil))))
0