如何快速计算特定日期对应的星期 Tomohiko Sakamoto 算法解析

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）来等价计算，从而实现少算一个量的目的

以下是等价性的证明：