Sstrhsrtj

ベルヌーイ数 (ベルヌーイすう、英: Bernoulli number) は数論における基本的な係数を与える数列であり、もともと、連続する整数のべき乗和を定式化する際の展開係数として1713年にヤコブ・ベルヌーイが著書 Ars Conjectandi (推測術) にて導入した^[1]ことからこの名称がついた。ベルヌーイ数は、べき乗和の展開係数にとどまらず、級数展開の係数や剰余項、リーマンゼータ関数においても登場する。また、ベルヌーイ数はすべてが有理数である。

定義

ベルヌーイ数 B_n は下に示すマクローリン展開 (テイラー展開) の展開係数として定義される。

f(x)=xex−1=∑n=0∞Bnn!xn.{displaystyle f(x)={frac {x}{e^{x}-1}}=sum _{n=0}^{infty }{frac {B_{n}}{n!}}x^{n}.}

この定義から、関数 x/(e^x − 1) を繰り返し微分していけばベルヌーイ数を得ることができるが、そのような手段でベルヌーイ数を得るのは容易ではない。 ベルヌーイ数を計算するには、マクローリン展開ではなく、次の漸化式を用いる。この漸化式から、ベルヌーイ数がすべて有理数であることがわかる。

B0=1,Bn=−1n+1∑k=0n−1(n+1k)Bk.{displaystyle B_{0}=1,quad {}B_{n}=-{1 over n+1}sum _{k=0}^{n-1}{n+1 choose k}B_{k}.}

ここで、(n+1k){displaystyle left({textstyle {n+1 atop k}}right)} は二項係数である。

上の漸化式を用いて、ベルヌーイ数を第 29 項までを算出すると下表のようになる。この表には、ベルヌーイ数が有理数であるので、分子と分母を記載している。

ベルヌーイ数の漸化式は、上記の関数 f(x) = x/(e^x − 1) の逆数をテイラー展開し、その 2 つの積が 1 になることから導出できる。その漸化式は厳密な計算には有用であるが、n が大きくなると途中の式の値が非常に大きくなるため、浮動小数点数を使って計算する場合、精度が著しく悪くなる計算として知られている。

奇数番目のベルヌーイ数は B₁ 以外はすべて 0 であり、偶数番目は B₀ を除いて正の数と負の数が交互に並ぶ。 ベルヌーイ数の第 3 項以降の奇数項が 0 となることは、正接関数 (tangent) のマクローリン展開から証明^[2]できる。

ベルヌーイ数の一般項

第 2 種スターリング数との関係から、次のようなベルヌーイ数の一般項を算出する公式が存在する。

Bn=∑j=0n(−1)jjn∑m=jn1m+1(mj).{displaystyle B_{n}=sum _{j=0}^{n}(-1)^{j},j^{n}sum _{m=j}^{n}{frac {1}{m+1}}{m choose j}.}

この公式は、総和記号が二重になっているため、上に示した漸化式ほど手軽にベルヌーイ数を計算する公式ではない。

漸近的性質

ベルヌーイ数とリーマンゼータ関数の関係から、

B2n=(−1)n+12(2n)!(2π)2n(1+122n+132n+142n+⋯){displaystyle B_{2n}=(-1)^{n+1}{frac {2(2n)!}{(2pi )^{2n}}}left(1+{frac {1}{2^{2n}}}+{frac {1}{3^{2n}}}+{frac {1}{4^{2n}}}+dotsb right)}

が成り立つ。従ってスターリングの公式から、n → ∞ のとき、

|B2n|∼4πn(nπe)2n{displaystyle |B_{2n}|sim 4{sqrt {pi n}}left({frac {n}{pi e}}right)^{2n}}

が成り立つ。

ベルヌーイ数を用いた級数展開

ベルヌーイ数は、いくつかの双曲線関数と三角関数の級数展開における展開係数となる。 ベルヌーイ数を展開係数とする関数とそのローラン級数による表現を挙げる。 まず、余接関数 (cotangent) のローラン級数展開は次のようになる。

