int get_yday(int year, int month, int day)
{
static int t[] = {0, 3, 2, 5, 0, 3, 5, 1, 4, 6, 2, 4};
year -= (month < 3);
return (year + year/4 - year/100 + year/400 + t[month-1] + day) % 7;
}
这个优雅的实现是Tomohiko Sakamoto 算法。其中比较难懂的有两点:
1.为什么year在month小于3的情况下要减一?
2.t数组的意义何在?
要弄懂这一点,首先我们要了解Disparate Gaussian公式,即
公式中的符号含义如下:
- w:星期(计算所得的数值对应的星期:0-星期日;1-星期一;2-星期二;3-星期三;4-星期四;5-星期五;6-星期六)
- c:年份前两位数
- y:年份后两位数
- m:月(m的取值范围为3至14,即在蔡勒公式中,某年的1、2月要看作上一年的13、14月来计算,比如2003年1月1日要看作2002年的13月1日来计算)(这里回答了为什么为什么year在month小于3的情况下要减一)
- d:日
- [ ]:称为高斯符号,代表向下取整,即,取不大于原数的最大整数。
- mod:同余(这里代表括号里的答案除以7后的余数)
把每一项offset除以7,余数刚好与t数组相同。我们可以推断出,t数组的意义就在计算这个month offsets
Tomohiko Sakamoto 算法对Disparate Gaussian公式进行了简化,他选择用一个不含有y(即年份后两位数)的形式(用c来表示y)来等价计算,从而实现少算一个量的目的
以下是等价性的证明: