BellB
予備知識
- BellBは,ベル数またはベル多項式返す数学関数である.中でも,BellB[n,x]は

番目のベル多項式
を返し,BellB[n]は 
番目のベル数
を返す.ベル多項式は指数母関数
で決定することができる.ベル数は漸化式
を満足する.最初の数個のベル多項式
は
で,最初の数個のベル数
は
である. - ベル多項式は,指数多項式,あるいはより明示的に,「完全指数ベル多項式」と呼ばれることがあり,
と表示されることがある.ベル多項式は,数学者で数学の解説者でもあり,1934年にこの多項式についての論文を発表したEric Temple Bellに因んで名付けられた. - 多項式
は,
を
個の部分に分割する方法が
通りあるなら
であると解釈される.さらに,
を分割する方法が全部で
通りなら
である.例えば,元の数が
の集合
を
の部分に分割する方法は
通り
,
つの部分に分割する方法は
通り (
),
つの部分に分割する方法は
通り(
,
と
),
つの部分に分割する方法は
通り(
)で,
となる.
を分割する方法は全部で5通りなので,
である. - ベル多項式とベル数は,BellY関数の特殊ケース(
および
)である.
でStirlingS2が返す第2種スターリング数を表すとすると,
となる.
例題
すべて開く すべて閉じるスコープ (5)
BellB[1000, 0.5]BellB[1000, 0.500000000000000000000000]BellBは要素単位でリストに並列的な関数の適用を行う:
BellB[{1, 2, 3, 4, 5, 6}]BellBは,ベキ級数に適用することができる:
BellB[n, Sin[x] + O[x] ^ 3]TraditionalFormによる表示:
BellB[n]//TraditionalFormBellB[n, x]//TraditionalFormアプリケーション (4)
BellB数とその漸近線との対比:
Show[ListPlot[Table[(Log[BellB[n]]/n), {n, 50}]], Plot[Log[n] + (Log[Log[n]] + 1) ((1/Log[n]) - 1) + (1/2) ((Log[Log[n]]/Log[n]))^2, {n, 5, 50}]]最初の10個のcomplementary Bell numbersを計算する:
Table[BellB[n, -1], {n, 10}]Table[Underoverscript[∑, k = 0, n](-1)^kStirlingS2[n, k], {n, 10}]ヘッセンベルグ行列式によるベル数についての式の最初の数例を検証する:
Table[BellB[n + 1] == Det[SparseArray[{Band[{2, 1}] -> 1, Band[{1, 1}] -> 2, {j_, k_} /; j < k :> (-1)^j + kBinomial[k - 1, j - 1]}, {n, n}]], {n, 2, 9}]ベル数BellB[n]は,これから作られた2つのハンケル行列式がどちらもBarnesG[n+2]と等しい一意的な集合であるという事実で特徴付けられる.最初のいくつかの例についてこれを検証する:
Table[BarnesG[n + 2] == Det[HankelMatrix[BellB[Range[0, n]], BellB[Range[n, 2n]]]] == Det[HankelMatrix[BellB[Range[n + 1]], BellB[Range[n + 1, 2n + 1]]]], {n, 9}]特性と関係 (7)
BellBの指数型母関数:
ExponentialGeneratingFunction[BellB[n], n, x]Sum[BellB[n]x^n / n!, {n, 0, Infinity}]Sum[k^n / k!, {k, 0, Infinity}, Assumptions -> n∈Integers && n ≥ 0]Sum[StirlingS2[n, k], {k, 0, n}]PoissonDistributionの
番目のモーメントは平均
の
番目のBellB多項式で与えられる:
Moment[PoissonDistribution[μ], n]FullSimplifyを使ってBellBを含む式を簡約する:
FullSimplify[Mod[BellB[k + p] + BellB[k] - BellB[k + 1], p], p ∈ Primes]Table[Sum[(n!/Underoverscript[∏, k = 1, n]k!^Subscript[m, k]Subscript[m, k]!)Boole[Subsuperscript[∑, k = 0, n]k Subscript[m, k] == n∧Subsuperscript[∑, k = 1, n]Subscript[m, k] == Subscript[m, 0]], Evaluate[Array[{Subscript[m, #], 0, n}&, n + 1, 0, Sequence]]], {n, 1, 5}]IntegerPartitionsを使って指標についての条件を満足する項の総和を直接計算する:
Table[Sum[(Multinomial@@p) / Product[k!, {k, Length /@ Split[p]}], {p, IntegerPartitions[n]}], {n, 10}]BellBの結果と比較する:
Table[BellB[k], {k, 1, 10}]Table[BellY[ConstantArray[1, {n, 2}]], {n, 10}]Table[BellB[n], {n, 10}]Table[BellY[ConstantArray[{1, x}, n]], {n, 5}]Table[BellB[n, x], {n, 5}]FindSequenceFunctionはBellBを認識する:
Table[BellB[n], {n, 10}]FindSequenceFunction[%, n]考えられる問題 (1)
テクニカルノート
-
▪
- 組合せ関数
関連リンク
テキスト
Wolfram Research (2007), BellB, Wolfram言語関数, https://reference.wolfram.com/language/ref/BellB.html.
CMS
Wolfram Language. 2007. "BellB." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/BellB.html.
APA
Wolfram Language. (2007). BellB. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/BellB.html
BibTeX
@misc{reference.wolfram_2026_bellb, author="Wolfram Research", title="{BellB}", year="2007", howpublished="\url{https://reference.wolfram.com/language/ref/BellB.html}", note=[Accessed: 15-August-2026]}
BibLaTeX
@online{reference.wolfram_2026_bellb, organization={Wolfram Research}, title={BellB}, year={2007}, url={https://reference.wolfram.com/language/ref/BellB.html}, note=[Accessed: 15-August-2026]}