首页 > 其他分享 >【毒瘤】儒略日 做题笔记

【毒瘤】儒略日 做题笔记

时间：2022-08-27 20:00:38浏览次数：82  
标签：int ans ++ 笔记 month 毒瘤 儒略 year lld

众所周知：

P7075 [CSP-S2020] 儒略日 是一道十分 \(EX\) 的大模拟

作为 \(CSP-S\) \(2020\) 的 \(T1\) 远难于 \(T2\)。

更是有巨佬将过这题的记录放在主页作为自豪的象征。

以前我也略有耳闻，但都没仔细看过题。

于是，今天，我下定决心，要切掉它！（还不是因为算法废了只能练模拟qwq）

于是，\(2022.8.27\) \(8:30 a.m.\)，我打开了 儒略日 !

0X00 做题前

首先花了好长时间读题，因为大模拟万一题目哪里理解错了直接寄掉。然后又花了好长时间来构思怎么解决。

首先，公元前 \(B.C\) 与 公元 \(A.D\) 的顺序就不一样，肯定分开来做。但公元 \(A.D\) 中有被删除的一段日期，而且前后对于闰年的定义不一样，所以也想到分开来做。

再考虑头疼的闰年问题。\(1583\) 年及以后的闰年就常规判。之前的话，公元前 \(B.C\) 一共只有四千多年， 公元 \(A.D\) 里只有一千多年，所以完全可以预处理。之后算年月就常规做。

但是最多年份有 \(10^9\)，而且询问有 \(10^5\)，所以一年一年排年份肯定超时。这时，很容易想到，公元 \(1582\) 年之前可以把 \(4\) 年打包起来做，取模 \(4\) 年总天数再处理零头即可，因为每 \(4\) 年的天数是一样的。而公元 \(1582\) 年之后就可以把 \(400\) 年打包起来做，取模 \(400\) 年总天数再处理零头即可，因为每 \(400\) 年的天数是一样的。

这样一想，也不是很难啊！于是，这个菜ji的脸上露出了蜜汁自信の微笑（请记住这个笑容）

0X01 计算

想到了分段做的思路，就要看怎么分段了。因为给出的是 公元前 4713 年 1 月 1 日正午 12 点到此后某一时刻间所经过的天数，所以对于分段标准也要做处理。

首先要算出的是 公元前 4713 年 1 月 1 日 到 公元前 1 年 12 月 31 日 一共有多少天。这个不难，只要找之前思路打好公元前 \(B.C\) 闰年表再加起来即可。因为公元前 4713 年 1 月 1 日不算，所以要 \(-1\)：

for(int i=1;i<=4713;i++) if(i%4==1) bc_run[i]=1;
for(int i=1;i<=4713;i++) t+=(bc_run[i])?366:365;
cout<<t-1;

于是，我们得出了 \(BC\) 段有 \(1721423\) 天。

其次要算出的是 公元 1 年 1 月 1 日 到 公元 1582 年 10 月 4 日 一共有多少天。这个也不难，只要找之前思路打好公元 \(A.D\) 闰年表再加起来即可。加到 \(1581\) 年，再单独加上 \(1582\) 年的 \(9\) 个月零 \(4\) 天：

for(int i=1;i<=1582;i++) if(i%4==1) ad_run[i]=1;
for(int i=1;i<=1581;i++) t+=(ad_run[i])?366:365;
cout<<t-1+31+28+31+30+31+30+31+31+30+4;

于是，我们得出了 \(AD\) 儒略历段有 \(577737\) 天。

然后要算出的是公元 \(1582\) 年之前 \(4\) 年的总天数 和 公元 \(1582\) 年之后 \(400\) 年的总天数，方便打包取模。

公元 \(1582\) 年之前 \(4\) 年的总天数 手算即可，有 \(1461\) 天。

公元 \(1582\) 年之后 \(400\) 年的总天数 常规判闰年算：

bool check_run(int x){
	if(x%400==0) return 1;
	if(x%100!=0&&x%4==0) return 1;
	return 0;
}
for(int i=1;i<=400;i++){
	if(check_run(i)) t++;
	t+=365;
}
cout<<t;

