基姆拉尔森计算公式 推导

需求：

给定一个xxxx-xx-xx日期，计算为星期几。

设定

int y;  //年
int m;  //月
int d;  //日
int w;  //周几
           

从 公元0年1月1日星期日 开始

推导

对于第一个月

w = (d-1) % 7  --------- 公式(1)
           

对于年

• 不考虑闰年

  在不考虑闰年的情况下，一年365天，365%7=1，就是说一年的第一天和最后一天是相同的。

  等价于，下一年的第一天星期几是会比这一年的最后一天+1的。

  完善公式(1)

w = (d-1 + y) % 7 --------- 公式(2)
           

• 考虑闰年

  因为闰年会多出来一天，所以相当于，计算当前年份前面有多少个闰年，将日期数w额外+1

  计算闰年的公式为:

y/4 - y/100 + y/400
           

结合之前的公式1，2

w = [d-1+y + (y-1)/4-(y-1)/100+(y-1)/400] % 7 -----公式(3)
           

对于其它月份

• 假设每个月都是28天

  因为28%7=0，也就是说每个月的w是相同的。

• 按正常月份计算

  一月是31天，比28多3天，也就是说，2月的w值，是应该比1月按28计算的往后推迟3天。

三月的值，因为二月刚好28天，不影响，相当于还是推后3天。

以此类推。

因为12月已是最后一个月，所以不用考虑12月的误差天数，同理，1月份的误差天数是0，因为前面没有月份影响它。

误差表

月	误差	累计	模7
1	3
2
4		6
5		8
		11
7		13
		16
9		19
10		21
		24
12	-	26

如果用一个数组记录就是

e[] = {0,3,3,6,1,4,6,2,5,0,3,5}
           

完善公式

w = [d-1+y + e[m-1] + (y-1)/4-(y-1)/100+(y-1)/400] % 7 --公式(4)
           

• 将闰年的情况考虑进去

  如果是闰年的话，2月之后的都会顺移一天

w = (d-1 + y + e[m-1] + (y-1)/4 - (y-1)/100 + (y-1)/400);
if(m>2 && (y%4==0 && y%100!=0 || y%400==0) && y!=0)
        ++w;
    w %= 7;
           

以上为基本推导过程

• 数学大佬对公式进行了优化
  • W= (d+2m+3(m+1)/5+y+y/4-y/100+y/400+1)%7