菜农星期公式揭秘及蔡勒公式推想

2万 · 2008-9-7 10:04

讨论及论证过程详见：https://bbs.21ic.com/club/bbs/list.asp?boardid=11&page=1&t=3063061

今日方敢自称菜农星期算法为“菜农星期公式”~~~

因为“菜农公式”竟然与“蔡勒公式”结果完全相同！
哈哈~~~不同菜农就栽了~~~
结果当然必须相同，但过程和推导方法虽然雷同，可公式的优化和“物理意义”及“路径”差距很大。
菜农不知蔡勒（Zeller）是何须人也，也不知是“什么朝代”人士，
也不想知道~~~估计也是“卖脑浆的主”~~~
首先对比2个算法：
蔡勒公式如下：
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后不会出现负数

菜农公式如下：
W = (C%4*5+Y+Y/4+(13*M+8)/5+D)%7
菜农的(13*M+8)/5是根据月星期表推导出的。
菜农依据表的万年某日星期算法：
   公式：某日星期几=(百年%4*5天+年+年/4+月星期表+日+2天)%7
   通用星期偏差表
   闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  29  31  30  31  30  31  31  30  31   30   31
   星期   3   6   0   3   5   1   3   6   2   4    0    2
   平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  28  31  30  31  30  31  31  30  31   30   31
   星期   4   0   0   3   5   1   3   6   2   4    0    2
或(1月2月为5月6月星期表)：
   公式：某日星期几=(百年%4*5天+年+年/4+月星期表+日-0天)%7
   通用星期偏差表(此表是菜农公式中的(13*M+8)/5)
   闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  29  31  30  31  30  31  31  30  31   30   31
   星期   5   1   2   5   0   3   5   1   4   6    2    4
   平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  28  31  30  31  30  31  31  30  31   30   31
   星期   6   2   2   5   0   3   5   1   4   6    2    4

故蔡勒用的星期表为(1月2月为13月14月星期表)：
   公式：某日星期几=(百年/4-2*百年+年+年/4+月星期表+日-1天)%7
   蔡勒星期偏差表(此表是蔡勒公式中的13*(M+1)/5)
   闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月 13月 14月
   天数   31  29  31  30  31  30  31  31  30  31   30   31   31   29
   星期   5   1   2   5   0   3   5   1   4   6    2    0    3    5
   平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月  13月 14月
   天数   31  28  31  30  31  30  31  31  30  31   30   31   31   28
   星期   6   2   2   5   0   3   5   1   4   6    2    0    3    5

根据21ic网友xwj的推导：

xwj 发表于 2008-9-6 19:56 侃单片机 ←返回版面

18楼：典型的1阶比例计算，公式可以这样简化：

蔡勒公式如下：
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=y+[y/4]+[c/4]-2c+[26(m+1)/10]+d-1
-->>
w=(5/4)y -(7/4)c +(26/10)m +26/10 +d -1
-->>
w=(5/4)y -(7/4)c  +d +(26/10)m +16/10
-->>
w=(5/4)y -(7/4)c  +d +(26m+16)/10

-->>
w=(5/4)y -(7/4)c  +d +(26m+16)/10 +N个7

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

注意反推(13*M+8)/5+D到13*（M+1)/5+D-1
(13*M+8)/5+D=(13*M+8+5-5)/5+D
=(13*M+13)/5-5/5+D
= 13*（M+1)/5+D-1

对比菜农公式：
W = (C%4*5+Y+Y/4+(13*M+8)/5+D)%7
难道(203+C/4-2*C)%7和(C%4*5)%7相等吗？？？
编过计算机程序的人一眼可以看出C%4*5的表述实际为：
(C&3)*5
或((C&3)<<2)|(C&3).
即把C的值取低2位(D1D0)复制到高2位(D3D2)
因为C&3=0~3或00B~11B
所以C%4*5只有4个取值0(0000B),5(0101B),10(1010B),15(1111B)
则((C%4)*5)%7对应0,5,3,1
可以看出C%4*5的表述和“物理意义”非常清晰：

hotpower 发表于 2008-9-6 09:34 侃单片机 ←返回版面

12楼：这个算法实际是把月星期表进行了代换

月星期表+日+2天等效为公式 (13 * m + 8)/ 5

再“转帖”我当时的“四百年历”的帖子：
推导过程：
1。平年365天（52周+1天），闰年366天（52周+2天）。平年2月28天，闰年2月29天。
   由于公元1月1日设为星期六，故3月1日为星期三。         ——注意这个“三”
   为使算法达到最简，故本算法以“星期”为计算单位。且选3月1日为基月。