得出的结果是 \(146097\) 天。

0X02 代码

对于公元前，首先我想到的是将 \(a\) 赋为 \(1721423-a\)，然后倒着做。

先抹到零头。因为公元前 \(1,5,9...\) 是闰年，所以抹完后第一年应该是闰年，判一下到不到 \(366\) 天，到的话加一年，然后看到不到 \(365\) 天，一直到不满一年（最多判三次就到了）。然后再 \(+1\) 就得到了最终年份（因为从 \(1\) 年开始）。根据答案年份对照公元前 \(B.C\) 闰年表判闰年，修改 \(2\) 月天数，之后一个月一个月倒着排，最终剩下的就是日了。于是满心欢喜写出了这段代码：

if(a<=1721423){
	a=1721423-a;
	int t=a/1461;   //判断几个4年
	a%=1461;        //抹零头
	ans_year=t*4;  
	if(a>366) a-=366,ans_year++;
	while(a>365) a-=365,ans_year++;
	ans_year++;          //得出答案年份
	if(bc_run[ans_year]) M[2]=29;
	else M[2]=28;        //改2月天数
	for(int i=12;i>=1;i--){ 
		if(a>M[i]){
			a-=M[i];
			ans_month=i-1;
		}
		else break;
	}
	ans_day=M[ans_month]-a;
	printf("%lld %lld %lld BC\n",ans_day,ans_month,ans_year);
	return ;
}

过样例 \(1\) 了！就这？（不屑

然鹅，其实全是错误 \(qwq\)

暂且把它放着，先看其他部分。照模照样写了一下，但过程中发现 公元 1582 年 10 月 15 日 到 公元 1582 年 12 月 31日 一段很难处理，就单独放了出来，写特判（反正只有 \(78\) 天）。

于是，一个小时多一点，我就写出了以下的代码：

#include<bits/stdc++.h>
using namespace std;
bool bc_run[5005],ad_run[2005];
long long a,ans_year,ans_month,ans_day,t;
int T,M[13]={0,31,28,31,30,31,30,31,31,30,31,30,31};

int check_run(int x){
	if(x%400==0) return 1;
	if(x%100!=0&&x%4==0) return 1;
	return 0;
}

void solve(){
	
	scanf("%lld",&a);
	
//================================ 公元前================================ 

	if(a<=1721423){
		a=1721423-a;
		int t=a/1461;
		a%=1461;
		ans_year=t*4;
		if(a>366) a-=366,ans_year++;
		while(a>365) a-=365,ans_year++;
		ans_year++;
		if(bc_run[ans_year]) M[2]=29;
		else M[2]=28;
		for(int i=12;i>=1;i--){
			if(a>M[i]){
				a-=M[i];
				ans_month=i-1;
			}
			else break;
		}
		ans_day=M[ans_month]-a;
		printf("%lld %lld %lld BC\n",ans_day,ans_month,ans_year);
		return ;
	}
	a-=1721423;
	
//======================== 公元1.1.1 ~ 公元1582.10.4 ======================== 
	
	if(a<=577737){
		int t=a/1461;
		a%=1461;
		ans_year=t*4;
		while(a>365) a-=365,ans_year++;
		ans_year++;
		if(ad_run[ans_year])  M[2]=29;
		else M[2]=28;
		for(int i=1;i<=12;i++){
			if(a>M[i]){
				a-=M[i];
				ans_month=i+1;
			}
			else break;
		}
		ans_day=a;
		printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
		return ;
	}
	a-=577737;
	
//====================== 公元1582.10.15 ~ 公元1582.12.31 ====================== 
	
	if(a<=78){
		ans_year=1582,ans_month=10;
		if(a>17) ans_month++,a-=17;
		if(a>30) ans_month++,a-=30;
		ans_day=a;
		if(ans_month==10) ans_day+=14;
		printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
		return ;
	}
	a-=78;
	
//============================== 公元1583.1.1 ~ ? ============================== 
	
	ans_year=1583;
	int t=a/146097;
	a%=146097;
	ans_year+=t*400;
	for(int i=1;i<=400;i++){
		int tmp=365+check_run(i);
		//cout<<tmp<<endl;
		if(a>tmp) a-=tmp,ans_year++;
	}
	if(check_run(ans_year))  M[2]=29;
	else M[2]=28;
	ans_month=1; 
	for(int i=1;i<=12;i++){
		if(a>M[i]){
			a-=M[i];
			ans_month=i+1;
		}
		else break;
	}
	ans_day=a;
	printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
	return ;
}

