調和数 (発散列)

n = ⌊x⌋ に対する調和数 H_n,1 のグラフ（赤）。これは γ + ln(x)（青）に漸近収斂する。

数学において、n-番目の調和数（ちょうわすう、テンプレート:Lang-en-short）は 1 から n までの自然数の逆数和

$H_n=1+\frac{1}{2}+\frac{1}{3}+\cdots +\frac{1}{n}=\sum _{k=1}^n\frac{1}{k}$

である。これは、1 から n までの自然数の調和平均の逆数の n-倍に等しい。

調和数は遥か昔から研究され、数論の各分野において重要である。調和数の極限は、調和級数と呼ばれ（しばしば調和数も含めて一口に調和級数と呼ぶこともある）、リーマンゼータ函数と近しい関係にあり、また種々の特殊函数のさまざまな表示に現れる。

十分大きな数の標本について、その出現頻度がジップの法則に従って分布するとき、全体の中で n-番目の頻度で現れる標本の総頻度は n-番目の調和数である。このことはロングテールおよびテンプレート:仮リンクの帰結の一種を導く。

調和数の計算法

調和数の積分表示

$H_n=\underset{0}{\overset{1}{\int }}\frac{1-x^n}{1-x}dx$

はオイラーによる。この等式は簡単な代数的等式

$\frac{1-x^n}{1-x}=1+x+\dots +x^{n-1}$

を使えば明らかである。また、積分の変数を単純に x = 1 − u と変換すれば、H_n のきれいな組合せ論的展開

$\begin{array}{cc}{H}_n& =\underset{0}{\overset{1}{\int }}\frac{1-(1-u{)}^n}{u}du=\underset{0}{\overset{1}{\int }}[\sum _{k=1}^n(-1{)}^{k-1}(\genfrac{}{}{0ex}{}{n}{k})u^{k-1}]du\\ & =\sum _{k=1}^n(-1{)}^{k-1}(\genfrac{}{}{0ex}{}{n}{k})\underset{0}{\overset{1}{\int }}{u}^{k-1}du=\sum _{k=1}^n(-1{)}^{k-1}\frac{1}{k}(\genfrac{}{}{0ex}{}{n}{k})\end{array}$

が得られる。同じ表現は、第三テンプレート:仮リンクで x₁ = 1, ..., x_n = 1 とおき、

${\mathrm{\Pi }}_k(1,\dots ,n)=(-1{)}^{n-k}(k-1)!(n-k)!$

なる事実を用いることでも得られる。すなわち

$H_n=H_{n,1}=\sum _{k=1}^n\frac{1}{k}=(-1{)}^{n-1}n!\sum _{k=1}^n\frac{1}{k^2{\mathrm{\Pi }}_k(1,\dots ,n)}=\sum _{k=1}^n(-1{)}^{k-1}\frac{1}{k}(\genfrac{}{}{0ex}{}{n}{k})$

が成り立つ。また、レトケシュ恒等式を x₁ = 1², ..., x_n = n² に対して用いれば、この場合

${\mathrm{\Pi }}_k(1^2,2^2,\dots ,n^2)=(-1{)}^{n-k}\frac{(n-k)!(n+k)!}{2k^2}$

となるので、ζ(2) の第 n-部分和についての類似の公式

$H_{n,2}=\sum _{k=1}^n\frac{1}{k^2}=2\sum _{k=1}^n(-1{)}^{k-1}\frac{1}{k^2}\frac{(\genfrac{}{}{0ex}{}{n}{k})}{(\genfrac{}{}{0ex}{}{n+k}{k})}$

を得る。H_n の増大度は n の自然対数 ln(n) と同程度の速さである。このことは、H_n を積分

$\underset{1}{\overset{n}{\int }}\frac{dx}{x}(=\ln (n))$

で近似することによって確認できる。数列 (H_n - ln(n)) は単調に減少して、

$\underset{n\to \mathrm{\infty }}{\lim }{H}_n-\ln (n)=\gamma$