coth⁡z=1z+∑k=1∞22kB2k(2k)!z2k−1,cot⁡z=1z+∑k=1∞(−1)k22kB2k(2k)!z2k−1.{displaystyle {begin{aligned}coth z&={frac {1}{z}}+sum _{k=1}^{infty }{frac {2^{2k}B_{2k}}{(2k)!}}z^{2k-1},\cot z&={frac {1}{z}}+sum _{k=1}^{infty }(-1)^{k}{frac {2^{2k}B_{2k}}{(2k)!}}z^{2k-1}.end{aligned}}}

第 1 の関係式は、ベルヌーイ数が x/(e^x − 1) の展開係数であることを利用して数式変形すれば得られる。 第 2 の関係式は cot z = -i coth (-iz) であることを利用すれば、第 1 の関係式から導き出される。これらの級数の収束半径は |z| < π である。 次に正接関数 (tangent) のローラン級数展開は次のようになる。

tan⁡z=∑k=1∞(−1)k(22k−42k)B2k(2k)!z2k−1.{displaystyle tan z=sum _{k=1}^{infty }{frac {(-1)^{k},(2^{2k}-4^{2k}),B_{2k}}{(2k)!}}z^{2k-1}.}

この関係式は、tan z = cot z − 2cot 2z を利用して余接関数のローラン級数展開を変形すれば導出できる。 なお、この級数の収束半径は |z| < π/2 である。この正接関数のローラン級数展開の展開係数による数列はタンジェント数と呼ばれる。 一方、余割関数 (cosecant) は次のようにローラン級数展開される。

csc⁡z=1sin⁡z=1z+∑k=1∞(−1)k(2−22k)B2k(2k)!z2k−1.{displaystyle csc z={frac {1}{sin z}}={frac {1}{z}}+sum _{k=1}^{infty }{frac {(-1)^{k},(2-2^{2k}),B_{2k}}{(2k)!}}z^{2k-1}.}

この関係式は、csc 2z = (tan z + cot z)/2 を利用すれば導出できる。 なお、この級数の収束半径は |z| < π である。

べき乗和による導入

ベルヌーイ数は、もともと、連続する整数のべき乗和を定式化する際に、展開係数として導入された。 現代の表記法によって書くならば、定式化するべき乗和とは、

Sk(n)≡∑j=0n−1jk=0k+1k+2k+⋯+(n−1)k{displaystyle S_{k}(n)equiv sum _{j=0}^{n-1}j^{k}=0^{k}+1^{k}+2^{k}+cdots +(n-1)^{k}}

なる総和である。この総和は、ベルヌーイ数を用いて、

Sk(n)=1k+1∑j=0k(k+1j)Bjnk−j+1{displaystyle S_{k}(n)={frac {1}{k+1}}sum _{j=0}^{k}{k+1 choose j}B_{j},n^{k-j+1}}

のように書くことができる。 ベルヌーイ数の漸化式は、べき乗和を定式化した際の考察から得られる。 さらに、ベルヌーイ数の指数型母関数が x/(e^x − 1) となることから、その母関数を現在ではベルヌーイ数の定義とする。

ヤコブ・ベルヌーイは彼の著書『推測術』でベルヌーイ数を導入した際、べき乗和を上に書いたような 0 から n − 1 にわたる和でなく、1 から n にわたる和:

S^k(n)≡∑j=1njk=1k+2k+3k+⋯+nk{displaystyle {hat {S}}_{k}(n)equiv sum _{j=1}^{n}j^{k}=1^{k}+2^{k}+3^{k}+cdots +n^{k}}

として扱っていた。 ベルヌーイは、その著書で整数のべき乗 n^c の和を計算する公式として、次の数式を記している^[3]。

∫nc=1c+1nc+1+12nc+c2Anc−1+c.c−1.c−22.3.4Bnc−3+c.c−1.c−2.c−3.c−42.3.4.5.6Cnc−5+c.c−1.c−2.c−3.c−4.c−5.c−62.3.4.5.6.7.8Dnc−7+……{displaystyle {begin{aligned}int n^{c}&={}{frac {1}{c+1}}n^{c+1}+{frac {1}{2}}n^{c}+{frac {c}{2}}An^{c-1}+{frac {c.c-1.c-2}{2.3.4}}Bn^{c-3}\&quad quad quad {}+{frac {c.c-1.c-2.c-3.c-4}{2.3.4.5.6}}Cn^{c-5}\&quad quad quad {}+{frac {c.c-1.c-2.c-3.c-4.c-5.c-6}{2.3.4.5.6.7.8}}Dn^{c-7}+ldots ldots end{aligned}}}

