二階與三階線性遞迴序列和多項式

二階與三階線性遞迴序列和多項式

孫維良

1. 引言

一個數列 {an}n≥0 滿足

a_n+2 = pa_n+1+ qa_n

其中 a0, a₁, p, q 為常數, q ̸= 0, 且 p, q 與 n 無關, 被稱為二階線性遞迴數列 (a number sequence of linear recurrence relation of order 2)。 最常見的例子就是 Fibonacci 數列 {Fn}n≥0 滿足

F_n+2= F_n+1+ F_n,F₀ = 1,F₁ = 1。

長久以來數學家們已發現在 Fibonacci 數列上的諸多性質, 讀者可參考 Koshy [9]。其中最令 人感興趣的是 Fibonacci 數列與 Pascal 三角形 (或稱 「楊輝三角形」, 「賈憲三角形」) 的關係 (請參考第 3 節)。 而 Pascal 三角形又和多項式 (1 + x)ⁿ 展開式的係數有著密不可分的連結, 這份連結就自然而然地延伸到 Fibonacci 數列了。 因此, 我們可以更廣泛地思考: 對於一般的 二階線性遞迴序列, 是否存在某種多項式, 能夠類似 Fibonacci 數列與 (1 + x)ⁿ 有著不可分 離的關係? 答案是肯定的, 也是本文的主要目的之一。 我們的目的有:

1. 我們將利用一般二階線性遞迴序列 (包含數列與多項式序列) 的生成函數 (請參考第 2 節) 與其一般表示式來定義序列的對應多項式, 且這樣的對應是一對一的。

2. 探討一般三階線性遞迴序列的對應多項式 (請參考第 3.3 節)。

2. 生成函數之基本性質與例子

數列 {an}n≥0 的生成函數 (the generating function) 定義為 A(t) =

∑∞ n=0

a_ntⁿ。

⋆注意, 這裡我們不考慮 A(t) 的斂散性。

52

根據這個定義我們能夠知道, 如果有兩個數列的生成函數相同, 則這兩個數列的每一項都 是相同的, 也就是說是相同的數列。 對於生成函數, 我們會有如下的基本特性:

定理 1: 令 {an}n≥0 為一個二階線性遞迴數列滿足 an+2 = pa_n+1+ qa_n, 其中 a0, a₁, p, q 為常數, q ̸= 0, 則其生成函數 A(t) 為

A(t) = a₀+ (a₁− pa0)t

1− pt − qt² 。 (1)

證明: 我們考慮三個關係式

A(t) = a₀+ a₁t + a₂t² + a₃t³+· · · ptA(t) = pa₀t + pa₁t²+ pa₂t³+· · · qt²A(t) = qa₀t²+ qa₁t³+· · · 因為 an+2 = pa_n+1+ qa_n, 所以我們有

(1− pt − qt²)A(t) = a₀ + (a₁− pa0)t。

將等號兩邊同除以 (1 − pt − qt²), 可得到我們所要的 A(t)。  特別地, 當我們考慮二階線性遞迴多項式序列 {an(x)}n≥0 亦會有相同結果。

定理 2: 令 {an(x)}n≥0 為一個二階線性遞迴多項式序列滿足 a_n+2(x) = p(x)a_n+1(x) + q(x)a_n(x),

其中 a0(x), a₁(x), p(x), q(x)為變數 x 的多項式,q(x) ̸= 0, 則其生成函數 A(t) =

∑∞ n=0

a_n(x)tⁿ 可寫成

A(t) = a0(x) + (a1(x)− p(x)a0(x)) t

1− p(x)t − q(x)t² 。 (2)

我們來看幾個例子, 以下是幾種有名的數列與多項式序列:

例子 1: 令 {Fn}n≥0, {Ln}n≥0 與 {Fn(x)}n≥0, {Ln(x)}n≥0 分別為 Fibonacci 和 Lucas 數列與 Fibonacci 和 Lucas 多項式序列 (請參考 Koshy [9]), 其中 F0 = F₁ = 1;L₀ = 2, L₁ = 1;F₀(x) = 1, F₁(x) = x;L₀(x) = 2, L₁(x) = x 且滿足

F_n+2= F_n+1+ F_n L_n+2= L_n+1+ L_n;

F_n+2(x) = xF_n+1(x) + F_n(x) L_n+2(x) = xL_n+1(x) + L_n(x)。

則根據定理 1 與 2 所求得的生成函數式子, 我們有

F_n L_n F_n(x) L_n(x)