2。每400年整一闰，或每4年且不为百年的一闰。
   百年%4=0闰或（年%4=0并且年<>0）闰。
3。每 4年（3个平年+1个闰年）共208周+5天                  ——注意这个“5天”
   每百年共100*（208周+5天）-1天=5217周+5天              ——注意这个“5天”(整
百年暂设为平年)
   每400年共4*（5217周+5天）+1天（整400年闰）=20871周+0天——注意这个“0天”和
“1天”(4个整百年只有一个闰年)
   即400年一轮回！（原来万年历400年前是一家）
4。万年3月1日星期算法(记住那个三)
   公式：某年3月1日星期几=(3天+百年%4*5天+年/4*5天+年%4+月星期表+日-1天)%7
         某年3月1日星期几=(百年%4*5天+年/4*5天+年%4+月星期表+日+2天)%7
   或    某年3月1日星期几=(百年%4*5天+年+年/4+月星期表+日+2天)%7

   闰4百年3月1日星期算法（百年%4=0）
   公式：某年3月1日星期几=(年+年/4+月星期表+日+2天)%7
   例：0000年3月1日星期几=(0+0/4+0+1+2)%7=3%7=星期三
       1600年3月1日星期几=(0+0/4+0+1+2)%7=3%7=星期三
       2000年3月1日星期几=(0+0/4+0+1+2)%7=3%7=星期三
       2001年3月1日星期几=(1+1/4+0+1+2)%7=4%7=星期四
       2004年3月1日星期几=(4+4/4+0+1+2)%7=8%7=星期一
       2008年3月1日星期几=(8+8/4+0+1+2)%7=13%7=星期六
       2042年3月1日星期几=(42+42/4+0+1+2)%7=55%7=星期六

   平4百年3月1日星期算法（百年%4<>0）
   公式：某年3月1日星期几=(百年%4*5天+年+年/4+月星期表+日+2天)%7
   例：1700年3月1日星期几=（17%4*5+0+0/4+0+1+2）%7=8%7=星期一(注意：1700年是平
年)
       1800年3月1日星期几=（18%4*5+0+0/4+0+1+2）%7=13%7=星期六(注意：1800年是平
年)
       1900年3月1日星期几=（19%4*5+0+0/4+0+1+2）%7=18%7=星期四(注意：1900年是平
年)
       1901年3月1日星期几=（19%4*5+1+1/3+0+1+2）%7=19%7=星期五
       1918年3月1日星期几=（19%4*5+18+18/4+0+1+2）%7=(15+22+3)%7=40%7=星期五
       1958年3月1日星期几=（19%4*5+58/4*5+58%4+3）%7=(15+70+2+3)%7=90%7=星期六
       1988年3月1日星期几=（19%4*5+88/4*5+88%4+3）%7=(15+110+0+3)%7=128%7=星期二
       1999年3月1日星期几=（19%4*5+99/4*5+99%4+3）%7=(15+120+3+3)%7=141%7=星期一
       2100年3月1日星期几=（21%4*5+0/4*5+0%4+3）%7=(5+0+0+3)%7=8%7=星期一(注意：
2100年是平年)
       2101年3月1日星期几=（21%4*5+1/4*5+1%4+3）%7=(5+0+1+3)%7=9%7=星期二
       2102年3月1日星期几=（21%4*5+2/4*5+2%4+3）%7=(5+0+2+3)%7=10%7=星期三
       2103年3月1日星期几=（21%4*5+3/4*5+3%4+3）%7=(5+0+3+3)%7=11%7=星期四
       2104年3月1日星期几=（21%4*5+4/4*5+4%4+3）%7=(5+1+0+3)%7=9%7=星期二(注意：
2104年是闰年)
       9999年3月1日星期几=（99%4*5+99/4*5+99%4+3）%7=(120+15+3+3)%7=141%7=星期一
   注：按400年一轮回！（400年前是一家）的说法
       1600年，2000年是一样的
       1700年，2100年是一样的
       1800年，2200年是一样的
       1900年，2300年是一样的
       所以，万年历实际上是骗人的，应该叫“四百年历”才是。万只是多的意思罢了。
5. 万年某日星期算法
   公式：某日星期几=(百年%4*5天+年+年/4+月星期表+日+2天)%7
   通同星期偏差表
   闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  29  31  30  31  30  31  31  30  31   30   31
   星期   3   6   0   3   5   1   3   6   2   4    0    2
   平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
   天数   31  28  31  30  31  30  31  31  30  31   30   31
   星期   4   0   0   3   5   1   3   6   2   4    0    2

