表記法 本稿で使用されている上付き文字± は、ベルヌーイ数の符号規則を区別するために用いられています。影響を受けるのはn =1の 項のみです。
以下の式では、次の関係式を用いて、ある符号規則から別の符号規則に切り替えることができます。B n + = ( − 1 ) n B n − {\displaystyle B_{n}^{+}=(-1)^{n}B_{n}^{-}} または、整数n = 2 以上の場合は、単純に無視します。
すべての奇数n > 1に対して B n = 0 であり、多くの公式は偶数インデックスのベルヌーイ数のみを扱うため、一部の著者はB 2 n の代わりに「 B n 」と表記します。この記事ではその表記法には従いません。
歴史
初期の歴史 ベルヌーイ数は、古代から数学者たちの関心を集めてきた整数のべき乗の和の計算という、初期の歴史に根ざしている。
関隆一 の『勝陽散法』 (1712年)の一ページで、二項係数とベルヌーイ数を表にまとめたもの。最初のn個の正の整数の和、最初の n 個の正の整数の平方和、および立方和を計算する方法は知られていましたが、実際の「公式」はなく、言葉による説明のみでした。この問題を考察した古代の偉大な数学者には、ピタゴラス (紀元前 572 年頃 - 497 年、ギリシャ)、アルキメデス (紀元前 287 年 - 212 年 、イタリア)、アーリヤバタ (紀元前 476 年、インド)、アル・カラジ (1019 年、ペルシャ)、イブン・アル・ハイサム (965 年 - 1039 年、イラク) などがいます。
16世紀後半から17世紀初頭にかけて、数学者たちは目覚ましい進歩を遂げた。西洋では、イギリスのトーマス・ハリオット (1560年~1621年)、ドイツのヨハン・ファウルハーバー (1580年~1635年)、フランスのピエール・ド・フェルマー (1601年~1665年)、そして同じくフランスの数学者ブレーズ・パスカル (1623年~1662年)が重要な役割を果たした。
トーマス・ハリオットは、記号表記を用いてべき乗の和の公式を導出し、書き記した最初の人物であると思われるが、彼でさえ4乗までの和までしか計算していなかった。ヨハン・ファウルハーバーは、1631年の著書『アカデミア・アルゲブラエ』 の中で、それまでの誰よりもはるかに高い17乗までのべき乗の和の公式を示したが、一般的な公式は示さなかった。
1654年にブレーズ・パスカルは、( n +1) k +1を p =0,1,2,..., k の最初のn 個の正の整数のp 乗の和に関連付けるパスカルの恒等式 を証明しました。
スイスの数学者ヤコブ・ベルヌーイ(1654-1705)は、すべてのべき乗の和に対して統一的な公式を提供する定数列B 0 、B 1 、B 2 、...の存在を最初に認識した人物である。
ベルヌーイが、任意の正の整数c に対するc 乗の和の公式の係数を迅速かつ容易に計算するために必要なパターンを発見したときの喜びは、彼のコメントから見て取れる。彼は次のように書いている。
「この表のおかげで、最初の1000個の数の10乗を足し合わせると、91,409,924,241,424,243,424,241,924,242,500になるということが、15分半もかからずに分かりました。」 ベルヌーイの結果は、1713年に 『アルス・コンジェクタンディ』 に死後出版された。関隆和は ベルヌーイ数を独自に発見し、その結果は1年前の1712年に、やはり死後出版された。[ 2 ] しかし、関は定数列に基づく公式として自分の方法を提示しなかった。
ベルヌーイのべき乗和の公式は、現在までに最も有用で汎用性の高い定式化である。ベルヌーイの公式における係数は、アブラハム・ド・モアブル の提案に基づき、現在ではベルヌーイ数と呼ばれている。
ベルヌーイの公式は、べき乗の和を計算する驚くべき方法を発見したものの、ベルヌーイの公式を明言しなかったヨハン・ファウルハーバーにちなんで、ファウルハーバーの公式と呼ばれることもある。クヌース によると、ファウルハーバーの公式の厳密な証明は、 1834年にカール・ヤコビ によって初めて発表された[ 10 ]。 クヌースによるファウルハーバーの公式の詳細な研究は、次のように結論付けている。
「ファウルハーバーはベルヌーイ数を発見しなかった。つまり、定数の単一の列 B 0 、B 1 、B 2 、... が均一な ∑ n m = 1 m + 1 ( B 0 n m + 1 − ( m + 1 1 ) B 1 n m + ( m + 1 2 ) B 2 n m − 1 − ⋯ + ( − 1 ) m ( m + 1 m ) B m n ) {\textstyle \sum n^{m}={\frac {1}{m+1}}\left(B_{0}n^{m+1}-{\binom {m+1}{1}}B_{1}n^{m}+{\binom {m+1}{2}}B_{2}n^{m-1}-\cdots +(-1)^{m}{\binom {m+1}{m}}B_{m}n\right)} すべてのべき乗の和について。例えば、 Σ n m の公式をN の多項式から n の多項式に 変換した後、係数のほぼ半分がゼロになったという事実について、彼は一度も言及しなかった 。 上記でクヌートは、B 1 − {\displaystyle B_{1}^{-}} ; 代わりにB 1 + {\displaystyle B_{1}^{+}} この式は引き算を回避します。
∑ n m = 1 m + 1 ( B 0 n m + 1 + ( m + 1 1 ) B 1 + n m + ( m + 1 2 ) B 2 n m − 1 + ⋯ + ( m + 1 m ) B m n ) 。 {\textstyle \sum n^{m}={\frac {1}{m+1}}\left(B_{0}n^{m+1}+{\binom {m+1}{1}}B_{1}^{+}n^{m}+{\binom {m+1}{2}}B_{2}n^{m-1}+\cdots +{\binom {m+1}{m}}B_{m}n\right).}
定義 過去300年の間に、ベルヌーイ数の多くの特徴付けが発見されており、それぞれがこれらの数を導入するために使用できます。ここでは、最も有用な4つの特徴付けのみを紹介します。
4つのアプローチの等価性 の証明については、 Ireland & Rosen (1990) またはConway & Guy (1996) を参照してください。
再帰的定義 ベルヌーイ数は和の公式に従う[ 1 ]
∑ k = 0 m ( m + 1 k ) B k − = δ m 、 0 ∑ k = 0 m ( m + 1 k ) B k + = m + 1 {\displaystyle {\begin{aligned}\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{-{}}&=\delta _{m,0}\\\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{+{}}&=m+1\end{aligned}}} どこm = 0 、 1 、 2... {\displaystyle m=0,1,2...} δ はクロネッカーデルタ を表す。
これらのうち最初のものは、式(m > 1の場合)として [ 16 ]と表記されることがある。 ( B + 1 ) m − B m = 0 、 {\displaystyle (B+1)^{m}-B_{m}=0,} ここで、べき乗は二項定理を用いて正式に展開され、B k {\displaystyle B^{k}} に置き換えられますB k {\displaystyle B_{k}} 。
解決するB m ∓ {\displaystyle B_{m}^{\mp {}}} 再帰式[ 17 ]を与える
B m − = δ m 、 0 − ∑ k = 0 m − 1 ( m k ) B k − m − k + 1 B m + = 1 − ∑ k = 0 m − 1 ( m k ) B k + m − k + 1 。 {\displaystyle {\begin{aligned}B_{m}^{-{}}&=\delta _{m,0}-\sum _{k=0}^{m-1}{\binom {m}{k}}{\frac {B_{k}^{-{}}}{m-k+1}}\\B_{m}^{+}&=1-\sum _{k=0}^{m-1}{\binom {m}{k}}{\frac {B_{k}^{+}}{m-k+1}}.\end{aligned}}}
明示的な定義 1893年にルイ・ザールシュッツは ベルヌーイ数に関する明示的な公式を合計38個挙げたが、[ 18 ] 通常は古い文献を参照している。そのうちの1つは(m ≥ 1 {\displaystyle m\geq 1} ):
B m − = ∑ k = 0 m 1 k + 1 ∑ j = 0 k ( k j ) ( − 1 ) j j m B m + = ∑ k = 0 m 1 k + 1 ∑ j = 0 k ( k j ) ( − 1 ) j ( j + 1 ) m 。 {\displaystyle {\begin{aligned}B_{m}^{-}&=\sum _{k=0}^{m}{\frac {1}{k+1}}\sum _{j=0}^{k}{\binom {k}{j}}(-1)^{j}j^{m}\\B_{m}^{+}&=\sum _{k=0}^{m}{\frac {1}{k+1}}\sum _{j=0}^{k}{\binom {k}{j}}(-1)^{j}(j+1)^{m}.\end{aligned}}}
生成関数 指数生成関数 は
t e t − 1 = t 2 ( 布 t 2 − 1 ) = ∑ m = 0 ∞ B m − t m m ! t e t e t − 1 = t 1 − e − t = t 2 ( 布 t 2 + 1 ) = ∑ m = 0 ∞ B m + t m m ! 。 {\displaystyle {\begin{alignedat}{3}{\frac {t}{e^{t}-1}}&={\frac {t}{2}}\left(\operatorname {coth} {\frac {t}{2}}-1\right)&&=\sum _{m=0}^{\infty }{\frac {B_{m}^{-{}}t^{m}}{m!}}\\{\frac {te^{t}}{e^{t}-1}}={\frac {t}{1-e^{-t}}}&={\frac {t}{2}}\left(\operatorname {coth} {\frac {t}{2}}+1\right)&&=\sum _{m=0}^{\infty }{\frac {B_{m}^{+}t^{m}}{m!}}.\end{alignedat}}} ここでの置換はt → − t {\displaystyle t\to -t} 生成関数間の算術差はB m + {\displaystyle B_{m}^{+}} そしてB m − {\displaystyle B_{m}^{-}} t です。
(通常の)生成関数
z − 1 ψ 1 ( z − 1 ) = ∑ m = 0 ∞ B m + z m {\displaystyle z^{-1}\psi _{1}(z^{-1})=\sum _{m=0}^{\infty }B_{m}^{+}z^{m}} これは漸近級数 です。これはトリガンマ関数 ψ 1 を含みます。
積分表現 上記の母関数から、偶数ベルヌーイ数に対する以下の積分公式が得られる。
B 2 n = 4 n ( − 1 ) n + 1 ∫ 0 ∞ t 2 n − 1 e 2 π t − 1 d t {\displaystyle B_{2n}=4n(-1)^{n+1}\int _{0}^{\infty }{\frac {t^{2n-1}}{e^{2\pi t}-1}}\mathrm {d} t}
ベルヌーイ数とリーマンゼータ関数 ベルヌーイ数(B 1 には 1/2 を使用)は、負の実数のリーマンゼータ関数に関連しています。 偶数インデックスのベルヌーイ数の絶対値とリーマンゼータ関数との関係 ベルヌーイ数は、リーマンゼータ関数 を用いて次のように表すことができる。
B n + = − n ζ ( 1 − n ) {\displaystyle B_{n}^{+}=-n\,\zeta (1-n)\quad } n ≥ 1 の 場合 。ここで、ゼータ関数の引数は0 または負の値になります。ζ ( k ) {\displaystyle \zeta (k)} 負の偶数の場合はゼロ(自明なゼロ )、n>1 が奇数の場合は、ζ ( 1 − n ) {\displaystyle \zeta (1-n)} ゼロです。
ゼータ関数方程式 とガンマ反射式 を用いると、次の関係が得られる。
B 2 n = ( − 1 ) n + 1 2 ( 2 n ) ! ( 2 π ) 2 n ζ ( 2 n ) {\displaystyle B_{2n}={\frac {(-1)^{n+1}2(2n)!}{(2\pi )^{2n}}}\zeta (2n)\quad } n ≥ 1 の 整数の場合 。これで、ゼータ関数の引数は正になります。
すると、ζ → 1 ( n → ∞ ) とスターリングの公式 から、
| B 2 n | ~ 4 π n ( n π e ) 2 n {\displaystyle |B_{2n}|\sim 4{\sqrt {\pi n}}\left({\frac {n}{\pi e}}\right)^{2n}\quad } n → ∞ の場合 。
ベルヌーイ数の効率的な計算 ベルヌーイ数を計算するための秋山・谷川アルゴリズムの可視化 秋山と谷川は、ベルヌーイ数を計算するための単純な「三角形アルゴリズム」(パスカルの三角形 に類似)を提示している。 [ 20 ] まず、b 0 、 m = 1 m + 1 {\textstyle b_{0,m}={\frac {1}{m+1}}} (m ≥ 0 の場合)。すると、三角形の連続する項は、漸化式を用いて計算できる。
b n + 1 、 m = ( m + 1 ) ( b n 、 m − b n 、 m + 1 ) {\displaystyle b_{n+1,m}=(m+1)(b_{n,m}-b_{n,m+1})}
条件b n 、 0 {\textstyle b_{n,0}} n 番目のベルヌーイ数B n に対応する。[ 21 ] [ 22 ]
ブレント とハーベイは、ベルヌーイ数を計算するためのいくつかのアルゴリズムを提示しており、その中には秋山・谷川アルゴリズムよりも高速でメモリ使用量が少ない単純なアルゴリズムも含まれている。このアルゴリズムは漸化式を用いて接線数 T n を計算し、
T n = ( − 1 ) n − 1 2 2 n ( 2 2 n − 1 ) B 2 n 2 n {\displaystyle T_{n}=(-1)^{n-1}2^{2n}(2^{2n}-1){\frac {B_{2n}}{2n}}}
ベルヌーイ数を計算する。これは、KnuthとBuckholtzによって与えられた方法に関連している。[ 23 ] [ 24 ]
一部のアプリケーションでは、 p が 素数である場合、ベルヌーイ数B 0 からB p − 3 を法pで計算できると便利です。たとえば、 pに対して Vandiver 予想が 成り立つかどうかをテストしたり、 p が不規則な素数 かどうかを判定したりする場合などです。上記の再帰式を使用してこのような計算を実行することは現実的ではありません。少なくともp 2 回 の算術演算が必要になるためです。幸いなことに、より高速な方法が開発されており[ 25 ] 、 O ( p (log p ) 2 ) 回の演算のみで済みます(ビッグ O 表記を 参照)。
David Harvey [ 26 ] は、多数の小さな素数pに対して B n を 法p で計算し、中国剰余定理を 使用してB n を 再構築することによってベルヌーイ数を計算するアルゴリズムについて説明しています。Harvey は、このアルゴリズムの漸近的な 時間計算量は O ( n 2 log( n ) 2 + ε ) であると述べており、この実装は 他の方法に基づく実装よりも大幅に高速であると主張しています。この実装を使用して、Harvey はn = 10 8 の場合のB n を計算しました。Harvey の実装は、バージョン 3.1 以降SageMath に含まれています。それ以前には、Bernd Kellner [ 27 ] が 2002 年 12 月にn = 10 6の場合の B n を 完全な精度で計算し、Oleksandr Pavlyk [ 28 ] が 2008 年 4 月にMathematica を使用してn = 10 7 の場合のB n を計算しました。
* 桁数は、 B n が 正規化された科学的記数法 で実数として書かれている場合の 10 の指数として理解されます。
ベルヌーイ数の応用
漸近解析 数学におけるベルヌーイ数の最も重要な応用は、オイラー・マクローリンの公式での使用であると言えるでしょう 。f が十分に頻繁に微分可能な関数であると仮定すると、オイラー・マクローリンの公式は次のように書くことができます
∑ k = 1 b − 1 f ( k ) = ∫ 1 b f ( x ) d x + ∑ k = 1 m B k − k ! ( f ( k − 1 ) ( b ) − f ( k − 1 ) ( 1 ) ) + R − ( f 、 m ) 。 {\displaystyle \sum _{k=a}^{b-1}f(k)=\int _{a}^{b}f(x)\,dx+\sum _{k=1}^{m}{\frac {B_{k}^{-}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R_{-}(f,m).} この式は、 B − 1 = − 1 / 2 という慣例を前提としています。B + 1 = + 1 / 2 という慣例を用いると、式は次のようになります。
∑ k = 1 + 1 b f ( k ) = ∫ 1 b f ( x ) d x + ∑ k = 1 m B k + k ! ( f ( k − 1 ) ( b ) − f ( k − 1 ) ( 1 ) ) + R + ( f 、 m ) 。 {\displaystyle \sum _{k=a+1}^{b}f(k)=\int _{a}^{b}f(x)\,dx+\sum _{k=1}^{m}{\frac {B_{k}^{+}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R_{+}(f,m).} ここf ( 0 ) = f {\displaystyle f^{(0)}=f} (つまり、ゼロ次導関数)f {\displaystyle f} はただf {\displaystyle f} さらに、f ( − 1 ) {\displaystyle f^{(-1)}} の原始関数 を表すf {\displaystyle f} 微積分学の基本定理 により、
∫ 1 b f ( x ) d x = f ( − 1 ) ( b ) − f ( − 1 ) ( 1 ) 。 {\displaystyle \int _{a}^{b}f(x)\,dx=f^{(-1)}(b)-f^{(-1)}(a).} したがって、最後の式はさらに簡略化して、オイラー・マクローリンの公式の次の簡潔な形式にすることができる。
∑ k = 1 + 1 b f ( k ) = ∑ k = 0 m B k k ! ( f ( k − 1 ) ( b ) − f ( k − 1 ) ( 1 ) ) + R ( f 、 m ) 。 \displaystyle \sum _{k=a+1}^{b}f(k)=\sum _{k=0}^{m}{\frac {B_{k}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R(f,m).} この形式は、例えばゼータ関数の重要なオイラー・マクローリン展開の源泉となっている。
ζ ( s ) = ∑ k = 0 m B k + k ! s k − 1 ¯ + R ( s 、 m ) = B 0 0 ! s − 1 ¯ + B 1 + 1 ! s 0 ¯ + B 2 2 ! s 1 ¯ + ⋯ + R ( s 、 m ) = 1 s − 1 + 1 2 + 1 12 s + ⋯ + R ( s 、 m ) 。 {\displaystyle {\begin{aligned}\zeta (s)&=\sum _{k=0}^{m}{\frac {B_{k}^{+}}{k!}}s^{\overline {k-1}}+R(s,m)\\&={\frac {B_{0}}{0!}}s^{\overline {-1}}+{\frac {B_{1}^{+}}{1!}}s^{\overline {0}}+{\frac {B_{2}}{2!}}s^{\overline {1}}+\cdots +R(s,m)\\&={\frac {1}{s-1}}+{\frac {1}{2}}+{\frac {1}{12}}s+\cdots +R(s,m).\end{aligned}}} ここで、s k は 上昇階乗 を表します。
ベルヌーイ数は、他の種類の漸近展開 でも頻繁に使用されます。次の例は、ディガンマ関数 ψ の古典的なポアンカレ型の漸近展開です。
ψ ( z ) ~ ln z − ∑ k = 1 ∞ B k + k z k {\displaystyle \psi (z)\sim \ln z-\sum _{k=1}^{\infty }{\frac {B_{k}^{+}}{kz^{k}}}}
べき乗の合計 ベルヌーイ数は、最初の n 個の正の整数の m 乗の和の閉形式表現において重要な役割を果たします。m 、n ≥ 0 に対して、次のように定義 し ます 。
S m ( n ) = ∑ k = 1 n k m = 1 m + 2 m + ⋯ + n m 。 {\displaystyle S_{m}(n)=\sum _{k=1}^{n}k^{m}=1^{m}+2^{m}+\cdots +n^{m}.} この式は常に、 n に関する次数m + 1の 多項式 として書き換えることができます。これらの多項式の係数は 、ベルヌーイの公式 によってベルヌーイ数と関連付けられます。
S m ( n ) = 1 m + 1 ∑ k = 0 m ( m + 1 k ) B k + n m + 1 − k = m ! ∑ k = 0 m B k + n m + 1 − k k ! ( m + 1 − k ) ! 、 {\displaystyle S_{m}(n)={\frac {1}{m+1}}\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{+}n^{m+1-k}=m!\sum _{k=0}^{m}{\frac {B_{k}^{+}n^{m+1-k}}{k!(m+1-k)!}},} ここで( m + 1 k ) は二項係数 を表します。
例えば、m を 1 とすると、三角数は 0, 1, 3, 6, ...になります 。OEIS : A000217 。
1 + 2 + ⋯ + n = 1 2 ( B 0 n 2 + 2 B 1 + n 1 ) = 1 2 ( n 2 + n ) 。 {\displaystyle 1+2+\cdots +n={\frac {1}{2}}(B_{0}n^{2}+2B_{1}^{+}n^{1})={\tfrac {1}{2}}(n^{2}+n).} m を 2とすると、正四角錐数 0, 1, 5, 14, ... OEIS : A000330 が得られます。
1 2 + 2 2 + ⋯ + n 2 = 1 3 ( B 0 n 3 + 3 B 1 + n 2 + 3 B 2 n 1 ) = 1 3 ( n 3 + 3 2 n 2 + 1 2 n ) 。 {\displaystyle 1^{2}+2^{2}+\cdots +n^{2}={\frac {1}{3}}(B_{0}n^{3}+3B_{1}^{+}n^{2}+3B_{2}n^{1})={\tfrac {1}{3}}\left(n^{3}+{\tfrac {3}{2}}n^{2}+{\tfrac {1}{2}}n\right).} ベルヌーイ数に関して別の表記法を用いる著者もおり、ベルヌーイの公式を次のように表している。
S m ( n ) = 1 m + 1 ∑ k = 0 m ( − 1 ) k ( m + 1 k ) B k − n m + 1 − k 。 {\displaystyle S_{m}(n)={\frac {1}{m+1}}\sum _{k=0}^{m}(-1)^{k}{\binom {m+1}{k}}B_{k}^{-{}}n^{m+1-k}.} ベルヌーイの公式は、べき乗の和 を計算する驚くべき方法も発見したヨハン・ファウルハーバー にちなんで、ファウルハーバーの公式 と呼ばれることもある。
ファウルハーバーの公式は、V. GuoとJ. Zengによってq- 類似の公式 に一般化された。[ 31 ]
テイラーシリーズ ベルヌーイ数は、多くの三角関数 や双曲線関数 のテイラー級数 展開に現れる。
タン x = 1 x ∑ n = 1 ∞ ( − 1 ) n − 1 2 2 n ( 2 2 n − 1 ) B 2 n ( 2 n ) ! x 2 n − 1 、 | x | < π 2 。 ベビーベッド x = 1 x ∑ n = 0 ∞ ( − 1 ) n B 2 n ( 2 x ) 2 n ( 2 n ) ! 、 0 < | x | < π 。 タン x = 1 x ∑ n = 1 ∞ 2 2 n ( 2 2 n − 1 ) B 2 n ( 2 n ) ! x 2 n − 1 、 | x | < π 2 。 布 x = 1 x ∑ n = 0 ∞ B 2 n ( 2 x ) 2 n ( 2 n ) ! 、 0 < | x | < π 。 {\displaystyle {\begin{aligned}\tan x&={\hphantom {1 \over x}}\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}2^{2n}(2^{2n}-1)B_{2n}}{(2n)!}}\;x^{2n-1},&&\left|x\right|<{\frac {\pi }{2}}.\\\cot x&={1 \over x}\sum _{n=0}^{\infty }{\frac {(-1)^{n}B_{2n}(2x)^{2n}}{(2n)!}},&0<&|x|<\pi .\\\tanh x&={\hphantom {1 \over x}}\sum _{n=1}^{\infty }{\frac {2^{2n}(2^{2n}-1)B_{2n}}{(2n)!}}\;x^{2n-1},&&|x|<{\frac {\pi }{2}}.\\\coth x&={1 \over x}\sum _{n=0}^{\infty }{\frac {B_{2n}(2x)^{2n}}{(2n)!}},&0<&|x|<\pi .\end{aligned}}}
ローランシリーズ ベルヌーイ数は次のローラン級数 に現れる:
ディガンマ関数 :ψ ( z ) = ln z − ∑ k = 1 ∞ B k + k z k {\displaystyle \psi (z)=\ln z-\sum _{k=1}^{\infty }{\frac {B_{k}^{+{}}}{kz^{k}}}}
トポロジーでの使用 平行化可能な多様体 を境界とするエキゾチックな (4 n − 1) 球面 の微分同相類の巡回群の位数に関する Kervaire–Milnor 公式には、ベルヌーイ数 が含まれる。ES n を n ≥ 2 の場合のそのようなエキゾチックな球面の数とすると、
ES n = ( 2 2 n − 2 − 2 4 n − 3 ) 分子 ( B 4 n 4 n ) 。 {\displaystyle {\textit {ES}}_{n}=(2^{2n-2}-2^{4n-3})\operatorname {Numerator} \left({\frac {B_{4n}}{4n}}\right).} 次元 4nの滑らかな 向き付けられた 閉多様体 のL 種数 に関するヒルツェブルフの 署名定理 もベルヌーイ数に関係している。
組み合わせ数との関連性 ベルヌーイ数と様々な組み合わせ数との関連性は、有限差分の古典理論と、ベルヌーイ数を基本的な組み合わせ原理である包含排除原理 の一例として組み合わせ論的に解釈することに基づいている。
ウォルピツキー数との関連性 定義は1883年にユリウス・ヴォルピツキーによって考案された。基本的な算術演算の他に、階乗関数n ! とべき乗関数k m のみが用いられる。符号なしヴォルピツキー数は次のように定義される。
W n 、 k = ∑ v = 0 k ( − 1 ) v + k ( v + 1 ) n ( k v ) 。 {\displaystyle W_{n,k}=\sum _{v=0}^{k}(-1)^{v+k}(v+1)^{n}{k \choose v}\ .} これらは第二種スターリング数 によっても表現できる。
W n 、 k = k ! { n + 1 k + 1 } 。 {\displaystyle W_{n,k}=k!\left\{{n+1 \atop k+1}\right\}.} ベルヌーイ数は 、 調和数列 1, 1/2 , 1/3 , ... で重み付けされたウォルピツキー数 の包含排除和として導入される 。
B n = ∑ k = 0 n ( − 1 ) k W n 、 k k + 1 = ∑ k = 0 n 1 k + 1 ∑ v = 0 k ( − 1 ) v ( v + 1 ) n ( k v ) 。 {\displaystyle B_{n}=\sum _{k=0}^{n}(-1)^{k}{\frac {W_{n,k}}{k+1}}\ =\ \sum _{k=0}^{n}{\frac {1}{k+1}}\sum _{v=0}^{k}(-1)^{v}(v+1)^{n}{k \choose v}\ .} B 0 = 1B 1 = 1 − 1 / 2 B 2 = 1 − 3 / 2 + 2 / 3 B 3 = 1 − 7 / 2 + 12 / 3 − 6 / 4 B 4 = 1 − 15 / 2 + 50 / 3 − 60 / 4 + 24 / 5 B 5 = 1 − 31 / 2 + 180 / 3 − 390 / 4 + 360 / 5 − 120 / 6 B 6 = 1 − 63 / 2 + 602 / 3 − 2100 / 4 + 3360 / 5 − 2520 / 6 + 720 / 7 この表現では、B + 1 = + 1 / 2 と なります。
数列s n , n ≥ 0 を考えます。Worpitzky の数OEIS : A028246 、OEIS : A163626 をs 0 、s 0 、s 1 、s 0 、s 1 、s 2 、s 0 、s 1 、s 2 、s 3 、...に適用すると、 s n に適用した秋山-谷川変換と同一になります( 第 1 種のスターリング数との関連を 参照)。これは次の表から確認できます。
最初の行はs 0 、s 1 、s 2 、s 3 、s 4 を表します。
したがって、第2の分数オイラー数OEIS : A198631 ( n ) / OEIS : A006519 ( n + 1 ) については、次のようになります。
E 0 = 1E 1 = 1 − 1 / 2 E 2 = 1 − 3 / 2 + 2 / 4 E 3 = 1 − 7 / 2 + 12 / 4 − 6 / 8 E 4 = 1 − 15 / 2 + 50 / 4 − 60 / 8 + 24 / 16 E 5 = 1 − 31 / 2 + 180 / 4 − 390 / 8 + 360 / 16 − 120 / 32 E 6 = 1 − 63 / 2 + 602 / 4 − 2100 / 8 + 3360 / 16 − 2520 / 32 + 720 / 64 ベルヌーイ数をウォルピツキー数で表す2番目の式は、n ≥ 1 の場合である。
B n = n 2 n + 1 − 2 ∑ k = 0 n − 1 ( − 2 ) − k W n − 1 、 k 。 {\displaystyle B_{n}={\frac {n}{2^{n+1}-2}}\sum _{k=0}^{n-1}(-2)^{-k}\,W_{n-1,k}.} 簡略化された第2ベルヌーイ数のウォルピツキー表現は次のとおりです。
OEIS : A164555 ( n + 1 ) / OEIS : A027642 ( n + 1 ) = n + 1 / 2 n + 2 − 2 × OEIS : A198631 ( n ) / OEIS : A006519 ( n + 1 )
これは、第2ベルヌーイ数と第2分数オイラー数を関連付けるものです。その始まりは次のとおりです。
1 / 2 , 1 / 6 , 0, − 1 / 30 , 0, 1 / 42 , ... = ( 1 / 2 , 1 / 3 , 3 / 14 , 2 / 15 , 5 / 62 , 1 / 21 , ...) × (1, 1 / 2 , 0, − 1 / 4 , 0, 1 / 2 , ...)最初の括弧内の分子はOEIS : A111701 です(第 1 種のスターリング数との関連を 参照)。
第二種スターリング数との関連性 第2種スターリング数 S ( k , m ) は、 [ 33 ] のような性質を持つ。
x k = ∑ m = 0 k x m _ S ( k 、 m ) {\displaystyle x^{k}=\sum _{m=0}^{k}{x^{\underline {m}}}S(k,m)} ここで、x m は下降階乗関数 を表す。
ベルヌーイ多項式 B k ( x ) は次のように書くことができます。[ 34 ]
B k ( x ) = k ∑ m = 0 k − 1 ( x m + 1 ) S ( k − 1 、 m ) m ! + B k {\displaystyle B_{k}(x)=k\sum _{m=0}^{k-1}{\binom {x}{m+1}}S(k-1,m)m!+B_{k}} ここで、k = 0, 1, 2,... のB k はベルヌーイ数である。
二項係数 の以下の性質:
( x m ) = ( x + 1 m + 1 ) − ( x m + 1 ) {\displaystyle {\binom {x}{m}}={\binom {x+1}{m+1}}-{\binom {x}{m+1}}} したがって、
x k = B k + 1 ( x + 1 ) − B k + 1 ( x ) k + 1 。 {\displaystyle x^{k}={\frac {B_{k+1}(x+1)-B_{k+1}(x)}{k+1}}.} ベルヌーイ多項式については、次のことも成り立つ。[ 34 ]
B k ( x ) = ∑ n = 0 k ( k n ) B n x k − n 。 {\displaystyle B_{k}(x)=\sum _{n=0}^{k}{\binom {k}{n}}B_{n}x^{k-n}.} ( x m + 1 ) におけるx の係数は (−1) m / m + 1 です。
最初の式のx の係数はB k ( x ) {\displaystyle B_{k}(x)} はk ∑ m = 0 k − 1 ( − 1 ) m m + 1 S ( k − 1 、 m ) m ! {\displaystyle k\sum _{m=0}^{k-1}{\frac {(-1)^{m}}{m+1}}S(k-1,m)m!} 一方、2番目の表現ではk B k − 1 。 {\displaystyle kB_{k-1}.} 交換するk − 1 {\displaystyle k-1} によるk {\displaystyle k} これにより、以下が得られます。
B k = ∑ m = 0 k ( − 1 ) m m ! m + 1 S ( k 、 m ) {\displaystyle B_{k}=\sum _{m=0}^{k}(-1)^{m}{\frac {m!}{m+1}}S(k,m)} (結果としてB 1 = + 1 / 2 となる) これはベルヌーイ数の明示的な公式であり、フォン・シュタウト・クラウゼン定理を 証明 するために使用できる。[ 35 ] [ 36 ] [ 37 ]
第一種スターリング数との関連性 符号なし第1種スターリング数 [ nm ] とベルヌーイ数(B1 = +1/2 ) を 関連付ける2つの主要な公式は以下 の とおり で ある。
1 m ! ∑ k = 0 m ( − 1 ) k [ m + 1 k + 1 ] B k = 1 m + 1 、 {\displaystyle {\frac {1}{m!}}\sum _{k=0}^{m}(-1)^{k}\left[{m+1 \atop k+1}\right]B_{k}={\frac {1}{m+1}},} そしてこの和の逆数(n ≥ 0 、m ≥ 0 の場合)
1 m ! ∑ k = 0 m ( − 1 ) k [ m + 1 k + 1 ] B n + k = A n 、 m 。 {\displaystyle {\frac {1}{m!}}\sum _{k=0}^{m}(-1)^{k}\left[{m+1 \atop k+1}\right]B_{n+k}=A_{n,m}.} ここで、A n 、m は有理数である秋山-谷川数であり、最初の数個が次の表に示されています。
秋山・谷川数は単純な漸化式を満たし、これを利用してベルヌーイ数を反復的に計算することができる。OEIS : A051714 / OEIS : A051715 を 参照のこと。
オートシーケンス とは、逆二項変換が符号付きシーケンスと等しくなるシーケンスのことです。主対角線が 0 列 = OEIS : A000004 の場合、オートシーケンスは第 1 種です。例: OEIS : A000045 、フィボナッチ数列。主対角線が最初の上側対角線に 2 を掛けたものの場合、それは第 2 種です。例: OEIS : A164555 / OEIS : A027642 、第 2 ベルヌーイ数列 ( OEIS : A190339を参照)。2 − n = 1 / OEIS : A000079に 秋山 -谷川変換を適用すると、 OEIS : A198631 ( n ) / OEIS : A06519 ( n + 1 )になります。したがって:
OEIS : A209308 およびOEIS : A227577 を参照してください。OEIS : A198631 ( n ) / OEIS : A006519 ( n + 1 ) は、 第 2 種 (分数) オイラー数であり、第 2 種のオートシーケンスです。
( OEIS : A164555 ( n + 2 ) / OEIS : A027642 ( n + 2 ) = 1 / 6 , 0, − 1 / 30 , 0, 1 / 42 , ... ) × ( 2 n + 3 − 2 / n + 2 = 3, 14 / 3 , 15 / 2 , 62 / 5 , 21, ... ) = OEIS : A198631 ( n + 1 ) / OEIS : A006519 ( n + 2 ) = 1 / 2 , 0, − 1 / 4 , 0 , 1 / 2 , ... 。 OEIS : A027641 / OEIS : A027642 にも有用です(Worpitzky 番号との関連を 参照)。
パスカルの三角形との関連性パスカルの三角形とベルヌーイ数を結びつける公式がある[ b ]
B n + = | A n | ( n + 1 ) ! {\displaystyle B_{n}^{+}={\frac {|A_{n}|}{(n+1)!}}~~~} 行と列の定義付き[ A n ] 私 、 k := { 0 もし k > 1 + 私 ( 私 + 1 k − 1 ) さもないと {\displaystyle [A_{n}]_{i,k}:={\begin{cases}0&{\text{if }}k>1+i\\{i+1 \choose k-1}&{\text{otherwise}}\end{cases}}} どこ| A n | {\displaystyle |A_{n}|} n×nヘッセンベルグ行列の行列式は パスカルの三角形 の一部である。
例:
B 6 + = | 1 2 0 0 0 0 1 3 3 0 0 0 1 4 6 4 0 0 1 5 10 10 5 0 1 6 15 20 15 6 1 7 21 35 35 21 | 7 ! = 120 5040 = 1 42 {\displaystyle B_{6}^{+}={\frac {\begin{vmatrix}1&2&0&0&0&0\\1&3&3&0&0&0\\1&4&6&4&0&0\\1&5&10&10&5&0\\1&6&15&20&15&6\\1&7&21&35&35&21\end{vmatrix}}{7!}}={\frac {120}{5040}}={\frac {1}{42}}} [ 38 ]
オイラー数との関連性 オイラー数 ⟨n m⟩ と ベルヌーイ数を結びつける 公式が存在する。
∑ m = 0 n ( − 1 ) m ⟨ n m ⟩ = 2 n + 1 ( 2 n + 1 − 1 ) B n + 1 n + 1 、 ∑ m = 0 n ( − 1 ) m ⟨ n m ⟩ ( n m ) − 1 = ( n + 1 ) B n 。 {\displaystyle {\begin{aligned}\sum _{m=0}^{n}(-1)^{m}\left\langle {n \atop m}\right\rangle &=2^{n+1}(2^{n+1}-1){\frac {B_{n+1}}{n+1}},\\\sum _{m=0}^{n}(-1)^{m}\left\langle {n \atop m}\right\rangle {\binom {n}{m}}^{-1}&=(n+1)B_{n}.\end{aligned}}} B 1 を 1 / 2 に設定した場合、両方の式 は n ≥ 0 に対して有効です。B 1 を− 1 / 2 に 設定した場合、 それぞれn ≥ 1 およびn ≥ 2 に対してのみ有効です。
二分木表現 スターリング多項式σ n ( x ) は 、ベルヌーイ数とB n = n ! σ n (1)の関係にある。SC Woon は、 σ n (1) を二分木として計算するアルゴリズムを説明した。 [ 39 ]
Woon の再帰アルゴリズム ( n ≥ 1 の場合) は、ルートノードN = [1,2] に を割り当てることから始まります。ツリーのノードN = [ a 1 , a 2 , ..., a k ]が与えられた場合、ノードの左の子はL ( N ) = [− a 1 , a 2 + 1, a 3 , ..., a k ] であり、右の子はR ( N ) = [ a 1 , 2, a 2 , ..., a k ] です。ノードN = [ a 1 , a 2 , ..., a k ]は、上記のツリーの初期部分では±[ a 2 , ..., a k ] と書かれ、± はa 1 の符号を表します。
ノードNが与えられた場合、 N の階乗は次のように定義されます。
N ! = 1 1 ∏ k = 2 長さ ( N ) 1 k ! 。 {\displaystyle N!=a_{1}\prod _{k=2}^{\operatorname {length} (N)}a_{k}!.} 固定されたツリーレベル n のノード N に限定すると、 1 / N ! の 合計はσ n (1) となり、したがって
B n = ∑ 樹木レベル n N ノード n ! N ! 。 {\displaystyle B_{n}=\sum _{\stackrel {N{\text{ node of}}}{{\text{ tree-level }}n}}{\frac {n!}{N!}}.} 例えば:
B 1 = 1!( 1 / 2 ! ) B 2 = 2!(− 1 / 3! + 1 / 2 ! 2 ! ) B 3 = 3!( 1 / 4! − 1 / 2!3! − 1 / 3!2! + 1 / 2!2!2! )
積分表現と積分継続 積分
b ( s ) = 2 e s 私 π / 2 ∫ 0 ∞ s t s 1 − e 2 π t d t t = s ! 2 s − 1 ζ ( s ) π s ( − 私 ) s = 2 s ! ζ ( s ) ( 2 π 私 ) s {\displaystyle b(s)=2e^{si\pi /2}\int _{0}^{\infty }{\frac {st^{s}}{1-e^{2\pi t}}}{\frac {dt}{t}}={\frac {s!}{2^{s-1}}}{\frac {\zeta (s)}{{}\pi ^{s}{}}}(-i)^{s}={\frac {2s!\zeta (s)}{(2\pi i)^{s}}}} n > 0 の場合、 b (2 n ) = B 2 n という特別な値を持つ。
例えば、b (3) = 3/2 ζ (3) π - 3i 、 b ( 5) = -15/2 ζ ( 5 ) π - 5iで ある 。 ここで、ζ はリーマンゼータ関数 、iは 虚数 単位 である。レオンハルト・オイラー(『 全集』 第1シリーズ、第10巻、351ページ)は これらの数値を考察し、計算を行った。
p = 3 2 π 3 ( 1 + 1 2 3 + 1 3 3 + ⋯ ) = 0.0581522 … q = 15 2 π 5 ( 1 + 1 2 5 + 1 3 5 + ⋯ ) = 0.0254132 … {\displaystyle {\begin{aligned}p&={\frac {3}{2\pi ^{3}}}\left(1+{\frac {1}{2^{3}}}+{\frac {1}{3^{3}}}+\cdots \right)=0.0581522\ldots \\q&={\frac {15}{2\pi ^{5}}}\left(1+{\frac {1}{2^{5}}}+{\frac {1}{3^{5}}}+\cdots \right)=0.0254132\ldots \end{aligned}}} 別の同様の積分 表現は
b ( s ) = − e s 私 π / 2 2 s − 1 ∫ 0 ∞ s t s シン π t d t t = 2 e s 私 π / 2 2 s − 1 ∫ 0 ∞ e π t s t s 1 − e 2 π t d t t 。 {\displaystyle b(s)=-{\frac {e^{si\pi /2}}{2^{s}-1}}\int _{0}^{\infty }{\frac {st^{s}}{\sinh \pi t}}{\frac {dt}{t}}={\frac {2e^{si\pi /2}}{2^{s}-1}}\int _{0}^{\infty }{\frac {e^{\pi t}st^{s}}{1-e^{2\pi t}}}{\frac {dt}{t}}.}
オイラー数とπ の関係オイラー数は 、 ベルヌーイ数と密接に関係する整数列です。ベルヌーイ数とオイラー数の漸近展開を比較すると、オイラー数E 2 n はベルヌーイ数B 2 n よりも約2 / π (4 2 n − 2 2 n ) 倍大きいことがわかります。したがって、次のようになります。
π ~ 2 ( 2 2 n − 4 2 n ) B 2 n E 2 n 。 {\displaystyle \pi \sim 2(2^{2n}-4^{2n}){\frac {B_{2n}}{E_{2n}}}.} この漸近方程式は、πが ベルヌーイ数とオイラー数の共通根に含まれることを示している。実際、πは これらの有理数近似から計算することができる。
ベルヌーイ数はオイラー数で表すことができ、またその逆も可能です。nが奇数の場合、B n = E n = 0 (B 1 を除く) となるため、 nが 偶数 の場合 のみ を考慮 すれ ば十分 です。
B n = ∑ k = 0 n − 1 ( n − 1 k ) n 4 n − 2 n E k n = 2 、 4 、 6 、 … E n = ∑ k = 1 n ( n k − 1 ) 2 k − 4 k k B k n = 2 、 4 、 6 、 … {\displaystyle {\begin{aligned}B_{n}&=\sum _{k=0}^{n-1}{\binom {n-1}{k}}{\frac {n}{4^{n}-2^{n}}}E_{k}&n&=2,4,6,\ldots \\[6pt]E_{n}&=\sum _{k=1}^{n}{\binom {n}{k-1}}{\frac {2^{k}-4^{k}}{k}}B_{k}&n&=2,4,6,\ldots \end{aligned}}} これらの変換式は、ベルヌーイ数とオイラー数の間の関係を表しています。しかし、より重要なのは、両方の種類の数に共通する深い算術的根源があり、それはπにも密接に関係する、より基本的な数列を通して表現できます。これらの数は、 n ≥ 1 に対して次のように定義されます[ 40 ] [ 41 ]
S n = 2 ( 2 π ) n ∑ k = 0 ∞ ( − 1 ) k n ( 2 k + 1 ) n = 2 ( 2 π ) n リム K → ∞ ∑ k = − K K ( 4 k + 1 ) − n 。 {\displaystyle S_{n}=2\left({\frac {2}{\pi }}\right)^{n}\sum _{k=0}^{\infty }{\frac {(-1)^{kn}}{(2k+1)^{n}}}=2\left({\frac {2}{\pi }}\right)^{n}\lim _{K\to \infty }\sum _{k=-K}^{K}(4k+1)^{-n}.} これらの数の魔法は、それらが有理数になるという事実にある。これは、レオンハルト・オイラー が画期的な論文「逆数の級数の和について」( De summis serierum reciprocarum )で初めて証明し、それ以来数学者たちを魅了してきた。[ 42 ] これらの数の最初のいくつかは
S n = 1 、 1 、 1 2 、 1 3 、 5 24 、 2 15 、 61 720 、 17 315 、 277 8064 、 62 2835 、 … {\displaystyle S_{n}=1,1,{\frac {1}{2}},{\frac {1}{3}},{\frac {5}{24}},{\frac {2}{15}},{\frac {61}{720}},{\frac {17}{315}},{\frac {277}{8064}},{\frac {62}{2835}},\ldots } ( OEIS : A099612 / OEIS : A099617 )これらはsec x + tan x の展開における係数です。
ベルヌーイ数とオイラー数は、これらの数の特別な見方 として理解することができ、数列S n から選択され、特別な用途で使用するためにスケーリングされています。
B n = ( − 1 ) ⌊ n 2 ⌋ [ n 平 ] n ! 2 n − 4 n S n 、 n = 2 、 3 、 … E n = ( − 1 ) ⌊ n 2 ⌋ [ n 平 ] n ! S n + 1 n = 0 、 1 、 … {\displaystyle {\begin{aligned}B_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]{\frac {n!}{2^{n}-4^{n}}}\,S_{n}\ ,&n&=2,3,\ldots \\E_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]n!\,S_{n+1}&n&=0,1,\ldots \end{aligned}}} 式 [ n even] は、 n が偶数の場合は 1、それ以外の場合は 0 の値をとります(アイバーソン括弧 )。
これらの恒等式は、 この節の冒頭にあるベルヌーイ数とオイラー数の商が、n が偶数の場合のR n = 2 S n / S n + 1 の 特殊なケースであることを示しています。R n は π の有理数近似であり、連続する 2 つの項は常にπ の真の値を囲みます。n = 1 から始まる数列 は ( OEIS : A132049 / OEIS : A132050 ) です。
2 、 4 、 3 、 16 5 、 25 8 、 192 61 、 427 136 、 4352 1385 、 12465 3968 、 158720 50521 、 … ⟶ π 。 {\displaystyle 2,4,3,{\frac {16}{5}},{\frac {25}{8}},{\frac {192}{61}},{\frac {427}{136}},{\frac {4352}{1385}},{\frac {12465}{3968}},{\frac {158720}{50521}},\ldots \quad \longrightarrow \pi .} これらの有理数は、前述のオイラーの論文の最後の段落にも登場する。
シーケンスOEIS : A046978 ( n + 2 ) / OEIS : A016116 ( n + 1 ) に対する秋山-谷川変換を考えます。
2番目から、1列目の分子はオイラーの公式の分母になります。1列目は − 1 / 2 × OEIS : A163982 です。
アルゴリズム的視点:ザイデル三角形 数列S n には、もう一つ予期せぬ重要な性質があります。それは、 S n +1 の分母が階乗n ! を割り切るということです。言い換えれば、T n = S n + 1 n !という数 (オイラーのジグザグ数 と呼ばれることもあります)は整数です。
T n = 1 、 1 、 1 、 2 、 5 、 16 、 61 、 272 、 1385 、 7936 、 50521 、 353792 、 … n = 0 、 1 、 2 、 3 、 … {\displaystyle T_{n}=1,\,1,\,1,\,2,\,5,\,16,\,61,\,272,\,1385,\,7936,\,50521,\,353792,\ldots \quad n=0,1,2,3,\ldots } ( OEIS : A000111 )。( OEIS : A253671 )を参照。それらの指数生成関数は、 正割関数 と正接 関数の和である。
∑ n = 0 ∞ T n x n n ! = タン ( π 4 + x 2 ) = 秒 x + タン x {\displaystyle \sum _{n=0}^{\infty }T_{n}{\frac {x^{n}}{n!}}=\tan \left({\frac {\pi }{4}}+{\frac {x}{2}}\right)=\sec x+\tan x} 。したがって、上記のベルヌーイ数とオイラー数の表現は、この数列を用いて次のように書き換えることができる。
B n = ( − 1 ) ⌊ n 2 ⌋ [ n 平 ] n 2 n − 4 n T n − 1 n ≥ 2 E n = ( − 1 ) ⌊ n 2 ⌋ [ n 平 ] T n n ≥ 0 {\displaystyle {\begin{aligned}B_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]{\frac {n}{2^{n}-4^{n}}}\,T_{n-1}\ &n&\geq 2\\E_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]T_{n}&n&\geq 0\end{aligned}}} これらの恒等式により、ベルヌーイ数とオイラー数を簡単に計算できます。オイラー数E 2 n はT 2 n からすぐに得られ、ベルヌーイ数B 2 n は、簡単なシフト操作によってT 2 n - 1 から得られる分数であり、有理数演算は不要です。
残る課題は、数値T n を計算する便利な方法を見つけることである。しかし、1877 年にフィリップ・ルートヴィヒ・フォン・ザイデルは、 T n を簡単に計算できる独創的なアルゴリズムを発表した。[ 43 ]
1 → 1 1 2 2 1 ← → 2 4 5 5 16 16 14 10 5 ← {\displaystyle {\begin{array}{crrrcc}{}&{}&{\color {red}1}&{}&{}&{}\\{}&{\rightarrow }&{\color {blue}1}&{\color {red}1}&{}\\{}&{\color {red}2}&{\color {blue}2}&{\color {blue}1}&{\leftarrow }\\{\rightarrow }&{\color {blue}2}&{\color {blue}4}&{\color {blue}5}&{\color {red}5}\\{\color {red}16}&{\color {blue}16}&{\color {blue}14}&{\color {blue}10}&{\color {blue}5}&{\leftarrow }\end{array}}} T n に対するセイデルのアルゴリズム
まず、0行目に1を記入し、k を現在記入中の行番号とします。 k が奇数の場合、行k -1 の左端の数字を行k の最初の位置に置き、左から右へ行を埋めていき、各エントリは左の数字と上の数字の合計となる。行の最後に、最後の数字を複製します。 k が偶数の場合は、反対方向にも同様の手順で進めてください。セイデルのアルゴリズムは実際にははるかに一般的であり(ドミニク・デュモン[ 44 ] の解説を参照)、その後何度か再発見された。
ザイデルのアプローチと同様に、DE クヌートと TJ バックホルツは数T 2 n の漸化式を与え、この方法を「整数に対する単純な演算のみを使用する電子計算機上で」 B 2 n とE 2 n を計算するために推奨した。[ 24 ]
VI アーノルド [ 45 ] はセイデルのアルゴリズムを再発見し、後にミラー、スローン、ヤングはブストロフェドン変換 という名前でセイデルのアルゴリズムを普及させた。
三角形の形状:
OEIS には、1 が 1 つ含まれる OEIS : A000657 と1 が 2 つ含まれるOEIS : A214267のみが含まれています。
以下の行に、補助的な1と1つの0を付加した分布を示します。
これはOEIS : A239005で、 OEIS : A008280 の符号付きバージョンです。主対角線はOEIS : A122045 です。主対角線はOEIS : A155585 です。中央の列はOEIS : A099023 です。行の合計: 1、1、-2、-5、16、61...。OEIS : A163747 を参照してください。以下の 1、1、0、-2、0、16、0 で始まる配列を参照してください。
秋山・谷川アルゴリズムをOEIS : A046978 ( n + 1 ) / OEIS : A016116 ( n ) に適用すると、次の結果が得られます。
1. 最初の列はOEIS : A122045 です。その二項変換は次のようになります。
この配列の最初の行はOEIS : A155585 です。増加する対角線の絶対値はOEIS : A008280 です。対角線の合計は− OEIS : A163747 ( n + 1 ) です。
2. 2列目は1 1 −1 −5 5 61 −61 −1385 1385... です。その二項変換は次のようになります。
この配列の最初の行は1 2 2 −4 −16 32 272 544 −7936 15872 353792 −707584... です。2 番目の二分法の絶対値は、1 番目の二分法の絶対値の 2 倍です。
OEIS : A046978 ( n ) / ( OEIS : A158780 ( n + 1 ) = abs( OEIS : A117575 ( n )) + 1 = 1, 2, 2, 3 / 2 , 1, 3 / 4 、に適用される秋山 ・谷川アルゴリズムを考えてみ ましょ う。 3 / 4 、 7 / 8 、1、 17 / 16 、 17 / 16 、 33 / 32 ... 。
絶対値がOEIS : A000111 である最初の列は、三角関数の分子である可能性があります。
OEIS : A163747 は第一種オートシーケンスです (主対角線はOEIS : A000004 です)。対応する配列は次のとおりです。
最初の 2 つの上側対角線は−1 3 −24 402... = (−1) n + 1 × OEIS : A002832 です。反対角線の合計は0 −2 0 10... = 2 × OEIS : A122045 ( n + 1) です。
− OEIS : A163982 は、例えばOEIS : A164555 / OEIS : A027642 のような第 2 種のオートシーケンスです。したがって、配列は次のようになります。
主対角線(ここでは2 −2 8 −92...) は、最初の上部の対角線(ここではOEIS : A099023 )の 2 倍です。反対角線の合計は2 0 −4 0... = 2 × OEIS : A155585 ( n + 1) です。OEIS : A163747 − OEIS : A163982 = 2 × OEIS : A122045 です。
組み合わせ論的視点:交互順列 1880年頃、ザイデルのアルゴリズムの発表から3年後、デジレ・アンドレは 、現在では古典的な組合せ解析の結果を証明した。[ 46 ] [ 47 ] アンドレは、三角関数 tan x とsec x のテイラー展開の最初の項を見て、ある発見をした。
タン x = x + 2 x 3 3 ! + 16 x 5 5 ! + 272 x 7 7 ! + 7936 x 9 9 ! + ⋯ 秒 x = 1 + x 2 2 ! + 5 x 4 4 ! + 61 x 6 6 ! + 1385 x 8 8 ! + 50521 x 10 10 ! + ⋯ {\displaystyle {\begin{aligned}\tan x&=x+{\frac {2x^{3}}{3!}}+{\frac {16x^{5}}{5!}}+{\frac {272x^{7}}{7!}}+{\frac {7936x^{9}}{9!}}+\cdots \\[6pt]\sec x&=1+{\frac {x^{2}}{2!}}+{\frac {5x^{4}}{4!}}+{\frac {61x^{6}}{6!}}+{\frac {1385x^{8}}{8!}}+{\frac {50521x^{10}}{10!}}+\cdots \end{aligned}}} 係数は、それぞれ奇数と偶数のインデックスのオイラージグザグ数 です。したがって、tan x + sec x の通常の展開の係数は有理数S n になります。
タン x + 秒 x = 1 + x + 1 2 x 2 + 1 3 x 3 + 5 24 x 4 + 2 15 x 5 + 61 720 x 6 + ⋯ {\displaystyle \tan x+\sec x=1+x+{\tfrac {1}{2}}x^{2}+{\tfrac {1}{3}}x^{3}+{\tfrac {5}{24}}x^{4}+{\tfrac {2}{15}}x^{5}+{\tfrac {61}{720}}x^{6}+\cdots } アンドレはその後、漸化式を用いて、奇数サイズの交代順列 は奇数インデックスのオイラー数(接線数とも呼ばれる)で列挙され、偶数サイズの交代順列は偶数インデックスのオイラー数(割線数とも呼ばれる)で列挙されることを示すことに成功した。
第1および第2ベルヌーイ数の算術平均は、対応するベルヌーイ数である。B 0 = 1 、B 1 = 0 、B 2 = 1 / 6 、 B 3 = 0 、 B 4 = − 1 / 30 、 OEIS : A176327 / OEIS : A027642 。 その 逆 秋山 -谷川 変換 の2 行 目 OEIS : A177427 を介して 、 バルマー 級数 OEIS : A061037 / OEIS : A061038 が 得 られる。
OEIS : A060819 ( n + 4 ) / OEIS : A145979 ( n )に秋山-谷川アルゴリズムを適用すると、B 1 を含まないベルヌーイ数OEIS : A027641 / OEIS : A027642 、OEIS : A164555 / OEIS : A027642 、またはOEIS : A176327 OEIS : A176289 が得られ、これらは固有ベルヌーイ数B i ( n ) と呼ばれます。
したがって、 OEIS : A145979 ( n )を介して、固有ベルヌーイ数とバルマー級数の間にもう 1 つのリンクがあります。
OEIS : A145979 ( n − 2 ) = 0, 2, 1, 6,... は非負数の順列です。
最初の行の項は f(n) = 1 / 2 + 1 / n + 2 です。 2、f(n) は第 2 種の自己数列です。 3/2、f(n) は逆二項変換により 3/2 −1/2 1/3 −1/4 1/5 ... = 1/2 + log 2 となります。
g(n) = 1/2 – 1 / (n+2) = 0, 1/6, 1/4, 3/10, 1/3 を考えます。秋山-タナギワ変換は次のようになります。
0, g(n) は第 2 種の自己数列である。
オイラーOEIS : A198631 ( n ) / OEIS : A006519 ( n + 1 ) は、 第 2 項 ( 1 / 2 ) を除いた分数 固有 オイラー 数 E i ( n ) = 1, 0, − 1 / 4 , 0, 1 / 2 , 0 , − 17 / 8 , 0, ...です 。 対応する秋山変換は次のとおりです。
最初の行はEu ( n ) です。Eu ( n ) の前にゼロが付くと、第一種のオートシーケンスになります。これはオレーム数と関連しています。2 行目の分子はOEIS : A069834の 前に 0 が付きます。差分表は次のとおりです。
ベルヌーイ数の算術的性質 ベルヌーイ数は、リーマンゼータ関数を用いて、整数n ≥ 0に対して B n = − nζ (1 − n )と表すことができる。ただし、 n = 0 の場合、式− nζ (1 − n ) は極限値として理解され、慣例としてB 1 = 1 / 2 が 使用される。これは、負の整数におけるゼータ関数の値と密接に関連している。そのため、ベルヌーイ数は深い算術的性質を持つことが期待され、実際に持っている。例えば、Agoh –Giuga 予想では、 p が素数であるのは、 pB p − 1 が法p で −1 と合同である場合のみであると仮定している。ベルヌーイ数の整除性は、Kummer の定理とHerbrand-Ribet の定理による強化によって 円分 体のイデアル類群と関連付けられ、 Ankeny–Artin–Chowla によって実二次体の類数と関連付けられている。
クンマーの定理 ベルヌーイ数は、クンマー の定理[ 48 ] によってフェルマーの最終定理 (FLT)と関連付けられています。クンマーの定理は次のように述べています。
奇素数p がベルヌーイ数 B 2 、B 4 、 ...、B p − 3 の分子のいずれも割り切れない場合、 x p + y p + z p = 0 は非ゼロ整数の解を持たない。 この性質を持つ素数は正則素数 と呼ばれます。クンマーのもう一つの古典的な結果は、次の合同式 です。[ 49 ]
p を 奇素数、b を p − 1 で割り切れない偶数とする。このとき、任意の非負整数kに対して、 B k ( p − 1 ) + b k ( p − 1 ) + b ≡ B b b ( モジュール p ) 。 {\displaystyle {\frac {B_{k(p-1)+b}}{k(p-1)+b}}\equiv {\frac {B_{b}}{b}}{\pmod {p}}.} これらの合同関係を一般化したものが、 p 進連続性と呼ばれる。
p 進連続性b 、m 、n が正の整数であり、m と nがp - 1で割り切れ ず、m≡n (mod p b -1 ( p -1)) である場合、
( 1 − p m − 1 ) B m m ≡ ( 1 − p n − 1 ) B n n ( モジュール p b ) 。 {\displaystyle (1-p^{m-1}){\frac {B_{m}}{m}}\equiv (1-p^{n-1}){\frac {B_{n}}{n}}{\pmod {p^{b}}}.} B n = − nζ (1 − n ) なので、これは次のように書くこともできます。
( 1 − p − u ) ζ ( u ) ≡ ( 1 − p − v ) ζ ( v ) ( モジュール p b ) 、 {\displaystyle \left(1-p^{-u}\right)\zeta (u)\equiv \left(1-p^{-v}\right)\zeta (v){\pmod {p^{b}}},} ここでu = 1 − m およびv = 1 − n であるため、u とv は 非正であり、p − 1 を法として 1 と合同ではありません。このことから、オイラー積の公式から1 − p − s を除いたリーマンゼータ関数は、p − 1 を法として特定のa ≢ 1 mod ( p − 1) と合同な奇数の負の整数上のp 進数で連続であり、したがってすべてのp 進整数に対して連続関数ζ p ( s ) に拡張できることがわかります。Z p 、 {\displaystyle \mathbb {Z} _{p},} p進 ゼータ関数 。
ラマヌジャンの合同式ラマヌジャン による以下の関係式は、ベルヌーイ数を計算するための、元の再帰的な定義よりも効率的な方法を提供する。
( m + 3 m ) B m = { m + 3 3 − ∑ j = 1 m 6 ( m + 3 m − 6 j ) B m − 6 j 、 もし m ≡ 0 ( モジュール 6 ) ; m + 3 3 − ∑ j = 1 m − 2 6 ( m + 3 m − 6 j ) B m − 6 j 、 もし m ≡ 2 ( モジュール 6 ) ; − m + 3 6 − ∑ j = 1 m − 4 6 ( m + 3 m − 6 j ) B m − 6 j 、 もし m ≡ 4 ( モジュール 6 ) 。 {\displaystyle {\binom {m+3}{m}}B_{m}={\begin{cases}{\frac {m+3}{3}}-\sum \limits _{j=1}^{\frac {m}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 0{\pmod {6}};\\{\frac {m+3}{3}}-\sum \limits _{j=1}^{\frac {m-2}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 2{\pmod {6}};\\-{\frac {m+3}{6}}-\sum \limits _{j=1}^{\frac {m-4}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 4{\pmod {6}}.\end{cases}}}
フォン・シュタウト=クラウゼン定理フォン シュタウト-クラウゼンの定理は、 1840 年にカール ゲオルク クリスチャン フォン シュタウト [ 50 ] とトーマス クラウゼン [ 51 ] によって独立して与えられました。この定理は、すべてのn > 0 について、
B 2 n + ∑ ( p − 1 ) ∣ 2 n 1 p {\displaystyle B_{2n}+\sum _{(p-1)\,\mid \,2n}{\frac {1}{p}}} は整数です。和は、p − 1が 2 n を割り切るすべての素数 p について取られます。
このことから、 B 2 n の分母は、p − 1 が 2 n を割り切るすべての素数p の積で与えられることがわかる。特に、これらの分母は平方因子を含まず 、6 で割り切れる。
奇数ベルヌーイ数はなぜ消滅するのか?合計
φ k ( n ) = ∑ 私 = 0 n 私 k − n k 2 {\displaystyle \varphi _{k}(n)=\sum _{i=0}^{n}i^{k}-{\frac {n^{k}}{2}}} 負のインデックスn の値に対して評価できます。そうすると、k が偶数の場合に奇関数で あることがわかり、これは和が奇インデックスの項のみを持つことを意味します。これとベルヌーイ和の公式から、m が偶数で2 k + 1 − m > 1の場合、 B 2 k + 1 − m は 0 となり、 B 1 の項は減算によって打ち消されることがわかります。フォン・シュタウト・クラウゼンの定理とウォルピツキーの表現を組み合わせると、この問題に対する組み合わせ論的な答えも得られます ( n > 1 の場合に有効)。
フォン・シュタウト=クラウゼンの定理から、奇数n > 1 の場合、数2 B n は整数であることがわかっている。問題の整数がゼロであることが事前にわかっている場合は、これは自明に思える。しかし、ヴォルピツキーの表現を適用すると、次のようになる。
2 B n = ∑ m = 0 n ( − 1 ) m 2 m + 1 m ! { n + 1 m + 1 } = 0 ( n > 1 奇妙だ ) {\displaystyle 2B_{n}=\sum _{m=0}^{n}(-1)^{m}{\frac {2}{m+1}}m!\left\{{n+1 \atop m+1}\right\}=0\quad (n>1{\text{ is odd}})} 整数の和 として表すのは自明ではない。ここで、奇数インデックスでベルヌーイ数が消滅することを説明する組み合わせ論的事実が明らかになる。S n , m を { 1, 2, ..., n }から{1, 2, ..., m } への全射写像の数とすると、S n , m = m ! { n m } となる。最後の等式は、次の場合にのみ成り立つ。
∑ 奇数 m = 1 n − 1 2 m 2 S n 、 m = ∑ 平 m = 2 n 2 m 2 S n 、 m ( n > 2 さえ ) 。 {\displaystyle \sum _{{\text{odd }}m=1}^{n-1}{\frac {2}{m^{2}}}S_{n,m}=\sum _{{\text{even }}m=2}^{n}{\frac {2}{m^{2}}}S_{n,m}\quad (n>2{\text{ is even}}).} この等式は帰納法で証明できます。この等式の最初の2つの例は次のとおりです。
n = 4: 2 + 8 = 7 + 3 、n = 6: 2 + 120 + 144 = 31 + 195 + 40 。このように、ベルヌーイ数は奇数インデックスでゼロになる。これは、ベルヌーイ数の中にいくつかの自明でない組み合わせ的恒等式が内包されているためである。
リーマン予想の再確認 ベルヌーイ数とリーマンゼータ関数の間の関係は、ベルヌーイ数のみを使用するリーマン予想 (RH)の別の定式化を提供するのに十分なほど強い。実際、マルセル・リースは 、RHが次の主張と同等であることを証明した。[ 52 ]
任意のε > 1/4 に対して、 定数C ε > 0 ( ε に依存する)が存在し、 x → ∞ の とき| R ( x ) | < C ε x ε と なる。ここで、R ( x )はリース関数 である。
R ( x ) = 2 ∑ k = 1 ∞ k k ¯ x k ( 2 π ) 2 k ( B 2 k 2 k ) = 2 ∑ k = 1 ∞ k k ¯ x k ( 2 π ) 2 k β 2 k 。 {\displaystyle R(x)=2\sum _{k=1}^{\infty }{\frac {k^{\overline {k}}x^{k}}{(2\pi )^{2k}\left({\frac {B_{2k}}{2k}}\right)}}=2\sum _{k=1}^{\infty }{\frac {k^{\overline {k}}x^{k}}{(2\pi )^{2k}\beta _{2k}}}.} n k は DE Knuth の表記法における上昇階乗 を表します。 β n = B n / n という数はゼータ関数の研究で頻繁に現れ、 β n はp − 1が n を割り切らない素数pに対して p 整数であるため。 β n は分割ベルヌーイ数 と呼ばれます。
一般化ベルヌーイ数一般化ベルヌーイ数は 、ベルヌーイ数と同様に定義されるある種の代数的数であり、 ベルヌーイ数がリーマンゼータ関数の特殊値と関連付けられているのと同様に、ディリクレL関数の特殊 値と関連付けられています 。
χを f を法とするディリクレ 指標とする。χに付随する一般化ベルヌーイ数は次のように定義される。
∑ 1 = 1 f χ ( 1 ) t e 1 t e f t − 1 = ∑ k = 0 ∞ B k 、 χ t k k ! 。 {\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 、 χ ) = − B k 、 χ k 、 {\displaystyle L(1-k,\chi )=-{\frac {B_{k,\chi }}{k}},} ここで、 L ( s , χ )はχ のディリクレL 関数である。[ 53 ]
注記 ↑ テキストの翻訳: [ 12 ] " ...このようにして、段階的に、より高いべき乗に到達し、わずかな努力で次の表を作成できます: べき乗の合計 ∫ n = ∑ k = 1 n k = 1 2 n 2 + 1 2 n {\displaystyle \textstyle \int n=\sum _{k=1}^{n}k={\frac {1}{2}}n^{2}+{\frac {1}{2}}n} ⋮ ∫ n 10 = ∑ k = 1 n k 10 = 1 11 n 11 + 1 2 n 10 + 5 6 n 9 − 1 n 7 + 1 n 5 − 1 2 n 3 + 5 66 n {\displaystyle \textstyle \int n^{10}=\sum _{k=1}^{n}k^{10}={\frac {1}{11}}n^{11}+{\frac {1}{2}}n^{10}+{\frac {5}{6}}n^{9}-1n^{7}+1n^{5}-{\frac {1}{2}}n^{3}+{\frac {5}{66}}n} 実際、もしあなたがそこで算術数列の法則を熱心に調べていれば、このような回りくどい計算をせずに同じことを続けることもできるでしょう。なぜなら、c {\displaystyle \textstyle c} は任意のべき乗の指数、すべての合計として扱われます。n c {\displaystyle \textstyle n^{c}} 生産されるか、 ∫ n c = ∑ k = 1 n k c = 1 c + 1 n c + 1 + 1 2 n c + c 2 A n c − 1 + c ( c − 1 ) ( c − 2 ) 2 ⋅ 3 ⋅ 4 B n c − 3 + c ( c − 1 ) ( c − 2 ) ( c − 3 ) ( c − 4 ) 2 ⋅ 3 ⋅ 4 ⋅ 5 ⋅ 6 C n c − 5 + c ( c − 1 ) ( c − 2 ) ( c − 3 ) ( c − 4 ) ( c − 5 ) ( c − 6 ) 2 ⋅ 3 ⋅ 4 ⋅ 5 ⋅ 6 ⋅ 7 ⋅ 8 D n c − 7 + ⋯ {\displaystyle \textstyle \int n^{c}=\sum _{k=1}^{n}k^{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\cdot 3\cdot 4}}Bn^{c-3}+{\frac {c(c-1)(c-2)(c-3)(c-4)}{2\cdot 3\cdot 4\cdot 5\cdot 6}}Cn^{c-5}+{\frac {c(c-1)(c-2)(c-3)(c-4)(c-5)(c-6)}{2\cdot 3\cdot 4\cdot 5\cdot 6\cdot 7\cdot 8}}Dn^{c-7}+\cdots } そして、その力の指数n {\displaystyle n} 2ずつ減り続け、n {\displaystyle n} またはn 2 {\displaystyle n^{2}} 大文字A 、 B 、 C 、 D 、 {\displaystyle \textstyle A,B,C,D,} など、最後の項の係数を順に表す。∫ n 2 、 ∫ n 4 、 ∫ n 6 、 ∫ n 8 {\displaystyle \textstyle \int n^{2},\int n^{4},\int n^{6},\int n^{8}} など、すなわち A = 1 6 、 B = − 1 30 、 C = 1 42 、 D = − 1 30 {\displaystyle \textstyle A={\frac {1}{6}},B=-{\frac {1}{30}},C={\frac {1}{42}},D=-{\frac {1}{30}}} 。 ↑ この公式はジョルジオ・ピエトロコラによって発見(あるいは再発見)されました。彼の証明はイタリア語で入手可能です(ピエトロコラ 2008 )。
参考文献 1 2 3 Weisstein, Eric W. 、「ベルヌーイ数」、MathWorld 1 2 セリン、ヘレイン 編 (1997)、『 非西洋文化における科学、技術、医学の歴史百科事典』 、シュプリンガー、p. 819 (p. 891)、 Bibcode : 2008ehst.book.....S 、 ISBN 0-7923-4066-3 ↑ スミス、デイビッド・ユージン;三上義雄(1914)『 日本数学史』 オープンコート出版、 108ページ 再版、ドーバー出版、2005年、ISBN 9780486434827 ↑ 北川智子 (2021-07-23)、「ベルヌーイ数の起源:18世紀初頭のバーゼルと江戸の数学」、 The Mathematical Intelligencer 、 44 : 46–56 、 doi : 10.1007/s00283-021-10072-y 、 ISSN 0343-6993 ↑ メナブレア、LF(1842)、 「チャールズ・バベッジが発明した解析機関の概略、翻訳者エイダ・オーガスタ・ラブレース伯爵夫人の回想録に関する注釈付き」 、 ジュネーブ大学図書館 、 82 、 注G参照 ↑ ドナルド・クヌース (2022)、最近のニュース (2022): 具体的な数学とベルヌーイ。しかし昨年、私はピーター・ルシュニーのベルヌーイ宣言をじっくりと読みました。そこで彼は、$B_1$の値が実際にはプラス1/2であるべき12以上の正当な理由を挙げています。彼は、20世紀初頭の数学者の中には、負の値を使うと数式が少し見栄えが良くなるという理由で、一方的に慣習を変更した者がいたと説明しています。彼らの善意から出たものの、結局は誤った選択が、私が1950年代に教えられたものにつながったのです。[…] 残念ながら、今では「マイナス1/2」の慣習を用いた書籍が何百冊も出版されています。さらに悪いことに、主要な数式処理ソフトウェアシステムのすべてに、この20世紀の異常な慣習が深く根付いています。しかし、ルシュニーは、私たち全員が間違っていたこと、そして状況がさらに悪化する前に正しい定義に戻す時が来たことを私に確信させてくれました。 ↑ ピーター・ルシュニー (2013)、『ベルヌーイ宣言』 ↑ Jacobi、CGJ ( 1834)、 「De usu Legitimo Formulas summatoriae Maclaurinianae」 、 Journal für die reine und angewandte Mathematik 、 12 : 263–272 ↑ スミス、デイヴィッド・ユージン (1929)。 「ジャック(I)ベルヌーイ:ベルヌーイ数について」 「 .数学のソースブック . ニューヨーク: マグロウヒル. pp. 85–90 . 1 2 ベルヌーイ、ジェイコブ(1713)。 Ars Conjectandi (ラテン語)。バーゼル:インペンシス・トゥルニシオルム、フラトラム。 pp. 97–98 . doi : 10.5479/sil.262971.39088000323931 。 ↑ ミラー、ジェフ(2017年6月23日) 「微積分記号の初期の使用例」 、 2026年1月3日 取得 ↑ ジョーダン(1950)p 233 ↑ アイルランドとローゼン(1990)p 229 ↑ ザールシュッツ、ルイ (1893)、 Vorlesungen über die Bernoullischen Zahlen、ihren Zusammenhang mit den Secanten-Coefficienten und ihre wichtigeren Anwendungen 、ベルリン: Julius Springer 。↑ 秋山茂樹。谷川良雄(2001年12月1日)。 「非正整数における複数のゼータ値」 (PDF) 。 ラマヌジャンジャーナル 。 5 (4): 327–351 . 土井 : 10.1023/A:1013981102941 。 2026 年 2 月 11 日 に取得 。 ↑ 金子正信(2000年12月12日)。 「ベルヌーイ数の秋山・谷川アルゴリズム」 (PDF) 。 整数シーケンスのジャーナル 。 3 (1) . 2026 年 2 月 11 日 に取得 。 ↑ 川崎菜穂、大野康夫 (2023年6月12日)。 「ベルヌーイ多項式の三角形アルゴリズム」 (PDF) . Integers . 23 . doi : 10.5281/zenodo.8028914 . 2026年 2月11日 取得 . ↑ Brent, Richard P.; Harvey, David (2013). "ベルヌーイ数、正接数、割線の高速計算" . Computational and Analytical Mathematics . 50 : 127– 142. arXiv : 1108.0286 . doi : 10.1007/978-1-4614-7621-4_8 . 2026年 2月11日 取得. 1 2 Knuth, DE ; Buckholtz, TJ (1967)、「接線数、オイラー数、ベルヌーイ数の計算」、 Mathematics of Computation 、 21 (100)、American Mathematical Society: 663–688 、 doi : 10.2307/2005010 、 JSTOR 2005010 ↑ Buhler, J.; Crandall, R.; Ernvall, R.; Metsankyla, T.; Shokrollahi, M. (2001), "Irregular Primes and Cyclotomic Invariants to 12 Million", Journal of Symbolic Computation , 31 ( 1–2 ): 89–96 , doi : 10.1006/jsco.1999.1011 ↑ Harvey, David (2010), "ベルヌーイ数を計算するためのマルチモジュラーアルゴリズム", Math. Comput. , 79 (272): 2361– 2370, arXiv : 0807.1347 , doi : 10.1090/S0025-5718-2010-02367-1 , S2CID 11329343 , Zbl 1215.11016 ↑ ケルナー、ベルント (2002)、 プログラム Calcbn – ベルヌーイ数を計算するプログラム 。↑ パヴリク、オレクサンドル(2008年4月29日) 「今日、ベルヌーイの記録を破った:解析エンジンからMathematicaへ」 、 Wolfram News 。↑ Guo, Victor JW; Zeng, Jiang (2005年8月30日)、「べき乗和に関するファウルハーバーの公式のq-類似」、 電子組合せ論ジャーナル 、 11 (2) R19、 arXiv : math/0501441 、 Bibcode : 2005math......1441G 、 doi : 10.37236/1876 、 S2CID 10467873 ↑ Comtet, L. (1974). 高度な組み合わせ論。有限および無限展開の技法 (改訂増補 版)。ドルトレヒト-ボストン:D. Reidel Publ. 1 2 ラデマッハー、H. (1973)、 解析的数論 、ニューヨーク市:シュプリンガー・フェルラーク 。↑ ブール、G.(1880) 『有限差分法の論文』 (第3 版)、ロンドン:マクミラン 。↑ グールド、ヘンリー W. (1972)、「ベルヌーイ数の明示的な公式」、 アメリカ数学月報 、 79 (1): 44–51 、 doi : 10.2307/2978125 、 JSTOR 2978125 ↑ アポストル、トム M. (2010)、 『解析的整数論入門』 、シュプリンガー・フェルラーク、 197 ページ ↑ ピエトロコラ、ジョルジオ (2017)、「 連続する整数のべき乗の和とパスカルの三角形から導出されるベルヌーイ数の計算のための多項式について」 、 Academia.edu 。 ↑ Woon, SC (1997)、「ベルヌーイ数を生成するためのツリー」、 Math. Mag. 、 70 (1): 51–56 、 doi : 10.2307/2691054 、 JSTOR 2691054 ↑ Stanley, Richard P. (2010), "A survey of alternating permutations", Combinatorics and graphs , Contemporary Mathematics, vol. 531, Providence, RI: American Mathematical Society, pp. 165–196 , arXiv : 0912.4240 , doi : 10.1090/conm/531/10466 , ISBN 978-0-8218-4865-4 MR 2757798、S2CID 14619581 ↑ Elkies, ND (2003), "On the sums Sum_(k=-infinity...infinity) (4k+1)^(-n)", Amer. Math. Monthly , 110 (7): 561– 573, arXiv : math.CA/0101168 , doi : 10.2307/3647742 , JSTOR 3647742 ↑ Euler, Leonhard (1735)、「De summis serierum reciprocarum」、 Opera Omnia 、I.14、E 41: 73–86 、 arXiv : math/0506415 、 Bibcode : 2005math....6415E ↑ Seidel, L. (1877)、「Über eine einfache Entstehungsweise der Bernoullischen Zahlen und einiger verwandten Reihen」、 Sitzungsber。ミュンシュ。アカド。 、 4 : 157~ 187 ↑ Dumont, D. (1981)、 「Matrices d'Euler-seidel」 、 Séminaire Lotharingien de Combinatoire 、 B05c ↑ Arnold, VI (1991), "関数特異点に関連するベルヌーイ・オイラーのアップダウン数、その組み合わせ論と算術", Duke Math. J. , 63 (2): 537– 555, doi : 10.1215/s0012-7094-91-06323-4 ↑ André, D. (1879)、「Développements de sec x et Tan x」、 CR Acad.科学。 、 88 : 965–967 ↑ André, D. ( 1881)、「Mémoire sur les permutations alternées」、 Journal de Mathématiques Pures et Appliquées 、 7 : 167–184 ↑ Kummer, EE (1850)、 「Allgemeiner Beweis des Fermat'schen Satzes, dass die Gleichung x λ + y λ = z λ durch ganze Zahlen unlösbar ist, für alle diejenigen Potenz-Exponenten λ, welche ungerade Primzahlen sind und in den Zählern der」 ersten (λ-3)/2 Bernoulli'schen Zahlen als Factoren nicht vorkommen」 、 J. Reine Angew。数学。 、 40 : 131~ 138 ↑ Kummer、EE (1851)、 「Über eine allgemeine Eigenschaft derrationen Entwicklungscoefficienten einer bestimmten Gattung Analytischer Functionen」 、 J. Reine Angew。数学。 、 1851 ( 41): 368–372 ↑ フォン・シュタウト、KG Ch. (1840)、「Beweis eines Lehrsatzes, die Bernoullischen Zahlen betreffend」、 Journal für die reine und angewandte Mathematik 、 21 : 372– 374 ↑ Clausen、Thomas (1840)、「Lehrsatz aus einer Abhandlung über die Bernoullischen Zahlen」、 アストロン。ナクル。 、 17 (22): 351–352 、 土井 : 10.1002/asna.18400172205 ↑ Riesz, M. (1916)、「Sur l'hypothèse de Riemann」、 Acta Mathematica 、 40 : 185–90 、 doi : 10.1007/BF02418544 ↑ Neukirch、Jürgen (1999)、 Algebraische Zahlentheorie 、 Grundlehren der mathematischen Wissenschaften 、vol. 322、ベルリン: Springer-Verlag 、 ISBN 978-3-540-65399-8 MR 1697859、Zbl 0956.11021 §VII.2.↑ Charollois, Pierre; Sczech, Robert (2016), "Elliptic Functions According to Eisenstein and Kronecker: An Update", EMS Newsletter , 2016– 9 (101): 8– 14, doi : 10.4171/NEWS/101/4 , ISSN 1027-488X , S2CID 54504376 1 2 伴内 健一;小林 伸一 (2010)、 「代数シータ関数とアイゼンシュタイン-クロネッカー数の p 進補間」 、 Duke Mathematical Journal 、 153 (2)、 arXiv : math/0610163 、 doi : 10.1215/00127094-2010-024 、 ISSN 0012-7094 、 S2CID 9262012 1 2 Malenfant, Jerome (2011), "Finite, closed-form expressions for the partition function and for Euler, Bernoulli, and Stirling numbers", arXiv : 1103.1585 [ math.NT ] ↑ Euler、E41、Inventio summae cuiusque seriei ex dato termino generali ↑ von Ettingshausen, A. (1827)、 Vorlesungen über die höhere Mathematik 、vol. 1、ウィーン: カール・ゲロルト ↑ Carlitz, L. (1968), "Bernoulli Numbers", Fibonacci Quarterly , 6 (3): 71–85 , doi : 10.1080/00150517.1968.12431229 ↑ Agoh, Takashi; Dilcher, Karl (2008), "ベルヌーイ数の相互関係", American Mathematical Monthly , 115 (3): 237– 244, doi : 10.1080/00029890.2008.11920520 , JSTOR 27642447 , S2CID 43614118
参考文献 Abramowitz, M.; Stegun, IA (1972)、「§23.1: ベルヌーイ多項式とオイラー多項式およびオイラー・マクローリンの公式」、数式、グラフ、数表を含む数学関数ハンドブック (第9版 )、ニューヨーク:ドーバー出版、pp. 804–806 。Arfken, George (1970), Mathematical methods for physicists (2nd ed.), Academic Press, Bibcode : 1970mmp..book.....A , ISBN 978-0120598519 Arlettaz、D. (1998)、「Die Bernoulli-Zahlen: eine Beziehung zwischen Topologie und Gruppentheorie」、Math。学期 、45 : 61– 75、土井 : 10.1007/s005910050037、S2CID 121753654 。Ayoub, A. (1981)、「オイラーとゼータ関数」、Amer. Math. Monthly 、74 (2): 1067–1086 、doi : 10.2307/2319041、JSTOR 2319041 。コンウェイ、ジョン ;ガイ、リチャード (1996)『数字の書 』シュプリンガー・フェルラーク 。ディルチャー、K.、スクラ、L.、スラヴツキー、I. Sh. (1991)、「ベルヌーイ数。文献目録(1713~1990)」、クイーンズ純粋応用数学論文集 (87)、オンタリオ州キングストン 。Dumont, D.; Viennot, G. (1980), "ジェノッキ数のザイデル生成の組み合わせ論的解釈", Ann. Discrete Math. , Annals of Discrete Mathematics, vol. 6, pp. 77– 87, doi : 10.1016/S0167-5060(08)70696-4 , ISBN 978-0-444-86048-4 。Entringer、RC (1966)、「オイラー数とベルヌーイ数の組み合わせ解釈」、Nieuw.アーチ。 V.ヴィスクンデ 、14 : 241–6 。Fee, G.; Plouffe, S. (2007)、「ベルヌーイ数の計算のための効率的なアルゴリズム」、arXiv : math/0702300 。Graham, R. ; Knuth, DE ; Patashnik, O. (1989), Concrete Mathematics (2nd ed.), Addison-Wesley, ISBN 0-201-55802-5 アイルランド、ケネス、ローゼン、マイケル (1990)、現代数論への古典的入門 (第 2 版)、シュプリンガー・フェルラーク、ISBN 0-387-97329-X ジョーダン、チャールズ (1950)、有限差分法 、ニューヨーク:チェルシー出版。 。金子正人(2000)「ベルヌーイ数に対する秋山・谷川アルゴリズム」『整数列ジャーナル 』 12巻 29ページ、Bibcode :2000JIntS...3...29K 。Knuth, DE (1993)、「ヨハン・ファウルハーバーとべき乗の和」、Mathematics of Computation 、61 (203)、American Mathematical Society: 277–294 、arXiv : math/9207222 、doi : 10.2307/2152953、JSTOR 2152953 Kouba, Omran (2016). "講義ノート、ベルヌーイ多項式とその応用". arXiv : 1309.7560v2 [ math.CA ]. ルシュニー、ピーター(2007)ベルヌーイ数の包含 。Luschny, Peter (2011年10月8日)、「TheLostBernoulliNumbers」、OeisWiki 、 2019年 5月11日 取得 。数学系譜プロジェクト 、ファーゴ:ノースダコタ州立大学数学科、日付不明、2019年5月10日にオリジナルからアーカイブ、2019年5月11日に取得 。ミルナー、ジョン・W. 、スタシェフ、ジェームズ・D. (1974)、「付録B:ベルヌーイ数」、特性クラス 、数学研究年報 、第76巻 、プリンストン大学出版局および東京大学出版会、pp. 281–287 、ISBN 978-0-691-08122-9 MR 0440554 。Pietrocola、Giorgio (2008 年 10 月 31 日)、「Esplorando un antico Sentiero: teoremi sulla somma di potenze di interi successivi (Corollario 2b)」、Maecla (イタリア語) 、2017 年4 月 8 日取得 。スラヴツキー、イリヤ・Sh. (1995)、「シュタウトとベルヌーイ数の算術的性質」、Historia Scientiarum 、2 : 69–74 。フォン・シュタウト、KG Ch. (1845)、「De numeris Bernoullianis, commentationem alteram」、エアランゲン 。孫志偉 (2005–2006)、「ベルヌーイ多項式とオイラー多項式に関するいくつかの興味深い結果」 、2001年10月31日にオリジナルからアーカイブ済み 。Woon, SC (1998)、「リーマンゼータ関数とベルヌーイ数の関係の一般化」、arXiv : math.NT/9812143 。Worpitzky, J. ( 1883)、「Studien über die Bernoullischen und Eulerschen Zahlen」、Journal für die reine und angewandte Mathematik 、94 : 203–232 。
外部リンク 「ベルヌーイ数」、数学百科事典 、EMS Press 、2001年 [1994年] プロジェクト・グーテンベルク からの最初の498個のベルヌーイ数 ベルヌーイ数を計算するためのマルチモジュラーアルゴリズム ベルヌーイ数のページ LiterateProgramsのベルヌーイ数プログラム P. ルシュニー著『不規則素数の計算』 P. ルシュニー著『ベルヌーイ数の計算と漸近解析』 ゴットフリート・ヘルムズ、「パスカル(二項)行列の文脈におけるベルヌーイ数 (PDF)」 、2022年10月9日にオリジナルからアーカイブ(PDF) ゴットフリート・ヘルムス、「パスカル/ベルヌーイ行列における同類べき乗の加算 (PDF)」 、2022年10月9日にオリジナルからアーカイブ(PDF) ゴットフリート・ヘルムス、「ベルヌーイ数とその関連数の和に関するいくつかの特別な性質 (PDF)」 、2022年10月9日にオリジナルからアーカイブ(PDF)