生成函數 1

1− t − t²

2− t 1− t − t²

1 1− xt − t²

2− xt 1− xt − t²

例子 2: 令 {T^n(k)(x)}n≥0 為第 k 類 Chebyshev 多項式序列, k = 1, 2, 3 或 4 (請參考 Chen 和 Louck [4])。 對任意的 k 皆滿足 T0^(k)(x) = 1 且

T_n+2^(k)(x) = 2x· T_n+1^(k)(x)− T_n^(k)(x)

其中 T1⁽¹⁾(x) = x,T₁⁽²⁾(x) = 2x,T₁⁽³⁾(x) = 2x− 1 和 T1⁽⁴⁾(x) = 2x + 1。 利用定理 2 的式 子, 我們有生成函數分別為

Tn^(k)(x) Tn⁽¹⁾(x) Tn⁽²⁾(x) Tn⁽³⁾(x) Tn⁽⁴⁾(x) 生成函數 1− xt

1− 2xt + t²

1 1− 2xt + t²

1− t 1− 2xt + t²

1 + t 1− 2xt + t²

例子 3: 令 {D⁽¹⁾n (x, a)}n≥0 與 {D⁽²⁾n (x, a)}n≥0 分別為第1類與第2類的 Dickson 多項式序 列, a ∈ R (請參考 Lidl, Mullen 和 Turnwald [10])。 二者皆滿足遞迴關係

D⁽ⁱ⁾_n+2(x, a) = xD_n+1⁽ⁱ⁾ (x, a)− aDn⁽ⁱ⁾(x, a), i = 1, 2

其中 D⁽¹⁾0 (x, a) = 2, D₁⁽¹⁾(x, a) = x 和 D⁽²⁾0 (x, a) = 1, D⁽²⁾₁ (x, a) = x。 因此根據定理 2 所求得的式子, 我們有生成函數

Dn⁽ⁱ⁾(t) Dn⁽¹⁾(t) D⁽²⁾n (t) 生成函數 2− xt

1− xt + at²

1 1− xt + at²

3. 對應多項式

現在考慮 Fibonacci 多項式 {Fn(x)}n≥0。Bicknell 在 [2] 中指出, 當我們把 Fn(x) 列 出來 (見表 4)。再把係數按照次序排好 (見表 5), 不考慮係數為 0 的數, 則會發現, 若將所有斜 率為 1 的對角線順時針旋轉 45^◦, 會形成 Pascal 三角形 (表 5 粗體部分與表 6)。

而 Pascal 三角形其實就與 (x + 1)ⁿ 的各項係數是一樣的, 我們將 (x + 1)ⁿ 展開並降 冪排好 (見表 7), 一樣考慮所有斜率為 1 的對角線, 每條對角線左端點都剛好對應到一個 n, 將該對角線的所有元素累加恰好就是 Fn(x),例如表 7 中, n = 4 對應到右邊的粗體部分即為

F4(x)。 事實上, 我們可以從 [2] 中所給的式子

F_n(x) =

⌊ⁿ₂⌋

∑

i=0

(n− i i

) xⁿ⁻²ⁱ

來看出這件事。 其中 (_n−i

i

)xⁿ⁻²ⁱ 剛好是 (x + 1)ⁿ⁻ⁱ 的某一項。

表 4:

F₀(x) = 1 F₁(x) = x F₂(x) = x²+ 1 F₃(x) = x³+ 2x F₄(x) = x⁴+ 3x²+ 1 F₅(x) = x⁵+ 4x³+ 3x F₆(x) = x⁶+ 5x⁴+ 6x²+ 1

...

表 5:

n x⁰ x¹ x² x³ x⁴ x⁵ x⁶ 0 1

1 0 1 2 1 0 1 3 0 2 0 1 4 1 0 3 0 1 5 0 3 0 4 0 1 6 1 0 6 0 5 0 1

... ...

表 6:

1

1 1

1 2 1

1 3 3 1

1 4 6 4 1

...

表 7:

n (x + 1)ⁿ 0 1

1 x 1 2 x² 2x 1 3 x³ 3x² 3x 1 4 x⁴ 4x³ 6x² 4x 1

... ...

但是這種多項式序列 {Fn(x)}n≥0與多項式 (x+1)ⁿ之間的對應關係並不是足夠好的, 舉 例來說: 考慮 Jacobsthal 多項式序列 {Jn(x)}n≥0 (請參考 [9]), 其中 J0(x) = J1(x) = 1,並 滿足

