關於 「 Fibonacci 與 Padovan 的對話 及 F-P 卷積恆等式」 之迴響
廖信傑 · 薛昭雄
1. 前言
Fibonacci 數 Fn 和 Padovan 數 Pn 分別滿足遞迴關係式 Fn+2 = Fn+1+ Fn 及初始 條件 F0 = 0, F1 = 1 和遞迴式 Pn+3 = Pn+1+ Pn 及初始條件 P0 = P1 = P2 = 1 (見廖信 傑[9])。
在 「Fibonacci 與 Padovan 的對話(上)」[6] 及 「Fibonacci 與 Padovan 的對話(下):
F-P 卷積恆等式」[7] 與 [8] 中陳建燁先生證明了 Fn和 Pn 可以表成變數取特別值的完全齊次 對稱多項式
Fn= hn−1(α, β), Pn = hn+2(a, b, c), (1) 其中 α, β 及 a, b, c 分別為 x2− x − 1 = 0 及 x3− x − 1 = 0 的根。 並且利用 (1) 及其稱作
「自由分解重組恆等式」[5] 之完全齊次對稱多項式 hn 的性質, 得到 Fibonacci 數和 Padovan 數有一個漂亮的卷積(convolution) 公式:
n i=0
FiPn−i = Fn+3− Pn+3. (2) 本文的目的是從生成函數 (generating function) 的角度來看 Fn 和 Pn 可以用完全齊次對稱 多項式來表示這件事。 第二節我們先很概略的介紹兩類對稱多項式及之後會用到的性質, 接著 利用生成函數的方法證明陳建燁先生 [6] 的兩個結果; 第三節, 我們利用生成函數得到更一般 的 k 階線性遞迴數列同樣也可以用完全對稱多項式來表示; 最後一節, 我們用生成函數方法證 明卷積公式 (2), 並且也提供另一個有類似卷積公式的例子 : Fibonacci 數和 Tribonacci 數。
2. F
n, P
n用h
n表示
這一節我們會用生成函數方法證明 Fibonacci 數 Fn 和 Padovan 數 Pn 可完全齊次對 稱多項式 hn 表示, 為此我們先回顧何謂完全齊次對稱多項式, 有關一些基本的生成函數操作 讀者可以參考 Wilf [4]。
50
給定正整數n, 一個m個變數的完全齊次多項式定義作 hn(x1, x2, . . . , xm) =
1≤i1≤i2≤...≤in≤m
xi1xi2. . . xin,
例如 h2(x1, x2, x3) = x21+ x22+ x23+ x1x2+ x1x3+ x2x3。 由 hn 的定義我們很容易可以得 到其生成函數為
∞ n=0
hn(x1, . . . , xm)tn= 1
(1− x1t)(1− x2t) . . . (1− xmt), (3) 只要把等號右邊看成 m 個級數 ∞
k=0xkitk (i = 1, 2, . . . , m) 的乘積後比較等式兩邊 tn 係數就可以得到。 另一類很常出現的對稱多項式是基本對稱多項式 (elementary symmetric polynomial) en, 定義為
en(x1, x2, . . . , xm) =
1≤i1<i2<...<in≤m
xi1xi2. . . xin,
例如 e1(x1, x2, x3) = x1+ x2+ x3, e2(x1, x2, x3) = x1x2+ x1x3+ x2x3, 注意到當 n > m 時 en(x1, . . . , xm) = 0。 並且其生成函數為
∞ n=0
en(x1, . . . , xm)tn= (1 + x1t)(1 + x2t) . . . (1 + xnt), (4) 類似地只要比較等號兩邊 tn 係數馬上可以得到。 在對稱函數的理論中 en 和 hn 是關係較密切 的兩類, 這某種程度也反映在他們的生成函數上, 觀察 (3) 和 (4) 注意到若將 (4) 中 t 用 −t 代換後與 (3) 相乘可得
1 =
m i=1
(1− xit)
m i=1
1 1− xit =
∞
n=0
(−1)nentn
∞
n=0
hntn
= 1 +
∞ n=1
n
k=0
(−1)n−khken−k
tn 比較等號兩邊 tn 係數後得到對任意 n ≥ 1 都有
n k=0
(−1)kekhn−k = 0 或等價地 hn =
n i=1
(−1)k−1ekhn−k (5) 這個公式讓我們利用遞迴的方式得到 en 和 hn 之間的轉換, 以上是本文關於對稱多項式具備 的知識, 實際上以上內容並沒有觸及對稱函數真正有趣的部份, 對這方面內容有興趣的讀者可 以參考 Macdonald [2] 或 Stanley 與 Fomin [3]的第七章。
有了這些預備知識後, 現在回到我們關心的 Fibonacci 數和 Padovan 數上, 我們已經知 道 {Fn} 和 {Pn} 的遞迴式, 對這類遞迴式一般大學數學系離散數學或組合數學課程中有一套 標準的方法可以讓我們得到 {Fn} 和 {Pn} 的生成函數
n≥0
Fntn = t
1− t − t2,
n≥0
Pntn = 1 + t 1− t2− t3.
定理 2.1 ([6]). 設 α1, α2 為 {Fn} 之遞迴式 Fn+2 = Fn+1+ Fn 特徵多項式的兩根, 則對 正整數 n ≥ 1 有
Fn = hn(α1, α2).
證明: 假設 α1, α2 為 {Fn} 的特徵方程 t2 − t − 1 = 0 之兩根, 觀察 {Fn} 的生成函數注意 到分母可改寫 1 − t − t2 = t2(t12 −1t− 1) = (1 − α1t)(1− α2t), 將其與 hn 的生成函數 (3) 比較
n≥0
Fntn= t
(1− α1t)(1− α2t) = t
∞ n=0
hn(α1, α2)tn =
n≥0
hn(α1, α2)tn+1, 可以得到當 n ≥ 1 時有 Fn = hn−1(α1, α2)。 定理 2.2. 設 β1, β2, β3 為 {Pn} 之遞迴式 Pn+3 = Pn+1+ Pn 特徵多項式的三根, 則對正整 數 n ≥ 1 有
Pn= hn+2(β1, β2, β3).
證明: 類似的方法可將 {Pn} 的生成函數改寫並與 (3) 比較得
n≥0
Pntn= 1 + t
(1− β1t)(1− β2t)(1− β3t) = (1 + t)
n≥0
hn(β1, β2, β3)tn
=
n≥0
hntn+
n≥0
hntn+1,
故當 n ≥ 1 時有 Pn = hn(β1, β2, β3) + hn−1(β1, β2, β3), 接下來同 [6] 中的做法利用 (5) 式 可進一步化簡, 當 n ≥ 1 時有 Pn= hn+2(β1, β2, β3)。
3. 一般的 k 階線性遞迴數列用 h
n表示
有了前一節的 Fn和 Pn的兩個例子, 現在我們來看更一般的情況。 給定初始值 b0, b1, . . ., bk−1, 考慮數列 {bn} 滿足 k 階線性遞迴方程式
bn+ a1bn−1+· · · + ak−1bn−k+1+ akbn−k = 0 (ak= 0).
用一般我們求數列生成函數的方法可以得到 {bn} 的生成函數為
n≥0bntn= p(t) q(t), 其中 q(t) = 1 + a1t +· · · + aktk,
p(t) = b0+ (b1 + a1b0)t + (b2+ a1b1 + a2b0)t2+· · · +
bk−1+
k−1 i=1
aibk−1−i
tk−1. 同時 {bn} 之特徵多項式為
r(t) = tk+ a1tk−1+· · · + ak−1t + ak,
可以發現 q(t) 和 r(t) 之間滿足 q(t) = tkr(1t),我們稱 q(t) 為 r(t)之reflected polynomial。
假設 α1, α2, . . . , αk 為 r(t) 之 k 個根, 這邊 αi 不一定要相異, 則特徵多項式可分解為 r(t) = (t− α1)(t− α2) . . . (t− αk).
由 q(t) = tkr(1t)可知
q(t) = (1− α1t)(1− α2t) . . . (1− αkt), 從而可以得到
n≥0
bntn =
b0+(b1+a1b0)t+(b2+a1b1+a2b0)t2+· · · +
bk−1+k−1
i=1 aibk−1−i
tk−1 (1− α1t)(1− α2t) . . . (1− αkt)
=
b0+(b1+a1b0)t+· · ·+
bk−1+
k−1 i=1
aibk−1−i
tk−1
n≥0
hn(α1, . . . , αk)tn, 比較等號兩邊 tn 係數可知, 當 n ≥ k − 1 時有
bn= b0hn(α1, . . . , αk) + (b1 + a1b0)hn−1(α1, . . . , αk) +· · · +
bk−1+
k−1 i=1
aibk−1−i
hn−k+1(α1, . . . , αk).
將上面的結果整理後我們有以下定理
定理 3.1. 令 {bn} 為一數列滿足 k 階線性遞迴式
bn+ a1bn−1+· · · + ak−1bn−k+1+ akbn−k = 0 (ak = 0) 並且給定 b0, b1, . . . , bk−1, 則其生成函數為
n≥0
bntn=
b0+ (b1+a1b0)t+(b2+a1b1+a2b0)t2+. . . +
bk−1+k−1
i=1aibk−1−i
tk−1 (1− α1t)(1− α2t) . . . (1− αkt) ,
其中 α1, . . . , αk 為 tk+ a1tk−1+· · · + ak−1t + ak = 0 的 k 個根, 並且當 n ≥ k − 1 時有 bn= b0hn(α1, . . . , αk) + (b1+ a1b0)hn−1(α1, . . . , αk) +· · ·
+
bk−1+
k−1 i=1
aibk−1−i
hn−k+1(α1, . . . , αk).
所以當 k = 2 時, 二階線性遞迴數列與完全齊次對稱多項式有以下關係
推論 3.2. 令 {an} 滿足二階線性遞迴式 an+2 = ban+1+ can 並給定初始值 a0, a1, 則其生
成函數為
n≥0
antn = a0+ (a1− ba0)t (1− αt)(1 − βt) 其中 α, β 為 t2− bt − c = 0 之兩根, 並且對 n ≥ 1 之 an 有
an = a0hn(α, β) + (a1− ba0)hn−1(α, β).
類似地, 若滿足遞迴式的是多項式數列則有
推論 3.3. 令多項式數列 {an(x)} 滿足二階線性遞迴式 an+2(x) = b(x)an+1(x) + c(x)an(x) 並給定初始值 a0(x), a1(x), 則其生成函數為
n≥0
an(x)tn= a0(x) + (a1(x)− b(x)a0(x))t (1− α(x)t)(1 − β(x)t) 其中 α(x), β(x) 為 t2− b(x)t − c(x) = 0 之兩根, 並且對 n ≥ 1 有
an(x) = a0(x)hn(α(x), β(x)) + (a1(x)− b(x)a0(x))hn−1(α(x), β(x)).
在 k = 3 時則有
推論 3.4. 令 {an} 滿足三階線性遞迴式 an+3= ban+2+ can+1+ dan 並給定初始值 a0, a1, a2, 則其生成函數為
n≥0
antn= a0+ (a1− ba0)t + (a2− ba1− ca0)t2 (1− αt)(1 − βt)(1 − γt) 其中 α, β, γ 為 t3− bt2− ct − d = 0 之三根, 並且對 n ≥ 1 之 an 有
an= a0hn(α, β, γ) + (a1− ba0)hn−1(α, β, γ) + (a2− ba1 − ca0)hn−2(α, β, γ).
若滿足遞迴式的是多項式則有
推論 3.5. 令 {an(x)} 滿足三階線性遞迴式 an+3(x) = b(x)an+2(x) + c(x)an+1(x) + d(x)an(x) 並給定初始值 a0(x), a1(x), a2(x), 則其生成函數為
n≥0
an(x)tn= a0(x) + (a1(x)− b(x)a0(x))t + (a2(x)− b(x)a1(x)− c(x)a0(x))t2 (1− α(x)t)(1 − β(x)t)(1 − γ(x)t) , 其中 α(x), β(x), γ(x) 為 t3−b(x)t2−c(x)t−d(x)=0 之三根, 並且對 n≥2 之 an(x) 有
an(x) = a0(x)hn(α, β, γ) + (a1(x)− b(x)a0(x))hn−1(α, β, γ) +(a2(x)− b(x)a1(x)− c(x)a0(x))hn−2(α, β, γ).
4. 生成函數方法證明卷積公式
這一節我們會利用生成函數的方法重新證明 (2)。 已知 {Fn} 和 {Pn} 的遞迴式, 對這類 遞迴式一般大學數學系離散數學或組合數學課程中有一套標準的方法可以讓我們得到 {Fn} 和 {Pn} 的生成函數
n≥0
Fnxn = x
1− x − x2,
n≥0
Pnxn = 1 + x 1− x2− x3. 並且我們所求卷積的數列 {n
i=0FiPn−i}n≥0 其對應生成函數為兩個原生成函數的乘積
n≥0
n
i=0
PiFn−i
xn= x
1− x − x2 · 1 + x
1− x2− x3. (6) 假設 (6) 式等號右邊可以部份分式如下
x
1− x − x2 · 1 + x
1− x2− x3 = Ax + B
1− x − x2 + ax2+ bx + c 1− x2− x3 , 於是等號兩邊同乘 (1 − x − x2)(1− x2− x3) 後可得
x(1 + x) = (Ax + B)(1− x2− x3) + (ax2+ bx + c)(1− x − x2) 比較等號兩邊係數後我們知道 a = −1, b = c = −2, A = 1, B = 2。 因此 (6) 式成為
n≥0
n
i=0
PiFn−i
xn= 2 + x
1− x − x2 − 2 + 2x + x2
1− x2− x3. (7) 觀察 (7) 式等號右邊第一式可改寫為 x
1− x − x2 + 2
1− x − x2, 注意到 1
1− x − x2 為 {Fn+1}n≥0 之生成函數因為
1
1− x − x2 = 1 x·
x
1− x − x2 − F0
=
n≥0
Fn+1xn,
因此 2 + x
1− x − x2 為 {Fn+ 2Fn+1 = Fn+2+ Fn+1= Fn+3}n≥0 之生成函數。
現在觀察 (7) 式等號右邊第二項, 注意到 {Pn+1}n≥0 之生成函數為
n≥0
Pn+1xn= 1 x·
1 + x
1− x2− x3 − P0
= 1 + x + x2 1− x2 − x3, 所以 2 + 2x + x2
1− x2− x3 可改寫為 1 + x
1− x2− x3 + 1 + x + x2
1− x2− x3 即為 {Pn+ Pn+1 = Pn+3}n≥0 之 生成函數。 因此我們得到
n i=0
FiPn−i = Fn+3− Pn+3. 以下我們給出另一個例子具有和Fn及Pn類似的卷積公式。
首先我們定義 Tribonacci 數 Tn [1] 滿足下列遞迴多項式 Tn+3= Tn+2+ Tn+1+ Tn
及初始條件 T0 = 0, T1 = T2 = 1, 容易可以得到 {Tn}n≥0 的生成函數為
n≥0
Tnxn= x
1− x − x2− x3. 考慮 {Fn}n≥0 和 {Tn}n≥0 的卷積 {n
i=0FiTn−i}n≥0 對應的生成函數, 一樣利用部分分式
法有 x
1− x − x2 · x
1− x − x2 − x3 =− 1 + x
1− x − x2 + 1 + x + x2 1− x − x2− x3 其中由之前的討論知道 1 + x
1− x − x2 = 1
1− x − x2+ x
1− x − x2 為 {Fn+Fn+1 = Fn+2}n≥0
之生成函數; 另一方面注意到 1 + x + x2
1− x − x2− x3 = 1 x2
x
1− x − x2− x3 − T0− T1x
=
n≥0
Tn+2xn,
所以 1 + x + x2
1− x − x2 − x3 為 {Tn+2}n≥0 之生成函數, 於是由這兩個結果我們得到
n i=0
FiTn−i = Tn+2− Fn+2.
參考文獻
1. T. Koshy, Fibonacci and Lucas Numbers with Applications, John Wiley & Sons, Inc.
2001, Ch. 46.
2. I. G. Macdonald, Symmetric functions and Hall polynomials, Oxford university press, 1998.
3. R. P. Stanley and S. Fomin, Enumerative Combinatorics: Volume 2 (Cambridge Studies in Advanced Mathematics Book 62), 1999.
4. H. S. Wilf, Generating-functionology. AK Peters/CRC Press, 2005.
5. 陳建燁。 完全齊次對稱多項式(起) : 自由分解重組恆等式。 高中數學學科中心電子報第113期, 2016 年8月。
6. 陳建燁。 Fibonacci 與 Padovan 的對話 (上): 將 Padovan 數列用 「完全齊次對稱多項式」 表示。
數學傳播季刊, 42(1), 71-79, 2018。
7. 陳建燁, Fibonacci 與 Padovan 的對話 (下): F-P 卷積恆等式。 數學傳播季刊, 42(3), 66-73, 2018。
8. 陳建燁。 F-P 卷積恆等式的一頁證明。 數學傳播季刊, 43(3), 60-62, 2019。
9. 廖信傑。 用矩陣方法探討三階遞迴數列。 數學傳播季刊, 38(1), 36-55, 2014。
—本文作者廖信傑投稿時為美國 University of Miami 博士生, 薛昭雄任教於美國 University of Nevada, Las Vegas—