結果

問題 No.188 HAPPY DAY
ユーザー sansaquasansaqua
提出日時 2020-11-22 13:43:54
言語 Common Lisp
(sbcl 2.3.8)
結果
AC  
実行時間 10 ms / 1,000 ms
コード長 5,491 bytes
コンパイル時間 77 ms
コンパイル使用メモリ 42,436 KB
実行使用メモリ 28,488 KB
最終ジャッジ日時 2024-07-23 16:27:37
合計ジャッジ時間 452 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 10 ms
28,488 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
; compiling file "/home/judge/data/code/Main.lisp" (written 23 JUL 2024 04:27:36 PM):

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

ソースコード

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*
          (if (typep obj 'double-float) 'double-float *read-default-float-format*)))
    (prog1 (princ obj stream) (terpri stream))))

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

(declaim ((simple-array (unsigned-byte 8) (*)) *month-days*))
(defparameter *month-days*
  #.(coerce #(31 28 31 30 31 30 31 31 30 31 30 31)
            '(simple-array (unsigned-byte 8) (*))))

(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 -32044) 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 ()
  (println
   (loop for m from 1 to 12
         for max-d = (aref *month-days* (- m 1))
         sum (loop for d from 1 to max-d
                   count (= m (+ (floor d 10)
                                 (mod d 10)))))))

#-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*))
  (uiop:chdir *default-pathname-defaults*)
  (defparameter *dat-pathname* (uiop:merge-pathnames* "test.dat" *lisp-file-pathname*))
  (defparameter *problem-url* "https://yukicoder.me/problems/no/188"))

#+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)))

#+(and sbcl (not 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)
0