J_n+2(x) = J_n+1(x) + xJ_n(x)。

則我們觀察 {Jn(x)}n≥0 與多項式 (1 + x)ⁿ 展開後升冪排列之關係, 如下

表 8:

J0(x) = 1 J₁(x) = 1 J₂(x) = 1 + x J₃(x) = 1 + 2x

J₄(x) = 1 + 3x + x² ...

表 9:

n (1 + x)ⁿ 0 1

1 1 x 2 1 2x x² 3 1 3x 3x² x³ 4 1 4x 6x² 4x³ x⁴

... ... 事實上,

J_n(x) =

⌊ⁿ₂⌋

∑

i=0

(n− i i

) xⁱ。

因此, {Fn(x)}n≥0 與 {Jn(x)}n≥0 對應至同一個多項式。也就是說在這樣的對應之下, 若給定 一個與序列 {an(x)}n≥0 對應的多項式 Pn(x),則我們無法判斷原先序列 {an(x)}n≥0 的樣子 為何。 因此我們接下來的目的就是要改善這樣的情形, 我們將利用二階線性遞迴的一般項表示 式, 來定義一個可以與序列一一對應的多項式, 也就是說不同的序列會對應至不同的多項式。

3.1. 二階線性遞迴的一般項表示式

定理 3: 令 {an}n≥0 為一個二階線性遞迴數列, 且 a0, a₁, p, q 為已知, q ̸= 0, 滿足 an+2 = pa_n+1+ qa_n。 則

a_n = a₀

⌊ⁿ₂⌋

∑

k=0

(n− k k

)

p^n−2kq^k+ (a₁ − pa0)

⌊∑ⁿ⁻¹₂ ⌋ k=0

(n− 1 − k k

)

p^n−1−2kq^k, n≥ 1。

證明: 首先令

f (t) = 1 1− pt − qt², 則我們有

f (t) = 1

1− pt − qt² = 1 1− (pt + qt²)

=

∑∞ n=0

(pt + qt²)n

=

∑∞ n=0

( _n

∑

k=0

(n k

)

(pt)^n−k( qt²)k

)

=

∑∞ n=0

( _n

∑

k=0

(n k

)

p^n−kq^kt^n+k )

,

我們改變足標即可得到

f (t) =

∑∞ n=0



^⌊

n 2⌋

∑

k=0

(n− k k

)

p^n−2kq^k



 tⁿ。 (3)

現在, 令 A(t) 為 {an}n≥0 的生成函數。 利用定理 1 與式子 (3), 則

A(t) = a₀+ (a₁− pa0)t 1− pt − qt²

= a₀f (t) + (a₁− pa0)tf (t)

= a₀

∑∞ n=0



^⌊

n 2⌋

∑

k=0

(n−k k

)

p^n−2kq^k



 tⁿ+(a₁−pa0)

∑∞ n=0



^⌊

n 2⌋

∑

k=0

(n−k k

)

p^n−2kq^k



 tⁿ⁺¹。

tⁿ 的係數即是我們所要的 an, 所以

a_n = a₀

⌊ⁿ₂⌋

∑

k=0

(n− k k

)

p^n−2kq^k+ (a₁− pa0)

⌊∑ⁿ⁻¹₂ ⌋ k=0

(n− 1 − k k

)

p^n−1−2kq^k, n≥ 1。 

同樣地, 利用相同的證明, 對於二階線性遞迴多項式序列我們也會有相似的結果。¹

定理 4: 令 {an(x)}n≥0 為一個二階線性遞迴多項式序列, 且 a0(x), a₁(x), p(x), q(x) 為已 知, q(x) ̸= 0, 滿足 an+2(x) = p(x)a_n+1(x) + q(x)a_n(x)。 則

an(x) = a0(x)

⌊ⁿ₂⌋

∑

k=0

(n− k k

)

p(x)^n−2kq(x)^k

+ (a1(x)− p(x)a0(x))

⌊∑ⁿ⁻¹₂ ⌋ k=0

(n− 1 − k k

)

p(x)^n−1−2kq(x)^k, n≥ 1。

觀察定理 3 的式子中之

⌊ⁿ₂⌋

∑

k=0

(n− k k

)

p^n−2kq^k,其實是將 (p + q)ⁿ展開後斜率為 1 之加 總, 即

1an的一般表示式並不唯一, 一個強而有力的方法是利用二次方程式 X²− pX − q = 0 之根式解, 進而得到另一個一階線性遞迴關 係式以求得 an(請參考 [8])。