int main(){
	for(int i=1;i<=4713;i++) if(i%4==1) bc_run[i]=1;
	for(int i=1;i<=1582;i++) if(i%4==1) ad_run[i]=1;
	scanf("%d",&T);
	while(T--) solve();
	return 0;
}

同时获得了0分的好成绩

当时感觉可能是公元前倒着做法假了。于是打了一边正着做的，获得了 \(10\) 分的高分。手造几组数据后发现， 闰年的 12 月 31 日 被吃了。

仔细检查代码后发现，这一段出了问题：

if(a>366) a-=366,ans_year++;
while(a>365) a-=365,ans_year++;

如果 \(a\) 刚好剩下 \(366\) 天，那么不会进第一重判断，但任会进 \(while\) 里的判断，所以就变成了“假闰年”，尽管 \(2\) 月有 \(29\) 天，但一年总共就 \(365\) 天，所以最后一天就被吃了 \(qwq\)。

改成这样就好了：

if(a>366){
	a-=366,ans_year++;
	while(a>365) a-=365,ans_year++;
}

于是，就取得了 30 分的高分！

当时自认为公元前部分没问题了，于是就去改后面了。

可就在一次乱试数据的过程中，我发现公元前有一个输出，日变成了 \(0\)。

所以就又写了一个特判，防止它被取模后是零。如果是 \(0\)，就输出前一年的最后一天。这下公元前真没问题了。

但这个菜ji当时脑抽，并没有想到后面也可能出现这个问题！

公元前完整代码：

//================================== 公元前================================== 

if(a<=1721423){
	a++;
	int t=a/1461;
	a%=1461;
	ans_year=t*4;
	if(a==0){
		printf("31 12 %lld BC\n",4713-ans_year+1);
		return ;
	}
	if(a>366){
		a-=366,ans_year++;
		while(a>365) a-=365,ans_year++;
	}
	if(bc_run[4713-ans_year]) M[2]=29;
	else M[2]=28;
	ans_month=1;
	for(int i=1;i<=12;i++){
		if(a>M[i]){
			a-=M[i];
			ans_month=i+1;
		}
		else break;
	}
	ans_day=a;
	printf("%lld %lld %lld BC\n",ans_day,ans_month,4713-ans_year);
	return ;
}

然后，经过漫长的找错过程。

中途瞎改，获得了 \(40、50\) 的奇怪分数。

在几个手搓数据中发现，有些年 1 月 1 日 被重复输出两次，同时 12 月 31 日 被吃掉。有些年输出 1 月 0 日。（这时候还没想到！）

吃完中饭回来，终于在 公元 1583 年 1 月 1日 及之后 的部分找到了严重的错误：

for(int i=1;i<=400;i++){
	int tmp=365+check_run(i);
	if(a>tmp) a-=tmp,ans_year++;
}

这一处理零头年份部分的代码是不是很眼熟？没错，是算 \(400\) 年有多少天那里复制过来的。

这就出问题了，\(check\_run(i)\) 干什么？判断 \(1-400\) 年里是不是闰年？应该\(check\_run(ans\_year)\) 啊！反正就是脑抽的表现。应该改成：

for(int i=1;i<=400;i++){
	int tmp=365+check_run(ans_year);
	if(a>tmp) a-=tmp,ans_year++;
	else break;
}

然后就获得了 70 分的好成绩！

之后实在找不出了，就去下载了错误数据，仔细核对后才发现了 1 月 0 日的问题就是 \(a\) 取模完变 \(0\) 了。于是在 公元 1 年 1 月 1 日 到 公元 1582 年 10 月 4 日 一段加上特判。

