蔡勒公式如下：
W = [C/4] - 2C + y + [y/4] + [13 * (M+1) / 5] + d - 1
或者是:w=y+[y/4]+[c/4]-2c+[26(m+1)/10]+d-1
菜農補充“完整的蔡勒公式”:
W = (203 + C/4 - 2C + y + y/4 + 13 * (M+1) / 5 + d - 1) % 7
或者是:w=(203+y+y/4+c/4-2c+26(m+1)/10+d-1)%7
其中:C<=99,w在加203后不會出現負數
xwj 發表于 2008-9-6 20:12 侃單片機 ←返回版面    |
19樓： 考慮取值范圍，為了避免溢出，機器最后可以這樣計算： |
w=y +y/4 -c*7/4 +d +(13*m+8)/5 +N個7
|
|
當C<=99時，2*C<=198,為使蔡勒公式不出現負數，取N=29即N*7=203
故有蔡勒公式的完整表述：
W = (203+C/4-2*C+Y+Y/4+(13*M+8)/5+D)%7
| HotC51 發表于 2009-3-11 05:49 裸奔式實時操作系統HotTask51 ←返回版面 |
樓主： 菜農自賀5年后數組匯編82字節雷翻了dengm同志90字節的萬年歷匯編 |
/*-----------------------------------------------------------------------
函數Week(macro_U16_arg1, macro_U8_arg2, macro_U8_arg3)
入口:
macro_U16_arg1 年(0x00~0x9999)BCD碼
macro_U8_arg2 月(0x01~0x12)BCD碼
macro_U8_arg3 日(0x01~0x31)BCD碼
出口:
macro_U8_Result 星期(0~6)
長度: 82字節
-----------------------------------------------------------------------*/
macro_Function_Arg3 (Week, macro_U16_arg1, macro_U8_arg2, macro_U8_arg3)
{
#define Week_LOOP -20
asm_MOV_R1(4),//將百年、年、月、日的BCD碼換成二進制數
//Week_LOOP:
asm_MOV_A(0x4f),//DataBuffer偏移,程序長度
//分別取出macro_U16H_arg1,macro_U16L_arg1macro_U8_arg2, macro_U8_arg3
asm_MOVC_A_aDPTR(),
asm_MOV_R0_A(),
asm_MOV_A_iR0(),
asm_ANL_A(0x0f),
asm_XCH_A_iR0(),
asm_SWAP_A(),
asm_ANL_A(0x0f),
asm_MOV_B(10),
asm_MUL_AB(),
asm_ADD_A_iR0(),
asm_MOV_iR0_A(),
asm_INC_DPTR(),
asm_DJNZ_R1(Week_LOOP),//四次未完繼續
//(星期=百年%4*5+年+年/4+(13*月+8)/5+日)%7,特別注意1月2月的年為"去年"
asm_CJNE_Rn(macro_U8_arg2, 0x03, 0),//判斷月
asm_JNC(size_MOV_A_Rn() + size_ORL_A() + size_MOV_Rn_A() \
+ size_DEC_Rn() + size_CJNE_Rn() + size_DEC_Rn()),//月
//月<=2
asm_MOV_A_Rn(macro_U8_arg2),
asm_ORL_A(0x04),//1月2月同5月六月表
asm_MOV_Rn_A(macro_U8_arg2),
asm_DEC_Rn(macro_U16L_arg1),//1月2月的年為"去年"
asm_CJNE_Rn(macro_U16L_arg1, 0xff, size_DEC_Rn()),//Week_Start
asm_DEC_Rn(macro_U16H_arg1),
//Week_Start:
//百年%4*5
asm_MOV_A_Rn(macro_U16H_arg1),//百年
asm_ANL_A(0x03),
asm_MOV_Rn_A(macro_U16H_arg1),//百年%4
asm_CLR_C(),
asm_RLC_A(),
asm_CLR_C(),
asm_RLC_A(),//百年%4*4
asm_ORL_A_Rn(macro_U16H_arg1),//百年%4*5
asm_MOV_Rn_A(macro_U16H_arg1),//保存(百年%4*5),最大值3*5=15
//年+年/4
asm_MOV_A_Rn(macro_U16L_arg1),//年
asm_CLR_C(),
asm_RRC_A(),
asm_CLR_C(),
asm_RRC_A(),
asm_ADD_A_Rn(macro_U16L_arg1),//年+年/4
asm_MOV_Rn_A(macro_U16L_arg1),//保存(年+年/4)，最大值99+25=124
asm_MOV_A_Rn(macro_U8_arg2),//月
asm_MOV_B(13),
asm_MUL_AB(),//13*月->A
asm_MOV_B(8),
asm_ADD_A_B(),//13*月+8
asm_MOV_B(5),
asm_DIV_AB(),//(13*月+8)/5 商->A
asm_MOV_Rn_A(macro_U8_arg2),//保存((13*月+8)/5)最大值33
//日
asm_ADD_A_Rn(macro_U8_arg3),//日,最大值31
asm_ADDC_A_Rn(macro_U16H_arg1),//(百年%4*5),最大值15
asm_ADDC_A_Rn(macro_U16L_arg1),//(年+年/4)，最大值99+25=124
/*
此時ACC最大為15+124+33+31=203，
注意，蔡勒公式內無203.-2*C可能會越界,菜農給出的蔡勒公式完整表述：
W = (203+C/4-2*C+Y+Y/4+(13*M+8)/5+D)%7
菜農星期公式:
W = ( C%4*5+Y+Y/4+(13*M+8)/5+D)%7
*/
asm_MOV_B(7),//取余數%7
asm_DIV_AB(),
asm_XCH_A_B(),
asm_MOV_Rn_A(macro_U8_Result),//送返回值(星期)
macro_return(),
//數據區DataBuffer
asm_DB(macro_U16H_arg1),
asm_DB(macro_U16L_arg1),
asm_DB(macro_U8_arg2),
asm_DB(macro_U8_arg3)
};
|
|
注意：“此時ACC最大為15+124+33+31=203，”
再次從另外一路方法印證了蔡勒公式應該加203來防止溢出。
當然122年前人算是不會出現問題的，但直接用于編程肯定有bug |
|