表 10:

n = 0:1 n = 1:p n = 2:p²+ q n = 3:p³+ 2pq

...

表 11:

n (p + q)ⁿ 0 1

1 p q 2 p² 2pq 3 p³ · · ·

... ... 因此,an 會與

a₀(p + q)ⁿ+ (a₁− pa0)(p + q)ⁿ⁻¹ = (p + q)ⁿ⁻¹(a₁+ qa₀) 有類似的關係。 根據這個觀察, 我們就可以定義出對應的規則。

3.2. 二階線性遞迴之對應多項式

定義 5: 沿用定理 3 的二階線性遞迴數列 {an}n≥0。 我們定義 CPm(X,{an}n≥0), m ≥ 0, 為 {an}n≥0 的對應多項式 (corresponding polynomial) 滿足:

1. CP0(X,{an}n≥0) = a₀;

2. CPm(X,{an}n≥0) = (p + qX)^m−1(a₁+ qa₀X),當 m ≥ 1。

若 {an}n≥0 為已知, 則我們將對應多項式簡記為 CPm(X)。 並且以符號 {an}n≥0 → CPm(X)

表示其對應。 我們也會用

{an}n≥0 → CPm(X),CP0(X) 來特別表示 m = 0 的情形。

特別地, 若 {an}n≥N, N ∈ N, 則 m ≥ N, 且 CPN(X) = aN 和 CPm(X) = (p + qX)^m−1−N(a1+ qa0X), 當 m ≥ N + 1。

接下來的引理代表其 「對應」 的意義:

引理 6: 令 {an}n≥0 → CPm(X), 則 a_n =∑

k≥0

[X^k] (CPn−k(X)) , n≥ 0, 其中 [X^k] (CPn−k(X)) 是指在 CPn−k(X)中 X^k 項的係數。

證明: 已知 a0 = [X⁰] (a0) = [X⁰] (CP0(X))。 當 n ≥ 1, CPn(X) = (p + qX)ⁿ⁻¹(a₁+ qa₀X)

= a0(p + qX)ⁿ+ (a1− pa0)(p + qX)ⁿ⁻¹

= a0

∑n l=0

(n l

)

p^n−lq^lX^l+ (a1− pa0)

n−1

∑

l=0

(n− 1 l

)

p^n−1−lq^lX^l。 所以

[X^k] (CPn−k(X)) = a₀

(n− k k

)

p^n−2kq^k+ (a₁ − pa0)

(n− k − 1 k

)

p^n−2k−1q^k。 因此

∑

k≥0

[X^k] (CPn−k(X)) =∑

k≥0

( a₀

(n− k k

)

p^n−2kq^k+ (a₁− pa0)

(n− k − 1 k

)

p^n−2k−1q^k )

= a₀

⌊ⁿ₂⌋

∑

k=0

(n−k k

)

p^n−2kq^k+(a₁−pa0)

⌊∑ⁿ⁻¹₂ ⌋ k=0

(n−k−1 k

)

p^n−2k−1q^k。 從定理 3 我們可以得到

an =∑

k≥0

[X^k] (CPn−k(X))。  因此, 在定義 5 中, 我們所謂 「對應」 的意思即為:an 就相當於以 CPn(X)為首, 將斜率 為 1 之各項係數累加,例如 a3 = [X⁰] (CP3(X)) + [X¹] (CP2(X)) = p²a₁+ (pqa₀ + qa₁) (見表 12 與13 粗體字部分)。

表 12:

a₀ = a₀ a₁ = a₁

a₂ = pa₁+ qa₀

a₃ = p²a₁+ (pqa₀+ qa₁) ...

表 13:

m X⁰ X¹ X² · · ·

CP0(X) a₀

CP1(X) a₁ qa₀

CP2(X) pa₁ pqa₀+ qa₁ q²a₀ CP3(X) p²a₁ · · ·

... ...

根據對應多項式的定義我們可以知道任何的二階線性遞迴數列皆可找到其對應多項式。 接 著我們要說明, 從對應多項式可以對應回唯一的一個二階線性遞迴數列, 也就是說這樣的對應 是一個好的對應。

引理 7: 令 {an}n≥0 與 {bn}n≥0 為兩個二階線性遞迴數列滿足

a_n+2= pa_n+1+ qa_n 與 bn+2 = p^′b_n+1+ q^′b_n, 其中 a0, a1, p, q, b0, b1, p^′, q^′ 已知, q ̸= 0, q^′ ̸= 0。 令

