ツェラーの公式

テンプレート:出典の明記 ツェラーの公式（ツェラーのこうしき、テンプレート:Lang-en-short、テンプレート:Lang-de-short）とは西暦（グレゴリオ暦またはユリウス暦）の年・月・日から、その日の曜日を算出する公式である。ドイツの数学者であるクリスティアン・ツェラー (テンプレート:Interlang) が1882年に考案した。

特定の日からの通算日数（ユリウス通日やRata Dieなど）を求め、その剰余から曜日を求める計算と本質は同じである。

呼称および表記

Zeller's congruence は直訳すると「Zeller の合同式」といった意味になるが、日本では主に「Zeller の公式」と呼ばれる。またカタカナ表記はドイツ語読みで「ツェラーの公式」と記すのが一般的だが、英語読みで「ゼラーの公式」と表記されることもある。

公式

テンプレート:Mvar 年 テンプレート:Mvar 月 テンプレート:Mvar 日の曜日を求める。

ただし、1月と2月は、前年のそれぞれ13月・14月として扱う。たとえば、2025年1月1日・2月1日は、2024年13月1日・14月1日とする。また、紀元前 $\hat{y}$ 年は西暦 $1-\hat{y}$ 年として扱う（天文学的紀年法）。たとえば、紀元前1年・前2年・前3年は、0年・－1年・－2年となる。

曜日	日	月	火	水	木	金	土
h	1	2	3	4	5	6	0

曜日 テンプレート:Mvar は次の式で求められる（ テンプレート:Mvar は、0 から 6 でそれぞれ土曜日から金曜日を表す）：

$h=(d+\lfloor \frac{26(m+1)}{10}\rfloor +Y+\lfloor \frac{Y}{4}\rfloor +\mathrm{\Gamma })mod7$

$\mathrm{\Gamma }$ はグレゴリオ暦 (テンプレート:En) かユリウス暦 (テンプレート:En) かで変わる項で、

