求C语言农历转公历算法？

[Claude]

本人逆着公历转农历写了下，得到不同的结果，有正确，有错误。弄得有点头大，结果还是没有找出问题。所以……直接散分，求算法。

[珍惜生命远离CPP]

查表法

农历每隔多少年修订一次,你用啥公式算?

[jznhljg]

http://www.delphibbs.com/delphibbs/dispq.asp?lid=182846

[Claude]

是用查表法，没有用什么公式啊！就是硬算，

[waizqfor]

够复杂啊 顶吧

[lpqss1]

//由于有闰月，所以函数形式：
int toSolar(int ly, int lm, int ld, int *sy, int *sm, int *sd, int *sy2, int *sm2, int *sd2);


引理：公历y-m-d一定在农历y-m-d之前

由引理，我们的探测可以从公历y-m-d开始，步骤如下：

1)toLunar(y, m, d, &ly0, &lm0, &ld0);
2)估算ly0-lm0-ld0和ly-lm-ld的日子差 N = (ly-ly0)*12*29 + (lm-lm0)*29 + (ld-ld0);
3)如果 N > 0，把探测公历y-m-d加上这个天数， dateAddDays(&y, &m, &d, N);
4)如果 N = 0，如果(ly==ly0)&&(lm==lm0)&&(ld==ld0)得到一个解sy=y;sm=m;sd=d; (如果此时sy已经被赋值过，则sy2=y;sm2=m;sd2=d;返回退出)
否则表示出现误差，修正N = 1;dateAddDays(&y, &m, &d, N);
从而开始探测闰月,令N = 30;dateAddDays(&y, &m, &d, N);
5)如果 N < 0，且N >= -30，肯定是上次探测N+29不够（该月大），应该+30，
所以令N = 1; dateAddDays(&y, &m, &d, N);



经过测试，最多需要3次可以命中，还是比较理想的算法了:)

PS:经过完全测试，公元1500-2100年没有问，最多尝试次数为3次

#define LUNAR_TO_SOLAR_STEPS 3 // 公历转公历搜索最长步数
// 农历返回公历
// 注意：这是我在掌心万年历中写的代码，版权属于我所有，仅供个人学习使用，严禁在任何商业软件中使用本代码
// 输入 农历年-月-日 ly-lm-ld (lm < 0表示闰月)
// 输出 公历年-月-日 y-m-d
/*
公元1500年前的情况太特殊,比如公元9-1月,缺少农历11月,直接导致计算ldiff错误...
所以toSolar只支持1500年以后的转换. --ZV
*/
int toSolar(int LunarYear, int LunarMonth, int LunarDay)
{
if (LunarYear < 1500) { // 非常复杂的问题
return CALENDAR_ERROR;
}

int ret;
int y = LunarYear;
int m = abs(LunarMonth);
int d = LunarDay;

int i, ldiff = 0, lunDate;
int ly, lm, ld;

validateDate(&y, &m, &d);
for (i = 0; i < LUNAR_TO_SOLAR_STEPS; i++) {
lunDate = toLunar(y, m, d);
DATE_TO_YMD(lunDate, ly, lm, ld);
if (abs(lm) == abs(LunarMonth) && ld == LunarDay) {
if (lm == LunarMonth) {
ret = YMD_TO_DATE(y, m, d);
return ret;
}
// 如果找到一个，则开始探测闰月
ldiff = 30; // 要么直接转到，要么多了一天而已
} else {
ldiff = (LunarYear - ly) * 12 * 29 + (abs(LunarMonth) - abs(lm)) * 29 + (LunarDay - ld);
if (ldiff == 0) { // 稍做修正
ldiff = 1;
}
}
dateAddDays(&y, &m, &d, ldiff); // 注意可能出现负数
}
return CALENDAR_ERROR;
}

[潇湘秦]

农历经常变，还经常有闰月，大进，小进什么啊？
没有办法算吧!

[lbh2001]

转载一个