{an}n≥0 → CPm(X,{an}n≥0) 與 {bn}n≥0 → CPm(X,{bn}n≥0)。

如果對所有 m ≥ 0 我們都有 CPm(X,{an}n≥0) = CPm(X,{bn}n≥0), 則 an = bn, 對所 有 n ≥ 0, 也就是說兩數列相同。

證明: 因為對所有 m ≥ 0, CPm(X,{an}n≥0) =CPm(X,{bn}n≥0), 根據引理 6 我們可以 得到

a_n =∑

k≥0

[X^k] (CPn−k(X),{an}n≥0) =∑

k≥0

[X^k] (CPn−k(X),{bn}n≥0) = b_n。  引理 7 說明了數列與對應多項式可以一一對應, 所以我們可以將對應符號改寫為

{an}n≥0 ↔ CPm(X) 表示其對應。 我們來看一個與 Pascal 三角形相關的例子:

例子 4: 令 a0 = 1,a₁ = p,且 an+2= pa_n+1+ qa_n 則從定義 5 我們有 {an}n≥0 ↔ (p + qX)^m。

表 14:

a0 = 1 a₁ = p a₂ = p²+ q a₃ = p³+ 2pq

...

表 15:

m X⁰ X¹ X² · · · (p + qX)⁰ 1

(p + qX)¹ p q (p + qX)² p² 2pq q² (p + qX)³ p³ · · ·

... ... 我們將表 15 的對應改寫成三角形的樣子

1

↙ ↘

p q

↙ ↘ ↙ ↘

p² 2pq q²

往左邊的方向是乘以一個 p 往下累加,往右邊則是乘以一個 q 。當 p = q = 1 時, 即為 Fi- bonacci 數列與 Pascal 三角形的對應關係,而 an 就恰好為 Fn。

事實上, 依相同的方法我們亦可定義出一般二次線性遞迴多項式序列的對應多項式:

定義 8: 沿用定理 4 的二階線性遞迴多項式序列 {an(x)}n≥0。 我們定義 CPm(X,{an(x)}n≥0), m≥ 0, 為 {an(x)}n≥0 的對應多項式滿足:

1. CP0(X,{an(x)}n≥0) = a0(x);

2. CPm(X,{an(x)}n≥0) = (p(x) + q(x)X)^m−1(a₁(x) + q(x)a₀(x)X), 當 m ≥ 1。

我們可以利用對應多項式做個簡單的應用:

例子 5: 考慮多項式序列 {wn(x, y) = ^xn+1_x^−y_−yⁿ⁺¹}n≥0, 不難驗證 {wn(x, y)} 為一個二階線性 遞迴滿足

w_n+2(x, y) = (x + y)w_n+1(x, y)− xywn(x, y) 且 w0(x, y) = 1,w₁(x, y) = x + y,則根據定義 8 我們可以得到

w_n(x, y)↔ ((x + y) − xyX)^m。 令 x + y = 1,−xy = 1, 則可解得 x = ¹⁺₂^√5,y = ¹⁻₂^√5, 所以得到

w_n (

1 +√ 5

2 ,1−√ 5 2

)

↔ (1 + X)^m。

注意到 {Fn}n≥0 ↔ (1 + X)^m (請參考例子 6), 因此根據對應多項式的唯一性我們就會有 w_n

(1+√ 5 2 ,¹⁻₂^√5

)

= F_n,也就是說

F_n = 1

√5



 (

1 +√ 5 2

)n+1

− (

1−√ 5 2

)n+1

 。

我們在以下列出一些常見的數列與多項式序列所的對應多項式。

例子 6: 令 {an}n≥0 為數列或多項式序列滿足 an+2 = pan+1 + qan。 從例子 1, 2, 3 和 Jacobsthal 多項式我們可以得到

表 16:

{an}n≥0 (a0, a1, p, q) ↔ CPm(X),CP0(X) {Fn}n≥0 (1, 1, 1, 1) ↔ (1 + X)^m, 1

{Ln}n≥0 (2, 1, 1, 1) ↔ (1 + X)^m−1(1 + 2X), 2 {Pn}n≥0 (0, 1, 2, 1) ↔ (2 + X)^m−1, 0

{Fn(x)}n≥0 (1, x, x, 1) ↔ (x + X)^m, 1

{Ln(x)}n≥0 (2, x, x, 1) ↔ (x + X)^m−1(x + 2X)), 2 {Jn(x)}n≥0 (1, 1, 1, x) ↔ (1 + xX)^m, 1

