「题解」CSP-S 2020 儒略日
谨以此篇题解，纪念我那炸掉的T1。。
基本思路
大模拟，我写了100多行。。
总之就是按照时间依次向后推进
先判断使用哪一历法，当r大于等于2299162时，使用的是格里高利历，反之，使用儒略历。之后分类讨论。
儒略历（r<2299162）
儒略历的话，因为前4713年正好是一个闰年，我就将四年分为了一组，也就是每4年有1461天，算出经过多少年，再算出减去这些年之后还剩的天数（实际上后面都是这个思路）
因为着四年中，第一年就是一个闰年，所以如果剩下的日子小于366，那就按照闰年算出月和日，反之减去366，算出多了多少年，再算日期。
格里高利历（r>=2299162）
真正麻烦的是格里高利历。
格里高利历又多了“百年不闰，四百年再闰”的规则，于是我将400年作为一个周期，每400年97闰，那么每400年就要经过146097天。
总之，为了方便计算，先把r减去2299162。
因为格里高利历开始的时间是1582年的10月15日，如果在这个日期后，但在1582年里，就要分类讨论。而1582年10月15日距离1583年正好77天，这就有了77这个常数。
1583年前
这部分的处理简单粗暴，分成了三类，分别是10月，11月，12月。每类单独算日期。
1583年后（r>77）
我当时想着是直接快进到1600年，因为这一年正好是400的整倍数，可以作为一个周期的起点，但1583年并不是一个闰年，如果不把1600年之前的每4年组合在一起，将会十分难算。于是我打算将1584年作为一个节点。于是我将在1584年之后的日期和在这之前的日期分开计算。
而当时间快进到1584年时，就可以将剩下的16年以4年一组的形式划分，这样算到1600年。
1600年后
当到了1600后，剩下的事情就相对简单了。
首先将每146097年分成一组，这就是一个400年，首先计算出有多少个400年。
然后判断剩下的日子是否属于第一个百年，因为第一个一百年有25个闰年，而剩下的三个百年中，每一百年只有24个闰年，这也是常数36525的由来。
然后将4年确定为一个周期，确定在这100或300年中，这个日期属于哪个4年。
最后，在确定这个日期是平年或是闰年，分别算出月和日。
代码（看到最后有惊喜）
#include <cstdio>
#include <iostream>
long long q, r, day, month, year;
const int months[15]={31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
const int monthr[15]={31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
//const int julian=2297591;
const int julian=2299162;
int Month(), Monthr();
int main(){
freopen("julian.in", "r", stdin);
freopen("julian.out", "w", stdout);
scanf("%lld", &q);
for(int asdf=0; asdf<q; asdf++){
scanf("%lld", &r);
r++;
if(r<julian){
year=(r/1461)*4;
r%=1461;
if(r<=59){
month=r>31?2:1;
day=r>31?r-31:r;
}
else{
if(r<=366){
month=Monthr();
day=r;
}
else{
year+=(r-366)/365+1;
r-=366; r%=365;
month=Month();
day=r;
}
}
if(year<4713)
if(day==0 && month==1)
printf("31 12 %lld BC\n", 4713-year+1);
else
printf("%lld %lld %lld BC\n", day, month, 4713-year);
else
if(day==0 && month==1)
printf("31 12 %lld\n", year-4713);
else
printf("%lld %lld %lld\n", day, month, year-4713+1);
}
else{
r-=julian;
year=1582;
if(r<=77){
if(r<=16){ month=10; day=15+r; }
else if(r<=46){ month=11; day=r-16; }
else { month=12; day=r-46; }
}
else {
r-=77;
if(r<=365){
year++;
month=Month();
day=r;
}
else{
r-=365; year=1584;
if(r<=5844) {
year+=(r/1461)*4;
r%=1461;
if(r<=366){
month=Monthr();
day=r;
}
else {
r-=366; year++;
year+=r/365;
r%=365;
month=Month();
day=r;
}
}
else{
r-=5844; year=1600;
year+=(r/146097)*400;
r%=146097;
if(r<=36525){
year+=(r/1461)*4;
r%=1461;
if(r<=366){ month=Monthr(); day=r; }
else {
r-=366; year++;
year+=r/365; r%=365;
month=Month(); day=r;
}
}
else{
r-=36525; year+=100;
year+=(r/36524)*100;
r%=36524;
if(r<=1460){
year+=r/365;
r%=365;
month=Month();
day=r;
}
else{
r-=1460; year+=4;
year+=(r/1461)*4;
r%=1461;
if(r<=366) { month=Monthr(); day=r; }
else{
r-=366; year++;
year+=r/365;
r%=365;
month=Month();
day=r;
}
}
}
}
}
}
if(day==0 && month==1)
printf("31 12 %lld\n", year-1);
else
printf("%lld %lld %lld\n", day, month, year);
}
// printf("%d %d %d\n", day, month, year);
}
return 0;
}
int Month(){
int i;
for(i=0; i<12; i++){
if(r>months[i]) r-=months[i];
else break;
}
return i+1;
}
int Monthr(){
int i;
for(i=0; i<12; i++){
if(r>monthr[i]) r-=monthr[i];
else break;
}
return i+1;
}
/*
int Month(int x){
if(0<x && x<=31) return 1;
if(31<x && x<=59) return 2;
if(59<x && x<=90) return 3;
if(90<x && x<=120) return 4;
if(120<x && x<=151) return 5;
if(151<x && x<=181) return 6;
if(181<x && x<=212) return 7;
if(212<x && x<=243) return 8;
if(243<x && x<=273) return 9;
if(273<x && x<=304) return 10;
if(304<x && x<=334) return 11;
if(334<x && x<=365) return 12;
if(0<x && x<=31) return 1;
}
int Monthr(int x){
if(0<x && x<=31) return 1;
if(31<x && x<=60) return 2;
if(60<x && x<=91) return 3;
if(91<x && x<=121) return 4;
if(121<x && x<=152) return 5;
if(152<x && x<=182) return 6;
if(182<x && x<=213) return 7;
if(213<x && x<=244) return 8;
if(244<x && x<=274) return 9;
if(274<x && x<=305) return 10;
if(305<x && x<=335) return 11;
if(335<x && x<=366) return 12;
if(0<x && x<=31) return 1;
}*/
最后附上我算这些常数的过程（python控制台信息）：
# ===== console log in CSP-S 2020 =====
noilinux@ubuntu:~/Desktop/SX-00019/julian$ python
Python 2.7.6 (default, Mar 22 2014, 22:59:38)
[GCC 4.8.2] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> quit()
noilinux@ubuntu:~/Desktop/SX-00019/julian$ python3
Python 3.4.0 (default, Apr 11 2014, 13:05:18)
[GCC 4.8.2] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> 365*400+97
146097
>>> 4713*365
1720245
>>> 4713*365+1581*365
2297310
>>> 31+28+31+30+31+30+31+31+30+4
277
>>> 4713*365+1581*365+277+(4713%4+1)+(1582%4)
2297591
>>> 365*4+1
1461
>>> 365*3+31+28
1154
>>> 31+28
59
>>> 59+31
90
>>> 17+30+31
78
>>> 3000000-2297591
702409
>>> 702409-(16+30+31)
702332
>>> 4*(702332/1461)
1922.880219028063
>>> 4*(702332//1461)
1920
>>> 702332-4*1461
696488
>>> 702332-4*1920
694652
>>> 702332%1461
1052
>>> 1583+1920
3503
>>> 3501-1582
1919
>>> 3501-1583
1918
>>> 198/4
49.5
>>> 1918/4
479.5
>>> 1918*365+479
700549
>>> 700549+17+30+31
700627
>>> 3000000-700627
2299373
>>> 2299373-(31+28+31+30+31+30+31+15)
2299146
>>> 31+28+31+30+31+30+31+31+30+4+(1582+4713-1)*365
2297587
>>> (1582+4713-1)/4
1573.5
>>> 2297587+1573
2299160
>>> 3000000-2299160
700840
>>> 700840-77
700763
>>> 4*(700763//1461)
1916
>>> 1916*1461
2799276
>>> 700763%1461
944
>>> 700763-365
700398
>>> 700398//1461
479
>>> 4*479
1916
>>> 700398%1461
579
>>> 579-366
213
>>> 1916+1584+1
3501
>>> 365*4+1
1461
>>> 100*365+97
36597
>>> 1600-1584
16
>>> 16*365+4
5844
>>> 365*400+97
146097
>>> 100*365+25
36525
>>> 100*365+24
36524
>>> 365*4
1460
>>>
