Transformation of a mathematical sequence
組合せ論 において 、 二項変換は、 順方向差分 を計算する 列変換(すなわち、 列 の変換) である。これは、 通常の生成関数 に関連付けられた列に二項変換を適用した結果である オイラー変換 と密接に関連している 。
意味 シーケンス { a n } の二 項変換 T は、 シーケンス { s n } が次のように定義されるもの
である。
s n = ∑ k = 0 n ( − 1 ) k ( n k ) a k . {\displaystyle s_{n}=\sum _{k=0}^{n}(-1)^{k}{\binom {n}{k}}a_{k}.}
正式には次のように書くことができる。
s n = ( T a ) n = ∑ k = 0 n T n k a k {\displaystyle s_{n}=(Ta)_{n}=\sum _{k=0}^{n}T_{nk}a_{k}}
変換は、 Tが 行列要素 T nk を持つ無限次元 演算子 であるときに行われる。この変換は 反転 であり、つまり
T T = 1 {\displaystyle TT=1}
または、インデックス表記を使用すると、
∑ k = 0 ∞ T n k T k m = δ n m {\displaystyle \sum _{k=0}^{\infty }T_{nk}T_{km}=\delta _{nm}}
クロネッカーのデルタ はどこに ありますか ?元の級数は次のようにして復元できます。 δ n m {\displaystyle \delta _{nm}}
a n = ∑ k = 0 n ( − 1 ) k ( n k ) s k . {\displaystyle a_{n}=\sum _{k=0}^{n}(-1)^{k}{\binom {n}{k}}s_{k}.}
数列の二項変換は 数列の n 番目の 前方差分であり、奇数差分には負の符号が付きます。つまり、次のようになります。
s 0 = a 0 s 1 = − ( Δ a ) 0 = − a 1 + a 0 s 2 = ( Δ 2 a ) 0 = − ( − a 2 + a 1 ) + ( − a 1 + a 0 ) = a 2 − 2 a 1 + a 0 ⋮ s n = ( − 1 ) n ( Δ n a ) 0 {\displaystyle {\begin{aligned}s_{0}&=a_{0}\\s_{1}&=-(\Delta a)_{0}=-a_{1}+a_{0}\\s_{2}&=(\Delta ^{2}a)_{0}=-(-a_{2}+a_{1})+(-a_{1}+a_{0})=a_{2}-2a_{1}+a_{0}\\&\;\;\vdots \\s_{n}&=(-1)^{n}(\Delta ^{n}a)_{0}\end{aligned}}}
ここで Δは 前進差分演算子 です 。
一部の著者は、二項変換を余分な符号で定義し、それが自己逆にならないようにしています。
t n = ∑ k = 0 n ( − 1 ) n − k ( n k ) a k {\displaystyle t_{n}=\sum _{k=0}^{n}(-1)^{n-k}{\binom {n}{k}}a_{k}}
その逆は
a n = ∑ k = 0 n ( n k ) t k . {\displaystyle a_{n}=\sum _{k=0}^{n}{\binom {n}{k}}t_{k}.}
この場合、前者の変換は 逆二項変換 と呼ばれ、後者は単に 二項変換 と呼ばれます。これは、例えば 「On-Line Encyclopedia of Integer Sequences」 などにおける標準的な用法です。
例 二項変換の両方のバージョンは差分表に現れます。次の差分表を考えてみましょう。
0 1 10 63 324 1485 1 9 53 261 1161 8 44 208 900 36 164 692 128 528 400
各行は前の行の差分です。( m 行目 の n番目の数は a m , n = 3 n −2 (2 m +1 n 2 + 2 m (1+6 m ) n + 2 m -1 9 m 2 )であり、差分方程式 a m +1, n = a m , n +1 - a m , n が成り立ちます。)
一番上の行を左から右に読むと、{ a n } = 0、1、10、63、324、1485、... となります。同じ開始点 0 を持つ対角線は、{ t n } = 0、1、8、36、128、400、... となります。{ t n } は、{ a n }の非反転二項変換です 。
一番上の行を右から左に読むと、{ b n } = 1485, 324, 63, 10, 1, 0, ... となります。同じ開始点 1485 を持つ交差対角線は、{ s n } = 1485, 1161, 900, 692, 528, 400, ... となります。{ s n } は、{ b n }の逆二項変換です 。
通常の生成関数 変換は 級数に関連付けられた 生成関数を結び付ける。 通常の生成関数 については、
f ( x ) = ∑ n = 0 ∞ a n x n {\displaystyle f(x)=\sum _{n=0}^{\infty }a_{n}x^{n}}
そして
g ( x ) = ∑ n = 0 ∞ s n x n {\displaystyle g(x)=\sum _{n=0}^{\infty }s_{n}x^{n}}
それから
g ( x ) = ( T f ) ( x ) = 1 1 − x f ( − x 1 − x ) . {\displaystyle g(x)=(Tf)(x)={\frac {1}{1-x}}f{\left({\frac {-x}{1-x}}\right)}.}
通常の母関数間の関係は、 オイラー変換 と呼ばれることもあります。これは通常、2つの異なる方法で現れます。1つの形は、 交代級数 の 収束を加速する ために使用されます。つまり、次の恒等式が成り立ちます。
∑ n = 0 ∞ ( − 1 ) n a n = ∑ n = 0 ∞ ( − 1 ) n ( Δ n a ) 0 2 n + 1 {\displaystyle \sum _{n=0}^{\infty }{\left(-1\right)}^{n}a_{n}=\sum _{n=0}^{\infty }{\left(-1\right)}^{n}{\frac {(\Delta ^{n}a)_{0}}{2^{n+1}}}}
これは、上記の最後の式にx = 1/2 を 代入することで得られます 。右辺の項は通常、はるかに小さくなり、はるかに急速に減少するため、数値積分を迅速に行うことができます。
オイラー変換は一般化できます (Borisov B. および Shkodrov V.、2007)。
∑ n = 0 ∞ ( − 1 ) n ( n + p n ) a n = ∑ n = 0 ∞ ( − 1 ) n ( n + p n ) ( Δ n a ) 0 2 n + p + 1 , {\displaystyle \sum _{n=0}^{\infty }{\left(-1\right)}^{n}{\binom {n+p}{n}}a_{n}=\sum _{n=0}^{\infty }{\left(-1\right)}^{n}{\binom {n+p}{n}}{\frac {(\Delta ^{n}a)_{0}}{2^{n+p+1}}},}
ここで p = 0, 1, 2,... 。
オイラー変換は オイラー超幾何積分 にも頻繁に適用されます。ここで、オイラー変換は次の形をとります。 2 F 1 {\displaystyle \,_{2}F_{1}}
2 F 1 ( a , b ; c ; z ) = ( 1 − z ) − b 2 F 1 ( c − a , b ; c ; z z − 1 ) . {\displaystyle \,_{2}F_{1}(a,b;c;z)=(1-z)^{-b}\,_{2}F_{1}\left(c-a,b;c;{\frac {z}{z-1}}\right).}
[他の超幾何級数への一般化については [1] を参照]
二項変換、およびそのオイラー変換としての変形は、 数の 連分数表現 との関連で注目される。連分数表現が 0 < x < 1 {\displaystyle 0<x<1}
x = [ 0 ; a 1 , a 2 , a 3 , ⋯ ] {\displaystyle x=[0;a_{1},a_{2},a_{3},\cdots ]}
それから
x 1 − x = [ 0 ; a 1 − 1 , a 2 , a 3 , ⋯ ] {\displaystyle {\frac {x}{1-x}}=[0;a_{1}-1,a_{2},a_{3},\cdots ]}
そして
x 1 + x = [ 0 ; a 1 + 1 , a 2 , a 3 , ⋯ ] . {\displaystyle {\frac {x}{1+x}}=[0;a_{1}+1,a_{2},a_{3},\cdots ].}
指数生成関数 指数生成関数 については 、
f ¯ ( x ) = ∑ n = 0 ∞ a n x n n ! {\displaystyle {\overline {f}}(x)=\sum _{n=0}^{\infty }a_{n}{\frac {x^{n}}{n!}}}
そして
g ¯ ( x ) = ∑ n = 0 ∞ s n x n n ! {\displaystyle {\overline {g}}(x)=\sum _{n=0}^{\infty }s_{n}{\frac {x^{n}}{n!}}}
それから
g ¯ ( x ) = ( T f ¯ ) ( x ) = e x f ¯ ( − x ) . {\displaystyle {\overline {g}}(x)=(T{\overline {f}})(x)=e^{x}{\overline {f}}(-x).}
ボレル 変換は 、通常の生成関数を指数生成関数に変換します。
二項畳み込み と を 複素数列とします。これらの二項畳み込みは で定義されます。 この畳み込みは、RL Graham、DE Knuth、O. Patashnik 共著『 Concrete Mathematics : A Foundation for Computer Science』(Addison-Wesley、1989年)に記載されています。二項畳み込みが結合法則と可換法則を満たし、と で 定義される 列が 二項畳み込みのもとで恒等式となることは容易に理解できます。さらに、 の列 には 逆数が存在することも容易に理解できます。したがって、 の列の集合は二項畳み込みのもと で アーベル群 を形成します 。 ( a n ) n ∈ N {\displaystyle (a_{n})_{n\in \mathbb {N} }} ( b n ) n ∈ N {\displaystyle (b_{n})_{n\in \mathbb {N} }} ( a ∘ b ) n = ∑ k = 0 n ( n k ) a k b n − k , n = 0 , 1 , 2 , … {\displaystyle (a\circ b)_{n}=\sum _{k=0}^{n}{\binom {n}{k}}a_{k}b_{n-k},\ \ n=0,1,2,\ldots } { e n } {\displaystyle \{e_{n}\}} e 0 = 1 {\displaystyle e_{0}=1} e n = 0 {\displaystyle e_{n}=0} n = 1 , 2 , … , {\displaystyle n=1,2,\ldots ,} { a n } {\displaystyle \{a_{n}\}} a 0 ≠ 0 {\displaystyle a_{0}\neq 0} { a n } {\displaystyle \{a_{n}\}} a 0 ≠ 0 {\displaystyle a_{0}\neq 0}
二項畳み込みは指数関数の積から自然に生じます。実際、 ( ∑ n = 0 ∞ a n x n n ! ) ( ∑ n = 0 ∞ b n x n n ! ) = ∑ n = 0 ∞ ( a ∘ b ) n x n n ! . {\displaystyle \left(\sum _{n=0}^{\infty }a_{n}{\frac {x^{n}}{n!}}\right)\left(\sum _{n=0}^{\infty }b_{n}{\frac {x^{n}}{n!}}\right)=\sum _{n=0}^{\infty }(a\circ b)_{n}{\frac {x^{n}}{n!}}.}
二項変換は二項畳み込みを用いて表すことができます。 すべての に対してと とします 。すると、
は二項畳み込みの下で の逆数である ため、 この式は メビウス反転型の式として解釈できます 。 λ n = ( − 1 ) n {\displaystyle \lambda _{n}=(-1)^{n}} 1 n = 1 {\displaystyle 1_{n}=1} n {\displaystyle n} ( T a ) n = ( λ a ∘ 1 ) n . {\displaystyle (Ta)_{n}=(\lambda a\circ 1)_{n}.} t n = ∑ k = 0 n ( − 1 ) n − k ( n k ) a k ⟺ a n = ∑ k = 0 n ( n k ) t k {\displaystyle t_{n}=\sum _{k=0}^{n}{\left(-1\right)}^{n-k}{\binom {n}{k}}a_{k}\iff a_{n}=\sum _{k=0}^{n}{\binom {n}{k}}t_{k}} t n = ( a ∘ λ ) n ⟺ a n = ( t ∘ 1 ) n {\displaystyle t_{n}=(a\circ \lambda )_{n}\iff a_{n}=(t\circ 1)_{n}} λ n {\displaystyle \lambda _{n}} 1 n {\displaystyle 1_{n}}
数学文献には、別の二項畳み込みも存在します。算術関数 と の二項畳み込みは 、 と定義されます。 ここで 、 は正の整数 の標準因数分解 、は 二項係数 です 。この畳み込みは、PJ McCarthy (1986) の著書に掲載され、L. Toth と P. Haukkanen (2009) によってさらに研究されました。 f {\displaystyle f} g {\displaystyle g} ( f ∘ B g ) ( n ) = ∑ d ∣ n ( ∏ p ( ν p ( n ) ν p ( d ) ) ) f ( d ) g ( n / d ) , {\displaystyle (f\circ _{B}g)(n)=\sum _{d\mid n}\left(\prod _{p}{\binom {\nu _{p}(n)}{\nu _{p}(d)}}\right)f(d)g(n/d),} n = ∏ p p ν p ( n ) {\displaystyle n=\prod _{p}p^{\nu _{p}(n)}} n {\displaystyle n} ( ν p ( n ) ν p ( d ) ) {\displaystyle {\binom {\nu _{p}(n)}{\nu _{p}(d)}}}
積分表現 シーケンスを 複素解析 関数で補間できる場合、シーケンスの二項変換は補間関数上の Nörlund–Rice 積分 によって表すことができます。
一般化 プロディンガーは、関連する モジュラー的な 変換を提示している。
u n = ∑ k = 0 n ( n k ) a k ( − c ) n − k b k {\displaystyle u_{n}=\sum _{k=0}^{n}{\binom {n}{k}}a^{k}{\left(-c\right)}^{n-k}b_{k}}
与える
U ( x ) = 1 c x + 1 B ( a x c x + 1 ) {\displaystyle U(x)={\frac {1}{cx+1}}B{\left({\frac {ax}{cx+1}}\right)}}
ここで 、U と B はそれぞれ、 級数とに関連付けられた通常の生成関数です 。 { u n } {\displaystyle \{u_{n}\}} { b n } {\displaystyle \{b_{n}\}}
上昇 k 二項変換は次のように定義されることがある。
∑ j = 0 n ( n j ) j k a j . {\displaystyle \sum _{j=0}^{n}{\binom {n}{j}}j^{k}a_{j}.}
下降 k 二項変換は
∑ j = 0 n ( n j ) j n − k a j . {\displaystyle \sum _{j=0}^{n}{\binom {n}{j}}j^{n-k}a_{j}.}
どちらも、級数のハンケル変換の 核 の 準同型です 。
二項変換が次のように定義されている場合
∑ i = 0 n ( − 1 ) n − i ( n i ) a i = b n . {\displaystyle \sum _{i=0}^{n}{\left(-1\right)}^{n-i}{\binom {n}{i}}a_{i}=b_{n}.}
これを関数と等しくする J ( a ) n = b n . {\displaystyle {\mathfrak {J}}(a)_{n}=b_{n}.}
新しい 順方向差分 表を作成し、この表の各行の最初の要素を取って新しいシーケンスを形成すると 、元のシーケンスの2番目の二項変換は次のようになります。 { b n } {\displaystyle \{b_{n}\}}
J 2 ( a ) n = ∑ i = 0 n ( − 2 ) n − i ( n i ) a i . {\displaystyle {\mathfrak {J}}^{2}(a)_{n}=\sum _{i=0}^{n}(-2)^{n-i}{\binom {n}{i}}a_{i}.}
同じプロセスを k 回繰り返すと、
J k ( a ) n = b n = ∑ i = 0 n ( − k ) n − i ( n i ) a i . {\displaystyle {\mathfrak {J}}^{k}(a)_{n}=b_{n}=\sum _{i=0}^{n}(-k)^{n-i}{\binom {n}{i}}a_{i}.}
その逆は、
J − k ( b ) n = a n = ∑ i = 0 n k n − i ( n i ) b i . {\displaystyle {\mathfrak {J}}^{-k}(b)_{n}=a_{n}=\sum _{i=0}^{n}k^{n-i}{\binom {n}{i}}b_{i}.}
これを一般化すると、
J k ( a ) n = b n = ( E − k ) n a 0 {\displaystyle {\mathfrak {J}}^{k}(a)_{n}=b_{n}=(\mathbf {E} -k)^{n}a_{0}}
ここでは シフト演算子 です 。 E {\displaystyle \mathbf {E} }
その逆は
J − k ( b ) n = a n = ( E + k ) n b 0 . {\displaystyle {\mathfrak {J}}^{-k}(b)_{n}=a_{n}=(\mathbf {E} +k)^{n}b_{0}.}
参照
参考文献 ^ Miller, Allen R.; Paris, RB (2010). 「一般化超幾何関数のオイラー型変換」. Z. Angew. Math. Phys . 62 (1): 31– 45. doi :10.1007/s00033-010-0085-0. S2CID 30484300. ジョン・H・コンウェイとリチャード・K・ガイ、1996年、 『民数記』 Donald E. Knuth, The Art of Computer Programming Vol. 3 、(1973) Addison-Wesley、Reading、MA。 Helmut Prodinger、「二項変換についてのいくつかの情報」、The Fibonacci Quarterly 32 (1994)、412-415。 Spivey, Michael Z.; Steil, Laura L. (2006). 「k-二項変換とハンケル変換」. Journal of Integer Sequences . 9 : 06.1.1.書誌 コード :2006JIntS...9...11S. Borisov, B.; Shkodrov, V. (2007). 「一般化二項変換における発散級数」. 上級研究継続数学 . 14 (1): 77– 82. Khristo N. Boyadzhiev、 「二項変換に関する注記 、理論と表、およびスターリング変換の付録」(2018)、World Scientific。 RL Graham、DE Knuth、O. Patashnik:「具体的な数学:コンピュータサイエンスの基礎」、Addison-Wesley(1989年)。 PJ McCarthy, 算術関数入門、Springer-Verlag、1986年。 P. Haukkanen、算術関数の二項畳み込みについて、Nieuw Arch。ウィスク。 (IV) 14 (1996)、no. 2、209--216。 L. Toth と P. Haukkanen、「算術関数の二項畳み込みについて」、J. Combinatorics and Number Theory 1(2009)、31-48。 P. Haukkanen, 通常の母関数を用いたいくつかの二項反転. Publ. Math. Debr. 47, No. 1-2, 181-191 (1995).
外部リンク Wolfram MathWorldでの二項変換 OEIS wiki の二項変換