/*------------农历转换函数-----------*/
char *GetDayOf(PSYSTEMTIME pSt)
{ 
 /*天干名称*/
 const char *cTianGan[] = {"甲","乙","丙","丁","戊","己","庚","辛","壬","癸"};
 /*地支名称*/
 const char *cDiZhi[] = {"子","丑","寅","卯","辰","巳","午","未","申","酉","戌","亥"};
 /*属相名称*/
 const char *cShuXiang[] = {"鼠","牛","虎","兔","龙","蛇","马","羊","猴","鸡","狗","猪"};
 /*农历日期名*/
 const char *cDayName[] = {"*","初一","初二","初三","初四","初五",
       "初六","初七","初八","初九","初十","十一","十二","十三","十四","十五","十六","十","十八","十九","二十","廿一","廿二","廿三","廿四","廿五","廿六","廿七","廿八","廿九","三十"};
 /*农历月份名*/
 const char *cMonName[] = {"*","正","二","三","四","五","六","七","八","九","十","十一","腊"};

 /*公历每月前面的天数*/
 const int wMonthAdd[12] = {0,31,59,90,120,151,181,212,243,273,304,334};
 /*农历数据*/
 const int wNongliData[100] = {2635,333387,1701,1748,267701,694,2391,133423,1175,396438
       ,3402,3749,331177,1453,694,201326,2350,465197,3221,3402
       ,400202,2901,1386,267611,605,2349,137515,2709,464533,1738
       ,2901,330421,1242,2651,199255,1323,529706,3733,1706,398762
       ,2741,1206,267438,2647,1318,204070,3477,461653,1386,2413
       ,330077,1197,2637,268877,3365,531109,2900,2922,398042,2395
       ,1179,267415,2635,661067,1701,1748,398772,2742,2391,330031
       ,1175,1611,200010,3749,527717,1452,2742,332397,2350,3222
       ,268949,3402,3493,133973,1386,464219,605,2349,334123,2709
       ,2890,267946,2773,592565,1210,2651,395863,1323,2707,265877};
 static int wCurYear,wCurMonth,wCurDay;
 static int nTheDate,nIsEnd,m,k,n,i,nBit;
 TCHAR szNongli[30], szNongliDay[10],szShuXiang[10];
 /*---取当前公历年、月、日---*/
 wCurYear = pSt->wYear;
 wCurMonth = pSt->wMonth;
 wCurDay = pSt->wDay;
 /*---计算到初始时间1921年2月8日的天数：1921-2-8(正月初一)---*/
 nTheDate = (wCurYear - 1921) * 365 + (wCurYear - 1921) / 4 + wCurDay + wMonthAdd[wCurMonth - 1] - 38;
 if((!(wCurYear % 4)) && (wCurMonth > 2))
  nTheDate = nTheDate + 1;



 /*--计算农历天干、地支、月、日---*/
 nIsEnd = 0;
 m = 0;
 while(nIsEnd != 1)
 {
  if(wNongliData[m] < 4095)
   k = 11;
  else
   k = 12;
  n = k;
  while(n>=0)
  {
   //获取wNongliData(m)的第n个二进制位的值
   nBit = wNongliData[m];
   for(i=1;i<n+1;i++)
    nBit = nBit/2;

   nBit = nBit % 2;

   if (nTheDate <= (29 + nBit))
   {
    nIsEnd = 1;
    break;
   }

   nTheDate = nTheDate - 29 - nBit;
   n = n - 1;
  }
  if(nIsEnd)
   break;
  m = m + 1;
 }
 wCurYear = 1921 + m;
 wCurMonth = k - n + 1;
 wCurDay = nTheDate;
 if (k == 12)
 {
  if (wCurMonth == wNongliData[m] / 65536 + 1)
   wCurMonth = 1 - wCurMonth;
  else if (wCurMonth > wNongliData[m] / 65536 + 1)
   wCurMonth = wCurMonth - 1;
 }



 /*--生成农历天干、地支、属相 ==> wNongli--*/
 wsprintf(szShuXiang,"%s",cShuXiang[((wCurYear - 4) % 60) % 12]);
 wsprintf(szNongli,"%s(%s%s)年",szShuXiang,cTianGan[((wCurYear - 4) % 60) % 10],cDiZhi[((wCurYear - 4) % 60) % 12]);



 /*--生成农历月、日 ==> wNongliDay--*/
 if (wCurMonth < 1)
  wsprintf(szNongliDay,"闰%s",cMonName[-1 * wCurMonth]);
 else
  strcpy(szNongliDay,cMonName[wCurMonth]);



 strcat(szNongliDay,"月");
 strcat(szNongliDay,cDayName[wCurDay]);
 return strcat(szNongli,szNongliDay);



}

[hjzwl1018]

顶。

[潇湘秦]

在中国农历论坛 看到的一些算法的讨论 转过来 看看