{Tn⁽¹⁾(x)}n≥0 (1, x, 2x,−1) ↔ (2x − X)^m−1(x− X), 1 {Tn⁽²⁾(x)}n≥0 (1, 2x, 2x,−1) ↔ (2x − X)^m, 1

{Tn⁽³⁾(x)}n≥0 (1, 2x− 1, 2x, −1) ↔ (2x − X)^m−1(2x− 1 − X), 1 {Tn⁽⁴⁾(x)}n≥0 (1, 2x + 1, 2x,−1) ↔ (2x − X)^m−1(2x + 1− X), 1 {Dn⁽¹⁾(x, a)}n≥0 (2, x, x,−a) ↔ (x − aX)^m−1(x− 2aX), 2 {Dn⁽²⁾(x, a)}n≥0 (1, x, x,−a) ↔ (x − aX)^m, 1

{Bn(x, r)}n≥0 (1, x + r + 1, x + 2,−1) ↔ ((x + 2) − X)^m−1((x + r + 1)− X) , 1 {Bn(x)}n≥0 (1, x + 2, x + 2,−1) ↔ ((x + 2) − X)^m, 1

{bn(x)}n≥0 (1, x + 1, x + 2,−1) ↔ ((x + 2) − X)^m−1((x + 1)− X) , 1 {Wn(x, y)}n≥0 (2, x + y, x + y,−xy) ↔ ((x + y) − xyX)^m−1((x + y− 2xyX), 2 {wn(x, y)}n≥0 (1, x + y, x + y,−xy) ↔ ((x + y) − xyX)^m, 1

其中

• {Pn}n≥0 為 Pell 數列 (請參考 [3])。

• {Bn(x, r)}n≥0 為廣義的 Morgan-Voyce 多項式序列, r ∈ R (請參考 Andr´e-Jeannin [1]), {Bn(x) = Bn(x, 1)}n≥0 與 {bn(x) = Bn(x, 0)}n≥0 為古典的 Morgan-Voyce 多 項式序列。

• Wn(x, y) = xⁿ+ yⁿ 和 wn(x, y) = ^xn+1_x^−y_−yⁿ⁺¹,不難驗證兩者滿足遞迴關係 W_n+2(x, y) = (x + y)W_n+1(x, y)− xyWn(x, y);

w_n+2(x, y) = (x + y)w_n+1(x, y)− xywn(x, y)。

事實上, 定理 4 與定義 8 皆可應用至多變數多項式函數, 而不改其性質。 因此我們可以仿 照定理 4 的一般項表示來求得

xⁿ+ yⁿ=

⌊ⁿ₂⌋

∑

k=0

(−1)^k n n− k

(n− k k

)

(x + y)^n−2k(xy)^k;

xⁿ⁺¹− yⁿ⁺¹ x− y =

⌊ⁿ2⌋

∑

k=0

(−1)^k

(n− k k

)

(x + y)^n−2k(xy)^k。 這就是著名的 Girard-Waring 公式 (請參考 [6])。

3.3. 三階線性遞迴之對應多項式

現在, 我們將利用先前的概念來定義一般三階線性遞迴的對應多項式。

定義 9: 令 {an}n≥0 為三階線性遞迴數列滿足

a_n+3 = pa_n+2+ qa_n+1+ ra_n

其中 a0, a1, a2, p, q, r 為已知, r ̸= 0。 定義 {an}n≥0 的對應多項式 CPm(X) ( 或 CPm(X,{an}n≥0))

滿足:

1. CP0(X) = a₀;

2. CP1(X) = a₁+ qa₀X + ra₀X²;

3. CPm(X) = a₀(p + qX + rX²)^m+ (a₁− pa0)(p + qX + rX²)^m−1 + (a₂− pa1− qa0)(p + qX + rX²)^m−2, m≥ 2。

記為 {an}n≥0 → CPm(X) 。

接著我們要說明兩者的對應關係, 相當於要證明 an=∑

k≥0

[X^k] (CPn−k(X)) , ∀n ≥ 0。

在證明這件事之前, 我們需要引用多項式定理 (multinomial theorem) 與 an 的一般表示式。

定理 10:[多項式定理 [11]] 對於任意的正整數 m 與非負整數 n, 我們皆有 (x₁+ x₂ +· · · + xm)ⁿ= ∑

k1+k2+···+k^m=n k1,...,km≥0

( n

k₁, k₂, . . . , k_m )

x₁^k1x₂^k2· · · xmkm

其中 ( _n

k1,k2,...,km

)= _k ^n!

1!k2!···km! 。

引理 11:[Chou, Hsu 和 Shiue [5]] 令 {an}n≥0 滿足定義 9 的遞迴關係式且令 c0 = a₀,c₁ = (a₁ − pa0),c₂ = (a₂− pa1− qa0), 則可以得到

a_n =

∑2 i=0

c_i ∑

k1+2k2+3k3=n−i k1+k2+k3=k

k1,k2,k3≥0

( k k₁, k₂, k₃

)

p^k1q^k2r^k3, 對所有可能的 k。 (4)

當 n − i < 0 時, 右邊的和定義為 0。

注意到在引理 11 裡, k1+ k2+ k3 的上界為 n − i, 也就是說 0 ≤ k ≤ n − i, 所以我們 能夠將 k → n − i − k 以得到 k1+ k₂+ k₃ = n− i − k 而不改變 k 值的所有可能性, 因此 a_n 的一般表示式 (4) 可改寫為

a_n =

∑2 i=0

c_i ∑

k1+2k2+3k3=n−i k1+k2+k3=n−i−k

k1,k2,k3≥0

(n− i − k k1, k2, k3

)

p^k1q^k2r^k3, 對所有可能的 k。 (5)

接著我們就會有最重要的對應關係式:

引理 12: 令 {an}n≥0 滿足定義 9 的遞迴關係式,則 an=∑

k≥0

[X^k] (CPn−k(X)) , ∀n ≥ 0。

證明: 當 0 ≤ n < 2 時, 從定義 9 我們可以得到 a0 = [X⁰] (a0) = [X⁰] (CP0(X)) 和

a₁ = [X⁰](a₁+ qa₀X + ra₀X²) + [X¹](a₀) = [X⁰] (CP1(X)) + [X¹] (CP0(X))。 當 n ≥ 2 時,令 c0 = a₀,c₁ = (a₁ − pa0),c₂ = (a₂− pa1− qa0),則從定義 9 和多項式定理 我們有

CPn−k(X)

= c₀(p + qX + rX²)^n−k+ c₁(p + qX + rX²)^n−k−1+ c₂(p + qX + rX²)^n−k−2

=

∑2 i=0

c_i(p + qX + rX²)^n−k−i

=

∑2 i=0

c_i ∑

k1+k2+k3=n−k−i k1,k2,k3≥0

(n− k − i k₁, k₂, k₃

)

p^k1q^k2r^k3X^k2+2k3。

因此 ∑

k≥0

[X^k] (CPn−k(X))

=∑

k≥0

∑2 i=0

c_i ∑

k1+k2+k3=n−k−i k2+2k3=k k1,k2,k3≥0

(n− k − i k₁, k₂, k₃

)

p^k1q^k2r^k3

=

∑2 i=0

ci

∑

k1+2k2+3k3=n−i k1+k2+k3=n−k−i

k1,k2,k3≥0

(n− k − i k₁, k₂, k₃

)

p^k1q^k2r^k3, 對所有可能的 k。

所以從式子 (5) 我們就有

a_n=∑

k≥0

[X^k] (CPn−k(X)) , ∀n ≥ 0。 

同樣地, 我們也會有對應多項式對應至數列的唯一性:

引理 13: 令 {an}n≥0 與 {bn}n≥0 為兩個三階線性遞迴數列滿足 a_n+3= pa_n+2+ qa_n+1+ ra_n 和 b_n+3= p^′b_n+2+ q^′b_n+1+ r^′b_n 且令

{an}n≥0 → CPm(X,{an}n≥0) 與 {bn}n≥0 → CPm(X,{bn}n≥0)。

如果對所有 m ≥ 0 我們都有 CPm(X,{an}n≥0) = CPm(X,{bn}n≥0), 則 an = b_n, 對所 有 n ≥ 0, 也就是說兩序列相同。

證明: 同引理 7 的證明。 

實際上我們可效仿定義 8, 將定義 9 推廣至一般三階線性遞迴的多項式序列 {an(x)}n≥0。 最後, 我們列出一些常見的三階遞迴關係之對應多項式。

例子 7: 令 {an}n≥0為一般三階線性遞迴數列或多項式序列, 滿足 an+3= pa_n+2+qa_n+1ra_n。 則從定義 9 有

表 17:

{an}n≥0 (a₀, a₁, a₂, p, q, r) ↔ CPm(X),CP0(X),CP1(X) {Tn}n≥0 (0, 1, 1, 1, 1, 1) ↔ (1+X +X²)^m−1, 0, 1

{(P a)n}n≥0 (1, 1, 1, 0, 1, 1) ↔ (X +X²)^m−1(1+X +X²), 1, (1+X +X²)

{(P e)n}n≥0 (3, 0, 2, 0, 1, 1) ↔ (X +X²)^m−2(−1+3X²+6X³+3X⁴), 3, (3X +3X²) {FnFn+1}n≥0 (1, 2, 6, 2, 2,−1) ↔ (2+2X −X²)^m, 1, (2+2X−X²)

{Tn(x)}n≥0 (0, 1, x², x², x, 1) ↔ (x²+xX +X²)^m−1, 0, 1 其中

• {Tn}n≥0 與 {Tn(x)}n≥0 為 tribonacci 數列與多項式序列 (請參考 [9])。

• {(P a)n}n≥0 為 Padovan 數列 (請參考 [7] 和 [12])。

• {(P e)n}n≥0 為 Perrin 數列 (請參考 [13])。

• {Fn}n≥0 為 Fibonacci 數列。更一般來說, 若二階線性遞迴數列 {an}n≥0 滿足 an+2 = pa_n+1+ qa_n, 則 {bn= a_na_n+1}n≥0 會形成一個三階線性遞迴數列, 滿足

b_n+3 = (p²+ q)b_n+2+ (p²q + q²)b_n+1− q³b_n。

備註: 四階以上的線性遞迴關係式亦有其相同性質, 可定義對應多項式。 但較為繁瑣與複雜, 我 們將以另文討論。

致謝

本文是參加中央研究院數學研究所的暑期研習課程 「組合數學與圖論專題」 時撰寫的, 由 衷地感謝美國內華達大學數學系薛昭雄教授的親自指導, 老師除了訓練我們做研究與找資料的 能力之外, 更不厭其煩地幫我們修改文章中的錯誤與提供建議, 使本文得以順利完成。 非常感 謝。

參考資料

1. R. Andr´e-Jeannin, A generalization of Morgan-Voyce polynomials, Fibonacci Quart. 32 (1994), no. 3, 228-231.

2. M. Bicknell, A primer for the Fibonacci numbers: part VII, Fibonacci Quart. 8 (1970), no. 4, 407-420.

3. M. Bicknell, A primer on the Pell sequence and related sequences, Fibonacci Quart. 13 (1975), no. 4, 345-349.

4. W. Y. C. Chen and J. D. Louck, The combinatorial power of the companion matrix, Linear Algebra Appl. 232 (1996), no. 1, 261-278.

5. W.-S. Chou, Leetsch C. Hsu, and Peter J.-S. Shiue, Application of Fa`a di Bruno’s formula in characterization of inverse relations, J. Comput. Appl. Math. 190 (2006), 151-169.

6. H. W. Gould, The Girard-Waring power sum formulas for symmetric functions and Fibonacci sequences, Fibonacci Quart. 37 (1999), no. 2, 135-140.

7. Thomas M. Green, Recurrent sequences and Pascal’s triangle, Math. Mag. 41 (1968), no. 1, 13-21.

8. T.-X. He and Peter J.-S. Shiue, On sequences of numbers and polynomials defined by linear recurrence relations of order 2, Int. J. Math. Math. Sci. 2009 (2009), Art. ID 709386, 21 pp.

9. T. Koshy, Fibonacci and Lucas Numbers with Applications, Wiley-Interscience, 2001.

10. R. Lidl, G. L. Mullen, and G. Turnwald, Dickson polynomials, Longman Scientific &

Technical, 1993.

11. National Institute of Standards and Technology, NIST Digital Library of Mathematical Functions, Section 26.4. http://dlmf.nist.gov/26.4.

12. I. Stewart, Mathematical recreations, Sci. Amer. 274 (1996), no. 6, 102-103.

13. F. Yilmaz and D. Bozkurt, Hessenberg matrices and the Pell and Perrin numbers, J.

Number Theory 131 (2011), no. 8, 1390-1396.

—本文作者2012年畢業於國立成功大學應用數學所碩士班^—

第 37 _卷第 2 _期 (146 _號 ) _勘誤

第35頁第8行

s

d = h

y− h − s− s^′ d− d^′

為 s

d = h

y− h = s− s^′ d− d^′ 之誤。