形式的な冪級数。係数は自然数でインデックスされたシーケンスに関する情報をエンコードします。
数学 において 、 生成関数は、 無限の数列を 形式的な冪級数 の 係数 として 表現したものです。生成関数は、形式的な級数に対する演算を含む何らかの式によって、級数としてではなく 閉じた形式 で表現されることがよくあります 。
生成関数には、 通常の生成関数 、 指数生成関数 、 ランバート級数 、 ベル級数 、 ディリクレ級数 など、さまざまな種類があります。すべてのシーケンスには、原則として各タイプの生成関数がありますが (ランバート級数とディリクレ級数では、インデックスが 0 ではなく 1 から始まる必要がある点が異なります)、処理のしやすさは大きく異なる場合があります。特定のコンテキストで最も役立つ特定の生成関数 (ある場合) は、シーケンスの性質と対処する問題の詳細によって異なります。
生成関数は、 項の級数がその項係数の列の生成元であると言えることから
、 生成級数と呼ばれることもあります [1] 。
歴史
生成関数は、 一般線形回帰問題を解くために、1730年に アブラハム・ド・モアブルによって初めて導入されました。 [2]
ジョージ・ポリアは 『数学ともっともらしい推論』 の中でこう書いています 。
「生成関数」という名前は ラプラス に由来します。しかし、名前を付けることなく、 オイラーはラプラスよりずっと前に生成関数という装置を使用していました[...]。彼はこの数学的なツールを組合せ解析と 数論の いくつかの問題に適用しました 。
意味
生成関数は、バッグに似た装置です。たくさんの小さな物をバラバラに持ち運ぶのは恥ずかしいかもしれませんが、代わりに、それらをすべてバッグに入れれば、持ち運ぶ物はバッグだけになります。
生成関数は、表示用に一連の数字を吊るす物干しロープです。
収束
通常の級数とは異なり、 形式的な 冪級数は 収束する 必要はありません 。実際、生成関数は実際には関数とは見なされず 、 「変数」は 不定値のままです。複数の不定値における形式的な冪級数に一般化して、無限の多次元の数値配列に関する情報をエンコードすることができます。したがって、生成関数は、 ドメインから コドメイン への マッピングという形式的な意味での関数ではありません 。
不定値x に関するこれらの式には、 算術演算、 x に関する微分 、および他の生成関数との合成 (つまり、代入) が含まれる場合があります。これらの演算は関数に対しても定義されているため、結果は xの関数のように見えます。実際、閉じた形式の式は、 xの (十分に小さい) 具体的な値で評価でき、その 級数展開 として形式級数を持つ関数 として解釈できることがよくあります。これが「生成関数」という呼称の理由です。ただし、このような解釈は可能である必要はありません。なぜなら、形式級数は、 x にゼロ以外の数値を代入したときに 収束級数 を与える必要がないためです 。
制限事項
x の関数として意味のあるすべての式が 、形式級数を指定する式としても意味があるわけではありません。たとえば、 x の負の累乗と分数の累乗は、対応する形式級数を持たない関数の例です。
種類
通常生成関数 (OGF)
生成関数という 用語 が修飾語なしで使用される場合、それは通常、通常の生成関数を意味すると解釈されます。 シーケンス a n の通常の生成関数は 次のとおりです。a
n が 離散ランダム変数 の 確率質量関数 である
場合 、その通常の生成関数は 確率生成関数 と呼ばれます。
グ
(
1つの
ん
;
x
)
=
∑
ん
=
0
∞
1つの
ん
x
ん
。
{\displaystyle G(a_{n};x)=\sum _{n=0}^{\infty }a_{n}x^{n}.}
指数生成関数 (EGF)
数列 a n の指数生成 関数 は
例えば
(
1つの
ん
;
x
)
=
∑
ん
=
0
∞
1つの
ん
x
ん
ん
!
。
{\displaystyle \operatorname {EG} (a_{n};x)=\sum _{n=0}^{\infty }a_{n}{\frac {x^{n}}{n!}}.}
一般に、ラベル付きオブジェクトを含む組み合わせ列挙 問題 では、指数生成関数は通常の生成関数よりも便利です。 [3]
指数生成関数のもう一つの利点は、線形 再帰関係を 微分方程式 の領域に移植する際に役立つことです 。たとえば、線形再帰関係 f n +2 = f n +1 + f n を満たす フィボナッチ数列 { f n } を考えます。対応する指数生成関数は次の形式になります。
EF
(
x
)
=
∑
ん
=
0
∞
ふ
ん
ん
!
x
ん
{\displaystyle \operatorname {EF} (x)=\sum _{n=0}^{\infty }{\frac {f_{n}}{n!}}x^{n}}
そして、その導関数は、上記の再帰関係と直接類似して、微分方程式 EF″( x ) = EF ′ ( x ) + EF( x ) を満たすことが容易に示されます。この見方では、階乗項 n ! は、 x n に作用する微分演算子を正規化するための単なるカウンター項です 。
ポアソン生成関数
数列 a n のポアソン生成 関数 は
PG
(
a
n
;
x
)
=
∑
n
=
0
∞
a
n
e
−
x
x
n
n
!
=
e
−
x
EG
(
a
n
;
x
)
.
{\displaystyle \operatorname {PG} (a_{n};x)=\sum _{n=0}^{\infty }a_{n}e^{-x}{\frac {x^{n}}{n!}}=e^{-x}\,\operatorname {EG} (a_{n};x).}
ランバートシリーズ
数列 a n のランバート級数は 、 ランバート級数ではインデックス n は 0 ではなく 1 から始まることに注意してください
。そうしないと、最初の項が未定義になります。
LG
(
a
n
;
x
)
=
∑
n
=
1
∞
a
n
x
n
1
−
x
n
.
{\displaystyle \operatorname {LG} (a_{n};x)=\sum _{n=1}^{\infty }a_{n}{\frac {x^{n}}{1-x^{n}}}.}
整数n ≥ 1 のべき 級数展開におけるランバート級数の係数は、
約数の和 によって関連付けられる。 本稿で は、 数論 における特殊 算術関数 に関連する、より古典的な、あるいは少なくともよく知られた例をいくつか提供する 。本稿で示されていないランバート級数の恒等式の例として、 | x |、| xq | < 1 に対して、次が成り立つことを示すことができる [4]。
b
n
:=
[
x
n
]
LG
(
a
n
;
x
)
{\displaystyle b_{n}:=[x^{n}]\operatorname {LG} (a_{n};x)}
b
n
=
∑
d
|
n
a
d
.
{\displaystyle b_{n}=\sum _{d|n}a_{d}.}
∑
n
=
1
∞
q
n
x
n
1
−
x
n
=
∑
n
=
1
∞
q
n
x
n
2
1
−
q
x
n
+
∑
n
=
1
∞
q
n
x
n
(
n
+
1
)
1
−
x
n
,
{\displaystyle \sum _{n=1}^{\infty }{\frac {q^{n}x^{n}}{1-x^{n}}}=\sum _{n=1}^{\infty }{\frac {q^{n}x^{n^{2}}}{1-qx^{n}}}+\sum _{n=1}^{\infty }{\frac {q^{n}x^{n(n+1)}}{1-x^{n}}},}
ここで、除数関数 の生成関数の特別な場合の恒等式は 、 d ( n ) ≡ σ 0 ( n ) で、次のように与えられる。
∑
n
=
1
∞
x
n
1
−
x
n
=
∑
n
=
1
∞
x
n
2
(
1
+
x
n
)
1
−
x
n
.
{\displaystyle \sum _{n=1}^{\infty }{\frac {x^{n}}{1-x^{n}}}=\sum _{n=1}^{\infty }{\frac {x^{n^{2}}\left(1+x^{n}\right)}{1-x^{n}}}.}
ベルシリーズ
数列 a n のベル 級数は不定数 x と素数 p の両方を用いた表現であり 、次のように与えられる: [5]
BG
p
(
a
n
;
x
)
=
∑
n
=
0
∞
a
p
n
x
n
.
{\displaystyle \operatorname {BG} _{p}(a_{n};x)=\sum _{n=0}^{\infty }a_{p^{n}}x^{n}.}
ディリクレ級数生成関数 (DGF)
形式的ディリクレ級数は 厳密には形式的冪級数ではないが、しばしば生成関数として分類される。 数列 a n のディリクレ級数生成関数は 次のようになる: [6]
DG
(
a
n
;
s
)
=
∑
n
=
1
∞
a
n
n
s
.
{\displaystyle \operatorname {DG} (a_{n};s)=\sum _{n=1}^{\infty }{\frac {a_{n}}{n^{s}}}.}
ディリクレ級数生成関数は、 n が 乗法関数 である場合に特に有用であり 、その場合には 関数のベル級数に関して
オイラー積 表現 [7]を持つ。
DG
(
a
n
;
s
)
=
∏
p
BG
p
(
a
n
;
p
−
s
)
.
{\displaystyle \operatorname {DG} (a_{n};s)=\prod _{p}\operatorname {BG} _{p}(a_{n};p^{-s})\,.}
もし n が ディリクレ指標 ならば、そのディリクレ級数生成関数は ディリクレ L 級数 と呼ばれる。また、上記の ランベルト級数 展開の係数のペア とそのDGFの間には関係がある。すなわち、次の式が成り立つことを証明できる。
もし成り立つならば、そしてその場合に限る。
ここで ζ ( s )は リーマンゼータ関数 である 。 [8]
[
x
n
]
LG
(
a
n
;
x
)
=
b
n
{\displaystyle [x^{n}]\operatorname {LG} (a_{n};x)=b_{n}}
DG
(
a
n
;
s
)
ζ
(
s
)
=
DG
(
b
n
;
s
)
,
{\displaystyle \operatorname {DG} (a_{n};s)\zeta (s)=\operatorname {DG} (b_{n};s),}
ディリクレ級数 生成関数 (DGF) によって生成される シーケンス a k は、次の通常の生成関数を持ちます。
DG
(
a
k
;
s
)
=
ζ
(
s
)
m
{\displaystyle \operatorname {DG} (a_{k};s)=\zeta (s)^{m}}
∑
k
=
1
k
=
n
a
k
x
k
=
x
+
(
m
1
)
∑
2
≤
a
≤
n
x
a
+
(
m
2
)
∑
a
=
2
∞
∑
b
=
2
∞
a
b
≤
n
x
a
b
+
(
m
3
)
∑
a
=
2
∞
∑
c
=
2
∞
∑
b
=
2
∞
a
b
c
≤
n
x
a
b
c
+
(
m
4
)
∑
a
=
2
∞
∑
b
=
2
∞
∑
c
=
2
∞
∑
d
=
2
∞
a
b
c
d
≤
n
x
a
b
c
d
+
⋯
{\displaystyle \sum _{k=1}^{k=n}a_{k}x^{k}=x+{\binom {m}{1}}\sum _{2\leq a\leq n}x^{a}+{\binom {m}{2}}{\underset {ab\leq n}{\sum _{a=2}^{\infty }\sum _{b=2}^{\infty }}}x^{ab}+{\binom {m}{3}}{\underset {abc\leq n}{\sum _{a=2}^{\infty }\sum _{c=2}^{\infty }\sum _{b=2}^{\infty }}}x^{abc}+{\binom {m}{4}}{\underset {abcd\leq n}{\sum _{a=2}^{\infty }\sum _{b=2}^{\infty }\sum _{c=2}^{\infty }\sum _{d=2}^{\infty }}}x^{abcd}+\cdots }
多項式列生成関数
生成関数の考え方は、他のオブジェクトのシーケンスに拡張できます。したがって、たとえば、 二項式 の多項式シーケンスは次のように生成されます。
ここで、 p n ( x ) は多項式のシーケンス、 f ( t ) は特定の形式の関数です。 シェファー シーケンスも 同様の方法で生成されます。詳細については、
メインの記事 「一般化された Appell 多項式」を参照してください。
e
x
f
(
t
)
=
∑
n
=
0
∞
p
n
(
x
)
n
!
t
n
{\displaystyle e^{xf(t)}=\sum _{n=0}^{\infty }{\frac {p_{n}(x)}{n!}}t^{n}}
より複雑な生成関数によって生成される多項式シーケンス の例には 次のものがあります。
その他の生成関数
より複雑な生成関数によって生成される他のシーケンスには次のものがあります。
二重指数関数生成関数。例: Aitken の配列: 数の三角形
生成関数と対角生成関数のアダマール積、およびそれに対応する 積分変換
畳み込み多項式
クヌースの論文「 畳み込み多項式 」 [9]は、 F (0)=1 となるべき級数展開を持つ
解析関数 F の形式の特殊生成関数によって
畳み込み多項式 列の一般化されたクラスを定義している 。
F
(
z
)
x
=
exp
(
x
log
F
(
z
)
)
=
∑
n
=
0
∞
f
n
(
x
)
z
n
,
{\displaystyle F(z)^{x}=\exp {\bigl (}x\log F(z){\bigr )}=\sum _{n=0}^{\infty }f_{n}(x)z^{n},}
deg f n ≤ n であり、すべての x 、 y およびすべての n ≥ 0 に対して
次の畳み込み条件が成り立つ場合、 多項式族 f 0 、 f 1 、 f 2 、...は 畳み込み族を 形成すると言います。
f
n
(
x
+
y
)
=
f
n
(
x
)
f
0
(
y
)
+
f
n
−
1
(
x
)
f
1
(
y
)
+
⋯
+
f
1
(
x
)
f
n
−
1
(
y
)
+
f
0
(
x
)
f
n
(
y
)
.
{\displaystyle f_{n}(x+y)=f_{n}(x)f_{0}(y)+f_{n-1}(x)f_{1}(y)+\cdots +f_{1}(x)f_{n-1}(y)+f_{0}(x)f_{n}(y).}
非同一ゼロ畳み込み族の場合、この定義は、シーケンスが上記で示した最初の形式の通常の生成関数を持つことを要求することと同等であることがわかります。
上記の表記法で定義された畳み込み多項式のシーケンスには、次の特性があります。
数列 n ! · f n ( x )は 二項式型 である
この数列の特別な値には、 f n (1) = [ z n ] F ( z ) と f n (0) = δ n ,0 があり、
任意の(固定された) に対して 、これらの多項式は次のような畳み込み公式を満たす。
x
,
y
,
t
∈
C
{\displaystyle x,y,t\in \mathbb {C} }
f
n
(
x
+
y
)
=
∑
k
=
0
n
f
k
(
x
)
f
n
−
k
(
y
)
f
n
(
2
x
)
=
∑
k
=
0
n
f
k
(
x
)
f
n
−
k
(
x
)
x
n
f
n
(
x
+
y
)
=
(
x
+
y
)
∑
k
=
0
n
k
f
k
(
x
)
f
n
−
k
(
y
)
(
x
+
y
)
f
n
(
x
+
y
+
t
n
)
x
+
y
+
t
n
=
∑
k
=
0
n
x
f
k
(
x
+
t
k
)
x
+
t
k
y
f
n
−
k
(
y
+
t
(
n
−
k
)
)
y
+
t
(
n
−
k
)
.
{\displaystyle {\begin{aligned}f_{n}(x+y)&=\sum _{k=0}^{n}f_{k}(x)f_{n-k}(y)\\f_{n}(2x)&=\sum _{k=0}^{n}f_{k}(x)f_{n-k}(x)\\xnf_{n}(x+y)&=(x+y)\sum _{k=0}^{n}kf_{k}(x)f_{n-k}(y)\\{\frac {(x+y)f_{n}(x+y+tn)}{x+y+tn}}&=\sum _{k=0}^{n}{\frac {xf_{k}(x+tk)}{x+tk}}{\frac {yf_{n-k}(y+t(n-k))}{y+t(n-k)}}.\end{aligned}}}
固定された非ゼロパラメータ に対して 、これらの畳み込み多項式列に対する生成関数は次のように修正される。
ここで 𝓕 t ( z )は、 𝓕 t ( z ) = F ( x 𝓕 t ( z ) t ) という形式の 関数方程式 によって暗黙的に定義される 。さらに、行列法(参考文献にあるように)を使用して、2つの畳み込み多項式列 ⟨ f n ( x ) ⟩ と ⟨ g n ( x ) ⟩ が与えられ、それぞれ対応する生成関数 F ( z ) x と G ( z ) x が ある場合、任意のt に対して次の 恒等式が成り立つことを証明できる。
t
∈
C
{\displaystyle t\in \mathbb {C} }
z
F
n
(
x
+
t
n
)
(
x
+
t
n
)
=
[
z
n
]
F
t
(
z
)
x
,
{\displaystyle {\frac {zF_{n}(x+tn)}{(x+tn)}}=\left[z^{n}\right]{\mathcal {F}}_{t}(z)^{x},}
[
z
n
]
(
G
(
z
)
F
(
z
G
(
z
)
t
)
)
x
=
∑
k
=
0
n
F
k
(
x
)
G
n
−
k
(
x
+
t
k
)
.
{\displaystyle \left[z^{n}\right]\left(G(z)F\left(zG(z)^{t}\right)\right)^{x}=\sum _{k=0}^{n}F_{k}(x)G_{n-k}(x+tk).}
畳み込み多項式シーケンスの例には、 二項べき級数 、 𝓑 t ( z ) = 1 + z 𝓑 t ( z ) t 、いわゆる ツリー多項式 、 ベル数 、 B ( n ) 、 ラゲール多項式 、 スターリング畳み込み多項式 などがあります。
通常の生成関数
単純なシーケンスの例
多項式は、通常の生成関数の特殊なケースであり、有限シーケンス、または同等に、特定のポイントの後に消えるシーケンスに対応します。これらは、多くの有限シーケンスが ポアンカレ多項式 などの生成関数として有用に解釈できるという点で重要です。
基本的な生成関数は定数列 1, 1, 1, 1, 1, 1, 1, 1, 1, ... の関数であり、その通常の生成関数は 等比級数である。
∑
n
=
0
∞
x
n
=
1
1
−
x
.
{\displaystyle \sum _{n=0}^{\infty }x^{n}={\frac {1}{1-x}}.}
左辺は 右辺の マクローリン級数展開です。あるいは、左辺のべき級数に 1 − x を掛けて、結果が定数べき級数 1 であることを確認することで等式を正当化できます (言い換えると、 x 0 の係数を除くすべての係数が 0 に等しい)。さらに、この特性を持つ他のべき級数は存在しません。したがって、左辺は、 べき級数の環における
1 − x の 逆乗を表します。
他の数列の通常の生成関数の式は、この式から簡単に導き出せます。例えば、 x → ax と代入すると、定数 aに対して 等比数列 1、 a 、 a 2 、 a 3 、... の生成関数が得られます 。
∑
n
=
0
∞
(
a
x
)
n
=
1
1
−
a
x
.
{\displaystyle \sum _{n=0}^{\infty }(ax)^{n}={\frac {1}{1-ax}}.}
(この等式は、左辺が右辺のマクローリン級数展開であるという事実からも直接導かれる。)特に、
∑
n
=
0
∞
(
−
1
)
n
x
n
=
1
1
+
x
.
{\displaystyle \sum _{n=0}^{\infty }(-1)^{n}x^{n}={\frac {1}{1+x}}.}
x を x のべき乗に 置き換えることで、シーケンスに規則的なギャップを導入することもできます 。たとえば、シーケンス 1, 0, 1, 0, 1, 0, 1, 0, ... ( x 、 x 3 、 x 5 、 ... をスキップします ) の場合、生成関数は次のようになります。
∑
n
=
0
∞
x
2
n
=
1
1
−
x
2
.
{\displaystyle \sum _{n=0}^{\infty }x^{2n}={\frac {1}{1-x^{2}}}.}
初期生成関数を2乗するか、両辺のx に関する微分を求め、実行変数 nを n + 1 に 変更すると、係数が 1、2、3、4、5、...という 数列を形成することがわかります 。したがって、
∑
n
=
0
∞
(
n
+
1
)
x
n
=
1
(
1
−
x
)
2
,
{\displaystyle \sum _{n=0}^{\infty }(n+1)x^{n}={\frac {1}{(1-x)^{2}}},}
そして、3乗は 1、3、6、10、15、21、… の三角数 を係数とし、その項 n は 二項係数 ( 2 + 2 2 ) 、 となることによって
∑
n
=
0
∞
(
n
+
2
2
)
x
n
=
1
(
1
−
x
)
3
.
{\displaystyle \sum _{n=0}^{\infty }{\binom {n+2}{2}}x^{n}={\frac {1}{(1-x)^{3}}}.}
より一般的には、任意の非負整数 k と非ゼロ実数値 a に対して、
∑
n
=
0
∞
a
n
(
n
+
k
k
)
x
n
=
1
(
1
−
a
x
)
k
+
1
.
{\displaystyle \sum _{n=0}^{\infty }a^{n}{\binom {n+k}{k}}x^{n}={\frac {1}{(1-ax)^{k+1}}}\,.}
以来
2
(
n
+
2
2
)
−
3
(
n
+
1
1
)
+
(
n
0
)
=
2
(
n
+
1
)
(
n
+
2
)
2
−
3
(
n
+
1
)
+
1
=
n
2
,
{\displaystyle 2{\binom {n+2}{2}}-3{\binom {n+1}{1}}+{\binom {n}{0}}=2{\frac {(n+1)(n+2)}{2}}-3(n+1)+1=n^{2},}
二項係数生成列の線形結合によって、
平方 数列 0、1、4、9、16、... の通常の生成関数を見つけることができます。
G
(
n
2
;
x
)
=
∑
n
=
0
∞
n
2
x
n
=
2
(
1
−
x
)
3
−
3
(
1
−
x
)
2
+
1
1
−
x
=
x
(
x
+
1
)
(
1
−
x
)
3
.
{\displaystyle G(n^{2};x)=\sum _{n=0}^{\infty }n^{2}x^{n}={\frac {2}{(1-x)^{3}}}-{\frac {3}{(1-x)^{2}}}+{\frac {1}{1-x}}={\frac {x(x+1)}{(1-x)^{3}}}.}
交互に拡張して、同じ平方数のシーケンスを次の形式の
等比級数 の導関数の合計として生成することもできます。
G
(
n
2
;
x
)
=
∑
n
=
0
∞
n
2
x
n
=
∑
n
=
0
∞
n
(
n
−
1
)
x
n
+
∑
n
=
0
∞
n
x
n
=
x
2
D
2
[
1
1
−
x
]
+
x
D
[
1
1
−
x
]
=
2
x
2
(
1
−
x
)
3
+
x
(
1
−
x
)
2
=
x
(
x
+
1
)
(
1
−
x
)
3
.
{\displaystyle {\begin{aligned}G(n^{2};x)&=\sum _{n=0}^{\infty }n^{2}x^{n}\\[4px]&=\sum _{n=0}^{\infty }n(n-1)x^{n}+\sum _{n=0}^{\infty }nx^{n}\\[4px]&=x^{2}D^{2}\left[{\frac {1}{1-x}}\right]+xD\left[{\frac {1}{1-x}}\right]\\[4px]&={\frac {2x^{2}}{(1-x)^{3}}}+{\frac {x}{(1-x)^{2}}}={\frac {x(x+1)}{(1-x)^{3}}}.\end{aligned}}}
帰納法により、正の整数m ≥ 1 について も同様に証明できる [10] [11]
n
m
=
∑
j
=
0
m
{
m
j
}
n
!
(
n
−
j
)
!
,
{\displaystyle n^{m}=\sum _{j=0}^{m}{\begin{Bmatrix}m\\j\end{Bmatrix}}{\frac {n!}{(n-j)!}},}
どこ { nk } は 第二種スターリング数 を表し 、生成関数は
∑
n
=
0
∞
n
!
(
n
−
j
)
!
z
n
=
j
!
⋅
z
j
(
1
−
z
)
j
+
1
,
{\displaystyle \sum _{n=0}^{\infty }{\frac {n!}{(n-j)!}}\,z^{n}={\frac {j!\cdot z^{j}}{(1-z)^{j+1}}},}
すると、上記の平方根の場合の結果を一般化した、
積分 m乗上の類似の生成関数を形成することができる。特に、次のように書くことができる。
z
k
(
1
−
z
)
k
+
1
=
∑
i
=
0
k
(
k
i
)
(
−
1
)
k
−
i
(
1
−
z
)
i
+
1
,
{\displaystyle {\frac {z^{k}}{(1-z)^{k+1}}}=\sum _{i=0}^{k}{\binom {k}{i}}{\frac {(-1)^{k-i}}{(1-z)^{i+1}}},}
スターリング数 を含むよく知られた有限和恒等式を適用して 、 [12]
∑
n
=
0
∞
n
m
z
n
=
∑
j
=
0
m
{
m
+
1
j
+
1
}
(
−
1
)
m
−
j
j
!
(
1
−
z
)
j
+
1
.
{\displaystyle \sum _{n=0}^{\infty }n^{m}z^{n}=\sum _{j=0}^{m}{\begin{Bmatrix}m+1\\j+1\end{Bmatrix}}{\frac {(-1)^{m-j}j!}{(1-z)^{j+1}}}.}
有理関数
数列の通常の生成関数は、 数列が定数係数 の線形再帰数列である場合に限り、 有理関数 (2 つの有限次多項式の比)として表現できます。これは上記の例を一般化したものです。逆に、多項式の分数によって生成されるすべての数列は、定数係数の線形再帰を満たします。これらの係数は、分数分母多項式の係数と同一です(したがって、直接読み取ることができます)。この観察は、定数係数の線形 差分方程式 によって定義される数列の生成関数を解くことが容易であり、したがって、これらの生成関数の係数の明示的な閉じた形式の式を解くことが容易であることを示しています。ここでの典型的な例は、生成関数の手法を使用して
フィボナッチ数列 の ビネの公式 を導出することです。
また、有理数生成関数のクラスは、次の形式の 準多項式 列を列挙する生成関数と正確に対応していることにも気づく [13]。
f
n
=
p
1
(
n
)
ρ
1
n
+
⋯
+
p
ℓ
(
n
)
ρ
ℓ
n
,
{\displaystyle f_{n}=p_{1}(n)\rho _{1}^{n}+\cdots +p_{\ell }(n)\rho _{\ell }^{n},}
ここで、逆根 、は固定されたスカラーであり、 p i ( n ) はすべての 1 ≤ i ≤ ℓ に対してn の多項式です 。
ρ
i
∈
C
{\displaystyle \rho _{i}\in \mathbb {C} }
一般に、有理関数の アダマール積は 有理関数を生成する。同様に、
F
(
s
,
t
)
:=
∑
m
,
n
≥
0
f
(
m
,
n
)
w
m
z
n
{\displaystyle F(s,t):=\sum _{m,n\geq 0}f(m,n)w^{m}z^{n}}
が2変量有理数生成関数である場合、対応する対 角生成関数 は、
diag
(
F
)
:=
∑
n
=
0
∞
f
(
n
,
n
)
z
n
,
{\displaystyle \operatorname {diag} (F):=\sum _{n=0}^{\infty }f(n,n)z^{n},}
は代数的で ある 。例えば、 [14]
F
(
s
,
t
)
:=
∑
i
,
j
≥
0
(
i
+
j
i
)
s
i
t
j
=
1
1
−
s
−
t
,
{\displaystyle F(s,t):=\sum _{i,j\geq 0}{\binom {i+j}{i}}s^{i}t^{j}={\frac {1}{1-s-t}},}
この生成関数の対角係数生成関数は、よく知られているOGF式で与えられる。
diag
(
F
)
=
∑
n
=
0
∞
(
2
n
n
)
z
n
=
1
1
−
4
z
.
{\displaystyle \operatorname {diag} (F)=\sum _{n=0}^{\infty }{\binom {2n}{n}}z^{n}={\frac {1}{\sqrt {1-4z}}}.}
この結果は、コーシーの積分公式 や 輪郭積分 、複素 留数 の取得、 2 つの変数の 形式的なべき級数 の直接操作など、さまざまな方法で計算されます 。
生成関数に対する演算
掛け算は畳み込みをもたらす
通常の生成関数の乗算は、シーケンスの離散 畳み込み ( コーシー積 )を生成します。たとえば、通常の生成関数 G ( a n ; x )
を持つシーケンスの累積和のシーケンス(やや一般的な オイラー・マクローリンの公式 と比較してください)
には、 生成関数があります
。
これは 、
(
a
0
,
a
0
+
a
1
,
a
0
+
a
1
+
a
2
,
…
)
{\displaystyle (a_{0},a_{0}+a_{1},a_{0}+a_{1}+a_{2},\ldots )}
G
(
a
n
;
x
)
⋅
1
1
−
x
{\displaystyle G(a_{n};x)\cdot {\frac {1}{1-x}}}
1 / 1 − x は、 シーケンス (1, 1, ...) の通常の生成関数です。生成関数の畳み込みと解釈による問題解決のさらなる例については、この記事の以下のアプリケーションセクションの畳み込みに関するセクションも参照してください。
シーケンスインデックスのシフト
整数 m ≥ 1に対して、 ⟨ g n − m ⟩ と ⟨ g n + m ⟩ のシフトされたシーケンスバリアントを列挙する修正生成関数に対して 、それぞれ次の2つの類似した恒等式が存在します。
z
m
G
(
z
)
=
∑
n
=
m
∞
g
n
−
m
z
n
G
(
z
)
−
g
0
−
g
1
z
−
⋯
−
g
m
−
1
z
m
−
1
z
m
=
∑
n
=
0
∞
g
n
+
m
z
n
.
{\displaystyle {\begin{aligned}&z^{m}G(z)=\sum _{n=m}^{\infty }g_{n-m}z^{n}\\[4px]&{\frac {G(z)-g_{0}-g_{1}z-\cdots -g_{m-1}z^{m-1}}{z^{m}}}=\sum _{n=0}^{\infty }g_{n+m}z^{n}.\end{aligned}}}
生成関数の微分と積分
生成関数の一次導関数とその積分には、それぞれ次のべき級数展開があります。
G
′
(
z
)
=
∑
n
=
0
∞
(
n
+
1
)
g
n
+
1
z
n
z
⋅
G
′
(
z
)
=
∑
n
=
0
∞
n
g
n
z
n
∫
0
z
G
(
t
)
d
t
=
∑
n
=
1
∞
g
n
−
1
n
z
n
.
{\displaystyle {\begin{aligned}G'(z)&=\sum _{n=0}^{\infty }(n+1)g_{n+1}z^{n}\\[4px]z\cdot G'(z)&=\sum _{n=0}^{\infty }ng_{n}z^{n}\\[4px]\int _{0}^{z}G(t)\,dt&=\sum _{n=1}^{\infty }{\frac {g_{n-1}}{n}}z^{n}.\end{aligned}}}
2 番目の恒等式の微分乗算演算を k 回繰り返して、シーケンスを n k 倍にすることができますが、そのためには微分と乗算を交互に行う必要があります。代わりに k 回の微分を順番に実行すると、 k 番目の 下降階乗 を乗算することになります 。
z
k
G
(
k
)
(
z
)
=
∑
n
=
0
∞
n
k
_
g
n
z
n
=
∑
n
=
0
∞
n
(
n
−
1
)
⋯
(
n
−
k
+
1
)
g
n
z
n
for all
k
∈
N
.
{\displaystyle z^{k}G^{(k)}(z)=\sum _{n=0}^{\infty }n^{\underline {k}}g_{n}z^{n}=\sum _{n=0}^{\infty }n(n-1)\dotsb (n-k+1)g_{n}z^{n}\quad {\text{for all }}k\in \mathbb {N} .}
第二種スターリング数を 使用すると 、次のように を掛ける別の式に変換できます( 生成関数変換 に関するメイン記事を参照 )。
n
k
{\displaystyle n^{k}}
∑
j
=
0
k
{
k
j
}
z
j
F
(
j
)
(
z
)
=
∑
n
=
0
∞
n
k
f
n
z
n
for all
k
∈
N
.
{\displaystyle \sum _{j=0}^{k}{\begin{Bmatrix}k\\j\end{Bmatrix}}z^{j}F^{(j)}(z)=\sum _{n=0}^{\infty }n^{k}f_{n}z^{n}\quad {\text{for all }}k\in \mathbb {N} .}
繰り返し積分の操作に対応するこの数列のべき乗の公式の負の順序の反転は、 ゼータ級数変換と、 生成関数の 微分ベースの変換として定義されるその一般化によって定義されます。 または、代わりに項ごとに、および数列生成関数に対して 積分変換を実行することによって定義されます。数列生成関数に対して 分数積分を 実行する関連操作については、 ここで 説明します 。
数列の等差数列を列挙する
このセクションでは、通常の生成関数 F ( z ) が与えられたときに、シーケンス { f an + b } を列挙する生成関数の公式を示します。ここで、 a ≥ 2 、 0 ≤ b < a 、 a と b は整数です ( 変換に関するメインの記事 を参照)。 a = 2 の場合、これは単に関数を 偶数部分と奇数部分 (つまり、偶数と奇数の累乗) に分解するだけです。
∑
n
=
0
∞
f
2
n
z
2
n
=
F
(
z
)
+
F
(
−
z
)
2
∑
n
=
0
∞
f
2
n
+
1
z
2
n
+
1
=
F
(
z
)
−
F
(
−
z
)
2
.
{\displaystyle {\begin{aligned}\sum _{n=0}^{\infty }f_{2n}z^{2n}&={\frac {F(z)+F(-z)}{2}}\\[4px]\sum _{n=0}^{\infty }f_{2n+1}z^{2n+1}&={\frac {F(z)-F(-z)}{2}}.\end{aligned}}}
より一般的には、 a ≥ 3 かつ ω a = exp と仮定する。 2πi / 1つの は 1 のa 番目 の原始根 を表します 。離散フーリエ変換 の応用として 、次の式が得られます [15]
∑
n
=
0
∞
f
a
n
+
b
z
a
n
+
b
=
1
a
∑
m
=
0
a
−
1
ω
a
−
m
b
F
(
ω
a
m
z
)
.
{\displaystyle \sum _{n=0}^{\infty }f_{an+b}z^{an+b}={\frac {1}{a}}\sum _{m=0}^{a-1}\omega _{a}^{-mb}F\left(\omega _{a}^{m}z\right).}
整数 m ≥ 1 の場合、 逆の床付き等差数列(各係数を m 回繰り返す )を提供する別の便利な式は、次の恒等式によって生成される [16]。
∑
n
=
0
∞
f
⌊
n
m
⌋
z
n
=
1
−
z
m
1
−
z
F
(
z
m
)
=
(
1
+
z
+
⋯
+
z
m
−
2
+
z
m
−
1
)
F
(
z
m
)
.
{\displaystyle \sum _{n=0}^{\infty }f_{\left\lfloor {\frac {n}{m}}\right\rfloor }z^{n}={\frac {1-z^{m}}{1-z}}F(z^{m})=\left(1+z+\cdots +z^{m-2}+z^{m-1}\right)F(z^{m}).}
ポ -再帰的シーケンスとホロノミック生成関数
定義
形式的な冪級数(または関数) F ( z )は、 次の形式の線型微分方程式を満たすとき、 ホロノミック であると言われる [17]。
c
0
(
z
)
F
(
r
)
(
z
)
+
c
1
(
z
)
F
(
r
−
1
)
(
z
)
+
⋯
+
c
r
(
z
)
F
(
z
)
=
0
,
{\displaystyle c_{0}(z)F^{(r)}(z)+c_{1}(z)F^{(r-1)}(z)+\cdots +c_{r}(z)F(z)=0,}
ここで係数 c i ( z ) は 有理関数の体にあります 。同様に、その導関数全体の張る ベクトル空間が 有限次元である場合、 はホロノミックです。
C
(
z
)
{\displaystyle \mathbb {C} (z)}
F
(
z
)
{\displaystyle F(z)}
C
(
z
)
{\displaystyle \mathbb {C} (z)}
前の式では必要に応じて分母を消去できるので、関数 c i ( z ) は z の多項式であると仮定できます。したがって、係数が次の 形式の
P 再帰 を満たす場合、生成関数がホロノミックであるという同等の条件を見ることができます。
c
^
s
(
n
)
f
n
+
s
+
c
^
s
−
1
(
n
)
f
n
+
s
−
1
+
⋯
+
c
^
0
(
n
)
f
n
=
0
,
{\displaystyle {\widehat {c}}_{s}(n)f_{n+s}+{\widehat {c}}_{s-1}(n)f_{n+s-1}+\cdots +{\widehat {c}}_{0}(n)f_{n}=0,}
十分に大きい n ≥ n 0 に対して成り立ち、 ĉ i ( n )は n 内の固定された有限次多項式です 。言い換えれば、シーケンスが P 再帰的であることとホロノミック生成関数を持つことは同値です。ホロノミック関数は、生成関数上の アダマール積 演算 ⊙ に関して閉じています 。
例
関数 e z 、 log z 、 cos z 、 arcsin z 、 、 二重対数 関数 Li 2 ( z ) 、 一般化された超幾何関数 p F q (...; ...; z ) 、およびべき級数によって定義される関数
1
+
z
{\displaystyle {\sqrt {1+z}}}
∑
n
=
0
∞
z
n
(
n
!
)
2
{\displaystyle \sum _{n=0}^{\infty }{\frac {z^{n}}{(n!)^{2}}}}
そして非収束的な
∑
n
=
0
∞
n
!
⋅
z
n
{\displaystyle \sum _{n=0}^{\infty }n!\cdot z^{n}}
すべてホロノミックです。
ホロノミック生成関数を持つP 再帰列 の例としては、 f n ≔ などが挙げられる。 1 / 1 + 1です ( 2 n n ) かつ f n ≔ 2 位 / 2 + 1 です 、ここで、および log n などのシーケンスは、 対応する生成関数の特異点の性質により P 再帰的では ありません。同様に、 tan z 、 sec z 、および Γ( z ) などの無限個の特異点を持つ関数はホロノミック関数
は ありません。
n
{\displaystyle {\sqrt {n}}}
作業用ソフトウェア ポ -再帰的シーケンスとホロノミック生成関数
Mathematica で P 再帰シーケンスを処理して操作するためのツールに は、RISC Combinatorics Group アルゴリズム組合せ論ソフトウェアサイトで非商用利用向けに提供されているソフトウェアパッケージが含まれます。ほとんどがクローズドソースであるにもかかわらず、このソフトウェアスイートの特に強力なツールは、 任意の入力シーケンスの P 再帰を 推測するためのパッケージ( 実験数学 と探究に便利) と、多くの合計の P 再帰を見つけ、 一般化 調和数を含む P 再帰の閉じた形式の解を解くことができるパッケージによって提供されます 。 [18] この特定の RISC サイトにリストされている他のパッケージは、 特に
ホロノミック 生成関数を扱うことを目的としています。 Guess Sigma
級数が 絶対収束する 場合、
は数列 a 0 、 a 1 、 ... の離散時間フーリエ変換です。
G
(
a
n
;
e
−
i
ω
)
=
∑
n
=
0
∞
a
n
e
−
i
ω
n
{\displaystyle G\left(a_{n};e^{-i\omega }\right)=\sum _{n=0}^{\infty }a_{n}e^{-i\omega n}}
数列の漸近的成長
微積分学では、冪級数の係数の成長率を使用して、冪級数の 収束半径を推測できることがよくあります。逆もまた成り立ち、生成関数の収束半径を使用して、基礎となる数列の 漸近的成長 を推測できることがよくあります 。
例えば、 有限の収束半径 rを持つ通常の生成関数 G ( a n ; x ) は次のように表すことができます。
G
(
a
n
;
x
)
=
A
(
x
)
+
B
(
x
)
(
1
−
x
r
)
−
β
x
α
{\displaystyle G(a_{n};x)={\frac {A(x)+B(x)\left(1-{\frac {x}{r}}\right)^{-\beta }}{x^{\alpha }}}}
ここで、 A ( x ) と B ( x ) はそれぞれ、 収束 半径が r より大きい(または 完全な )関数であり 、 B ( r ) ≠ 0 の場合には
、 ガンマ関数 、 二項係数 、または 多重集合係数を
使用します。n が無限大に近づくにつれて、これらの式のいずれに対する a n の比の極限 は1 になることが保証されます。単に a n が これらの式に比例する
だけではありません。
a
n
∼
B
(
r
)
r
α
Γ
(
β
)
n
β
−
1
(
1
r
)
n
∼
B
(
r
)
r
α
(
n
+
β
−
1
n
)
(
1
r
)
n
=
B
(
r
)
r
α
(
(
β
n
)
)
(
1
r
)
n
,
{\displaystyle a_{n}\sim {\frac {B(r)}{r^{\alpha }\Gamma (\beta )}}\,n^{\beta -1}\left({\frac {1}{r}}\right)^{n}\sim {\frac {B(r)}{r^{\alpha }}}{\binom {n+\beta -1}{n}}\left({\frac {1}{r}}\right)^{n}={\frac {B(r)}{r^{\alpha }}}\left(\!\!{\binom {\beta }{n}}\!\!\right)\left({\frac {1}{r}}\right)^{n}\,,}
このアプローチは、多くの場合、 n の 漸近級数の複数の項を生成するために反復される 。特に、
G
(
a
n
−
B
(
r
)
r
α
(
n
+
β
−
1
n
)
(
1
r
)
n
;
x
)
=
G
(
a
n
;
x
)
−
B
(
r
)
r
α
(
1
−
x
r
)
−
β
.
{\displaystyle G\left(a_{n}-{\frac {B(r)}{r^{\alpha }}}{\binom {n+\beta -1}{n}}\left({\frac {1}{r}}\right)^{n};x\right)=G(a_{n};x)-{\frac {B(r)}{r^{\alpha }}}\left(1-{\frac {x}{r}}\right)^{-\beta }\,.}
この生成関数の係数の漸近的増加は、 上記のように生成関数を記述するための
A 、 B 、 α 、 β 、および r を見つけることによって求めることができます。
同様の漸近解析は指数生成関数に対しても可能であり、指数生成関数の場合、それは な / ん ! は、これらの漸近公式に従って成長します。一般に、1 つのシーケンスの生成関数から 2 番目のシーケンスの生成関数を引いた収束半径が、個々の生成関数の収束半径よりも大きい場合、2 つのシーケンスは同じ漸近成長を持ちます。
正方形の列の漸近的成長
上記で導出されたように、平方数列の通常の生成関数は次のようになります。
G
(
n
2
;
x
)
=
x
(
x
+
1
)
(
1
−
x
)
3
.
{\displaystyle G(n^{2};x)={\frac {x(x+1)}{(1-x)^{3}}}.}
r = 1 、 α = −1 、 β = 3 、 A ( x ) = 0 、 B ( x ) = x + 1 とすると 、正方形が予想どおりに大きくなることが確認できます。
a
n
∼
B
(
r
)
r
α
Γ
(
β
)
n
β
−
1
(
1
r
)
n
=
1
+
1
1
−
1
Γ
(
3
)
n
3
−
1
(
1
1
)
n
=
n
2
.
{\displaystyle a_{n}\sim {\frac {B(r)}{r^{\alpha }\Gamma (\beta )}}\,n^{\beta -1}\left({\frac {1}{r}}\right)^{n}={\frac {1+1}{1^{-1}\,\Gamma (3)}}\,n^{3-1}\left({\frac {1}{1}}\right)^{n}=n^{2}.}
カタラン数の漸近的増加
カタラン数 の通常の生成関数 は
G
(
C
n
;
x
)
=
1
−
1
−
4
x
2
x
.
{\displaystyle G(C_{n};x)={\frac {1-{\sqrt {1-4x}}}{2x}}.}
r = の場合 1 / 4 、 α = 1 、 β = − 1 / 2 、 A ( x ) = 1 / 2 、そして B ( x ) = − 1 / 2 、カタロニア語の数字については次のように結論付けることができます。
C
n
∼
B
(
r
)
r
α
Γ
(
β
)
n
β
−
1
(
1
r
)
n
=
−
1
2
(
1
4
)
1
Γ
(
−
1
2
)
n
−
1
2
−
1
(
1
1
4
)
n
=
4
n
n
3
2
π
.
{\displaystyle C_{n}\sim {\frac {B(r)}{r^{\alpha }\Gamma (\beta )}}\,n^{\beta -1}\left({\frac {1}{r}}\right)^{n}={\frac {-{\frac {1}{2}}}{\left({\frac {1}{4}}\right)^{1}\Gamma \left(-{\frac {1}{2}}\right)}}\,n^{-{\frac {1}{2}}-1}\left({\frac {1}{\,{\frac {1}{4}}\,}}\right)^{n}={\frac {4^{n}}{n^{\frac {3}{2}}{\sqrt {\pi }}}}.}
二変量および多変量生成関数
複数の変数における生成関数は、複数のインデックスを持つ配列に一般化できます。これらの非多項式の二重和の例は、 多変量生成関数 または スーパー生成関数 と呼ばれます。2 つの変数の場合、これらはしばしば 2 変量生成関数 と呼ばれます。
二変量の場合
2次元配列 a m , n (ここで n と m は自然数)の通常の生成関数は次のようになります。
たとえば、 (1 + x ) n は 固定された nに対する 二項係数 の通常の生成関数であるため 、二項係数を生成する2変量生成関数を求めることができます (
G
(
a
m
,
n
;
x
,
y
)
=
∑
m
,
n
=
0
∞
a
m
,
n
x
m
y
n
.
{\displaystyle G(a_{m,n};x,y)=\sum _{m,n=0}^{\infty }a_{m,n}x^{m}y^{n}.}
nk ) が すべてのk と n に対して成り立つ 。これを行うには、 (1 + x ) n 自体をn 内のシーケンスとして 考え、これらのシーケンス値を係数として持つy 内の生成関数を見つけます 。 n の生成関数は次のようになるため 、 二
項係数の生成関数は次のようになります。
その他の例としては、 二項係数 、 スターリング数 、 オイラー数 に対する次の 2 変数生成関数があります。ここで、 ω と z は 2 つの変数を表します。 [19]
1
1
−
a
y
,
{\displaystyle {\frac {1}{1-ay}},}
∑
n
,
k
(
n
k
)
x
k
y
n
=
1
1
−
(
1
+
x
)
y
=
1
1
−
y
−
x
y
.
{\displaystyle \sum _{n,k}{\binom {n}{k}}x^{k}y^{n}={\frac {1}{1-(1+x)y}}={\frac {1}{1-y-xy}}.}
e
z
+
w
z
=
∑
m
,
n
≥
0
(
n
m
)
w
m
z
n
n
!
e
w
(
e
z
−
1
)
=
∑
m
,
n
≥
0
{
n
m
}
w
m
z
n
n
!
1
(
1
−
z
)
w
=
∑
m
,
n
≥
0
[
n
m
]
w
m
z
n
n
!
1
−
w
e
(
w
−
1
)
z
−
w
=
∑
m
,
n
≥
0
⟨
n
m
⟩
w
m
z
n
n
!
e
w
−
e
z
w
e
z
−
z
e
w
=
∑
m
,
n
≥
0
⟨
m
+
n
+
1
m
⟩
w
m
z
n
(
m
+
n
+
1
)
!
.
{\displaystyle {\begin{aligned}e^{z+wz}&=\sum _{m,n\geq 0}{\binom {n}{m}}w^{m}{\frac {z^{n}}{n!}}\\[4px]e^{w(e^{z}-1)}&=\sum _{m,n\geq 0}{\begin{Bmatrix}n\\m\end{Bmatrix}}w^{m}{\frac {z^{n}}{n!}}\\[4px]{\frac {1}{(1-z)^{w}}}&=\sum _{m,n\geq 0}{\begin{bmatrix}n\\m\end{bmatrix}}w^{m}{\frac {z^{n}}{n!}}\\[4px]{\frac {1-w}{e^{(w-1)z}-w}}&=\sum _{m,n\geq 0}\left\langle {\begin{matrix}n\\m\end{matrix}}\right\rangle w^{m}{\frac {z^{n}}{n!}}\\[4px]{\frac {e^{w}-e^{z}}{we^{z}-ze^{w}}}&=\sum _{m,n\geq 0}\left\langle {\begin{matrix}m+n+1\\m\end{matrix}}\right\rangle {\frac {w^{m}z^{n}}{(m+n+1)!}}.\end{aligned}}}
多変量ケース
多変量生成関数は、実際には、指定された行と列の合計を持つ非負整数の 分割表 の数を計算するときに発生します。表に r 行と c 列があるとします。行の合計は t 1 、 t 2 ... t r で、列の合計は s 1 、 s 2 ... s c です。IJ Good によると 、 [20] このような表の数は、次の係数です
。
x
1
t
1
⋯
x
r
t
r
y
1
s
1
⋯
y
c
s
c
{\displaystyle x_{1}^{t_{1}}\cdots x_{r}^{t_{r}}y_{1}^{s_{1}}\cdots y_{c}^{s_{c}}}
∏
i
=
1
r
∏
j
=
1
c
1
1
−
x
i
y
j
.
{\displaystyle \prod _{i=1}^{r}\prod _{j=1}^{c}{\frac {1}{1-x_{i}y_{j}}}.}
連分数による表現(ヤコビ型) J -分数)
定義
h次の有理数収束が 2 h 次精度の べき級数 を表す (形式的な) ヤコビ型 連分数( それぞれ J 分数 と S 分数 )と スティルチェス型 連分数の展開は、多くの特殊な 1 変数および 2 変数シーケンスの典型的に発散する通常の生成関数を表現する別の方法です。ヤコビ型連分数 ( J 分数) の特定の形式は、次の式のように展開され、 特定のアプリケーション依存のコンポーネント シーケンス {ab i } と { c i }について、 z に関する次の対応するべき級数展開を持ちます。ここで、 z ≠ 0 は、以下に示す 2 番目のべき級数展開の形式的な変数を表します。 [21]
J
[
∞
]
(
z
)
=
1
1
−
c
1
z
−
ab
2
z
2
1
−
c
2
z
−
ab
3
z
2
⋱
=
1
+
c
1
z
+
(
ab
2
+
c
1
2
)
z
2
+
(
2
ab
2
c
1
+
c
1
3
+
ab
2
c
2
)
z
3
+
⋯
{\displaystyle {\begin{aligned}J^{[\infty ]}(z)&={\cfrac {1}{1-c_{1}z-{\cfrac {{\text{ab}}_{2}z^{2}}{1-c_{2}z-{\cfrac {{\text{ab}}_{3}z^{2}}{\ddots }}}}}}\\[4px]&=1+c_{1}z+\left({\text{ab}}_{2}+c_{1}^{2}\right)z^{2}+\left(2{\text{ab}}_{2}c_{1}+c_{1}^{3}+{\text{ab}}_{2}c_{2}\right)z^{3}+\cdots \end{aligned}}}
前の式でj n ≔ [ z n ] J [∞] ( z ) と略記された の係数は 、次の式の行列解に対応する。
z
n
{\displaystyle z^{n}}
[
k
0
,
1
k
1
,
1
0
0
⋯
k
0
,
2
k
1
,
2
k
2
,
2
0
⋯
k
0
,
3
k
1
,
3
k
2
,
3
k
3
,
3
⋯
⋮
⋮
⋮
⋮
]
=
[
k
0
,
0
0
0
0
⋯
k
0
,
1
k
1
,
1
0
0
⋯
k
0
,
2
k
1
,
2
k
2
,
2
0
⋯
⋮
⋮
⋮
⋮
]
⋅
[
c
1
1
0
0
⋯
ab
2
c
2
1
0
⋯
0
ab
3
c
3
1
⋯
⋮
⋮
⋮
⋮
]
,
{\displaystyle {\begin{bmatrix}k_{0,1}&k_{1,1}&0&0&\cdots \\k_{0,2}&k_{1,2}&k_{2,2}&0&\cdots \\k_{0,3}&k_{1,3}&k_{2,3}&k_{3,3}&\cdots \\\vdots &\vdots &\vdots &\vdots \end{bmatrix}}={\begin{bmatrix}k_{0,0}&0&0&0&\cdots \\k_{0,1}&k_{1,1}&0&0&\cdots \\k_{0,2}&k_{1,2}&k_{2,2}&0&\cdots \\\vdots &\vdots &\vdots &\vdots \end{bmatrix}}\cdot {\begin{bmatrix}c_{1}&1&0&0&\cdots \\{\text{ab}}_{2}&c_{2}&1&0&\cdots \\0&{\text{ab}}_{3}&c_{3}&1&\cdots \\\vdots &\vdots &\vdots &\vdots \end{bmatrix}},}
ここで、 j 0 ≡ k 0,0 = 1 、 j n = k 0, n ( n ≥ 1 ) 、 k r , s = 0 ( r > s の 場合 )、すべての整数 p 、 q ≥ 0に対して、次の 加法 関係式が成り立ちます 。
j
p
+
q
=
k
0
,
p
⋅
k
0
,
q
+
∑
i
=
1
min
(
p
,
q
)
ab
2
⋯
ab
i
+
1
×
k
i
,
p
⋅
k
i
,
q
.
{\displaystyle j_{p+q}=k_{0,p}\cdot k_{0,q}+\sum _{i=1}^{\min(p,q)}{\text{ab}}_{2}\cdots {\text{ab}}_{i+1}\times k_{i,p}\cdot k_{i,q}.}
の特性 h 収束関数
h ≥ 0 の場合(実際には h ≥ 2 の 場合)、 無限 J 分数 J [∞] ( z )に収束する有理数 h を 次のように定義できます 。
Conv
h
(
z
)
:=
P
h
(
z
)
Q
h
(
z
)
=
j
0
+
j
1
z
+
⋯
+
j
2
h
−
1
z
2
h
−
1
+
∑
n
=
2
h
∞
j
~
h
,
n
z
n
{\displaystyle \operatorname {Conv} _{h}(z):={\frac {P_{h}(z)}{Q_{h}(z)}}=j_{0}+j_{1}z+\cdots +j_{2h-1}z^{2h-1}+\sum _{n=2h}^{\infty }{\widetilde {j}}_{h,n}z^{n}}
次のように再帰的に定義されるシーケンス P h ( z ) と Q h ( z ) を成分ごとに調べる:
P
h
(
z
)
=
(
1
−
c
h
z
)
P
h
−
1
(
z
)
−
ab
h
z
2
P
h
−
2
(
z
)
+
δ
h
,
1
Q
h
(
z
)
=
(
1
−
c
h
z
)
Q
h
−
1
(
z
)
−
ab
h
z
2
Q
h
−
2
(
z
)
+
(
1
−
c
1
z
)
δ
h
,
1
+
δ
0
,
1
.
{\displaystyle {\begin{aligned}P_{h}(z)&=(1-c_{h}z)P_{h-1}(z)-{\text{ab}}_{h}z^{2}P_{h-2}(z)+\delta _{h,1}\\Q_{h}(z)&=(1-c_{h}z)Q_{h-1}(z)-{\text{ab}}_{h}z^{2}Q_{h-2}(z)+(1-c_{1}z)\delta _{h,1}+\delta _{0,1}.\end{aligned}}}
さらに、すべてのh ≥ 2 に対して 収束関数 Conv h ( z )が有理的であることは、 j n のシーケンスが満たす追加の差分方程式と合同性を意味し 、 M h ≔ ab 2 ⋯ ab h + 1 の 場合、 h ‖ M h であれば、合同性が得られます。
j
n
≡
[
z
n
]
Conv
h
(
z
)
(
mod
h
)
,
{\displaystyle j_{n}\equiv [z^{n}]\operatorname {Conv} _{h}(z){\pmod {h}},}
h ≥ 2 の場合 、つまり、以下の表の例のように、これらのシーケンスが q 、 x 、 R などの補助パラメータに暗黙的に依存しない場合、パラメータシーケンス {ab i } と { c i } の非記号的で確定的な選択。
例
次の表は、最初のサブセクションで定義した J 分数の一般的な展開によって生成された所定のシーケンス j n のいくつかの特殊なケースで計算によって見つかった(そして後に引用文献 [22]で正しいことが証明された)成分シーケンスの閉じた形式の式の例を示しています。ここでは、 0 < | a |、| b |、| q | < 1 と定義し 、パラメーター と x は これらの展開に関して不定であるとします。ここで、これらの J分数の展開によって列挙された所定のシーケンスは、 q -ポッホハマー記号 、 ポッホハマー記号 、および 二項係数 で定義されます 。
R
,
α
∈
Z
+
{\displaystyle R,\alpha \in \mathbb {Z} ^{+}}
上記のヤコビ型 J 分数の定義に対応するこれらの級数の収束半径は、一般に、これらの数列の通常の生成関数を定義する対応するべき級数展開の収束半径とは異なります。
例
平方数
平方数列 a n = n 2 の生成関数は次 のとおりです。
ここで ζ ( s)は リーマンゼータ関数 である 。
アプリケーション
生成関数は次の目的で使用されます。
フィボナッチ数列 などの再帰関係で与えられた数列の 閉じた式 を見つけます 。
シーケンスの再帰関係 を見つけます 。生成関数の形式から再帰式が推測される場合があります。
シーケンス間の関係を見つけます。2 つのシーケンスの生成関数が類似した形式である場合、シーケンス自体が関連している可能性があります。
シーケンスの漸近的な動作を調べます。
シーケンスを含む同一性を証明します。
組合せ論 における 列挙 問題を解き 、その解をエンコードします。 ルーク多項式は 組合せ論の応用例です。
無限和を評価します。
さまざまなテクニック: 和の評価と生成関数を使ったその他の問題への取り組み
生成関数は、合計を操作したり、合計間の同一性を確立したりするためのいくつかの方法を提供します。
最も単純なケースは s n = Σのときである。 nk = 0 です a k 。すると、 S ( z ) = あ ( ず ) / 1 − z 対応する通常の生成関数の場合。
例えば、
H k = 1 +
を操作できる。
s
n
=
∑
k
=
1
n
H
k
,
{\displaystyle s_{n}=\sum _{k=1}^{n}H_{k}\,,}
1 / 2 + ⋯ + 1 / け は 調和数 です 。
を調和数の通常の生成関数とします。すると
、
H
(
z
)
=
∑
n
=
1
∞
H
n
z
n
{\displaystyle H(z)=\sum _{n=1}^{\infty }{H_{n}z^{n}}}
H
(
z
)
=
1
1
−
z
∑
n
=
1
∞
z
n
n
,
{\displaystyle H(z)={\frac {1}{1-z}}\sum _{n=1}^{\infty }{\frac {z^{n}}{n}}\,,}
S
(
z
)
=
∑
n
=
1
∞
s
n
z
n
=
1
(
1
−
z
)
2
∑
n
=
1
∞
z
n
n
.
{\displaystyle S(z)=\sum _{n=1}^{\infty }{s_{n}z^{n}}={\frac {1}{(1-z)^{2}}}\sum _{n=1}^{\infty }{\frac {z^{n}}{n}}\,.}
分子との畳み込みを
使うと
、次のようにも書ける。
1
(
1
−
z
)
2
=
∑
n
=
0
∞
(
n
+
1
)
z
n
,
{\displaystyle {\frac {1}{(1-z)^{2}}}=\sum _{n=0}^{\infty }(n+1)z^{n}\,,}
s
n
=
∑
k
=
1
n
n
+
1
−
k
k
=
(
n
+
1
)
H
n
−
n
,
{\displaystyle s_{n}=\sum _{k=1}^{n}{\frac {n+1-k}{k}}=(n+1)H_{n}-n\,,}
∑
k
=
1
n
H
k
=
(
n
+
1
)
(
H
n
+
1
−
1
)
.
{\displaystyle \sum _{k=1}^{n}{H_{k}}=(n+1)(H_{n+1}-1)\,.}
生成関数を使用してシーケンスを関連付け、合計を操作する別の例として、任意のシーケンス ⟨ f n ⟩ に対して、
すべての n ≥ 0 に対して2 つの合計シーケンスを定義し
、2 番目の合計を最初の合計で表現しようとします。生成関数によるアプローチを提案します。
s
n
:=
∑
m
=
0
n
(
n
m
)
f
m
3
n
−
m
s
~
n
:=
∑
m
=
0
n
(
n
m
)
(
m
+
1
)
(
m
+
2
)
(
m
+
3
)
f
m
3
n
−
m
,
{\displaystyle {\begin{aligned}s_{n}&:=\sum _{m=0}^{n}{\binom {n}{m}}f_{m}3^{n-m}\\[4px]{\tilde {s}}_{n}&:=\sum _{m=0}^{n}{\binom {n}{m}}(m+1)(m+2)(m+3)f_{m}3^{n-m}\,,\end{aligned}}}
まず、 二項変換を 使用して、最初の合計の生成関数を次のように書きます。
S
(
z
)
=
1
1
−
3
z
F
(
z
1
−
3
z
)
.
{\displaystyle S(z)={\frac {1}{1-3z}}F\left({\frac {z}{1-3z}}\right).}
数列 ⟨ ( n + 1)( n + 2)( n + 3) f n ⟩ の生成関数は次のように表される
ので
、上で定義した2番目の和の生成関数は次の形式で表すことができる。
6
F
(
z
)
+
18
z
F
′
(
z
)
+
9
z
2
F
″
(
z
)
+
z
3
F
‴
(
z
)
{\displaystyle 6F(z)+18zF'(z)+9z^{2}F''(z)+z^{3}F'''(z)}
S
~
(
z
)
=
6
(
1
−
3
z
)
F
(
z
1
−
3
z
)
+
18
z
(
1
−
3
z
)
2
F
′
(
z
1
−
3
z
)
+
9
z
2
(
1
−
3
z
)
3
F
″
(
z
1
−
3
z
)
+
z
3
(
1
−
3
z
)
4
F
‴
(
z
1
−
3
z
)
.
{\displaystyle {\tilde {S}}(z)={\frac {6}{(1-3z)}}F\left({\frac {z}{1-3z}}\right)+{\frac {18z}{(1-3z)^{2}}}F'\left({\frac {z}{1-3z}}\right)+{\frac {9z^{2}}{(1-3z)^{3}}}F''\left({\frac {z}{1-3z}}\right)+{\frac {z^{3}}{(1-3z)^{4}}}F'''\left({\frac {z}{1-3z}}\right).}
特に、この修正和生成関数は、
a ( z ) = 6(1 − 3 z ) 3 、 b ( z ) = 18(1 − 3 z ) 3 、 c ( z ) = 9(1 − 3 z ) 3 、および d ( z ) = (1 − 3 z ) 3 (ただし、 (1 − 3 z ) 3 = 1 − 9 z + 27 z 2 − 27 z 3 )の形式で表記できます。
a
(
z
)
⋅
S
(
z
)
+
b
(
z
)
⋅
z
S
′
(
z
)
+
c
(
z
)
⋅
z
2
S
″
(
z
)
+
d
(
z
)
⋅
z
3
S
‴
(
z
)
,
{\displaystyle a(z)\cdot S(z)+b(z)\cdot zS'(z)+c(z)\cdot z^{2}S''(z)+d(z)\cdot z^{3}S'''(z),}
最後に、2 番目の合計から 1 番目の合計までを次の形式で表すことができます。
s
~
n
=
[
z
n
]
(
6
(
1
−
3
z
)
3
∑
n
=
0
∞
s
n
z
n
+
18
(
1
−
3
z
)
3
∑
n
=
0
∞
n
s
n
z
n
+
9
(
1
−
3
z
)
3
∑
n
=
0
∞
n
(
n
−
1
)
s
n
z
n
+
(
1
−
3
z
)
3
∑
n
=
0
∞
n
(
n
−
1
)
(
n
−
2
)
s
n
z
n
)
=
(
n
+
1
)
(
n
+
2
)
(
n
+
3
)
s
n
−
9
n
(
n
+
1
)
(
n
+
2
)
s
n
−
1
+
27
(
n
−
1
)
n
(
n
+
1
)
s
n
−
2
−
(
n
−
2
)
(
n
−
1
)
n
s
n
−
3
.
{\displaystyle {\begin{aligned}{\tilde {s}}_{n}&=[z^{n}]\left(6(1-3z)^{3}\sum _{n=0}^{\infty }s_{n}z^{n}+18(1-3z)^{3}\sum _{n=0}^{\infty }ns_{n}z^{n}+9(1-3z)^{3}\sum _{n=0}^{\infty }n(n-1)s_{n}z^{n}+(1-3z)^{3}\sum _{n=0}^{\infty }n(n-1)(n-2)s_{n}z^{n}\right)\\[4px]&=(n+1)(n+2)(n+3)s_{n}-9n(n+1)(n+2)s_{n-1}+27(n-1)n(n+1)s_{n-2}-(n-2)(n-1)ns_{n-3}.\end{aligned}}}
例3: 相互再帰シーケンスの生成関数
この例では、 Concrete Mathematics のセクション 7.3 に示されている生成関数の例を再定式化します (生成関数のシリーズのきれいな図については、同じ参考文献のセクション 7.1 も参照)。特に、 3 x n の 長方形にマークされていない 2 x 1 のドミノ ピースを並べる方法の総数 ( U n で表記) を求めるものとします。補助シーケンスV n を、完全な長方形の 3 x n の長方形から角を除いた部分を覆う方法の数として定義します 。これらの定義を使用して、垂直ドミノと水平ドミノの場合を処理するためにこの定義をさらに細分化することなく、 U n の 閉じた形式の 式を示します。2 つのシーケンスの通常の生成関数は、次のシリーズに対応していることに注意してください。
U
(
z
)
=
1
+
3
z
2
+
11
z
4
+
41
z
6
+
⋯
,
V
(
z
)
=
z
+
4
z
3
+
15
z
5
+
56
z
7
+
⋯
.
{\displaystyle {\begin{aligned}U(z)=1+3z^{2}+11z^{4}+41z^{6}+\cdots ,\\V(z)=z+4z^{3}+15z^{5}+56z^{7}+\cdots .\end{aligned}}}
3行n 列の長方形の左端から始まる可能な構成を考慮すると、U 0 = 1、U 1 = 0 、 V 0 = 0 、 V 1 = 1 の場合、上記のように定義された n ≥ 2の 2つのシーケンスについて、次の 相互 依存的、または 相互再帰的な 再帰関係を 表すことができます 。
U
n
=
2
V
n
−
1
+
U
n
−
2
V
n
=
U
n
−
1
+
V
n
−
2
.
{\displaystyle {\begin{aligned}U_{n}&=2V_{n-1}+U_{n-2}\\V_{n}&=U_{n-1}+V_{n-2}.\end{aligned}}}
すべての整数m ≥ 0 に対して 、インデックスシフトされた生成関数は [注1]
を満たすので
、上で指定した初期条件と前の2つの再帰関係を使用して、これらのシーケンスの生成関数を関連付ける次の2つの方程式が得られることがわかります。
これは、方程式系を解くことによって(これがここでの私たちの方法の特別なトリックです)、
z
m
G
(
z
)
=
∑
n
=
m
∞
g
n
−
m
z
n
,
{\displaystyle z^{m}G(z)=\sum _{n=m}^{\infty }g_{n-m}z^{n}\,,}
U
(
z
)
=
2
z
V
(
z
)
+
z
2
U
(
z
)
+
1
V
(
z
)
=
z
U
(
z
)
+
z
2
V
(
z
)
=
z
1
−
z
2
U
(
z
)
,
{\displaystyle {\begin{aligned}U(z)&=2zV(z)+z^{2}U(z)+1\\V(z)&=zU(z)+z^{2}V(z)={\frac {z}{1-z^{2}}}U(z),\end{aligned}}}
U
(
z
)
=
1
−
z
2
1
−
4
z
2
+
z
4
=
1
3
−
3
⋅
1
1
−
(
2
+
3
)
z
2
+
1
3
+
3
⋅
1
1
−
(
2
−
3
)
z
2
.
{\displaystyle U(z)={\frac {1-z^{2}}{1-4z^{2}+z^{4}}}={\frac {1}{3-{\sqrt {3}}}}\cdot {\frac {1}{1-\left(2+{\sqrt {3}}\right)z^{2}}}+{\frac {1}{3+{\sqrt {3}}}}\cdot {\frac {1}{1-\left(2-{\sqrt {3}}\right)z^{2}}}.}
したがって、前の式の生成関数の 2 番目の部分分数展開から得られるシーケンスに代数的簡略化を実行すると、 U 2 n + 1 ≡ 0 であり、
すべての整数 n ≥ 0に対してであることがわかります。また、 フィボナッチ数列 の 2 次 再帰に適用された同じシフトされた生成関数手法は、上記の 有理関数 のサブセクションですでに説明されている、または少なくともヒントが示されている、1 つの変数の再帰関係を解くために生成関数を使用する典型的な例であることにも注意してください 。
U
2
n
=
⌈
(
2
+
3
)
n
3
−
3
⌉
,
{\displaystyle U_{2n}=\left\lceil {\frac {\left(2+{\sqrt {3}}\right)^{n}}{3-{\sqrt {3}}}}\right\rceil \,,}
畳み込み(コーシー積)
2 つの形式的なべき級数の項の離散 畳み込みにより 、生成関数の積が、元の流れの項の畳み込み和を列挙する生成関数に変換されます ( コーシー積を 参照)。
A ( z ) と B ( z ) が通常の生成関数であるとします 。
C
(
z
)
=
A
(
z
)
B
(
z
)
⇔
[
z
n
]
C
(
z
)
=
∑
k
=
0
n
a
k
b
n
−
k
{\displaystyle C(z)=A(z)B(z)\Leftrightarrow [z^{n}]C(z)=\sum _{k=0}^{n}{a_{k}b_{n-k}}}
A ( z ) と B ( z ) が指数生成関数である とします。
C
(
z
)
=
A
(
z
)
B
(
z
)
⇔
[
z
n
n
!
]
C
(
z
)
=
∑
k
=
0
n
(
n
k
)
a
k
b
n
−
k
{\displaystyle C(z)=A(z)B(z)\Leftrightarrow \left[{\frac {z^{n}}{n!}}\right]C(z)=\sum _{k=0}^{n}{\binom {n}{k}}a_{k}b_{n-k}}
3つの通常の生成関数の積から得られる三重畳列を考える
C
(
z
)
=
F
(
z
)
G
(
z
)
H
(
z
)
⇔
[
z
n
]
C
(
z
)
=
∑
j
+
k
+
l
=
n
f
j
g
k
h
l
{\displaystyle C(z)=F(z)G(z)H(z)\Leftrightarrow [z^{n}]C(z)=\sum _{j+k+l=n}f_{j}g_{k}h_{l}}
ある正の整数m ≥ 1 に対して、シーケンスの m 倍畳み込みを考える (応用については以下の例を参照)。
C
(
z
)
=
G
(
z
)
m
⇔
[
z
n
]
C
(
z
)
=
∑
k
1
+
k
2
+
⋯
+
k
m
=
n
g
k
1
g
k
2
⋯
g
k
m
{\displaystyle C(z)=G(z)^{m}\Leftrightarrow [z^{n}]C(z)=\sum _{k_{1}+k_{2}+\cdots +k_{m}=n}g_{k_{1}}g_{k_{2}}\cdots g_{k_{m}}}
生成関数の乗算、またはそれらの基礎となるシーケンスの畳み込みは、特定のカウントおよび確率シナリオにおける独立したイベントの概念に対応することができます。たとえば、 ランダム変数 Zの 確率生成関数 、または pgfを G Z ( z ) と表記するという表記規則を採用すると、 X と Y が独立している
場合、 任意の 2 つのランダム変数 [23]に対して、次のことを示すことができます。同様に、 n ≥ 0 セントをセット {1、5、10、25、50} の値の硬貨
額面(それぞれ、ペニー、ニッケル、ダイム、クォーター、ハーフ ドル) で支払う方法の数は、積によって生成され、さらに、 n セントを任意の正の整数額面の硬貨で支払うことを
許可した場合、 分割関数 生成関数によって生成されるこのような変更の組み合わせの数の生成に到達します。この関数は 、無限の q -ポッホハマー記号 積によって拡張されます。
G
X
+
Y
(
z
)
=
G
X
(
z
)
G
Y
(
z
)
,
{\displaystyle G_{X+Y}(z)=G_{X}(z)G_{Y}(z)\,,}
C
(
z
)
=
1
1
−
z
1
1
−
z
5
1
1
−
z
10
1
1
−
z
25
1
1
−
z
50
,
{\displaystyle C(z)={\frac {1}{1-z}}{\frac {1}{1-z^{5}}}{\frac {1}{1-z^{10}}}{\frac {1}{1-z^{25}}}{\frac {1}{1-z^{50}}},}
∏
n
=
1
∞
(
1
−
z
n
)
−
1
.
{\displaystyle \prod _{n=1}^{\infty }\left(1-z^{n}\right)^{-1}\,.}
例: カタラン数の生成関数
生成関数の畳み込みが有用な例として、カタラン数 の通常の生成関数を表す特定の閉形式関数 C n を解くことができます 。 特に 、 この シーケンスは、積 x 0 · x 1 ·⋯· x n に括弧を挿入して乗算の順序を完全に指定する方法の数として組み合わせ論的に解釈されます。たとえば、 C 2 = 2 は、2 つの式 x 0 · ( x 1 · x 2 ) と ( x 0 · x 1 ) · x 2 に対応します。したがって、シーケンスは、次に示す再帰関係を満たし
、対応する畳み込み生成関数 C ( z ) が次を満たし
ます。
C
n
=
∑
k
=
0
n
−
1
C
k
C
n
−
1
−
k
+
δ
n
,
0
=
C
0
C
n
−
1
+
C
1
C
n
−
2
+
⋯
+
C
n
−
1
C
0
+
δ
n
,
0
,
n
≥
0
,
{\displaystyle C_{n}=\sum _{k=0}^{n-1}C_{k}C_{n-1-k}+\delta _{n,0}=C_{0}C_{n-1}+C_{1}C_{n-2}+\cdots +C_{n-1}C_{0}+\delta _{n,0}\,,\quad n\geq 0\,,}
C
(
z
)
=
z
⋅
C
(
z
)
2
+
1
.
{\displaystyle C(z)=z\cdot C(z)^{2}+1\,.}
C (0) = 1 ≠ ∞ なので 、この生成関数の式は次のようになる。
C
(
z
)
=
1
−
1
−
4
z
2
z
=
∑
n
=
0
∞
1
n
+
1
(
2
n
n
)
z
n
.
{\displaystyle C(z)={\frac {1-{\sqrt {1-4z}}}{2z}}=\sum _{n=0}^{\infty }{\frac {1}{n+1}}{\binom {2n}{n}}z^{n}\,.}
上記の
C ( z ) を 暗黙的に定義する最初の式は、
この生成関数の別の「単純な」(形式の) 連分数展開につながることを意味していることに注意してください。
C
(
z
)
=
1
1
−
z
⋅
C
(
z
)
,
{\displaystyle C(z)={\frac {1}{1-z\cdot C(z)}}\,,}
例: ファンの全域木と畳み込みの畳み込み
次数 n のファンとは、頂点 {0, 1, ..., n } 上のグラフ で、 2 n − 1 本の辺が次の規則に従って接続されているものと定義されます。頂点 0 は、他の n 個の頂点のそれぞれに 1 本の辺で接続されており、頂点は、 常に 1 ≤ k < nに対して、次の頂点 k + 1 に 1 本の辺で接続されています 。 [24] 次数 1 のファンは 1 つ、次数 2 のファンは 3 つ、次数 3 のファンは 8 つ、というように存在します。 全域木は 、元の頂点がすべて含まれ、この部分グラフを接続するのに十分な数の辺を含むグラフの部分グラフですが、部分グラフにサイクルが存在するほど多くの辺は含まれません。各 n ≥ 1に対して、次数 n のファンの全域木 f n がいくつ存在するかを調べます 。
k
{\displaystyle k}
観察として、隣接する頂点の集合を結合する方法の数を数えることで、この問題にアプローチすることができます。たとえば、 n = 4 のとき、 f 4 = 4 + 3 · 1 + 2 · 2 + 1 · 3 + 2 · 1 · 1 + 1 · 2 · 1 + 1 · 1 · 2 + 1 · 1 · 1 · 1 = 21 となり 、これは シーケンス g n = n = [ z n ] の m倍畳み込みの合計です 。 ず / (1− z ) 2 ただし、 m ≔ 1, 2, 3, 4 です 。より一般的には、この数列の式を次のように書くことができます。
この数列の通常の生成関数は、次の畳み込みの合計で与えられることがわかります。この畳み込みの合計から、
最後の生成関数の部分分数展開を
取ることで、数列の正確な式を抽出できます 。
f
n
=
∑
m
>
0
∑
k
1
+
k
2
+
⋯
+
k
m
=
n
k
1
,
k
2
,
…
,
k
m
>
0
g
k
1
g
k
2
⋯
g
k
m
,
{\displaystyle f_{n}=\sum _{m>0}\sum _{\scriptstyle k_{1}+k_{2}+\cdots +k_{m}=n \atop \scriptstyle k_{1},k_{2},\ldots ,k_{m}>0}g_{k_{1}}g_{k_{2}}\cdots g_{k_{m}}\,,}
F
(
z
)
=
G
(
z
)
+
G
(
z
)
2
+
G
(
z
)
3
+
⋯
=
G
(
z
)
1
−
G
(
z
)
=
z
(
1
−
z
)
2
−
z
=
z
1
−
3
z
+
z
2
,
{\displaystyle F(z)=G(z)+G(z)^{2}+G(z)^{3}+\cdots ={\frac {G(z)}{1-G(z)}}={\frac {z}{(1-z)^{2}-z}}={\frac {z}{1-3z+z^{2}}}\,,}
明示的な指定ではなく、関数方程式で指定される生成関数によく遭遇する。例えば、 n 個のノード(葉を含む)
上の二分木の数の生成関数 T(z)は、
T
(
z
)
=
z
(
1
+
T
(
z
)
2
)
{\displaystyle T(z)=z\left(1+T(z)^{2}\right)}
ラグランジュ の逆定理は、 このような方程式の解を明示的に評価するために使用されるツールです。
上記の定理を関数方程式に適用すると、次の式が得られます ( )。
ϕ
(
z
)
=
1
+
z
2
{\textstyle \phi (z)=1+z^{2}}
[
z
n
]
T
(
z
)
=
[
z
n
−
1
]
1
n
(
1
+
z
2
)
n
{\displaystyle [z^{n}]T(z)=[z^{n-1}]{\frac {1}{n}}(1+z^{2})^{n}}
二項定理展開により、 が偶数の場合 、式は を返します 。これは、二分木の葉の数が内部ノードの数より1つ多いことが証明されているため、合計は常に奇数になるはずであるため、予想どおりです。 ただし、 が奇数の場合、
n
{\displaystyle n}
0
{\displaystyle 0}
n
{\displaystyle n}
[
z
n
−
1
]
1
n
(
1
+
z
2
)
n
=
1
n
(
n
n
+
1
2
)
{\displaystyle [z^{n-1}]{\frac {1}{n}}(1+z^{2})^{n}={\frac {1}{n}}{\dbinom {n}{\frac {n+1}{2}}}}
を内部ノードの数と すると、式ははるかに簡潔になります。これで、式は単に 番目の カタロニア数になります。
n
{\displaystyle n}
n
{\displaystyle n}
フリーパラメータの導入(スネークオイル法)
合計 s n は 複雑な場合があり、必ずしも簡単に評価できるとは限りません。「フリー パラメーター」法は、これらの合計を評価する別の方法です (H. Wilf は「snake oil」と呼んでいます)。
これまで説明した両方の方法 では、合計の限界として nがあります。合計に n が明示的に表示されない場合は、 n を 「自由」パラメータと見なし、 s n を F ( z ) = Σ s n z n の係数として扱い、 n と k の合計の順序を変更して 、内部和を計算します。
例えば、計算したい場合、
nを 「自由」なパラメータとして
扱い、
s
n
=
∑
k
=
0
∞
(
n
+
k
m
+
2
k
)
(
2
k
k
)
(
−
1
)
k
k
+
1
,
m
,
n
∈
N
0
,
{\displaystyle s_{n}=\sum _{k=0}^{\infty }{{\binom {n+k}{m+2k}}{\binom {2k}{k}}{\frac {(-1)^{k}}{k+1}}}\,,\quad m,n\in \mathbb {N} _{0}\,,}
F
(
z
)
=
∑
n
=
0
∞
(
∑
k
=
0
∞
(
n
+
k
m
+
2
k
)
(
2
k
k
)
(
−
1
)
k
k
+
1
)
z
n
.
{\displaystyle F(z)=\sum _{n=0}^{\infty }{\left(\sum _{k=0}^{\infty }{{\binom {n+k}{m+2k}}{\binom {2k}{k}}{\frac {(-1)^{k}}{k+1}}}\right)}z^{n}\,.}
合計を交換(「スネークオイル」)すると、
F
(
z
)
=
∑
k
=
0
∞
(
2
k
k
)
(
−
1
)
k
k
+
1
z
−
k
∑
n
=
0
∞
(
n
+
k
m
+
2
k
)
z
n
+
k
.
{\displaystyle F(z)=\sum _{k=0}^{\infty }{{\binom {2k}{k}}{\frac {(-1)^{k}}{k+1}}z^{-k}}\sum _{n=0}^{\infty }{{\binom {n+k}{m+2k}}z^{n+k}}\,.}
内部の合計は 2m + 2k の / (1 − z ) m + 2 k + 1 . したがって
F
(
z
)
=
z
m
(
1
−
z
)
m
+
1
∑
k
=
0
∞
1
k
+
1
(
2
k
k
)
(
−
z
(
1
−
z
)
2
)
k
=
z
m
(
1
−
z
)
m
+
1
∑
k
=
0
∞
C
k
(
−
z
(
1
−
z
)
2
)
k
where
C
k
=
k
th Catalan number
=
z
m
(
1
−
z
)
m
+
1
1
−
1
+
4
z
(
1
−
z
)
2
−
2
z
(
1
−
z
)
2
=
−
z
m
−
1
2
(
1
−
z
)
m
−
1
(
1
−
1
+
z
1
−
z
)
=
z
m
(
1
−
z
)
m
=
z
z
m
−
1
(
1
−
z
)
m
.
{\displaystyle {\begin{aligned}F(z)&={\frac {z^{m}}{(1-z)^{m+1}}}\sum _{k=0}^{\infty }{{\frac {1}{k+1}}{\binom {2k}{k}}\left({\frac {-z}{(1-z)^{2}}}\right)^{k}}\\[4px]&={\frac {z^{m}}{(1-z)^{m+1}}}\sum _{k=0}^{\infty }{C_{k}\left({\frac {-z}{(1-z)^{2}}}\right)^{k}}&{\text{where }}C_{k}=k{\text{th Catalan number}}\\[4px]&={\frac {z^{m}}{(1-z)^{m+1}}}{\frac {1-{\sqrt {1+{\frac {4z}{(1-z)^{2}}}}}}{\frac {-2z}{(1-z)^{2}}}}\\[4px]&={\frac {-z^{m-1}}{2(1-z)^{m-1}}}\left(1-{\frac {1+z}{1-z}}\right)\\[4px]&={\frac {z^{m}}{(1-z)^{m}}}=z{\frac {z^{m-1}}{(1-z)^{m}}}\,.\end{aligned}}}
そして、
s
n
=
{
(
n
−
1
m
−
1
)
for
m
≥
1
,
[
n
=
0
]
for
m
=
0
.
{\displaystyle s_{n}={\begin{cases}\displaystyle {\binom {n-1}{m-1}}&{\text{for }}m\geq 1\,,\\{}[n=0]&{\text{for }}m=0\,.\end{cases}}}
同じ方法をもう一度合計に適用すると分かりやすいが、今回は nの代わりに mを 自由パラメータとして取る 。つまり、
G
(
z
)
=
∑
m
=
0
∞
(
∑
k
=
0
∞
(
n
+
k
m
+
2
k
)
(
2
k
k
)
(
−
1
)
k
k
+
1
)
z
m
.
{\displaystyle G(z)=\sum _{m=0}^{\infty }\left(\sum _{k=0}^{\infty }{\binom {n+k}{m+2k}}{\binom {2k}{k}}{\frac {(-1)^{k}}{k+1}}\right)z^{m}\,.}
合計を交換(「スネークオイル」)すると、
G
(
z
)
=
∑
k
=
0
∞
(
2
k
k
)
(
−
1
)
k
k
+
1
z
−
2
k
∑
m
=
0
∞
(
n
+
k
m
+
2
k
)
z
m
+
2
k
.
{\displaystyle G(z)=\sum _{k=0}^{\infty }{\binom {2k}{k}}{\frac {(-1)^{k}}{k+1}}z^{-2k}\sum _{m=0}^{\infty }{\binom {n+k}{m+2k}}z^{m+2k}\,.}
すると内部の合計は (1 + z ) n + k となる。したがって
G
(
z
)
=
(
1
+
z
)
n
∑
k
=
0
∞
1
k
+
1
(
2
k
k
)
(
−
(
1
+
z
)
z
2
)
k
=
(
1
+
z
)
n
∑
k
=
0
∞
C
k
(
−
(
1
+
z
)
z
2
)
k
where
C
k
=
k
th Catalan number
=
(
1
+
z
)
n
1
−
1
+
4
(
1
+
z
)
z
2
−
2
(
1
+
z
)
z
2
=
(
1
+
z
)
n
z
2
−
z
z
2
+
4
+
4
z
−
2
(
1
+
z
)
=
(
1
+
z
)
n
z
2
−
z
(
z
+
2
)
−
2
(
1
+
z
)
=
(
1
+
z
)
n
−
2
z
−
2
(
1
+
z
)
=
z
(
1
+
z
)
n
−
1
.
{\displaystyle {\begin{aligned}G(z)&=(1+z)^{n}\sum _{k=0}^{\infty }{\frac {1}{k+1}}{\binom {2k}{k}}\left({\frac {-(1+z)}{z^{2}}}\right)^{k}\\[4px]&=(1+z)^{n}\sum _{k=0}^{\infty }C_{k}\,\left({\frac {-(1+z)}{z^{2}}}\right)^{k}&{\text{where }}C_{k}=k{\text{th Catalan number}}\\[4px]&=(1+z)^{n}\,{\frac {1-{\sqrt {1+{\frac {4(1+z)}{z^{2}}}}}}{\frac {-2(1+z)}{z^{2}}}}\\[4px]&=(1+z)^{n}\,{\frac {z^{2}-z{\sqrt {z^{2}+4+4z}}}{-2(1+z)}}\\[4px]&=(1+z)^{n}\,{\frac {z^{2}-z(z+2)}{-2(1+z)}}\\[4px]&=(1+z)^{n}\,{\frac {-2z}{-2(1+z)}}=z(1+z)^{n-1}\,.\end{aligned}}}
したがって、前と同じように、 m ≥ 1 の
場合も 得られます
。
s
n
=
[
z
m
]
z
(
1
+
z
)
n
−
1
=
[
z
m
−
1
]
(
1
+
z
)
n
−
1
=
(
n
−
1
m
−
1
)
,
{\displaystyle s_{n}=\left[z^{m}\right]z(1+z)^{n-1}=\left[z^{m-1}\right](1+z)^{n-1}={\binom {n-1}{m-1}}\,,}
生成関数は合同性を証明する
2 つの生成関数 (べき級数) が m を 法として合同であるとは、 A ( z ) ≡ B ( z ) (mod m ) と書かれるとき、それらの係数が すべての n ≥ 0に対して m を 法として合同である、すなわち、 整数 nのすべての関連するケースに対して a n ≡ b n (mod m )である (ここで m が 整数である と仮定する必要はないことに注意。たとえば、不定の x において多項式値になる可能性は十分にある)。「より単純な」右側の生成関数 B ( z ) が z の有理関数である場合 、この数列の形は、整数値 m ≥ 2 の特定のケースを法として最終的に周期的になることを示唆している 。
たとえば 、 オイラー 数 が 3 を法として次の合同性を満たすことを証明できる :
[ 25]
⟨
E
n
⟩
=
⟨
1
,
1
,
5
,
61
,
1385
,
…
⟩
⟼
⟨
1
,
1
,
2
,
1
,
2
,
1
,
2
,
…
⟩
(
mod
3
)
,
{\displaystyle \langle E_{n}\rangle =\langle 1,1,5,61,1385,\ldots \rangle \longmapsto \langle 1,1,2,1,2,1,2,\ldots \rangle {\pmod {3}}\,,}
∑
n
=
0
∞
E
n
z
n
=
1
−
z
2
1
+
z
2
(
mod
3
)
.
{\displaystyle \sum _{n=0}^{\infty }E_{n}z^{n}={\frac {1-z^{2}}{1+z^{2}}}{\pmod {3}}\,.}
任意の整数(つまり、素数 p k だけではない)を法とする特殊な生成関数によって列挙される数列の合同性を得るための 1 つの便利な方法は、上記のJ 分数による(収束しないものも含む)通常の生成関数の連分数表現のセクションで示されています。ランドーの Lectures on Generating Functions から、連分数による表現を通じて展開される生成級数に関連する特定の結果を 次のように引用します。
定理: 連分数の展開によって生成される級数の合同性 — 生成関数 A ( z ) が形式の
無限 連分数 で表され
、 A p ( z ) が、すべての 0 ≤ n < 2 pに対して a n = [ z n ] A p ( z ) と定義されるこの連分数展開の p 次収束を表すとします 。この場合、次のようになります。
A
(
z
)
=
1
1
−
c
1
z
−
p
1
z
2
1
−
c
2
z
−
p
2
z
2
1
−
c
3
z
−
⋱
{\displaystyle A(z)={\cfrac {1}{1-c_{1}z-{\cfrac {p_{1}z^{2}}{1-c_{2}z-{\cfrac {p_{2}z^{2}}{1-c_{3}z-{\ddots }}}}}}}}
関数 A p ( z ) は、すべての p ≥ 2 に対して有理数である。ここで、 p | p 1 、 p 1 p 2 、 p 1 p 2 p 3 の割り切れる条件の1つが満たされていると仮定する 。 つまり、 ある k ≥ 1に対して p | p 1 p 2 ⋯ p k である。
整数 pが 積 p 1 p 2 ⋯ p kを 割り切る場合、 A ( z ) ≡ A k ( z ) (mod p ) となります 。
生成関数は、その係数の合同性を証明する際にも使用されます。次に、 第一種スターリング数 と 分割関数 p ( n ) の特殊なケースの合同性を導出する 2 つの具体的な例を示します。これは、整数列 を含む問題に取り組む際の生成関数の汎用性を示しています 。
小さな整数を法とするスターリング数
有限積によって生成されるスターリング数に関する
主要 記事
S
n
(
x
)
:=
∑
k
=
0
n
[
n
k
]
x
k
=
x
(
x
+
1
)
(
x
+
2
)
⋯
(
x
+
n
−
1
)
,
n
≥
1
,
{\displaystyle S_{n}(x):=\sum _{k=0}^{n}{\begin{bmatrix}n\\k\end{bmatrix}}x^{k}=x(x+1)(x+2)\cdots (x+n-1)\,,\quad n\geq 1\,,}
ウィルフの標準的な参考文献Generatingfunctionology のセクション4.6にあるように、これらの数の合同性の概要を、それらの生成関数の特性から厳密に導出する 。基本的な議論を繰り返すと、が2を法として減少する場合、これらの有限積生成関数はそれぞれを満たすことがわかります。
S
n
(
x
)
=
[
x
(
x
+
1
)
]
⋅
[
x
(
x
+
1
)
]
⋯
=
x
⌈
n
2
⌉
(
x
+
1
)
⌊
n
2
⌋
,
{\displaystyle S_{n}(x)=[x(x+1)]\cdot [x(x+1)]\cdots =x^{\left\lceil {\frac {n}{2}}\right\rceil }(x+1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }\,,}
これは、これらのスターリング数 の偶奇性が 二項係数の偶奇性と一致すること
を意味する。
[
n
k
]
≡
(
⌊
n
2
⌋
k
−
⌈
n
2
⌉
)
(
mod
2
)
,
{\displaystyle {\begin{bmatrix}n\\k\end{bmatrix}}\equiv {\binom {\left\lfloor {\frac {n}{2}}\right\rfloor }{k-\left\lceil {\frac {n}{2}}\right\rceil }}{\pmod {2}}\,,}
そしてその結果、 [ nk ] は k < ⌊ のときは常に偶数である ん / 2 ⌋ .
同様に、スターリング数生成関数を定義する右辺の積を3を法として簡約すると、少し複雑な式が得られます。
[
n
m
]
≡
[
x
m
]
(
x
⌈
n
3
⌉
(
x
+
1
)
⌈
n
−
1
3
⌉
(
x
+
2
)
⌊
n
3
⌋
)
(
mod
3
)
≡
∑
k
=
0
m
(
⌈
n
−
1
3
⌉
k
)
(
⌊
n
3
⌋
m
−
k
−
⌈
n
3
⌉
)
×
2
⌈
n
3
⌉
+
⌊
n
3
⌋
−
(
m
−
k
)
(
mod
3
)
.
{\displaystyle {\begin{aligned}{\begin{bmatrix}n\\m\end{bmatrix}}&\equiv [x^{m}]\left(x^{\left\lceil {\frac {n}{3}}\right\rceil }(x+1)^{\left\lceil {\frac {n-1}{3}}\right\rceil }(x+2)^{\left\lfloor {\frac {n}{3}}\right\rfloor }\right)&&{\pmod {3}}\\&\equiv \sum _{k=0}^{m}{\begin{pmatrix}\left\lceil {\frac {n-1}{3}}\right\rceil \\k\end{pmatrix}}{\begin{pmatrix}\left\lfloor {\frac {n}{3}}\right\rfloor \\m-k-\left\lceil {\frac {n}{3}}\right\rceil \end{pmatrix}}\times 2^{\left\lceil {\frac {n}{3}}\right\rceil +\left\lfloor {\frac {n}{3}}\right\rfloor -(m-k)}&&{\pmod {3}}\,.\end{aligned}}}
分割関数の合同性
この例では、無限積の冪級数展開が多くの特殊関数の展開を生成する仕組みの一部を取り入れ、分割関数を列挙する。特に、 分割 関数 p ( n )は 、次式で与えられる逆無限 q- ポッホハマー記号 積(または 場合によっては
z-ポッホハマー積)によって生成されることを思い出す。
∑
n
=
0
∞
p
(
n
)
z
n
=
1
(
1
−
z
)
(
1
−
z
2
)
(
1
−
z
3
)
⋯
=
1
+
z
+
2
z
2
+
3
z
3
+
5
z
4
+
7
z
5
+
11
z
6
+
⋯
.
{\displaystyle {\begin{aligned}\sum _{n=0}^{\infty }p(n)z^{n}&={\frac {1}{\left(1-z\right)\left(1-z^{2}\right)\left(1-z^{3}\right)\cdots }}\\[4pt]&=1+z+2z^{2}+3z^{3}+5z^{4}+7z^{5}+11z^{6}+\cdots .\end{aligned}}}
この分割関数は多くの既知の 合同性を 満たしており、その中には次のような結果も含まれるが、この関数の関連する整数合同の形式についてはまだ多くの未解決の疑問が残っている。 [26]
p
(
5
m
+
4
)
≡
0
(
mod
5
)
p
(
7
m
+
5
)
≡
0
(
mod
7
)
p
(
11
m
+
6
)
≡
0
(
mod
11
)
p
(
25
m
+
24
)
≡
0
(
mod
5
2
)
.
{\displaystyle {\begin{aligned}p(5m+4)&\equiv 0{\pmod {5}}\\p(7m+5)&\equiv 0{\pmod {7}}\\p(11m+6)&\equiv 0{\pmod {11}}\\p(25m+24)&\equiv 0{\pmod {5^{2}}}\,.\end{aligned}}}
生成関数と合同式の操作を形式的な冪級数に対して使用して、上に挙げた合同式の最初のものの非常に初歩的な証明を行う方法を示します。
まず、二項係数生成関数では、
1、 z 5 、 z 10 、…
のべき乗に対応するものを除いて、すべての係数が5で割り切れることに気づきます。 さらに、それらの場合、係数の余りは5を法とした1になります。したがって、
または同等に、
次の
ようになります。
1
(
1
−
z
)
5
=
∑
i
=
0
∞
(
4
+
i
4
)
z
i
,
{\displaystyle {\frac {1}{(1-z)^{5}}}=\sum _{i=0}^{\infty }{\binom {4+i}{4}}z^{i}\,,}
1
(
1
−
z
)
5
≡
1
1
−
z
5
(
mod
5
)
,
{\displaystyle {\frac {1}{(1-z)^{5}}}\equiv {\frac {1}{1-z^{5}}}{\pmod {5}}\,,}
1
−
z
5
(
1
−
z
)
5
≡
1
(
mod
5
)
.
{\displaystyle {\frac {1-z^{5}}{(1-z)^{5}}}\equiv 1{\pmod {5}}\,.}
(
1
−
z
5
)
(
1
−
z
10
)
(
1
−
z
15
)
⋯
(
(
1
−
z
)
(
1
−
z
2
)
(
1
−
z
3
)
⋯
)
5
≡
1
(
mod
5
)
.
{\displaystyle {\frac {\left(1-z^{5}\right)\left(1-z^{10}\right)\left(1-z^{15}\right)\cdots }{\left((1-z)\left(1-z^{2}\right)\left(1-z^{3}\right)\cdots \right)^{5}}}\equiv 1{\pmod {5}}\,.}
これを無限積展開すると、
z · ((1 − z )(1 − z 2 )⋯) 4 における z 5 m + 5 の係数は、すべてのm に対して 5 で割り切れること
がわかります 。 [27] 最後に、前の式の z 5 m + 5 の係数を等しくすることで
、 すべて の m
≥ 0 に対して p (5 m + 4) ≡ 0 (mod 5) で あることを証明できます 。
z
⋅
(
1
−
z
5
)
(
1
−
z
10
)
⋯
(
1
−
z
)
(
1
−
z
2
)
⋯
=
z
⋅
(
(
1
−
z
)
(
1
−
z
2
)
⋯
)
4
×
(
1
−
z
5
)
(
1
−
z
10
)
⋯
(
(
1
−
z
)
(
1
−
z
2
)
⋯
)
5
,
{\displaystyle z\cdot {\frac {\left(1-z^{5}\right)\left(1-z^{10}\right)\cdots }{\left(1-z\right)\left(1-z^{2}\right)\cdots }}=z\cdot \left((1-z)\left(1-z^{2}\right)\cdots \right)^{4}\times {\frac {\left(1-z^{5}\right)\left(1-z^{10}\right)\cdots }{\left(\left(1-z\right)\left(1-z^{2}\right)\cdots \right)^{5}}}\,,}
∑
n
=
1
∞
p
(
n
−
1
)
z
n
=
z
(
1
−
z
)
(
1
−
z
2
)
⋯
=
z
⋅
(
1
−
z
5
)
(
1
−
z
10
)
⋯
(
1
−
z
)
(
1
−
z
2
)
⋯
×
(
1
+
z
5
+
z
10
+
⋯
)
(
1
+
z
10
+
z
20
+
⋯
)
⋯
{\displaystyle {\begin{aligned}\sum _{n=1}^{\infty }p(n-1)z^{n}&={\frac {z}{(1-z)\left(1-z^{2}\right)\cdots }}\\[6px]&=z\cdot {\frac {\left(1-z^{5}\right)\left(1-z^{10}\right)\cdots }{(1-z)\left(1-z^{2}\right)\cdots }}\times \left(1+z^{5}+z^{10}+\cdots \right)\left(1+z^{10}+z^{20}+\cdots \right)\cdots \end{aligned}}}
生成関数の変換には、他の用途を提供するものがいくつかあります ( メイン記事 を参照)。シーケンスの 通常の生成関数 (OGF) の変換は、1 つのシーケンスの生成関数を別のシーケンスを列挙する生成関数に変換する方法を提供します。これらの変換には通常、シーケンス OGF を含む積分式 ( 積分変換 を参照) またはこれらの関数の高次導関数の加重和 ( 微分変換を 参照) が含まれます。
生成関数変換は、和の生成関数を表現しようとするときに役立ちます。
s
n
:=
∑
m
=
0
n
(
n
m
)
C
n
,
m
a
m
,
{\displaystyle s_{n}:=\sum _{m=0}^{n}{\binom {n}{m}}C_{n,m}a_{m},}
S ( z ) = g ( z ) A ( f ( z )) の形式で、 元の流れの生成関数を含みます。たとえば、合計がの場合
、修正された合計式の生成関数は [28]
で与えられます( 二項変換 と スターリング変換
も参照 )。
s
n
:=
∑
k
=
0
∞
(
n
+
k
m
+
2
k
)
a
k
{\displaystyle s_{n}:=\sum _{k=0}^{\infty }{\binom {n+k}{m+2k}}a_{k}\,}
S
(
z
)
=
z
m
(
1
−
z
)
m
+
1
A
(
z
(
1
−
z
)
2
)
{\displaystyle S(z)={\frac {z^{m}}{(1-z)^{m+1}}}A\left({\frac {z}{(1-z)^{2}}}\right)}
シーケンスのOGF、 F ( z ) とその指数生成関数、またはEGF、 F̂ ( z ) との間の変換を行う積分式もあり 、その逆も同様です。
F
(
z
)
=
∫
0
∞
F
^
(
t
z
)
e
−
t
d
t
,
F
^
(
z
)
=
1
2
π
∫
−
π
π
F
(
z
e
−
i
ϑ
)
e
e
i
ϑ
d
ϑ
,
{\displaystyle {\begin{aligned}F(z)&=\int _{0}^{\infty }{\hat {F}}(tz)e^{-t}\,dt\,,\\[4px]{\hat {F}}(z)&={\frac {1}{2\pi }}\int _{-\pi }^{\pi }F\left(ze^{-i\vartheta }\right)e^{e^{i\vartheta }}\,d\vartheta \,,\end{aligned}}}
ただし、これらの積分は適切な z の値に対して収束するものとする。
特殊生成関数の表
特殊な数学的級数の最初のリストは ここに あります。いくつかの有用で特殊な数列生成関数は、 Concrete Mathematicsのセクション5.4と7.4、およびWilfの Generatingfunctionology のセクション2.5に記載さ れています。その他の注目すべき特殊な生成関数には、次の表のエントリが含まれますが、これは決して完全ではありません。 [29]
参照
注記
^ちなみに、 m < 0 の
場合の対応する式は次のように表される。
∑
n
=
0
∞
g
n
+
m
z
n
=
G
(
z
)
−
g
0
−
g
1
z
−
⋯
−
g
m
−
1
z
m
−
1
z
m
.
{\displaystyle \sum _{n=0}^{\infty }g_{n+m}z^{n}={\frac {G(z)-g_{0}-g_{1}z-\cdots -g_{m-1}z^{m-1}}{z^{m}}}\,.}
参考文献
^ この代替用語は、EN Gilbert (1956)「Enumeration of Labeled graphs」、 Canadian Journal of Mathematics 3、p. 405–411 にすでに記載されていますが、2000 年以前には使用がまれであり、それ以降は使用が増えているようです。
^ Knuth, Donald E. (1997). 「§1.2.9 生成関数」. 基本的なアルゴリズム . コンピュータプログラミングの芸術 . 第 1 巻 (第 3 版). Addison-Wesley. ISBN 0-201-89683-4 。
^ フラジョレット&セジウィック 2009、95ページ
^ 「ランバート級数の恒等式」。Math Overflow 。2017年。
^ アポストル、トム・M. (1976)、 解析的数論入門 、数学の学部テキスト、ニューヨーク-ハイデルベルク:シュプリンガー・フェアラーク、 ISBN 978-0-387-90163-3 、 MR 0434929、 Zbl 0335.10001 42~43ページ
^ ウィルフ 1994、56 ページ
^ ウィルフ 1994、59 ページ
^ ハーディ、GH; ライト、EM; ヒースブラウン、DR; シルバーマン、JH (2008)。 数論入門 (第6版)。オックスフォード大学出版局。p. 339。ISBN 9780199219858 。
^ Knuth, DE (1992). 「畳み込み多項式」. Mathematica J. 2 : 67–78. arXiv : math/9207221 . Bibcode :1992math......7221K.
^ Spivey, Michael Z. (2007). 「組み合わせ和と有限差」. 離散数学 . 307 (24): 3130–3146. doi : 10.1016/j.disc.2007.03.052 . MR 2370116.
^ Mathar, RJ (2012). 「さらにもう一つの積分表」. arXiv : 1207.5845 [math.CA]. v4 等式 (0.4)
^ Graham、Knuth、Patashnik 1994、§6.1の表265、スターリング数三角形を含む有限和恒等式。
^ ランド 2003、§2.4
^ スタンレー、リチャード P.、フォミン、セルゲイ (1997)の例 。「§6.3」。列挙的組合せ論: 第 2 巻。ケンブリッジ高等数学研究。第 62 巻。ケンブリッジ大学出版局 。ISBN 978-0-521-78987-5 。
^ クヌース 1997, §1.2.9
^ Graham、Knuth、Patashnik 1994、p. 569、演習 7.36 の解答
^ フラジョレット&セジウィック 2009、§B.4
^ Schneider, C. (2007). 「記号的総和が組合せ論を支援する」 Sém. Lothar. Combin . 56 : 1–36.
^ これらの用語の使用法については、Graham、Knuth、Patashnik 1994、§7.4 の特殊シーケンス生成関数に関するセクションを参照してください。
^ Good, IJ (1986). 「対称ディリクレ分布とその混合分布の分割表への応用について」 Annals of Statistics . 4 (6): 1159–1189. doi : 10.1214/aos/1176343649 .
^ J 分数の特性に関するより詳しい情報については 以下を参照してください。
Flajolet, P. (1980). 「連分数の組合せ的側面」 (PDF) . 離散数学 . 32 (2): 125–161. doi :10.1016/0012-365X(80)90050-3.
ウォール、HS(2018)[1948]。連分数の解析理論。ドーバー 。ISBN 978-0-486-83044-5 。
^ 以下の記事を参照してください。
Schmidt, Maxie D. (2016). 「平方級数生成関数の連分数」 arXiv : 1612.02778 [math.NT].
— (2017). 「一般 化 階乗関数の通常生成関数のヤコビ型連分数」。 整数 シーケンスジャーナル 。20。arXiv : 1610.09691。17.3.4 。
— (2017). 「 h ≥ 2の整数を法とする二項係数のヤコビ型連分数と合同式 」. arXiv : 1702.01374 [math.CO].
^ グラハム、クヌース、パタシュニク、1994 年、§8.3
^ Graham、Knuth、Patashnik 1994、§7.3 の例 6 に、別の方法と、生成関数を使用したこの問題の完全なセットアップが記載されています。このより「複雑な」アプローチは、同じ参考文献のセクション 7.5 に記載されています。
^ ランド 2003、§5
^ ハーディ他 2008, §19.12
^ ハーディ、GH; ライト、EM 『数論入門』 。 p.288、Th.361
^ グラハム、クヌース、パタシュニック 1994、p. 535、演習 5.71
^ Plouffe, Simon (1992) にある 1031 の生成関数 も参照してください。 近似関数といくつかの推測 [ 母関数の近似といくつかの推測 ] (修士) (フランス語)。ケベック大学モントリオール。 arXiv : 0911.4975 。
引用
外部リンク