那一日，他们又想起了儒略日带来的恐惧

题目描述

为了简便计算，天文学家们使用儒略日（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）均不含前导零，中间用一个空格隔开。例如：公元 2024 年 11 月 7 日正午 12 点，输出为 7 11 2024。
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

【数据范围】

$Q \leq 10^5$，年份答案不超过 $10^9$

BUG 代码

#include<stdio.h>
#include<math.h>

#define ma 3000005
int Q;
int flag=0;
long long ri,year=-4713,month=1,day=1,p;
const int k=365*400+97;
int a[15]={0,31,28,31,30,31,30,31,31,30,31,30,31};
int ansy[ma],ansm[ma],ansd[ma];
int run(int year)
{
    if(year<0)
    {
        if(abs(year)%4==1) return 1;
        return 0;
    }
    if(year>1582)
    {
        if(year%400==0) return 1;
        if(year%4==0&&year%100!=0) return 1;
        return 0;
    }
    if(year>0)
        if(year%4==0) return 1;
    return 0;
}
int main()
{
    scanf("%d",&Q);
    for(int i=1;year<=1599;i++)
    {	
        day++;
        ansy[i]=year;ansm[i]=month;ansd[i]=day;
        if(month==2&&run(year)&&day==29)
        {
            month++,day=0;
            continue;
        }
        if(!run(year)&&day>=a[month])	month++,day=0;
        else if(month!=2&&run(year)&&day>=a[month])	month++,day=0;
        if(month==13) year++,month=1,day=0;
        if(year==1582&&month==10&&day==5) day=15;
    }
    ansy[0]=-4713,ansm[0]=1,ansd[0]=1;
    int tot;
    while(Q--)
    {
        scanf("%lld",&ri);
        if(ri<=2305447)
        {

            if(ansy[ri]<0) printf("%d %d %d BC\n",ansd[ri],ansm[ri],ansy[ri]);
            else printf("%d %d %d\n",ansd[ri],ansm[ri],ansy[ri]);
            continue;
        }
        ri-=2305447;
        year=1600,month=1,day=0;
        p=ri/k;
        year+=400*p;
        p=ri%k;
        while(1)
        {
            if(run(year)&&p>=366)
            {
                year++;
                p-=366;
                continue;
            }
            else if(!run(year)&&p>=365)
            {
                year++;
                p-=365;
                continue;
            }
            else break;
        }
        for(int i=1;i<=12;i++)
        {
            if(i==2&&run(year))
            {
                if(p>29)
                {
                    p-=29;
                    continue;	
                }
                printf("%d %d %lld\n",day+p,i,year);
                break;
            }
            else if(p>a[i])
            {
                p-=a[i];
                continue;
            }
            printf("%d %d %lld\n",day+p,i,year);
            break;
        }
    }
    return 0; 
}

Code: clay17