if(a==0){
	printf("31 12 %lld\n",ans_year);
	return ;
}

就取得了 90 分高分！

最后在 公元 1583 年 1 月 1日 及之后 也加了特判：

if(a==0){
	printf("31 12 %lld\n",ans_year-1);
	return ;
}

AC !

真的太激动了！

历时 \(3.5h +\)，切掉 [CSP-S2020] 儒略日

太不容易了 \(qwq\)

最后是完整代码：

Code:

#include<bits/stdc++.h>
using namespace std;
bool bc_run[5005],ad_run[2005];
long long a,ans_year,ans_month,ans_day,t;
int T,M[13]={0,31,28,31,30,31,30,31,31,30,31,30,31};

int check_run(int x){
	if(x%400==0) return 1;
	if(x%100!=0&&x%4==0) return 1;
	return 0;
}

void solve(){

//================================== 公元前================================== 

	if(a<=1721423){
		a++;
		int t=a/1461;
		a%=1461;
		ans_year=t*4;
		if(a==0){
			printf("31 12 %lld BC\n",4713-ans_year+1);
			return ;
		}
		if(a>366){
			a-=366,ans_year++;
			while(a>365) a-=365,ans_year++;
		}
		if(bc_run[4713-ans_year]) M[2]=29;
		else M[2]=28;
		ans_month=1;
		for(int i=1;i<=12;i++){
			if(a>M[i]){
				a-=M[i];
				ans_month=i+1;
			}
			else break;
		}
		ans_day=a;
		printf("%lld %lld %lld BC\n",ans_day,ans_month,4713-ans_year);
		return ;
	}
	a-=1721423;
	
//========================= 公元1.1.1 ~ 公元1582.10.4 =========================
	
	if(a<=577737){
		int t=a/1461;
		a%=1461;
		ans_year=t*4;
		if(a==0){
			printf("31 12 %lld\n",ans_year);
			return ;
		}
		ans_year++;
		while(a>365) a-=365,ans_year++;
		if(ad_run[ans_year])  M[2]=29;
		else M[2]=28;
		ans_month=1;
		for(int i=1;i<=12;i++){
			if(a>M[i]){
				a-=M[i];
				ans_month=i+1;
			}
			else break;
		}
		ans_day=a;
		printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
		return ;
	}
	a-=577737;
	
//====================== 公元1582.10.15 ~ 公元1582.12.31 =======================
	
	if(a<=78){
		ans_year=1582,ans_month=10;
		if(a>17) ans_month++,a-=17;
		if(a>30) ans_month++,a-=30;
		ans_day=a;
		if(ans_month==10) ans_day+=14;
		printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
		return ;
	}
	a-=78;
	
//============================== 公元1583.1.1 ~ ? ============================== 

	ans_year=1583;
	int t=a/146097;
	a%=146097;
	ans_year+=t*400;
	for(int i=1;i<=400;i++){
		int tmp=365+check_run(ans_year);
		if(a>tmp) a-=tmp,ans_year++;
		else break;
	}
	if(a==0){
		printf("31 12 %lld\n",ans_year-1);
		return ;
	}
	if(check_run(ans_year))  M[2]=29;
	else M[2]=28;
	ans_month=1; 
	for(int i=1;i<=12;i++){
		if(a>M[i]){
			a-=M[i];
			ans_month=i+1;
		}
		else break;
	}
	ans_day=a;
	printf("%lld %lld %lld\n",ans_day,ans_month,ans_year);
	return ;
	
//==============================================================================
	
}

int main(){
	//freopen("tt.txt","w",stdout);
	//freopen("r.in","r",stdin);
	//freopen("my.in","w",stdout);
	for(int i=1;i<=4713;i++) if(i%4==1) bc_run[i]=1;
	for(int i=1;i<=1582;i++) if(i%4==0) ad_run[i]=1;
	scanf("%d",&T);
	for(int i=1;i<=T;i++){
		scanf("%lld",&a);
		solve();
	}
	return 0;
}



