解析学

【和の公式の親玉】
ベルヌーイ数をわかりやすく説明します

「わかりやすいベルヌーイ数」と書かれた黒板

ベルヌーイ数の一つの理解の仕方は、
”和の公式の親玉”です。

1から100まで足した数
1の2乗から100の2乗まで足した数

を求める公式を高校数学で習いますが、

ベルヌーイ数を使うと
任意の自然数m乗の和の公式が求まります。

この記事では簡単な地道な計算のみで
ベルヌーイ数を説明したいと思います。

[mathjax]

ベルヌーイ数とは

ベルヌーイ数とは数列{Bn}のことです。

1から始まり以降、-1/2、1/6、と無限に続いて行きます。

B0B1B2B3B4B5B6
1-1/21/60-1/3001/42

奇数番目のベルヌーイ数はB1を除きすべて0
という特徴を持ちます。

    \[B_n = 0 \quad (n \in 2\mathbb{N}+1) \hspace{20cm}\]

正確な定義

ベルヌーイ数は正確には次の式で定義されます。

    \[B_n := \sum_{k=0}^n  \sum_{r=0}^k (-1)^r {}_k \mathrm{C}_r \frac{r^n}{k+1} \hspace{20cm}\]

ベルヌーイ数を与える式は沢山あり
これは中でも一番わかりやすい物です。

二重和

それでも難しいので見方を説明します。

上式は二重和で書かれています。

