将年月日转换为儒略日和约化儒略日的程序代码

儒略日（Julian Day Number, JD）是从公元前4713年1月1日中午（儒略历）开始算起的一天中的天数，包括闰年的额外日子。而约化儒略日（Reduced Julian Day, RJD）则是在儒略日的基础上去掉整百年份不是400倍数的那部分，简化了计算。

以下是Python中将年月日转换为儒略日的一个简单示例：

from datetime import datetime

def julian_day(year, month, day):
    a = (14 - month) // 12
    y = year + 4800 - a
    m = month + 12*a - 3
    jdn = day + (153*m+2)/5 + 365*y + y//4 - y//100 + y//400 - 32045
    return jdn

def reduced_julian_day(jdn):
    if jdn >= 2299161: # 跨过了公元元年的分界线
        rjd = jdn - 1524763.75
    else:
        rjd = jdn + 1524763.75
    return int(rjd)

year = int(input("请输入年份："))
month = int(input("请输入月份（1-12）："))
day = int(input("请输入日期："))

julian_date = julian_day(year, month, day)
rjd = reduced_julian_day(julian_date)

print(f"儒略日: {julian_date}")
print(f"约化儒略日: {rjd}")