/*
for(int i=1;i<=4713;i++) if(i%4==1) bc_run[i]=1;
for(int i=1;i<=4713;i++) t+=(bc_run[i])?366:365;
cout<<t-1;

BC: 1721423
*/


/*
for(int i=1;i<=1582;i++) if(i%4==1) ad_run[i]=1;
for(int i=1;i<=1581;i++) t+=(ad_run[i])?366:365;
cout<<t-1+31+28+31+30+31+30+31+31+30+4;

AD 儒略历: 577737
*/


/*
bool check_run(int x){
	if(x%400==0) return 1;
	if(x%100!=0&&x%4==0) return 1;
	return 0;
}
for(int i=1;i<=400;i++){
	if(check_run(i)) t++;
	t+=365;
}
cout<<t;

400年周期: 146097
*/

G L & H F ！

标签：int,ans,++,笔记,month,毒瘤,儒略,year,lld
From： https://www.cnblogs.com/binary1110011/p/16631329.html

相关文章

  • 2022-08-25 第二小组 张鑫 学习笔记
    实训四十七天元素操作BOM1.学习内容自定义属性设置元素属性<divhaha="abc"id="xyz"></div><script>letdiv=document.querySelector("div");//......
  • enote笔记法之附录1——“语法词”（即“关联词”）（ver0.2）
    章节：enote笔记法之附录1——“语法词”（即“关联词”）（ver0.2）上面的是截屏的完整版，分割线下面的是纯文字版本：  作者姓名（本人的真实姓名）：胡佳吉 居住地：上海作者网名：E......
  • 论文笔记 - ToAlign: Task-oriented Alignment for Unsupervised Domain Adaptation
    摘要文章提出了一个任务导向的无监督域自适应。认为不是所有的特征都需要对齐，而是根据任务的先验知识，将源域特征分解为与任务相关的要对齐的，和与任务无关的可忽略的。......
  • GNSS模块使用笔记
    目录目录GNSS芯片NMEA0183协议GNSS芯片ATGM336H-5N31（GPS+BDS双模）原理图NMEA0183协议《CASIC多模卫星导航接收机协议规范》前缀BD---北斗模式GP---GPS模......
  • 人物交互算法（HOI）学习笔记之 ——QPIC
    论文简介QPIC:Query-BasedPairwiseHuman-ObjectInteractionDetectionwithImage-WideContextualInformation[论文地址][https://arxiv.org/abs/2103.05399][代......
  • Rust学习笔记-初识
    很早以前，看到马斯克说：“我是Rust的粉丝，但是为了性能会选择C。”之后就对Rust产生了一些兴趣，不过由于时间关系没有及时深入了解。最近，又看到Meta选定四种编程语言作为首选......
  • js笔记整理2
    对象及日期定时器日期对象的定义（使用new关键词）1.获取当前的时间（本地的时间）vardate=newDate()//不传参就是获取当前时间2.获取指定的时间vardate=newDate(......
  • dragonfly 蜻蜓算法 学习笔记
    1、GettingStated1.1CommandLine使用方法：在pycharm中：cdexamplepython..\bin\dragonfly-script.py--configxxx.json--optionxxx.txt1)BasicUse全局优化......
  • 数据库学习笔记 （本数据库学习笔记以SQL sever 2019 为例进行学习） 20220823 第一节课
    教材及参考数据库课程讲什么？内容安排第一部分数据库原理部分第一章数据库系统概述为什么要学习数据库？数据库的发展改变了人们的工作和生活模式信息积累与运用......
  • 【笔记】线性代数
    矩阵乘法首先给出矩阵乘法的代数意义：结合一个具体的例子来理解：设答案矩阵为\(ans\)。根据公式：\(ans_{1,1}\)是由\(a\)矩阵的第一行与\(b\)矩阵的第一列逐位相......

赞助商

阅读排行

网站主页关于我们联系我们网站地图

本网站内容转载自其他媒体，侵权联系[admin##ips99.com]。

Copyright © 2020-2023 IPS99 版权所有 IPS99