可以看出：
月星期表实际是根据(闰年/平年,天数)得到星期数的。
它实际是个固定的数，只要分清闰年/平年即可。

新算法：星期=(百年%4*5+年+年/4+(13*月+8)/5+日)%7
旧算法：星期=(百年%4*5天+年+年/4+月星期表+日+2天)%7

昨日首次“拜读”了蔡勒公式，对-2C一直迷茫~~~
既然((C%4)*5)%7对应0,5,3,1
那么(203+C/4-2*C)也对应0,5,3,1？？？
C=0, (203+0/4-2*0)%7=203%7=0
C=1, (203+1/4-2*1)%7=(203-2)%7=5
C=2, (203+2/4-2*2)%7=(203-4)%7=3
C=3, (203+3/4-2*3)%7=(203-6)%7=1
C=4, (203+4/4-2*4)%7=(203+1-8)%7==(203+0)%7=0
C=5, (203+5/4-2*5)%7=(203+1-10)%7=(203-2)%7=(203+5)%7=5
C=6, (203+6/4-2*6)%7=(203+1-12)%7=(203-4)%7=(203+3)%7=3
C=7, (203+7/4-2*7)%7=(203+1-14)%7=(203+1-0)%7=(203+1)%7=1
C=8, (203+8/4-2*8)%7=(203+2-16)%7=(203-14)%7=(203+0)%7=0
晕呀！！！
蔡勒公式中的-2C竟然是想得到重复的0,5,3,1序列采取的“手段”

最后还有个1月2月变换问题，蔡勒公式取13月14月，菜农取5月6月。

因为菜农公式是在月星期表“进化”而来，又因为13月14月与5月6月的表值相同，
所以最后决定改为1月2月变换为5月6月。

这样做现在看来“物理意义”失去，但当时考虑“月表”的M大于12就不成“月”了.

总之把握：

1月2月为去年的以3月为基月表的表尾。

所以:

公式中出现1月2月的年是“去年”！！！

公式中出现1月2月的月是“明年”！！！(13月14月下一年的1月2月)

这样“去年”和“明年”就“能量守恒”了~~~

到此菜农“反汇编”蔡勒公式终成正果~~~
感觉蔡勒公式确实和菜农公式在思路上有很大的差距~~~
哈哈~~~看来蔡勒“同志”的脑浆还是“白颜色”的~~~

最后给出完整的菜农公式：

W = (C%4*5+Y+Y/4+(13*M+8)/5+D)%7

其中:

C=1~XX,Y=0~99,M=1~12,D=1~28/29/30/31

计算时要注意1月和2月的“去年”“明年”变换。

即：

当M=1或2时，M=M+12或M=M+4,Y=Y-1;当结果Y<0时，C=C-1.

哈哈~~刚搜索的：

蔡勒公式最早是由德国数学家克里斯蒂安 · 蔡勒（ Christian Zeller, 1822-1899 ）在 1886 年推导出的，因此通称为蔡勒公式（ Zeller’s Formula ）。为方便口算，式中的 [13 * (M+1) / 5] 也往往写成 [26 * (M+1) / 10] 。

哈哈~~~122年前就有蔡勒公式，菜农深感惭愧~~~

向蔡勒同志敬礼！！！

菜农HotPower@126.com 2008.9.6 23:18 于西安大雁塔菜地
相关链接:http://blog.**/hotpower/158297/message.aspx

2万 · 2008-9-7 17:13

首先声明菜农不喜欢数学，更不是数学家~~~

刚找到当年(13*M+8)/5)推导过程的“手稿”。它要从月星期表说起

通用星期偏差表(此表是菜农公式中的(13*M+8)/5)
分组           大  小  大  小  大  大  小   大   小   大   大   小  大  小  大
闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月 13月 14月 3月 4月 5月
天数   31  29  31  30  31  30  31  31  30  31   30   31   32   29  32  30   31
星期   5   1   2   5   0   3   5   1   4   6    2    4     0    3   5   1    3
星期差         3   2   3   2   3   3   2   3    2    3     3    2   3   2    3

分组           大  小  大  小  大  大  小   大   小   大   大   小  大  小  大
平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月 13月 14月 3月 4月 5月
天数   31  28  31  30  31  30  31  31  30  31   30   31   32   28  32  30   31
星期   6   2   2   5   0   3   5   1   4   6    2    4     0    3   5   1    3
星期差         3   2   3   2   3   3   2   3    2    3     3    2   3   2    3