なる定数（この定数 γ はオイラー・マスケローニ定数と呼ばれ、その値は 0.5772156649... である）を極限にもち、これに対応する漸近展開は

$H_n=\ln n+\gamma +\frac{1}{2}{n}^{-1}-\frac{1}{12}{n}^{-2}+\frac{1}{120}{n}^{-4}+𝒪(n^{-6})$

で与えられる。

分数パラメータに対する特殊値

調和数 H_n のパラメータ n を積分

$H_{\alpha }=\underset{0}{\overset{1}{\int }}\frac{1-x^{\alpha }}{1-x}dx$

によって拡張すれば、0 と 1 の間の分数値をもつパラメータ α に対する解析的な特殊値を定めることができる。あるいはさらに漸化式

$H_{\alpha }=H_{\alpha -1}+\frac{1}{\alpha }$

によって拡張することもでき、結局は任意の x > 0 に対して（x が整数のときもそうでないときも）

$H_x=x\sum _{k=1}^{\mathrm{\infty }}\frac{1}{k(x+k)}$

が成立する。いくつかの特殊値について計算すれば、以下のようになる。

$\begin{array}{cc}{H}_{3/4}& =\frac{4}{3}-3\ln 2+\frac{\pi }{2},\\ H_{2/3}& =\frac{3}{2}(1-\ln 3)+\frac{\pi }{6}\sqrt{3},\\ H_{1/2}& =2-2\ln 2,\\ H_{1/3}& =3-\frac{\pi }{2\sqrt{3}}-\frac{3}{2}\ln 3,\\ H_{1/4}& =4-\frac{\pi }{2}-3\ln 2,\\ H_{1/6}& =6-\frac{\pi }{2}\sqrt{3}-2\ln 2-\frac{3}{2}\ln (3),\\ H_{1/8}& =8-\frac{\pi }{2}-4\ln 2-\frac{1}{\sqrt{2}}\{\pi +\ln (2+\sqrt{2})-\ln (2-\sqrt{2})\},\\ H_{1/12}& =12-3(\ln 2+\frac{\ln 3}{2})-\pi (1+\frac{\sqrt{3}}{2})+2\sqrt{3}\ln \sqrt{2-\sqrt{3}}\end{array}$

調和数の母函数

調和数の列の母函数は

$\sum _{n=1}^{\mathrm{\infty }}{z}^n{H}_n=\frac{-\ln (1-z)}{1-z}$

で与えられる（ln(z) は自然対数）。また、冪指数型母函数は

$\sum _{n=1}^{\mathrm{\infty }}\frac{z^n}{n!}{H}_n=-e^z\sum _{k=1}^{\mathrm{\infty }}\frac{1}{k}\frac{(-z{)}^k}{k!}=e^z\text{Ein}(z)$

となる。ここで Ein(z) は整指数積分で、

$\text{Ein}(z)={\text{E}}_1(z)+\gamma +\ln z=\mathrm{\Gamma }(0,z)+\gamma +\ln z$

が成り立つものである（ただし、Γ(0, z) は不完全ガンマ函数）。

応用

調和数は、ディガンマ関数に対する

$\psi (n)=H_{n-1}-\gamma$

のような、いくつかの特殊函数に関する計算公式に現れる。このような関係式は、しばしば調和数のパラメータ n を整数以外に拡張するための定義式としても利用される。先の節で述べたような極限によって、調和数から定数 γ を定義することがよく行われるが、

$\gamma =\underset{n\to \mathrm{\infty }}{\lim }(H_n-\ln (n+\frac{1}{2}))$

としたほうが収斂が早い。

2002年にジェフリー・ラガリアスは、リーマン予想が「不等式

$\sigma (n)\le H_n+\ln (H_n)e^{H_n}$

が任意の自然数 n に対して成立し、かつ n > 1 のときは真の（等号無しの）不等式として成立する」という主張に等価であることを示した。ここで σ(n) は n の約数和である。

一般化

一般化調和数