$\mathrm{\Gamma }=\{\begin{array}{cc}-2C+\lfloor \frac{C}{4}\rfloor & \text{: Gregorian }(1582⪅y)\\ -C+5& \text{: Julian }(4⪅y⪅1582)\end{array}$

ただし

$C=\lfloor \frac{y}{100}\rfloor$

$Y=ymod100$

とする。

$\lfloor x\rfloor$は テンプレート:Mvar を超えない（テンプレート:Mvar以下）の最大の整数（床関数）であり、テンプレート:Math は テンプレート:Mvar を テンプレート:Mvar で割った剰余である。テンプレート:Mvar と テンプレート:Mvar は、（西暦が4桁の場合）西暦の上2桁と下2桁を表す中間変数で、たとえば2025年なら、それぞれ20と25になる。

適用範囲

以上の2つの計算式は、それぞれの暦法を過去または未来に単純に延長して使うのであれば、年数の有効範囲の限界はない（先発グレゴリオ暦、先発ユリウス暦を参照）。

しかし歴史上の実際の年月日についての曜日を知るには注意が必要である。

グレゴリウス13世は、ユリウス暦の1582年10月4日木曜日の翌日を、グレゴリオ暦の10月15日金曜日として実施した。

しかし、各国・各地方が実際にグレゴリオ暦を導入した時期は異なる場合がある。例えばイギリス帝国ではユリウス暦の1752年9月2日水曜日の翌日をグレゴリオ暦の9月14日木曜日として導入している。 テンプレート:1752年9月のカレンダー (英国) テンプレート:Clear

そのため、知りたい年月日がユリウス暦なのかグレゴリオ暦なのかによって使用すべき式を選ぶ必要がある。 テンプレート:Main

また、ガイウス・ユリウス・カエサルがユリウス暦を制定・実施したのは紀元前45年1年1月であるが、その直後の紀元前44年にカエサルが暗殺されたため、その後の閏年の設定が誤っていた時期と、それを補正した期間があった。これらの時期の閏年が何年に実施されたかについては、さまざまな主張があり議論になっている。少なくとも西暦8年以降は本来のユリウス暦のルールに則って閏年が実施されているとされている。 テンプレート:Main

ゆえに西暦8年以前の日付の実際の曜日を求めるのは困難である。

ISO 8601 との関係

曜日	日	月	火	水	木	金	土
h	1	2	3	4	5	6	0
D	7	1	2	3	4	5	6

日付の表記の国際規格である ISO 8601が定める テンプレート:Mvar（1 から 7 でそれぞれ月曜日から日曜日）を求めるには、次の式

$D=\left((d+\lfloor \frac{26(m+1)}{10}\rfloor +Y+\lfloor \frac{Y}{4}\rfloor +\mathrm{\Gamma }+5)mod7\right)+1$

を使う。

コンピュータでの計算

ツェラーの公式では年を2つの中間変数を使って計算しているがこれは筆算での計算の容易化のためで、コンピュータで計算する場合には中間変数を用いない式に変形して使用することがある。

変形したグレゴリオ暦の式の例は以下の通り（1月、2月を前年の13月、14月として扱うのは同様）。

$h^{\prime }=(y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{13m+8}{5}\rfloor +d)mod7$

同様に変形したユリウス暦の式の例は以下の通り。

$h^{\prime }=(y+\lfloor \frac{y}{4}\rfloor +\lfloor \frac{13m+8}{5}\rfloor +d+5)mod7$

曜日	日	月	火	水	木	金	土
h'	0	1	2	3	4	5	6

この場合の曜日 テンプレート:Mvar は、0 から 6 でそれぞれ日曜日から土曜日を表す。

実装するうえでは、いくつか注意すべき点があるテンプレート:Efn。

• 負数の剰余は、開発言語などによって異なる動作をすることがある（剰余演算を参照）。上記の式は数学的定義の通り、例えば テンプレート:Math であることが前提である。
• 負数の整数除算の丸めは、開発言語などによって異なる動作をすることがある（端数処理を参照）。整数除算を使って床関数と同じ結果を得るには、負の無限大への丸めである必要がある。上記の式を整数除算でそのまま記述するには、例えば テンプレート:Math が テンプレート:Math である必要がある。
• 開発言語やライブラリなどによっては、例えば1月を 0 とあらわすことがあるため、必要であれば変換を行う。

以下に Python3 でのグレゴリオ暦版の実装例を示す（以下の関数では月を 1 - 12 の値で受け取り、関数内に1月、2月を前年の13月、14月として扱うための処理を含む）。

def day_of_week(y, m, d):↲
    y, m = (y-1, m+12) if (m < 3) else (y, m)↲
    return (y + y//4 - y//100 + y//400 + (13*m + 8)//5 + d) % 7

同様の条件で、C言語での実装例を示すテンプレート:Efn。

int bmod(int n, int d)
{
    int r = n % d;
    return (r < 0) ? (r+d) : r;
}

int idiv(int n, int d)
{
    int q = n / d;
    return (n%d < 0) ? q-1 : q;
}

int day_of_week(int y, int m, int d)
{
    if (m < 3) {
        y--;
        m += 12;
    }
    return bmod(y + idiv(y, 4) - idiv(y, 100) + idiv(y, 400) + (13*m + 8)/5 + d, 7);
}

ただし、負数の年を考慮しなくてよい場合、つまり西暦0年（紀元前1年）3月1日以降のみを対象とすればよい場合には、C言語でもPython3と同様のシンプルな実装で問題ない。

int day_of_week(int y, int m, int d)
{
    if (m < 3) {
        y--;
        m += 12;
    }
    return (y + y/4 - y/100 + y/400 + (13*m + 8)/5 + d) % 7;
}

計算例

グレゴリオ暦テンプレート:Efnでの曜日の計算例を示す。

グレゴリオ暦の計算例

年月日	曜日	テンプレート:Mvar	テンプレート:Mvar
（紀元前44年3月15日） -43年3月15日	金曜日	6	5
（紀元前2年1月11日） -1年1月11日 （-2年13月11日として扱う）	月曜日	2	1
1年1月1日 （0年13月1日として扱う）	月曜日	2	1
1582年10月14日 （ユリウス暦1582年10月4日にあたる）	木曜日	5	4
1582年10月15日	金曜日	6	5
2000年2月29日 （1999年14月29日として扱う）	火曜日	3	2
2023年12月31日	日曜日	1	0

ユリウス暦テンプレート:Efnでの曜日の計算例を示す。

ユリウス暦の計算例

年月日	曜日	テンプレート:Mvar	テンプレート:Mvar
（紀元前44年3月15日） -43年3月15日	水曜日	4	3
（紀元前2年1月11日） -1年1月11日 （-2年13月11日として扱う）	土曜日	0	6
1年1月1日 （0年13月1日として扱う）	土曜日	0	6
1582年10月4日	木曜日	5	4
1582年10月5日 （グレゴリオ暦1582年10月15日にあたる）	金曜日	6	5
2000年2月29日 （1999年14月29日として扱う）	月曜日	2	1
2023年12月31日	土曜日	0	6

ツェラーの公式の導出

曜日は７日間で循環しているので、上記テンプレート:EquationNote式の 7 の剰余を求めることで曜日が判明する。即ち、

である。

曜日	日	月	火	水	木	金	土
h'	0	1	2	3	4	5	6

このとき テンプレート:Mvar のとり得る値は 0 から 6 で、それぞれ日曜日から土曜日を表す。

（現行のグレゴリオ暦は、1582年10月15日にこの日を金曜日であるとして施行されたが、この日を起点に遡ってグレゴリオ暦を適用すると1年1月1日は月曜日となるため テンプレート:Mvar の値が表す曜日がこのような並びになる）

テンプレート:EquationNote式が 7 の剰余である事を利用すると、以下の通り変形できる。

$\begin{array}{ccc}{h}^{\prime }& =& (7(52y-62)+y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{153(m+1)}{5}\rfloor +6+d)mod7\\ & =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{153(m+1)}{5}\rfloor +6+d)mod7\end{array}$

ここで、床関数の性質（ $\lfloor a\rfloor +b=\lfloor a+b\rfloor$ , ただし テンプレート:Mvar は整数）を利用すると、

$\begin{array}{ccc}{h}^{\prime }& =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{153(m+1)}{5}+6\rfloor +d)mod7\\ & =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{153m+153+30}{5}\rfloor +d)mod7\\ & =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{153m+183}{5}\rfloor +d)mod7\\ & =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{35(4m+5)+13m+8}{5}\rfloor +d)mod7\\ & =& (y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor 7(4m+5)+\frac{13m+8}{5}\rfloor +d)mod7\end{array}$