由于明年的1月(13月)和今年5月的月表值都是0，2月(14月)和今年6月的月表值都是3
故月表简化为：
通用星期偏差表(此表是菜农公式中的(13*M+8)/5)
分组           大  小  大  小  大  大  小   大   小   大
年           3月 4月 5月 6月 7月 8月 9月 10月 11月 12月
                      1月 2月
天数           31  30  31  30  31  31  30  31   30   31
                       31  28/29
星期值         2   5   0   3   5   1   4   6    2    4
星期差         3   2   3   2   3   3   2   3    2    3

可以看出"分组"为两组(大月小月大月小月大月)，每5(记住这个"5")个月1组(注意1月2月划归为5月6月):
3月 4月 5月 6月 7月和 8月 9月 10月 11月 12月
星期差同为：3   2   3   2   3

星期差求和=3+2+3+2+3=13(记住这个"13")

再对13个月的星期值求和=0+3+2+5+0+3+5+1+4+6+2+4=40
由于每5个月1组，则40/5=8(记住这个"8")

因为原查表算法为：
某日星期几=(百年%4*5+年+年/4+月星期表+日)%7

所以要琢磨出:
    月星期表 = ((13*M+8)/5)%7

因为每5个月1组，故肯定有((aM+b)/5)%7或((aM+b)%5)%7存在。
但是((aM+b)%5)%7可化简为(aM+b)%5，结果只有0,1,2,3,4.这肯定不过！！！

所以((aM+b)/5)%7将是推导的框架！！！

哈哈~~~"5"用到了。该"13"和"8"上阵了。

参数只有a,b.代入((aM+b)/5)%7几次即可搞定((13*M+8)/5)%7~~~

哈哈~~~菜农没文化让大家见笑了~~~搞笑~~~绝对的"非典"

菜农的"推导过程"绝对气晕数学教授，不过俺数学不好也难为菜农胡思乱想了~~~

只看该作者 · 2008-9-7 17:26

估计是知道公式.

2万 · 2008-9-7 19:31

没想搞什么"公式",自己玩玩而已.

本不再想上网,准备"养老",但看见有求助帖https://bbs.21ic.com/club/bbs/list.asp?boardid=11&page=1&t=3063061

本不想参与,看来楼主用的是老的随手写的算法,并说有错...

现在俺还没时间验证.因为有新的查表算法和"公式",所以不想再追究"错事"~~~

后来搜索了有关星期的算法,发现真有些"怪异",故就用心写了篇"作文"~~~

有喜欢研究星期算法的朋友可以点击http://blog.**/hotpower/26788/category.aspx

即可看到菜农从2004年探讨星期算法的全过程的帖子(算法推导及程序等)

又该到"养老"的时候了~~~

朋友们再见~~~

1万 · 2008-9-7 20:29

1万 · 2008-9-8 22:52

   1年只需12个字节，100年只需1200字节，如果用压缩BCD只需600字节，再找规律压缩一下，更省空间。那些复杂的算法可能不止600字节哦。
   月历表中1号前有几个空格（即偏差值，星期天不算），日期除以7的余数，加上空格数，就是星期数，太简单了。

    例：2008年9月8日  本月月历表中1日前面空格为0（星期天不算)，8%7=1，1+0=1，即是星期一；
   2032年9月8日  本月月历表中1日前面空格为2（星期天不算)，8%7=1，1+2=3，即是星期三；

   公式：   日期%7 + 偏差  (偏差来自表格)

   再变一下，取余数除法都不需要，大于等于7就减去7......

   偶读初中的时候就用这个方法“在同学面前耍聪明”


   修正一下用词： 1号的前一天是星期几，偏差就是几天。(星期天是0)

2万 · 2008-9-9 00:30

某日星期几=(百年%4*5+年+年/4+月星期表+日)%7

通用星期偏差表(此表是菜农公式中的(13*M+8)/5)
闰年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月 13月 14月 3月 4月 5月
天数   31  29  31  30  31  30  31  31  30  31   30   31   32   29  32  30   31
星期   5   1   2   5   0   3   5   1   4   6    2    4     0    3   5   1    3

平年  1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月 13月 14月 3月 4月 5月
天数   31  28  31  30  31  30  31  31  30  31   30   31   32   28  32  30   31
星期   6   2   2   5   0   3   5   1   4   6    2    4     0    3   5   1    3

例如2008年的“骗术”

根据星期表查表公式：

某日星期几=(百年%4*5+年+年/4+月星期表+日)%7
          =(20%4*5+08+08/4+月星期表+日)%7
          =(0+10+月星期表+日)%7
          =(3+月星期表+日)%7
查2008.9.8

