調和数

同義語:harmonic numberharmonic numbers

概要

調和数(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}$ をもち、ランダムな置換の巡回の個数やくじ集めの回数の平均を表す。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 調和級数, 対数関数

正の整数の逆数を $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. $\displaystyle\sum_{k=1}^nH_k=(n+1)H_n-n$。
  2. $\displaystyle\sum_{k=1}^nkH_k=\frac{n(n+1)}2H_n-\frac{n(n-1)}4$。
  3. $\displaystyle H_n=\int_0^1\frac{1-x^n}{1-x}\,dx$(被積分関数は $x=1$ で値 $n$ をとる多項式 $1+x+\cdots+x^{n-1}$ である)。
  4. $\displaystyle H_n=\sum_{k=1}^n(-1)^{k-1}\binom nk\frac1k$。
和の順序の交換と積分

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)。

調和数の差と $\log2$

正の整数 $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$ は整数でない。

2 で割り切れる回数が最大の項に着目する

整数 $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}$ を考える。

$H_{p-1}$ の分子

$p$ を奇素数とする。$H_{p-1}$ を既約分数 $\frac ND$ で表すと、$N$ は $p$ で割り切れる。

$k$ と $p-k$ を組にする

$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$ 個の元の部分集合ごとに数える

$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$ で割り切れない。

関連項目

参考文献

[2]
Philippe Flajolet and Robert Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009, §II.3 (32)(p. 117、くじ集めの回数の平均)、Example III.4 (11)(p. 160、置換の巡回の個数の平均と母関数)、Note IV.12(p. 238、二項係数による交代和)、Appendix A.2 (4)(p. 724、調和数の漸近展開)。頁は著者の e-version(2009-06-26)による

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する