公历转阴历算法

简述

阳历和阴历之间是不能换算的，二者采用的计算方法不同，没有特殊关系可以查询，只能分别计算。
公历是阳历的一种，阳历就是以地球绕太阳公转的运动周期为基础而制定的历法。以耶稣诞生之年作为纪年的开始。
阴历，又称太阴历，是以月球绕行地球一周为一个月，即以朔望月作为确定历月的基础的一种历法。

源代码

#include <bits/stdc++.h>
using namespace std;

/*
太阳公转周期，一年时间约365.2425天。
公历为了对齐公转周期，4年一闰，100年一停闰，400年加一闰。 400年共97闰， （365×400+97）/400 = 365.2425
农历，一个月相周期为一月，月相周期大约为29.53,   有大小月之分，大月30天，小月29天。
农历为了对齐公转周期，每19年加7个闰月，
因为，农历的推算比较复杂， 大多使用查表法进行计算.
*/

/*
使用比特位记录每年的情况
数据说明
0~4 共5bit 春节日份
5~6 共2bit 春节月份
7~19 共13bit 13个月的大小月情况(如果无闰月，最后位无效)，大月为1（30天）,小月为0（29天）
20~23 共4bit 记录闰月的月份，如果没有闰月为0
*/
int lunarYear[199] = {
    0x04AE53, 0x0A5748, 0x5526BD, 0x0D2650, 0x0D9544, 0x46AAB9, 0x056A4D, 0x09AD42, 0x24AEB6, 0x04AE4A, /*1901-1910*/
    0x6A4DBE, 0x0A4D52, 0x0D2546, 0x5D52BA, 0x0B544E, 0x0D6A43, 0x296D37, 0x095B4B, 0x749BC1, 0x049754, /*1911-1920*/
    0x0A4B48, 0x5B25BC, 0x06A550, 0x06D445, 0x4ADAB8, 0x02B64D, 0x095742, 0x2497B7, 0x04974A, 0x664B3E, /*1921-1930*/
    0x0D4A51, 0x0EA546, 0x56D4BA, 0x05AD4E, 0x02B644, 0x393738, 0x092E4B, 0x7C96BF, 0x0C9553, 0x0D4A48, /*1931-1940*/
    0x6DA53B, 0x0B554F, 0x056A45, 0x4AADB9, 0x025D4D, 0x092D42, 0x2C95B6, 0x0A954A, 0x7B4ABD, 0x06CA51, /*1941-1950*/
    0x0B5546, 0x555ABB, 0x04DA4E, 0x0A5B43, 0x352BB8, 0x052B4C, 0x8A953F, 0x0E9552, 0x06AA48, 0x6AD53C, /*1951-1960*/
    0x0AB54F, 0x04B645, 0x4A5739, 0x0A574D, 0x052642, 0x3E9335, 0x0D9549, 0x75AABE, 0x056A51, 0x096D46, /*1961-1970*/
    0x54AEBB, 0x04AD4F, 0x0A4D43, 0x4D26B7, 0x0D254B, 0x8D52BF, 0x0B5452, 0x0B6A47, 0x696D3C, 0x095B50, /*1971-1980*/
    0x049B45, 0x4A4BB9, 0x0A4B4D, 0xAB25C2, 0x06A554, 0x06D449, 0x6ADA3D, 0x0AB651, 0x093746, 0x5497BB, /*1981-1990*/
    0x04974F, 0x064B44, 0x36A537, 0x0EA54A, 0x86B2BF, 0x05AC53, 0x0AB647, 0x5936BC, 0x092E50, 0x0C9645, /*1991-2000*/
    0x4D4AB8, 0x0D4A4C, 0x0DA541, 0x25AAB6, 0x056A49, 0x7AADBD, 0x025D52, 0x092D47, 0x5C95BA, 0x0A954E, /*2001-2010*/
    0x0B4A43, 0x4B5537, 0x0AD54A, 0x955ABF, 0x04BA53, 0x0A5B48, 0x652BBC, 0x052B50, 0x0A9345, 0x474AB9, /*2011-2020*/
    0x06AA4C, 0x0AD541, 0x24DAB6, 0x04B64A, 0x69573D, 0x0A4E51, 0x0D2646, 0x5E933A, 0x0D534D, 0x05AA43, /*2021-2030*/
    0x36B537, 0x096D4B, 0xB4AEBF, 0x04AD53, 0x0A4D48, 0x6D25BC, 0x0D254F, 0x0D5244, 0x5DAA38, 0x0B5A4C, /*2031-2040*/
    0x056D41, 0x24ADB6, 0x049B4A, 0x7A4BBE, 0x0A4B51, 0x0AA546, 0x5B52BA, 0x06D24E, 0x0ADA42, 0x355B37, /*2041-2050*/
    0x09374B, 0x8497C1, 0x049753, 0x064B48, 0x66A53C, 0x0EA54F, 0x06B244, 0x4AB638, 0x0AAE4C, 0x092E42, /*2051-2060*/
    0x3C9735, 0x0C9649, 0x7D4ABD, 0x0D4A51, 0x0DA545, 0x55AABA, 0x056A4E, 0x0A6D43, 0x452EB7, 0x052D4B, /*2061-2070*/
    0x8A95BF, 0x0A9553, 0x0B4A47, 0x6B553B, 0x0AD54F, 0x055A45, 0x4A5D38, 0x0A5B4C, 0x052B42, 0x3A93B6, /*2071-2080*/
    0x069349, 0x7729BD, 0x06AA51, 0x0AD546, 0x54DABA, 0x04B64E, 0x0A5743, 0x452738, 0x0D264A, 0x8E933E, /*2081-2090*/
    0x0D5252, 0x0DAA47, 0x66B53B, 0x056D4F, 0x04AE45, 0x4A4EB9, 0x0A4D4C, 0x0D1541, 0x2D92B5            /*2091-2099*/
}; //阴历数据表