二重和は具体的には

    \[<span class="ql-right-eqno"> (1) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-5a8814063c79699e6a914ef6deed2113_l3.png" height="79" width="357" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*} \sum_{k=0}^2 \sum_{r=0}^k a_{r} &=& \sum_{r=0}^0 a_{r} +\sum_{r=0}^1 a_{r} +\sum_{r=0}^2 a_{r} \\ &=& a_0 +(a_0 +a_1) +(a_0 +a_1 +a_2) \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

の様に計算される
とても多くの項の和です。

B2の値

試しに定義に従いB2を求めてみます。

    \[B_2 := \sum_{k=0}^2  \sum_{r=0}^k (-1)^r {}_k \mathrm{C}_r \frac{r^2}{k+1} \hspace{20cm}\]

    \[\quad<span class="ql-right-eqno"> (2) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-73d7339955777d1fd4f42f387b4aad80_l3.png" height="114" width="345" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*} = \sum_{r=0}^0 (-1)^r {}_0 \mathrm{C}_r \frac{r^2}{0+1} &+&\sum_{r=0}^1 (-1)^r {}_1 \mathrm{C}_r \frac{r^2}{1+1} \\ &+&\sum_{r=0}^2 (-1)^r {}_2 \mathrm{C}_r \frac{r^2}{2+1} \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/> \hspace{20cm}\]

    \[<span class="ql-right-eqno"> (3) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-d822936d69dda3dc910a612591a347b5_l3.png" height="95" width="415" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*} = \{ (-1)^0 {}_0 \mathrm{C}_0 0^2 \} +\left\{ \frac{(-1)^0 {}_1 \mathrm{C}_0 0^2}{2} \right\} +\left\{ \frac{(-1)^1 {}_1 \mathrm{C}_1 1^2}{2} \right\} \\ +\left\{ \frac{(-1)^0 {}_2 \mathrm{C}_0 0^2}{3} \right\} +\left\{ \frac{(-1)^1 {}_2 \mathrm{C}_1 1^2}{3} \right\}  +\left\{ \frac{(-1)^2 {}_2 \mathrm{C}_2 2^2}{3} \right\}    \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

    \[\quad= 0 +0 -\frac{1}{2} +0 -\frac{2}{3} +\frac{4}{3}\hspace{20cm}\]

    \[\quad =\frac{2}{3} -\frac{1}{2}\hspace{20cm}\]

    \[\quad= \frac{1}{6} \hspace{20cm}\]

0乗と0個の組み合わせ

ベルヌーイ数の計算では整数aについて

    \[a^0 := 1 \quad {}_a \mathrm{C}_0 := 1\hspace{20cm}\]

とします。

なぜ大事?

何故こんな複雑な数を考えるかと言うと

一つの理由は、
べき乗の和の公式が求まるからです。

実は高校数学で習った公式

    \[\sum_{k=1}^n k \,\, = \frac{1}{2}n(n+1) \hspace{20cm}\]

    \[\sum_{k=1}^n k^2 =\frac{1}{6}n(n+1)(2n+1) \hspace{20cm}\]

はベルヌーイ数を使っても求める事ができます。

1乗の和の公式

    \[\sum_{k=1}^n k = \frac{1}{1+1} (B_0 n^{1+1} -{}_{1+1} \mathrm{C}_1 B_1 n   ) \hspace{20cm}\]

    \[\quad \quad \,\,= \frac{1}{2}\left\{1 \cdot n^2 -2 \cdot\left(-\frac{1}{2} \right) \cdot n   \right\} \hspace{20cm}\]

    \[\quad\quad \,\,= \frac{1}{2}(n^2 +n   )\hspace{20cm}\]

    \[\quad\quad \,\,=  \frac{1}{2}n(n+1) \hspace{20cm}\]

2乗の和の公式

    \[\sum_{k=1}^n k^2 = \frac{1}{2+1} (B_0 n^{2+1} -{}_{2+1} \mathrm{C}_1 B_1 n^2 +{}_{2+1} \mathrm{C}_2 B_2 n^{2-1}   )\hspace{20cm}\]

    \[\quad\quad \quad \!= \frac{1}{3}\left\{1 \cdot n^3 -3 \cdot\left(-\frac{1}{2} \right) \cdot n^2 +3 \cdot\left(\frac{1}{6} \right) \cdot n   \right\} \hspace{20cm}\]

    \[\quad\quad \quad \!=\frac{1}{3} \left\{ n^3 +\frac{3}{2}  n^2 +\frac{1}{2} n   \right\}\hspace{20cm}\]

    \[\quad\quad \quad \!=\frac{1}{6}n(2n^2 +3n +1  ) \hspace{20cm}\]

    \[\quad\quad \quad \!=\frac{1}{6}n(n+1)(2n+1)\hspace{20cm}\]

べき乗の和の公式

一般にm乗の和の公式は、

B0~Bmのベルヌーイ数で書ける
という次の定理が成立します。

ベルヌーイの公式

    \[\sum_{k=1}^n k^m = \frac{1}{m+1} \sum_{s=0}^m (-1)^s {}_{m+1} \mathrm{C}_s B_s n^{m+1-s}\hspace{20cm}\]

なのでベルヌーイ数を先の定義式で
コンピュータに計算させておけば

任意のべき乗の和の公式が手に入ります。

3乗の和の公式

    \[<span class="ql-right-eqno"> (4) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-cf468166357afb5251b1e4e058c94e64_l3.png" height="79" width="457" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*} \sum_{k=1}^n k^3 = \frac{1}{3+1} (B_0 n^{3+1} &-&{}_{3+1} \mathrm{C}_1 B_1 n^3\\ &+&{}_{3+1} \mathrm{C}_2 B_2 n^{3-1} -{}_{3+1} \mathrm{C}_3 B_3 n^{3-2}  ) \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

    \[\quad\quad \quad \! <span class="ql-right-eqno"> (5) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-d20554aba2110cc8c4f6ec567b39c15b_l3.png" height="92" width="316" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*}=\left. \frac{1}{4}\right\{1 \cdot n^4  &-&4 \cdot\left(-\frac{1}{2} \right) \cdot n^3  \\ &+& \left.6 \cdot\left(\frac{1}{6} \right) \cdot n^2 -4 \cdot 0 \cdot n   \right\} \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{4} (n^4 +2n^3 +n^2 -0   )\hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{4}n^2 (n^2 +2n +1   ) \hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{4} n^2 (n+1)^2\hspace{20cm}\]

    \[\quad\quad \quad =\left\{ \frac{1}{2}n(n+1) \right\}^2 \hspace{20cm}\]

4乗の和の公式

最後に、知っている人の方が少ない
4乗の和の公式を導出します。

    \[<span class="ql-right-eqno"> (6) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-9142bd6ef52889d61d215a47349ca33d_l3.png" height="79" width="450" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*} \sum_{k=1}^n k^4 &=& \frac{1}{4+1} (B_0 n^{4+1} -{}_{4+1} \mathrm{C}_1 B_1 n^4\\ &+&{}_{4+1} \mathrm{C}_2 B_2 n^{4-1} -{}_{4+1} \mathrm{C}_3 B_3 n^{4-2} +{}_{4+1} \mathrm{C}_4 B_4 n^{4-3}  ) \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

    \[\quad\quad \quad \!\! <span class="ql-right-eqno"> (7) </span><span class="ql-left-eqno">   </span><img src="https://shoma-life-blog.com/wp-content/ql-cache/quicklatex.com-8e363e44e13b2abdc317d17e60c51e63_l3.png" height="92" width="364" class="ql-img-displayed-equation quicklatex-auto-format" alt="\begin{eqnarray*}&=& \frac{1}{5} \left\{1 \cdot n^5 -5 \cdot\left(-\frac{1}{2} \right) \cdot n^4 \right. \\ &+& \left.10 \cdot\left(\frac{1}{6} \right) \cdot n^3 -10 \cdot 0 \cdot n^2 +5 \cdot \left(-\frac{1}{30} \right) \cdot n  \right\} \end{eqnarray*}" title="Rendered by QuickLaTeX.com"/>\hspace{20cm}\]

    \[\quad\quad \quad  = \frac{1}{5}\left\{ n^5 +\frac{5}{2} n^4+ \frac{5}{3} n^3 -0 -\frac{1}{6} n  \right\} \hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{30} n(6n^4 +15n^3 +10n^2 -1   )\hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{30} n(n+1)(6n^3 +9n^2 +n -1   ) \hspace{20cm}\]

    \[\quad\quad \quad =\frac{1}{30} n(n+1)(2n+1)(3n^2 +3n -1   ) \hspace{20cm}\]

ゴリゴリの数学になりますが
以下のサイトが勉強になります。

ベルヌーイ数の漸化式は

    \[B_n = -\sum_{k=0}^{n-1} \frac{1}{n+1} {}_{n+1} \mathrm{C}_k B_k\]

です。
(最後に両辺にnの階乗を掛けます)

-解析学