首先，我们要确定一个时刻，作为一天的起点（包括这个时刻）。然后我们以月亮通过朔望交界点的那一天作为每个月的第一天，以太阳通过雨水点的那个月作为正月，依次的，以太阳通过各中气点的那个月作为二月、三月等等。
经过长期的观察之后，我们就会发现每个月不是29天就是30天，每年不是12个月就是13个月，每个中气点唯一对应一个月，但是有的月却没有中气。
为了制定我们的历法，我们至少要知道足够精度的朔望月长度，记为ML，和两个中气的间隔时间记为YL，这两个值经过前人长期不懈的观测和计算，我们知道分别为29.5306日和30.4377日。接下来，我们要随便找到某一年，以太阳通过雨水点的时刻到那一年开始的那个时刻的时间，记为DY，以日为单位，还要找出这一年正月里，月亮通过朔望交界点的时刻到那一月开始的那个时刻的时间，记为DM，以日为单位。这两个值是要通过实地观测得出来的。
这个时候，我们就可以通过递推来制订我们的历法，以预测未来月球和太阳的运行情况，并进一步预测月相的变化和气候的变化。

在本算法中我们使用了三个近似处理：
每天的长度总是一样的；
每个朔望月的长度总是一样的；
每两个相邻中气的时间总是一样的。

接下来，我们就逐月推算每月的大小：
Select Case Fix(DM + ML)
Case is = 29
The month have 29 days.
The DM of next month is DM + ML - 29.
Case is = 30
The month have 30 days.
The DM of next month is DM + ML - 30.
End Select

每年中各月的月份、是否为闰月
For i = 1 to 12
This month has A days.
If DY>=A Then
This month is an Extra Month.
DY = DY - A
Else
This month is the i Month.
DY = DY + YL - A
i = i + 1
End If
Next

这就是我们的核心算法。通过它，我们可以从任何一个包含中气的月份开始向后推算各月各年的情况。
由于我们所用的DM、DY、ML、YL都是近似值，若干年后，我们的历法必有偏差，我们只要重新观测，更新DM、DY、ML、YL的值就可以了。事实上，诸朝气数最多不过400年，在这些值精确到万分位的时候，通行一朝是完全可以的。

下面对大小月的排列和闰月的出现作一些定性分析
显然0≤DM＜1，那么
DM＜0.4694时，当月为小月，下月为大月
DM≥0.4694时，当月为大月
DM＜0.4694 * 2时，下月为小月，下下月为大月
DM≥0.4694 * 2时，下月亦为大月，下下月为小月
所以，大小月一般是交替出现，有时会出现连续两个月是大月的情况，但不会出现连续三个月或更多的是大月的情况，也不会出现连续两个月更多的是小月的情况。也就是说，一般是小月大月相继出现一段时间后，接着就出现一个大月，然后又是小月大月相继出现，具体的间隔也不定，大家可以通过计算机推算后，找点规律出来。

一年中最少有十二个月，因为YL * 12 （代表十二个中气循环一次的时间）＞ ML * 12 （代表十二个朔望月的时间） + 1 （第一个月的DM总小于1），最多有十三个月，因为30（第一个月的DY总小于30） + YL * 12 ＜ ML * 14。两个相邻中气的间隔时间大于30，所以一个月不可能有两个中气，那么一年有十二个月的时候，12个中气必然刚好依次各居一个月，这一年就称为平年，一年有十三个月的时候，必然有且仅有一个月没有中气，那个月就时闰月，这一年就称为闰年。平年的时候，根据大小月的排列规则，一年最多有8个大月，最少有6个大月，因为ML * 12 + 1 ＜ 30 * 8 + 29 * （12 - 8），所以不可能有8个大月，最多为7个大月，那么一年就为354或355天。闰年的时候，根据大小月的排列规则，一年最多有9个大月，最少有6个大月，因为ML * 13 + 1 ＜ 30 * 8 + 29 * （13 - 8），所以不可能有8个或更多的大月，最多为7个大月，那么一年就为383或384天。

因为每月的天数不固定，所以闰月的出现规律也不好讨论，连续十二个月的天数也不固定，所以闰年的出现规律也不好讨论。不过可以肯定是，每年的第一个月肯定不是闰月，因为我们是以包含中气雨水的月份作为正月，也就是每年的第一个月的。所以，闰月必定出现在某个有名字的月份后面，出现在哪个后面就叫闰哪月。

[ARLENE_YANG]

up

[Claude]

谁有现成的C程序，可以发到我的邮箱[e-mail]，参考参考……

[bfhtian]

没有什么现成的公式吧

[lpqss1]

http://lpqss1.download.csdn.net/