P7075 [CSP-S2020] 儒略日题解

CSP-S2020儒略日题解

[CSP-S2020] 儒略日

题目描述

为了简便计算，天文学家们使用儒略日（Julian day）来表达时间。所谓儒略日，其定义为从公元前 4713 年 1 月 1 日正午 12 点到此后某一时刻间所经过的天数，不满一天者用小数表达。若利用这一天文学历法，则每一个时刻都将被均匀的映射到数轴上，从而得以很方便的计算它们的差值。

现在，给定一个不含小数部分的儒略日，请你帮忙计算出该儒略日（一定是某一天的中午 12 点）所对应的公历日期。

我们现行的公历为格里高利历（Gregorian calendar），它是在公元 1582 年由教皇格里高利十三世在原有的儒略历（Julian calendar）的基础上修改得到的（注：儒略历与儒略日并无直接关系）。具体而言，现行的公历日期按照以下规则计算：

1. 公元 1582 年 10 月 15 日（含）以后：适用格里高利历，每年一月 $31$ 天、 二月 $28$ 天或 $29$ 天、三月 $31$ 天、四月 $30$ 天、五月 $31$ 天、六月 $30$ 天、七月 $31$ 天、八月 $31$ 天、九月 $30$ 天、十月 $31$ 天、十一月 $30$ 天、十二月 $31$ 天。其中，闰年的二月为 $29$ 天，平年为 $28$ 天。当年份是 $400$ 的倍数，或日期年份是 $4$ 的倍数但不是 $100$ 的倍数时，该年为闰年。
2. 公元 1582 年 10 月 5 日（含）至 10 月 14 日（含）：不存在，这些日期被删除，该年 10 月 4 日之后为 10 月 15 日。
3. 公元 1582 年 10 月 4 日（含）以前：适用儒略历，每月天数与格里高利历相同，但只要年份是 $4$ 的倍数就是闰年。
4. 尽管儒略历于公元前 45 年才开始实行，且初期经过若干次调整，但今天人类习惯于按照儒略历最终的规则反推一切 1582 年 10 月 4 日之前的时间。注意，公元零年并不存在，即公元前 1 年的下一年是公元 1 年。因此公元前 1 年、前 5 年、前 9 年、前 13 年……以此类推的年份应视为闰年。

输入格式

第一行一个整数 $Q$，表示询问的组数。
接下来 $Q$ 行，每行一个非负整数 $r_i$，表示一个儒略日。

输出格式

对于每一个儒略日 $r_i$，输出一行表示日期的字符串 $s_i$。共计 $Q$ 行。 $s_i$ 的格式如下：

1. 若年份为公元后，输出格式为 Day Month Year。其中日（Day）、月（Month）、年（Year）均不含前导零，中间用一个空格隔开。例如：公元
   2020 年 11 月 7 日正午 12 点，输出为 7 11 2020。
2. 若年份为公元前，输出格式为 Day Month Year BC。其中年（Year）输出该年份的数值，其余与公元后相同。例如：公元前 841 年 2 月 1 日正午 12
   点，输出为 1 2 841 BC。

样例 #1

样例输入 #1

3
10
100
1000

样例输出 #1

11 1 4713 BC
10 4 4713 BC
27 9 4711 BC

样例 #2

样例输入 #2

3
2000000
3000000
4000000

样例输出 #2

14 9 763
15 8 3501
12 7 6239

样例 #3

样例输入 #3

见附件中的 julian/julian3.in

样例输出 #3

见附件中的 julian/julian3.ans

提示

【数据范围】

测试点编号	$Q=$	$r_i\le$
$1$	$1000$	$365$
$2$	$1000$	$10^4$
$3$	$1000$	$10^5$
$4$	$10000$	$3\times 10^5$
$5$	$10000$	$2.5\times 10^6$
$6$	$10^5$	$2.5\times 10^6$
$7$	$10^5$	$5\times 10^6$
$8$	$10^5$	$10^7$
$9$	$10^5$	$10^9$
$10$	$10^5$	年份答案不超过 $10^9$

题解

从题目给出的条件，我们不难发现将一个年份转化为儒略??将一个儒略日转化为年份要相对简单，于是我们可以写一个转化的函数，再使用二分答案找到年份
年份分公元前和公元，公元前用负数表示，公元用正数
注意：公元0年不存在，18xx年是19世纪所以我们往前加一年，从公元前4712年开始计算，如果答案是公元前就再减回来。而这里又要注意，公元前4712年恰好也是闰年，还要再加一天
可以得出以下代码

#include <iostream>
#define ull unsigned long long
using namespace std;
ull year_to_day(int year){
	int start = -4712;//公元元年不存在  
    if (year < start) return 0;
    ull res = (year - start + 1) * 365;
    if (year < 1582){
        res += (year - start) / 4 + 1;//儒略历 
    }else{
        res -= 10;//消失的10天
        res += (1581 - start)/4 + 1; //1582年以前用儒略历 
        res += (year-1580)/4 - (year-1500)/100 +(year-1200)/400; // 1582年后用格里高利历，容斥原理求闰年数 
    }
    return res;
}
//Written by Kevin

二分代码:

int l = -4712, r = 1e9+5, mid;
while(l < r){
	mid=(l+r)>>1;
	if(year_to_day(mid)>=d)r=mid;
	else l=mid+1;
}
//Written by Kevin

到现在，我们已经确定了年份
如何求剩下的天数呢？我们可以用总儒略日数减去截止去年的总儒略日，剩下的就是今年还剩多少天
用总天数依次减每个月的天数，直到不够减了，月份和日期就都出来了
记得，闰年2月有29天，1582年10月只有21天

AC代码

#include <iostream>
#define int long long
#define ull unsigned long long
using namespace std;
int T;
ull d,m;
int a[13]={0,31,28,31,30,31,30,31,31,30,31,30,31};//每月天数列表 
ull year_to_day(int year){
	int start = -4712;//公元元年不存在  
    if (year < start) return 0;
    ull res = (year - start + 1) * 365;
    if (year < 1582) {
        res += (year - start) / 4 + 1;//儒略历 
    } else {
        res -= 10;
        res += (1581 - start)/4 + 1; //1582年以前用儒略历 
        res += (year-1580)/4 - (year-1500)/100 +(year-1200)/400; // 1582年后用格里高利历，容斥原理求闰年数 
    }
    return res;
}
bool pd(int x){//判断是否是闰年 
    if (x>=1582)return(x%4==0&&x%100!=0)||(x%400==0); 
    else return x % 4 == 0;
}
signed main(){
	cin >> T;
	while(T--){
		cin >> d;
		d++;
		int l=-4712,r=1e9+5,mid;
		while(l<r){//二分答案 
			mid=(l+r)>>1;
			if(year_to_day(mid)>=d)r=mid;
			else l=mid+1;
		}
		d-=year_to_day(l-1);//减去去年 
		a[2]=28+pd(l);
		a[10] = l==1582 ? 21 : 31;//特判 
		for(int i=1;i<=12;i++){//确定月份和日期 
			if(d>a[i])d-=a[i];
			else{
				m=i;
				break;
			}
		}
		if(l==1582&&m==10&&d>=5)d+=10;//1582年10月没有5-14日，10.4明天就是10.15
		if(l>0)cout << d << " " << m << " " << l << endl;
		else cout << d << " " << m << " " << 1-l << " BC" << endl;//加的一年要减回去 
	}
	return 0;
}
//Written by Kevin