さらに、テンプレート:Mvar が 7 の剰余であることを利用して、

$h^{\prime }=(y+\lfloor \frac{y}{4}\rfloor -\lfloor \frac{y}{100}\rfloor +\lfloor \frac{y}{400}\rfloor +\lfloor \frac{13m+8}{5}\rfloor +d)mod7$

と変形できる。

テンプレート:Math ( テンプレート:Math, テンプレート:Mvar, テンプレート:Mvar は整数)と置くと、以下の通り変形される。

$\begin{array}{ccc}{h}^{\prime }& =& (100C+Y+\lfloor \frac{100C+Y}{4}\rfloor -\lfloor \frac{100C+Y}{100}\rfloor +\lfloor \frac{100C+Y}{400}\rfloor +\lfloor \frac{13m+8}{5}\rfloor +d)mod7\\ & =& (100C+Y+\lfloor 25C+\frac{Y}{4}\rfloor -\lfloor C+\frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}-1\rfloor +d)mod7\\ & =& (100C+Y+\lfloor 25C+\frac{Y}{4}\rfloor -\lfloor C+\frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}\rfloor -1+d)mod7\\ & =& (100C+Y+25C+\lfloor \frac{Y}{4}\rfloor -C-\lfloor \frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}\rfloor +d-1)mod7\\ & =& (124C+Y+\lfloor \frac{Y}{4}\rfloor -\lfloor \frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}\rfloor +d-1)mod7\\ & =& ((7\times 18)C-2C+Y+\lfloor \frac{Y}{4}\rfloor -\lfloor \frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}\rfloor +d-1)mod7\end{array}$

テンプレート:Mvar が 7 の剰余であることを利用して、

$h^{\prime }=(-2C+Y+\lfloor \frac{Y}{4}\rfloor -\lfloor \frac{Y}{100}\rfloor +\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor +\lfloor \frac{26m+26}{10}\rfloor +d-1)mod7$

ここで テンプレート:Mathより、テンプレート:Math であり、

$\lfloor \frac{Y}{100}\rfloor =0$

同様に テンプレート:Math であり、また テンプレート:Math の小数部分は、0, 0.25, 0.5, 0.75 の何れかの値を取る為、テンプレート:Math の小数部分は高々 テンプレート:Math であり、

$\lfloor \frac{C}{4}+\frac{Y}{400}\rfloor =\lfloor \frac{C}{4}\rfloor$

としてよい。よって、

$h^{\prime }=(-2C+Y+\lfloor \frac{Y}{4}\rfloor +\lfloor \frac{C}{4}\rfloor +\lfloor \frac{26m+26}{10}\rfloor +d-1)mod7$

曜日	日	月	火	水	木	金	土
h'	0	1	2	3	4	5	6
h	1	2	3	4	5	6	0

このとき テンプレート:Mvar のとり得る値は 0 から 6 で、それぞれ日曜日から土曜日を表す（テンプレート:EquationNote式の場合と同様）。

演算（項）数を減らすためこの式の被除数に 1 加算すると

$h=(-2C+Y+\lfloor \frac{Y}{4}\rfloor +\lfloor \frac{C}{4}\rfloor +\lfloor \frac{26(m+1)}{10}\rfloor +d)mod7$

となり、ツェラーの公式が導出できる。

このとき テンプレート:Mvar のとり得る値は 0 から 6 で、それぞれ土曜日から金曜日を表す（被除数に 1 加算しているためテンプレート:EquationNote式の場合と比べ、１日分のずれが生じた）。

日本での紹介

ツェラーの公式が日本で最初に紹介された時期は不明だが、以下の事例がある。

1921年（大正10年）の『日本中等教育数学会雑誌』にツェラーの公式が紹介されている。

ただし「ツェラーの公式を紹介している記事」を紹介しており、Zellerの名前は記載されていない。そのため、ツェラーの公式を Philip Franklin が1921年に考案したと間違えられることもある。

1938年（昭和13年）に『月刊数学』でZellerの名前で紹介されている。

どちらの事例でもグレゴリオ暦版の公式だけが紹介されている。

外部リンク

• テンプレート:Wikisourcelang-inline
• テンプレート:Webarchive
• テンプレート:IETF RFC - 付録Bに、ツェラーの公式の変形（中間変数は使うが、日曜日は0になるように）したものが記載されている。