この数式に記載されている展開係数 A,B,C,D,…{displaystyle A,B,C,D,ldots } がベルヌーイ数 (B2{displaystyle B_{2}} 以降) である。ベルヌーイが記した数式は、

S^k(n)=1k+1∑j=0k(k+1j)B^jnk−j+1{displaystyle {hat {S}}_{k}(n)={frac {1}{k+1}}sum _{j=0}^{k}{k+1 choose j}{hat {B}}_{j},n^{k-j+1}}

に相当する。この数式に用いた展開係数 B^j{displaystyle {hat {B}}_{j}} は、

B^1=1/2(=−B1),B^j=Bj(j≠1){displaystyle {hat {B}}_{1}=1/2,(=-B_{1}),quad {hat {B}}_{j}=B_{j}quad (jneq 1)}

のように、j≠1{displaystyle jneq 1} においてベルヌーイ数と一致する。一部の文献^[4][5]では Bj{displaystyle B_{j}} の代わりに B^j{displaystyle {hat {B}}_{j}} をベルヌーイ数と呼んでいる。

一方、日本ではベルヌーイとほぼ同時期に関孝和がべき乗和を定式化し、ベルヌーイ数を発見していた^[6]。 そのため、ベルヌーイ数を関・ベルヌーイ数と書いている文献^[7]もある。

一般ベルヌーイ数

一般ベルヌーイ数は代数的数で、ベルヌーイ数がリーマンゼータ函数の特殊値に関連する方法と同じ方法で、ディリクレの L-関数の特殊値に関連して定義される。

χ を mod f のディリクレ指標とすると、一般ベルヌーイ数 B_k,χ は、

∑a=1fχ(a)teateft−1=∑k=0∞Bk,χtkk!{displaystyle sum _{a=1}^{f}chi (a){frac {te^{at}}{e^{ft}-1}}=sum _{k=0}^{infty }B_{k,chi }{frac {t^{k}}{k!}}}

により定義される。B_1,1 = 1/2 を除き、任意のディリクレ指標 χ に対し、χ(−1) ≠ (−1)^k であれば、B_k,χ = 0 である。

正でない整数におけるリーマンゼータ関数の値とベルヌーイ数の間の関係を一般化し、全ての整数 k ≥ 1 に対し、

L(1−k,χ)=−Bk,χk{displaystyle L(1-k,chi )=-{frac {B_{k,chi }}{k}}}

が成り立つ。ここに L(s, χ) は χ のディリクレの L-関数である^[8]。

脚注

1. 
   ^ Julian Havil, "オイラーの定数 ガンマ," 新妻弘 訳, 共立出版, 初版, pp. 97–99, 2009.
   


2. 
   ^ 正接関数のマクローリン展開の結果において、実数変数を仮定した場合、ベルヌーイ数の第 3 項以降の奇数項は虚数項に対応する。 実数変数における正接関数が実数関数でなければならないので、そのマクローリン展開に虚数項に対応する項が存在してはならない。 よって、ベルヌーイ数の第3項以降の奇数項はゼロでなければならない。
   


3. 
   ^ E. Hairer, G. Wanner, "解析教程 上," 蟹江幸博 訳, シュプリンガー・ジャパン, 新装版, p. 18, 2006.
   


4. 
   ^ 例えば、 荒木恒男, 伊吹山知義, 金子昌信, "ベルヌーイ数とゼータ関数," 牧野書店, 2001.
   


5. 
   ^ Wikipedia ファウルハーバーの公式 もベルヌーイの記述に基づき、第 1 項を1/2とする記述で説明している。
   


6. 
   ^ 小川束, "関孝和によるベルヌーイ数の発見," 数理解析研究所講究録, 第1583巻, 2008.
   


7. 
   ^ 例えば、 桜井進, 中村義作, "天才たちが愛した美しい数式," PHP研究所, 第1版, p.205, 2008.
   


8. 
   ^ Neukirch 1999, §VII.2
   

関連項目

• ヤコブ・ベルヌーイ


• 関孝和


• ファウルハーバーの公式


• ベルヌーイ多項式


• ゼータ関数


• 三角関数