調和数(harmonic number)とは、正の整数 $n$ について $H_n=1+\frac12+\cdots+\frac1n$ で定まる有理数であり、調和級数の部分和である。$\log(n+1)< H_n\le1+\log n$ が成り立ち、$H_n-\log n$ は狭義単調減少して Euler の定数 $\gamma$ に収束し、その差は $\frac1n$ より小さい。連続する 2 個以上の正の整数の逆数の和は整数でないので、$n\ge2$ の $H_n$ は整数でなく、奇素数 $p$ について $H_{p-1}$ の分子は $p$ で割り切れる。和・積分表示・二項係数による表示・母関数 $\frac1{1-x}\log\frac1{1-x}$ をもち、ランダムな置換の巡回の個数やくじ集めの回数の平均を表す。
正の整数の逆数を $1$ から $\frac1n$ まで足した数 $H_n=1+\frac12+\cdots+\frac1n$ を 調和数(harmonic number)という。$H_n$ は 調和級数 $\sum\frac1k$ の部分和であり、$n$ とともに限りなく大きくなるが、その増え方は $\log n$ と同じくらい遅い。差 $H_n-\log n$ は Euler の定数 $\gamma=0.5772\ldots$ に収束する。
調和数は、和の計算・数論・確率のどこにも現れる。積分 $\int_1^n\frac{dx}x=\log n$ の和の版として和の公式に現れ、分母を調べると $n\ge2$ では整数でないことや、奇素数 $p$ について $H_{p-1}$ の分子が $p$ で割り切れることが分かる。確率では、ランダムに選んだ $n$ 文字の置換の巡回の個数の平均や、$n$ 種類のくじをすべて集めるまでの回数の平均が調和数で表される。
この記事では、調和数を定義して、基本の恒等式(和・積分表示・二項係数による表示)と母関数を証明する。次に、$\log n$ との比較 $\log(n+1)< H_n\le1+\log n$ と、差 $H_n-\log n$ の極限を示す。さらに、連続する整数の逆数の和が整数にならないことと、奇素数 $p$ について $H_{p-1}$ の分子が $p$ で割り切れることを証明し、置換の巡回の個数とくじ集めの問題への応用を述べる。
$0$ 以上の整数 $n$ について
$$
H_n:=\sum_{k=1}^{n}\frac1k=1+\frac12+\cdots+\frac1n\qquad(H_0:=0)
$$
を第 $n$ 調和数という。より一般に、実数 $r>0$ について $H_n^{(r)}:=\sum_{k=1}^nk^{-r}$ を 一般化調和数 という($H_n^{(1)}=H_n$)。
$H_n$ は有理数で、$H_{n}-H_{n-1}=\frac1n>0$ なので $(H_n)$ は狭義単調増加である。$r>1$ なら $H_n^{(r)}$ は $n\to\infty$ で収束し(調和級数 の命題「冪の逆数の級数」)、極限は Riemannゼータ関数 の値 $\zeta(r)$ である。たとえば $H_n^{(2)}\to\frac{\pi^2}6$(Basel問題)。$r=1$ では $H_n\to\infty$ である(調和級数 の定理「調和級数の発散」)。
$n=1,\dots,10$ の値は次のとおりである。
| $n$ | $H_n$ | 小数 |
|---|---|---|
| $1$ | $1$ | $1$ |
| $2$ | $\frac32$ | $1.5$ |
| $3$ | $\frac{11}6$ | $1.8333\ldots$ |
| $4$ | $\frac{25}{12}$ | $2.0833\ldots$ |
| $5$ | $\frac{137}{60}$ | $2.2833\ldots$ |
| $6$ | $\frac{49}{20}$ | $2.45$ |
| $7$ | $\frac{363}{140}$ | $2.5928\ldots$ |
| $8$ | $\frac{761}{280}$ | $2.7178\ldots$ |
| $9$ | $\frac{7129}{2520}$ | $2.8289\ldots$ |
| $10$ | $\frac{7381}{2520}$ | $2.9289\ldots$ |
既約分数の分母は $n=1$ を除いて偶数で、分子は奇数である(thm-hn-not-integer の証明がその理由を示す)。$H_4=\frac{25}{12}$、$H_6=\frac{49}{20}$、$H_{10}=\frac{7381}{2520}=\frac{11^2\cdot61}{2520}$ の分子は、それぞれ $5^2$、$7^2$、$11^2$ で割り切れる(prop-hn-prime-numerator)。
正の整数 $n$ について次が成り立つ。
1・2:$H_k=\sum_{j\le k}\frac1j$ を代入して $j$ と $k$ の和の順序を入れかえる。$j$ を固定すると $k$ は $j$ から $n$ まで動くので
$$
\sum_{k=1}^nH_k=\sum_{j=1}^n\frac{n-j+1}j=(n+1)H_n-n,
$$
$$
\sum_{k=1}^nkH_k=\sum_{j=1}^n\frac1j\sum_{k=j}^nk=\sum_{j=1}^n\frac1j\Bigl(\frac{n(n+1)}2-\frac{j(j-1)}2\Bigr)=\frac{n(n+1)}2H_n-\frac12\sum_{j=1}^n(j-1)
$$
であり、$\sum_{j=1}^n(j-1)=\frac{n(n-1)}2$ である。
3:$x\ne1$ で $\frac{1-x^n}{1-x}=\sum_{k=0}^{n-1}x^k$ であり、$\int_0^1x^k\,dx=\frac1{k+1}$ なので、積分は $\sum_{k=0}^{n-1}\frac1{k+1}=H_n$ である。
4:3 で $x=1-y$ と置き換えると $H_n=\int_0^1\frac{1-(1-y)^n}y\,dy$ である。二項定理から $1-(1-y)^n=\sum_{k=1}^n(-1)^{k-1}\binom nky^k$ なので、被積分関数は $\sum_{k=1}^n(-1)^{k-1}\binom nky^{k-1}$ であり、項ごとに積分して 4 を得る。$\square$
1 は、積分の公式 $\int_1^n\log x\,dx=n\log n-n+1$ に当たる和の公式である(和と積分の対応は 和と積分の対応表 で扱う)。4 は、たとえば $n=3$ で $3-\frac32+\frac13=\frac{11}6$ であり、$\sum_{k=1}^n\binom nk\frac{(-1)^k}k=-H_n$ の形で FS09 p. 238(Note IV.12)にある。
$|x|< 1$ について
$$
\sum_{n=1}^\infty H_nx^n=\frac1{1-x}\log\frac1{1-x}
$$
が成り立つ。
$1\le H_n\le n$ なので、左辺の冪級数は $|x|< 1$ で絶対収束する。その和を $G(x)$ とおくと、$H_0=0$ と $H_n-H_{n-1}=\frac1n$ から
$$
(1-x)G(x)=\sum_{n\ge1}H_nx^n-\sum_{n\ge1}H_{n-1}x^n=\sum_{n\ge1}(H_n-H_{n-1})x^n=\sum_{n\ge1}\frac{x^n}n
$$
である(2 つの級数はどちらも絶対収束するので、項ごとに引いてよい)。関数の近似とテイラー展開 の定理「対数関数の展開」を $-x$ に使うと、$|x|< 1$ で $\sum_{n\ge1}\frac{x^n}n=-\log(1-x)=\log\frac1{1-x}$ である。$1-x$ で割って主張を得る。$\square$
母関数は、$\frac1{1-x}=\sum x^m$ と $\log\frac1{1-x}=\sum\frac{x^k}k$ の積であり、その係数が $\sum_{k\le n}\frac1k=H_n$ になること(Cauchy積)とも読める。FS09 p. 160 は、この母関数を置換の巡回の個数の平均の計算に使っている(prop-hn-cycles)。
正の整数 $n$ について
$$
\log(n+1)< H_n\le1+\log n
$$
が成り立つ。右の不等号が等号になるのは $n=1$ のときに限る。
$k\ge1$ とする。$k< x< k+1$ で $\frac1{k+1}< \frac1x< \frac1k$ なので、連続関数を積分して
$$
\frac1{k+1}< \int_k^{k+1}\frac{dx}x=\log(k+1)-\log k< \frac1k
$$
である(対数関数 の定理「積分による表示」)。右の不等号を $k=1,\dots,n$ について加えると $\log(n+1)< H_n$ である。$n\ge2$ なら、左の不等号を $k=1,\dots,n-1$ について加えると $H_n-1< \log n$ である。$n=1$ では $H_1=1=1+\log1$ である。$\square$
$D_n:=H_n-\log n$($n\ge1$)とおく。$(D_n)$ は正の数からなる狭義単調減少列であり、極限 $\gamma:=\lim_{n\to\infty}D_n$ が存在して、すべての $n\ge1$ で
$$
0< D_n-\gamma< \frac1n
$$
が成り立つ。すなわち $H_n=\log n+\gamma+\frac{\theta_n}n$($0< \theta_n< 1$)と書ける。
thm-hn-log から $D_n>\log(n+1)-\log n>0$ である。prf-hn-log の評価から
$$
D_n-D_{n+1}=\bigl(\log(n+1)-\log n\bigr)-\frac1{n+1}
$$
は $0$ より大きく、$\frac1n-\frac1{n+1}$ より小さい。よって $(D_n)$ は下に有界な狭義単調減少列で、極限 $\gamma\ge0$ をもつ(単調収束定理)。$e_k:=D_k-D_{k+1}$、$d_k:=\frac1k-\frac1{k+1}$ とおくと $0< e_k< d_k$ であり、$N\to\infty$ とした望遠鏡和から
$$
D_n-\gamma=\sum_{k=n}^\infty e_k,\qquad\frac1n=\sum_{k=n}^\infty d_k
$$
である。各 $e_k$ は正なので $D_n-\gamma>0$ であり、$\frac1n-(D_n-\gamma)=\sum_{k\ge n}(d_k-e_k)\ge d_n-e_n>0$ なので $D_n-\gamma< \frac1n$ である。$\square$
$\gamma=0.5772156649\ldots$ は Euler の定数である。$n=1$ では $D_1=1$ なので、上の評価から $0\le\gamma< 1$ が分かる。より精密な評価 $\frac1{2(n+1)}< D_n-\gamma< \frac1{2n}$ と $\gamma$ のほかの表し方は Eulerの定数 で扱い、$H_n$ の漸近展開 $H_n=\log n+\gamma+\frac1{2n}-\frac1{12n^2}+O(n^{-4})$ は Euler–Maclaurinの公式 から得られる(FS09 p. 724)。
正の整数 $n$ について
$$
H_{2n}-H_n=\sum_{k=1}^{2n}\frac{(-1)^{k-1}}k
$$
であり、$\log2-\frac1n< H_{2n}-H_n< \log2$ である。とくに $H_{2n}-H_n\to\log2$ である。
右辺の和は、$H_{2n}$ から偶数番目の項を 2 回引いたものなので、$H_{2n}-2\sum_{j=1}^n\frac1{2j}=H_{2n}-H_n$ である。また $H_{2n}-H_n=\log2+D_{2n}-D_n$ であり、thm-hn-gamma から $D_{2n}< D_n$、$D_n-D_{2n}< D_n-\gamma< \frac1n$ なので、評価を得る。$\square$
$n\to\infty$ とすると、交代調和級数 $1-\frac12+\frac13-\cdots$ の偶数番目までの部分和が $\log2$ に収束することになる(調和級数 の系「交代調和級数の和」)。$n=10$ では $H_{20}-H_{10}=0.6687\ldots$、$\log2=0.6931\ldots$ である。
$H_n$ は $n\ge2$ で整数にならない。もっと一般に、連続する 2 個以上の整数の逆数の和は整数にならない。$2$ で割り切れる回数が最大の項がただ 1 つしかないことが理由である。
$1\le m< n$ を整数とすると、$\displaystyle\sum_{k=m}^n\frac1k$ は整数でない。とくに $n\ge2$ なら $H_n$ は整数でない。
整数 $k\ge1$ が $2$ で割り切れる回数を $v(k)$ と書き($k=2^{v(k)}\cdot(\text{奇数})$)、$s:=\max_{m\le k\le n}v(k)$ とおく。$m< n$ なので範囲には連続する 2 つの整数があり、その一方は偶数だから $s\ge1$ である。
段 1(最大の項はただ 1 つ)。範囲の中に $2^s$ の倍数が 2 つ、$u2^s< w2^s$ があるとすると、$u< u+1\le w$ より $(u+1)2^s$ も範囲にある。$u$ と $u+1$ の一方は偶数なので、$u2^s$ と $(u+1)2^s$ の一方は $2^{s+1}$ で割り切れ、$s$ の最大性に反する。よって $v(k)=s$ となる $k$ はただ 1 つで、それを $k_0$ とする。
段 2(偶奇を比べる)。$L$ を $m,m+1,\dots,n$ の最小公倍数とすると、$L$ が $2$ で割り切れる回数は $s$ である。$k\ne k_0$ なら $v(k)< s$ なので $L/k$ は偶数であり、$L/k_0$ は奇数である。したがって
$$
L\sum_{k=m}^n\frac1k=\sum_{k=m}^n\frac Lk
$$
は奇数である。和 $\sum_{k=m}^n\frac1k$ が整数なら、左辺は偶数 $L$($s\ge1$)の倍数なので偶数になり、矛盾する。$\square$
$m=1$ の場合が $H_n$ である。段 2 から、$H_n=\frac{(\text{奇数})}{L}$ と書けるので、既約分数にしても分子は奇数、分母は $2^s$ で割り切れる。たとえば $H_{10}=\frac{7381}{2520}$ では $s=3$($k_0=8$)で、$2520=2^3\cdot315$ である。2 進付値の言葉では $v_2(H_n)=-s$ であり、p進付値 の命題「調和数の 2 進付値」が同じことを述べている。
次に、分母に素数 $p$ が現れない調和数 $H_{p-1}$ を考える。
$p$ を奇素数とする。$H_{p-1}$ を既約分数 $\frac ND$ で表すと、$N$ は $p$ で割り切れる。
$k$ と $p-k$ を組にすると
$$
2H_{p-1}=\sum_{k=1}^{p-1}\Bigl(\frac1k+\frac1{p-k}\Bigr)=p\sum_{k=1}^{p-1}\frac1{k(p-k)}
$$
である。右の和を通分すると分母を $B:=((p-1)!)^2$ にとれるので、整数 $A$ を使って $H_{p-1}=\frac{pA}{2B}$ と書ける。$p$ は奇数で $(p-1)!$ の因数はどれも $p$ より小さいので、$2B$ は $p$ で割り切れない。$\frac ND=\frac{pA}{2B}$ から $2BN=pAD$ であり、$p$ は素数で $2B$ を割り切らないので $p\mid N$ である。$\square$
$p\ge5$ なら $N$ は $p^2$ で割り切れる。これが Wolstenholmeの定理 であり(HW60 Theorem 115、§7.8、p. 88)、ex-hn-values の $H_4$、$H_6$、$H_{10}$ がその例である。$p=3$ では $H_2=\frac32$ の分子 $3$ は $9$ で割り切れない。
調和数は、確率で「$k$ 番目の寄与が $\frac1k$」となる量の平均として現れる。
$\{1,\dots,n\}$ の置換 $\sigma$ を $n!$ 通りから等しい確率で選ぶ。$1\le k\le n$ について、$\sigma$ の長さ $k$ の巡回(不動点は長さ $1$ の巡回と数える)の個数の期待値は $\frac1k$ であり、巡回の個数の期待値は $H_n$ である。
$k$ 個の元からなる部分集合 $A$ を固定する。$A$ がちょうど $\sigma$ の 1 つの巡回をなす置換 $\sigma$ は、$A$ の元を 1 つの巡回に並べる方法 $(k-1)!$ 通りと、補集合の置換 $(n-k)!$ 通りを組み合わせて得られる。よって、$\sigma$ と、$\sigma$ の長さ $k$ の巡回の台 $A$ の組の個数は $\binom nk(k-1)!(n-k)!=\frac{n!}k$ である。$n!$ で割ると、長さ $k$ の巡回の個数の期待値は $\frac1k$ である。巡回の個数は長さごとの個数の和なので、期待値 の定理「期待値の線形性」から、その期待値は $\sum_{k=1}^n\frac1k=H_n$ である。$\square$
thm-hn-gamma と合わせると、ランダムな置換の巡回の個数は平均して $\log n+\gamma$ 程度である(FS09 p. 160)。相異なる $n$ 個の数をランダムに並べて左から見ていくとき、それまでの最大値を更新する回数の平均も $H_n$ である($k$ 番目の数がそれまでの $k$ 個の中で最大になる確率が $\frac1k$ であることによる。記録更新の回数)。
$n$ 種類の景品が等しい確率 $\frac1n$ で独立に出るくじを、全種類がそろうまで引くとき、引く回数の期待値は $nH_n$ である。$i-1$ 種類($1\le i\le n$)がそろった後、新しい種類が出る確率は $\frac{n-i+1}n$ なので、新しい種類が出るまでの回数の期待値は $\frac n{n-i+1}$ である(成功の確率が $q>0$ の試行を成功するまで繰り返すと、回数の期待値は $\sum_{j\ge1}j(1-q)^{j-1}q=\frac1q$)。これを $i=1,\dots,n$ について加えると $n\bigl(\frac1n+\frac1{n-1}+\cdots+1\bigr)=nH_n$ である。さいころの 6 つの目をすべて出すには、平均 $6H_6=6\cdot\frac{49}{20}=14.7$ 回振る。同じ値は FS09 p. 117 の (32) にあり、そこでは交代和 $n\sum_j\binom nj\frac{(-1)^{j-1}}j$ とも表されて、prop-hn-identities の 4 と結びつく。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| 逆数をとる整数が連続していること | $\frac12+\frac13+\frac16$ | 和が整数でない |
| 項が 2 個以上($m< n$) | $m=n=1$ | 和が整数でない |
| $p$ が素数 | $p=9$ の $H_8=\frac{761}{280}$ | $H_{p-1}$ の分子が $p$ で割り切れる |
| $p\ge5$ | $p=3$ | $H_{p-1}$ の分子が $p^2$ で割り切れる |
1 行目:$\frac12+\frac13+\frac16=\frac{3+2+1}6=1$ である。分母のうち 2 で割り切れる回数が最大(1 回)のものが $2$ と $6$ の 2 つあるので、prf-hn-not-integer の段 1 が成り立たない。
2 行目:$\frac11=1$ は整数である。
3 行目:$761=3\cdot253+2$ は $9$ で($3$ でも)割り切れない。合成数 $n$ では $k$ と $n-k$ の組の分母 $k(n-k)$ が $n$ の約数を含みうるので、prf-hn-prime-numerator の議論が使えない。
4 行目:$H_2=\frac32$ の分子 $3$ は $9$ で割り切れない。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する