n-番目の m-次一般化調和数 テンプレート:Lang は

$H_{n,m}=\sum _{k=1}^n\frac{1}{k^m}$

で与えられる。n を無限大に飛ばした極限が存在するのは m > 1 の時に限られることに注意。一般化調和数を表す記号としては

$H_{n,m}=H_n^{(m)}=H_m(n).$

なども使われることがある。なお、m = 1 の場合が通常の調和数であり、添字 m を落として

$H_n=\sum _{k=1}^n\frac{1}{k}$

と書く。また、n → ∞ の極限で一般化調和数はリーマンゼータ函数に収斂する。つまり

$\underset{n\to \mathrm{\infty }}{\lim }{H}_{n,m}=\zeta (m)$

が成り立つ。一般化調和数はベルヌーイ数を調べる際に現れ、またスターリング数を調べる際にも現れる。一般化調和数の母函数は

$\sum _{n=1}^{\mathrm{\infty }}{z}^n{H}_{n,m}=\frac{{\text{Li}}_m(z)}{1-z}$

である。ここで Li_m(z) は多重対数函数で |z| < 1 とする。この式で m = 1 としたものは、先に述べた調和数列の母函数に一致する。

複素平面への一般化

調和数についてのオイラーの積分公式は次の積分等式

$\underset{a}{\overset{1}{\int }}\frac{1-x^s}{1-x}dx=-\sum _{k=1}^{\mathrm{\infty }}\frac{1}{k}(\genfrac{}{}{0ex}{}{s}{k})(a-1{)}^k$

から従うが、この式は s を一般の複素数としても（二項係数を適切に拡張すれば）成り立つ。a = 0 とすれば、この公式から調和数を補間して複素平面へ拡張した函数の積分表示と級数表示が両方得られる。この積分等式自体はニュートン級数（ニュートンの一般二項定理）

$\sum _{k=0}^{\mathrm{\infty }}(\genfrac{}{}{0ex}{}{s}{k})(-x{)}^k=(1-x{)}^s$

から簡単な操作で得られる。調和数を補間する函数は、実はディガンマ関数 ψ(x) をつかって

$\psi (s+1)+\gamma =\underset{0}{\overset{1}{\int }}\frac{1-x^s}{1-x}dx$

と書ける（γ はオイラー-マスケローニ定数）。この積分の過程を繰り返せば

$H_{s,2}=-\sum _{k=1}^{\mathrm{\infty }}\frac{(-1{)}^k}{k}(\genfrac{}{}{0ex}{}{s}{k})H_k$

を得る。

関連項目

• テンプレート:仮リンク
• Tajima's D
• クーポン・コレクター問題
• テンプレート:仮リンク

参考文献

• Arthur T. Benjamin, Gregory O. Preston, Jennifer J. Quinn, A Stirling Encounter with Harmonic Numbers, (2002) Mathematics Magazine, 75 (2) pp 95-103.
• Donald Knuth. The Art of Computer Programming, Volume 1: Fundamental Algorithms, Third Edition. Addison-Wesley, 1997. ISBN 0-201-89683-4. Section 1.2.7: Harmonic Numbers, pp.75–79.
• Ed Sandifer, How Euler Did It -- Estimating the Basel problem (2003)
• テンプレート:Mathworld
• Peter Paule and Carsten Schneider, Computer Proofs of a New Family of Harmonic Number Identities, (2003) Adv. in Appl. Math. 31(2), pp. 359-378.
• Wenchang CHU, A Binomial Coefficient Identity Associated with Beukers' Conjecture on Apery Numbers, (2004) The Electronic Journal of Combinatorics, 11, #N15.
• Ayhan Dil and Istvan Mezo, A Symmetric Algorithm for Hyperharmonic and Fibonacci Numbers, (2008) Applied Mathematics and Computation 206, 942--951.
• Zoltán Retkes, "An extension of the Hermite–Hadamard Inequality", Acta Sci. Math. (Szeged), 74 (2008), pages 95–106.