const int beginYear = 1901; //计算阴历的开始年份

const int days[2][12] = {
    {1, 32, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335},
    {1, 32, 61, 92, 122, 153, 183, 214, 245, 275, 306, 336}};

struct lunar
{
    int year, month, day; // 年月日
    bool isLeap; //是否为闰月，即阴历一年中多月的那一个月
    lunar() {}
    lunar(int y, int m, int d, bool isL) : year(y),
                                           month(m), day(d), isLeap(isL){};
};

bool isleap(int year)
{
    if ((year % 4 == 0 && year % 4 != 0) || (year % 400 == 0))
    {
        return 1;
    }
    return 0;
}//判断是否闰年

lunar getLunar(int y, int m, int d)
{
    lunar nowlunar(y, 0, 0, false); //默认初始化
    int index = y - beginYear;                        //年份下标
    int springMonth = (lunarYear[index] & 0x60) >> 5; //春节的月份
    int springDay = lunarYear[index] & 0x1f;          // 春节的日历
    int todaySolar = days[isleap(y)][m - 1] + d - 1;  //公历第几天
    int springLunar = days[isleap(y)][springMonth - 1] + springDay - 1;
    // 春节是公历年的第几天
    int gap = todaySolar - springLunar + 1; //计算公历这天与春节在公历的日期的间隔
    if (gap < 0) //今天在春节的前面 重新计算
    {
        index--;
        nowlunar.year--;
        if (index < 0)
            return nowlunar;                          //少于1901，不在计算范围内
        springMonth = (lunarYear[index] & 0x60) >> 5; //春节的月份
        springDay = lunarYear[index] & 0x1f;          // 春节的日历
        springLunar = days[isleap(y)][springMonth - 1] + springDay - 1;
        int yearToday = days[isleap(y - 1)][11] + 30;
        gap = todaySolar + yearToday - springLunar + 1;
    }
    int lunarMonth = 1;
    for (; lunarMonth <= 13; ++lunarMonth)
    {
        int monthDay = 29;
        if ((lunarYear[index] >> (6 + lunarMonth)) & 0x1)
            monthDay = 30;
        if (gap <= monthDay)
            break;
        gap -= monthDay;
    }
    nowlunar.day = gap;
    int leapMonth = (lunarYear[index] >> 20) & 0xf;
    if (leapMonth > 0 && leapMonth < lunarMonth)
    {
        lunarMonth--; //如果当前月为闰月，设置闰月标志
        if (lunarMonth == leapMonth)
            nowlunar.isLeap = true;
    }
    nowlunar.month=lunarMonth;
    return nowlunar;
}

int main(){
    int y,m,d;
    while(cin>>y>>m>>d){
        lunar nowlunar  = getLunar(y,m,d);
        printf("%d-%d-%d ",nowlunar.year,nowlunar.month,nowlunar.day);
    }
}