因为9月表=4，故(3+月星期表+日)%7=(3+4+日)%7=(0+日)%7
即9月的日就是“星期”

所以你就会“聪明”地立即“说出”2008.9.8是星期一(8=7+1=0+1=1)

2008.1.8如何查表？？？

这里有2种查法。
注意顶部的通用星期偏差表的1月2月和13月14月,
1月2月需要知道今年是否闰年，13月14月不需要，但要将年-1。

1.不需知道闰年的查法
  因为是1月，需要变换为13月或5月，且年变换为2007年。
  即需要查2007年的13月的8日是星期几。
某日星期几=(百年%4*5+年+年/4+月星期表+日)%7
          =(20%4*5+07+07/4+月星期表+日)%7
          =(0+8+月星期表+日)%7
          =(1+月星期表+日)%7
          =(1+13月星期表+日)%7
          =(1+0+日)%7
故2008.1.8是(1+8)%7=(2+7)%7=2

2.需知道闰年的查法
因为2008是闰年，故需要查闰年表
某日星期几=(百年%4*5+年+年/4+月星期表+日)%7
          =(20%4*5+08+08/4+月星期表+日)%7
          =(0+10+月星期表+日)%7
          =(3+1月星期闰年表+日)%7
          =(3+5+日)%7
          =(1+日)%7
故2008.1.8是(1+8)%7=(2+7)%7=2

哈哈~~~这种“心算”骗人的把戏就是当知道年月时，日+偏差就是星期~~~

1万 · 2008-9-9 11:18

这种骗人术在一些低端CPU系统，很实用，仅需要查表、判断、加减法运算即可。且一般的产品(非日历类产品)，设置10年就可以了。

2万 · 2008-9-9 19:06

只看该作者 · 2008-9-9 22:29

哦！　头疼　这些没有物理意义的东东了　

2万 · 2008-9-9 22:37

只看该作者 · 2008-9-10 01:06

2万 · 2008-9-10 01:13

哈哈~~~欢迎大家到文潭漫游---菜农公式栏目

1万 · 2008-9-10 20:19

    15分钟前的、2天内的、某时间段的等等时间运算。原来直接用年月日时分秒来运算，真的很繁琐。你要知道，有时候增加一秒钟，年月日时分都得变。要找15分钟前的数据，有可能跨小时、跨天、跨月、跨年......
    后来，把年月日时分秒统一换算为秒(定义2000年1月1日00:00:00为第0秒)，然后用秒运算，简单多了。
    在没有硬件RTC的系统，如果需要RTC，直接用秒累加，非常方便，需要的时候，把秒化成年月日时分秒就可以了。

2万 · 2008-9-14 09:48

此时月末，年末特别重要。
故此时日期翻转后可能还需年翻转。

只看该作者 · 2011-3-22 10:55

看懂了。。。

sInitTime.u32cDayOfWeek =(sInitTime.u32Year/100%4*5+sInitTime.u32Year%100+sInitTime.u32Year%100/4+(13*sInitTime.u32cMonth+8)/5+sInitTime.u32cDay)%7;

只看该作者 · 2011-3-30 18:56

非常不错，另类思维

只看该作者 · 2011-3-30 19:03

只看该作者 · 2011-3-31 17:41

有意思

2万 · 2011-4-1 13:56

哈哈，挖墓？？？

菜农星期公式揭秘及蔡勒公式推想

相关下载

倒塌的(13*M+8)/5推导过程

以前看电视上采访一个天才小神童,专算这个的.

本来只是自己在GPS项目中和LPCARM的RTC模块上解决星期问题

哈哈~~~厉害……俺还是查表再计算吧~~~

实际上用查表法很简单

“在同学面前耍聪明”~~~实际在骗人~~~

哈哈，是骗人术

哈哈~~~日历本身已隐藏了%7~~~

！ 头疼

日历算法对RTC调整和GPS时区日翻转及跨月跨年等都有帮助

非典中的非典

哈哈~~~菜农从外星球带来的怪物21ic腐败的"斑竹"们是看不明

说到跨年跨月，我想起来当年开发汽车行驶记录仪的时候

伦敦时间下午4点整东八区日期就要翻转

永恒之光奖章

技术导师奖章

七世轮回

坚毅之洋流

核心会员奖章

沉静之湖泊

以坛为家

技术奇才奖章

无冕之王奖章

湍急之河流

缘定三生

经常做客

技术新星奖章

涓涓之细流

哈哈~厉害……俺还是查表再计算吧~

“在同学面前耍聪明”~实际在骗人~

哈哈~日历本身已隐藏了%7~

！　头疼