翻轉塗色—Thue-Morse 數列的同色 3-AP 個數
莊沅蓉
國立臺南女子高級中學
目錄
3 3
—AP3(Tn) = 4n−2− 2n−2 4
( ) #{Tn I 3-AP} = #{Tn II 3-AP} = AP3(Tn−1) . . . 6
( ) #{Tn III 3-AP} = #{Tn IV 3-AP} . . . . 6
( ) #{Tn III 3-AP} = 4n−3 . . . 9
( ) AP3(Tn) = 2· AP3(Tn−1) + 2· 4n−3 . . . 13
—AP3(Tn,p) 13 ( ) #{Tn iii 3-AP} = #{Tp+1 III 3-AP} . . . 15
( ) #{Tn iv 3-AP} . . . 15 17
一 問題
會考慮本問題是在我看完 [3] 之後所獲得的靈感。令 A = a1a2· · · an 為一長度為 n 的雙色字 串,也就是任一個 a1, a2, ..., an不是 0 就是 1。如果存在正整數 i 及 d 滿足 ai = ai+d= ai+2d, 則我們稱 (i, i + d, i + 2d) 為 A 的一個同色 3-AP。本文的目的在決定長度為 2n的 Thue-Morse 雙色字串 1 0 01 0110 01101001| {z · · ·}
2nbits
有幾個同色 3-AP。
參考文獻 [1],[2]
二 定義列表
定義 一. 一個長度為 n 的雙色字串定義為一個二進位表示的 n 位數字。
定義 二. 令 A = a1a2· · · an, B = b1b2· · · bm 為兩個長度分別為 n, m 的雙色字串,定義 AB = a1a2· · · anb1b2· · · bm.
定義 三. 若 A = a1a2· · · an 為一個長度為 n 的雙色字串,則定義 A 的互補雙色字串為 A = (1− a1)(1− a2)· · · (1 − an).
定義 四. 令 A = a1a2· · · an 為一長度為 n 的雙色字串,定義 A 的一個同色 3-AP 為一數組 (i, i + d, i + 2d),此數組滿足 ai = ai+d = ai+2d,其中 d 稱為此同色 3-AP 的公差。並將 A 的 所有同色 3-AP 的個數記為
AP3(A) = #{(i, i + d, i + 2d) | ai = ai+d = ai+2d}.
例子 1. 假設 A = 1001011001101001,(1, 4, 7) 為 A 的一個同色 3-AP,(2, 8, 14) 為 A 的另一 個同色 3-AP,
A = 1
a1
a02
a03
a14
a05
a16
a17
a08
a09
a110
a111
a012
a113
a014
a015
a116
定義 五. 令 T0 = 1,並對於所有的非負整數 n,定義長度為 2n 且以 1 開頭的 Thue-Morse 雙色字串為 Tn= Tn−1Tn−1。
例子 2.
T0 = 1,
T1 = 11 = 10, T2 = 10 10 = 1001,
T3 = 1001 1001 = 10010110,
T4 = 10010110 10010110 = 1001011001101001, ...
定義 . 令 M 為一元素為二進位數字的矩陣,定義 N(M ) 表示矩陣 M 中,全部為 0 的行 數加上全部為 1 的行數。
例子 3. 例如 N(
0 1 0 1 1 1 1 1 0 0 1 1 1 1 0 0 0 1
) = 3.
三 主要的結果—AP
3(T
n) = 4
n−2− 2
n−2引理 七. Tn
不會有公差
d = 1或
d = 2或
d > 2n−1− 1的同色
3-AP。
定義 . 給定一個固定的正整數 d,其中 3 ≤ d ≤ 2n−1− 1,我們將 Tn 按下圖的方式排列,
稱為 Tn 的 d-中心排列。
d bits
z }| {
Tn−1
...
(2n−1 mod d) bitsz }| {
a2n−1−2d+1 · · · · · · · a2n−1−2d+i · · · a2n−1−d
← MIII(d) a2n−1−d+1 · · · · a2n−1−d+i · · · a2n−1
Tn−1 a2n−1+1 · · · · a2n−1+i · · · a2n−1+d
a2n−1+d+1 · · · · · · · · · · · · · · a2n−1+2d
↑
.. .
MIV(d)
例子 4. 例如 T4 = 1001 0110| 0110 1001 的 d = 6-中心排列為
6 bits
z }| {
2 bits
z}|{
T3 1 0 ← MIII(6)
0 1 0 1 1 0
T3 0 1 1 0 1 0
0 1
↑ MIV(6)
例如 T4 = 1001 0110| 0110 1001 的 d = 5-中心排列為
5 bits
z }| { T3 z }| {3 bits
1 0 0 ← MIII(5) 1 0 1 1 0
T3 0 1 1 0 1 0 0 1
↑ MIV(5)
例如 T5 = 1001 0110 0110 1001| 0110 1001 1001 0110 的 d = 7-中心排列為
7 bits
z }| {
7 bits
z }| { MIII(7)
T4 1 0 ↙
0 1 0 1 1 0 0 1 1 0 1 0 0 1
T4 0 1 1 0 1 0 0
1 1 0 0 1 0 1
↗ 1 0
MIV(7)
觀察. 注意到,上圖中若三個 0 或三個 1 出現在同一行,就對應到 Tn 的一個公差為 d 的同 色 3-AP。例如上圖中 MIII(d) 矩陣中的一個全 0 行或全 1 行就對應到 Tn 的一個公差為 d 的同色 3-AP。利用這個 d-中心排列,我們可以把求 Tn 的同色 3-AP 的個數問題,轉化成求 矩陣中全 0 行及全 1 行的個數問題。
定義 九. 在 Thue-Morse 雙色字串 Tn= Tn−1Tn−1 = a1a2· · · a2n−1a2n−1+1· · · a2n 中,令 (i, i + d, i + 2d) 為 Tn 的一個同色 3-AP,
• 若 i < i + d < i + 2d≤ 2n−1,則稱 (i, i + d, i + 2d) 為 Tn 的第 I 類同色 3-AP。直覺來 看,表示這個 3-AP 是落在 d-中心排列的上半部。
• 若 2n−1 < i < i + d < i + 2d,則稱 (i, i + d, i + 2d) 為 Tn 的第 II 類同色 3-AP。直覺來 看,表示這個 3-AP 是落在 d-中心排列的下半部。
• 若 i < i + d≤ 2n−1 < i + 2d,則稱 (i, i + d, i + 2d) 為 Tn 的第 III 類同色 3-AP。直覺來 看,表示這個 3-AP 是落在 MIII(d)中的。
• 若 i≤ 2n−1 < i + d < i + 2d,則稱 (i, i + d, i + 2d) 為 Tn 的第 IV 類同色 3-AP。直覺來 看,表示這個 3-AP 是落在 MIV(d) 中的。
接下來我們分四個部份
• (一):證明 #{Tn 的第 I 類同色3-AP} = #{Tn 的第 II 類同色3-AP} = AP3(Tn−1),於 是我們可以得到
AP3(Tn) = #{Tn 的第 I 類同色3-AP} + #{Tn 的第 III 類同色3-AP} + #{Tn 的第 II 類同色3-AP} + #{Tn 的第 IV 類同色3-AP}
= 2· AP3(Tn−1) + #{Tn 的第 III 類同色3-AP} + #{Tn 的第 IV 類同色3-AP}
• (二):證明 #{Tn 的第 III 類同色3-AP} = #{Tn 的第 VI 類同色3-AP},於是我們可以 得到
AP3(Tn) = 2· AP3(Tn−1) + 2· #{Tn 的第 III 類同色3-AP}
• (三):證明 #{Tn 的第 III 類同色3-AP} = 4n−3,於是我們可以得到 AP3(Tn) = 2· AP3(Tn−1) + 2· 4n−3
• (四):求解遞迴方程式 AP3(Tn) = 2· AP3(Tn−1) + 2· 4n−3,於是我們可以得到最主要的 結果
AP3(Tn) = 4n−2− 2n−2.
(一) # {T
n的第 I 類同色3-AP } = #{T
n的第 II 類同色3-AP } = AP
3(T
n−1)
引理 十.
假設雙色字串
A = a1a2· · · an,則
AP3(A) = AP3(A)。
Proof. 因為 ai = ai+d= ai+2d ⇔ ai = ai+d= ai+2d,所以AP3(A) = #{(i, i + d, i + 2d) | ai = ai+d = ai+2d}
= #{(i, i + d, i + 2d) | ai = ai+d = ai+2d} = AP3(A).
引理 十一. #{Tn
的第
I類同色
3-AP} = #{Tn的第
II類同色
3-AP} = AP3(Tn−1) Proof. 假設 Thue-Morse 雙色字串 Tn = Tn−1Tn−1 = a1a2· · · a2n−1a2n−1+1· · · a2n,#{Tn 的第 I 類同色3-AP}
= #{(i, i + d, i + 2d) | ai = ai+d= ai+2d, i < i + d < i + 2d≤ 2n−1}
= AP3(Tn−1)
引理 十↓
= AP3(Tn−1)
Tn=Tn−1Tn−1
=↓ #{(i, i + d, i + 2d) | ai = ai+d= ai+2d, 2n−1 < i < i + d < i + 2d}
= #{Tn 的第 II 類同色3-AP}
(二) # {T
n的第 III 類同色3-AP } = #{T
n的第 IV 類同色3-AP }
由觀察 三,我們可以得到 引理 十二.
• #{Tn
的第
III類同色
3-AP} =2n∑−1−1 d=3
N(MIII(d))
,
• #{Tn
的第
IV類同色
3-AP} =2n−1∑−1 d=3
N(MIV(d))
。
觀察.T0 = 1 T1 = 10 T2 = 10| 01
T3 = 1001| 0110 ...
更一般地,想像一下在 Tn 中,Tn−1 與 Tn−1 中間有一條中心分隔線,然後考慮與中心分 隔線等距的對應字元。
Tn= a1a2· · · a2n−j· · · a2n−1a2n
| {z }
Tn−1
| a2n+1a2n+2· · · a2n+(j+1)· · · a2n
| {z }
Tn−1
由這個觀察我們可以得到
引理 十三.
假設
Tn= Tn−1Tn−1 = a1a2· · · a2n−1a2n−1+1· · · a2n,
•
如果
n是偶數,則有
a2n−1−j = a2n−1+(j+1)•
如果
n是奇數,則有
a2n−1−j = a2n−1+(j+1)。
觀察. 我們來看一下這些跨越 Tn 中心分隔線的對應字元在 d-中心排列中會如何展現。
• a2n−1−2d+i = a2n−1−(2d−i) 對應到 a2n−1+((2d−i)+1)。
• a2n−1−d+i = a2n−1−(d−i) 對應到 a2n−1+((d−i)+1)。
• a2n−1+i = a2n−1+((i−1)+1) 對應到 a2n−1−(i−1)(右到左)。
d bits
z }| {
...
a2n−1−2d+1 · · · · · · · · · · · · · a2n−1−2d+i · · · a2n−1−d
← MIII(d) a2n−1−d+1 · · · a2n−1−(i−1) · · · · a2n−1−d+i · · · a2n−1
a2n−1+1 · · · a2n−1+((d−i)+1) · · · · a2n−1+i · · · a2n−1+d
a2n−1+d+1 · · · a2n−1+((2d−i)+1) · · · · · · · · · · · · · a2n−1+2d
↑
...
MIV(d)
注意到
a2n−1−2d+i, a2n−1−d+i,
a2n−1+i
三個字元在同一行上,而其分別對應的
a2n−1+((2d−i)+1), a2n−1+((d−i)+1),
a2n−1−(i−1)
三個字元也在 同一行上,因為他們的位置公差亦為 d。
引理 十四.
如果
a2n−1−2d+i = a2n−1−d+i = a2n−1+i,則
a2n−1+((2d−i)+1) = a2n−1+((d−i)+1) = a2n−1−(i−1)。也就是
a2n−1−2d+i
||
a2n−1−d+i
||
a2n−1+i
⇒
a2n−1+((2d−i)+1)
||
a2n−1+((d−i)+1)
||
a2n−1−(i−1)
.
Proof. 當 n 是偶數時,
a2n−1−2d+i = a2n−1−d+i= a2n−1+i
⇒ a2n−1−(2d−i)= a2n−1−(d−i) = a2n−1+((i−1)+1)
引理 十三
⇒↓ a2n−1+((2d−i)+1) = a2n−1+((d−i)+1) = a2n−1−(i−1)
當 n 是奇數時,
a2n−1−2d+i = a2n−1−d+i= a2n−1+i
⇒ a2n−1−(2d−i)= a2n−1−(d−i) = a2n−1+((i−1)+1)
引理 十三
⇒↓ a2n−1+((2d−i)+1) = a2n−1+((d−i)+1) = a2n−1−(i−1)
⇒ a2n−1+((2d−i)+1) = a2n−1+((d−i)+1) = a2n−1−(i−1)
引理 十四表明了,MIII(d) 中全 0 行及全 1 行的個數與 MIV(d) 中全 0 行及全 1 行的個 數會是相等的,也就是 N(MIII(d)) = N(MIV(d))。於是再由引理 十二,可以得到
#{Tn 的第 III 類同色3-AP}
引理 十四
=↓
2n−1∑−1 d=3
N(MIII(d))
引理 十二
=↓
2n∑−1−1 d=3
N(MIV(d))
引理 十四
=↓ #{Tn 的第 IV 類同色3-AP}
(三) # {T
n的第 III 類同色3-AP } = 4
n−3由引理 十二,我們有 #{Tn 的第 III 類同色3-AP} = 2
n∑−1−1 d=3
N(MIII(d)),所以我們來計算 N(MIII(d))。
1 觀察 MIII(d) 並分組
定義 十五. 令 Tm(ai, aj) 表示 Thue-Morse 雙色字串 Tm 中,第 i 個字元到第 j 個字元組成 的雙色字串。
我們觀察 MIII(d),發現可將其按照 d 的值分成三類,
觀察.
MIII(d) =
Tn−1(a1, as) Tn−1(a2n−1−s+1, a2n−1) Tn−1(a2n−1−2s+1, a2n−1−s)
當d∈ S1={2n−1− s | s = 1, 2, ..., 2n−2};
Tn−2k(a2s+1, a2n−2k) Tn−2k(a1, as) Tn−2k(as+1, a2n−2k−s) Tn−2k(a2n−2k−s+1, a2n−2k)
Tn−2k(a1, a2n−2k−2s) Tn−2k(a2n−2k−2s+1, a2n−2k−s)
當d∈ S2=
{
2n−2k− s | k=1,2,...,⌊n−22 ⌋ s=1,2,...,2n−2k−1
}
;
Tn−2k+1(a2s+1, a2n−2k+1) Tn−2k+1(a1, as)
Tn−2k+1(as+1, a2n−2k+1−s) Tn−2k+1(a2n−2k+1−s+1, a2n−2k+1) Tn−2k+1(a1, a2n−2k+1−2s) Tn−2k+1(a2n−2k+1−2s+1, a2n−2k+1−s)
當d ∈ S3= {
2n−2k+1− s | k=2,...,⌊n−22 ⌋
s=1,2,...,2n−2k
} .
2 考如何計算 N(MIII(d))
定義 十 . 假設雙色字串 A = a1a2· · · an−1an,定義 A 的逆向雙色字串 ←−
A = anan−1· · · a2a1。 引理 十七.
假設
Thue-Morse雙色字串
Tn= a1a2· · · a2n−1a2n,且假設
←T−n = b1b2· · · b2n−1b2n = a2na2n−1· · · a2a1,則
bi = a2n−i+1,其中
i≤ 2n。 別地,我們有
Tm(ai, aj) =←−Tm(a2m−j+1, a2m−i+1)
。
Proof.Tm(ai, aj) = ←−
Tm(aj, ai)
= ←−
Tm(bj, bi)
= ←−
Tm(b2m−j+1, b2m−i+1)
重新命名←−−
Tm 的字元
=↓ ←−
Tm(a2m−j+1, a2m−i+1)
觀察. 因為我們是考慮 N(MIII(d)),所以我們可以對 MIII(d) 做下面兩種動作而不影響 N(MIII(d)) 的值,
• 同時將 MIII(d)的每一列取逆向或互補,
• 做矩陣列互換的動作,
所以
當 d ∈ S1 時,
N(MIII(d)) = N(
Tn−1(a1, as) Tn−1(a2n−1−s+1, a2n−1) Tn−1(a2n−1−2s+1, a2n−1−s)
)引理= N(↓
←−−Tn−1(a2n−1, a2n−1−s+1)
←−−Tn−1(as, a1)
←−−Tn−1(a2s, as+1)
)
標號逆向而已
=↓ N(
←−−Tn−1(a2n−1−s+1, a2n−1)
←−−Tn−1(a1, as)
←−−Tn−1(as+1, a2s)
)
1, 3 列互換
=↓ N(
←−−Tn−1(as+1, a2s)
←−−Tn−1(a1, as)
←−−Tn−1(a2n−1−s+1, a2n−1)
)
化簡↓
= N(
Tn−1(as+1, a2s) Tn−1(a1, as) Tn−1(a2n−1−s+1, a2n−1)
)我們定義=↓ N(ABB(n− 1, s)).
當 d ∈ S2 時,
N(MIII(d)) = N(
Tn−2k(a2s+1, a2s) Tn−2k(a1, as) Tn−2k(as+1, a2n−2k−s) Tn−2k(a2n−2k−s+1, a2n−2k)
Tn−2k(a1, a2n−2k−s) Tn−2k(a2n−2k−2s+1, a2n−2k−s)
)
= N(
Tn−2k(a2s+1, a2s) Tn−2k(as+1, a2n−2k−s)
Tn−2k(a1, a2n−2k−s)
) + N(
Tn−2k(a1, as) Tn−2k(a2n−2k−s+1, a2n−2k) Tn−2k(a2n−2k−2s+1, a2n−2k−s)
)
我們定義↓
= N(γ(n− 2k, s)) + N(AAA(n − 2k, s)).
當 d ∈ S3 時,類似地,
N(MIII(d)) = N(
Tn−2k+1(a2s+1, a2n−2k+1) Tn−2k+1(a1, as)
Tn−2k+1(as+1, a2n−2k+1−s) Tn−2k+1(a2n−2k+1−s+1, a2n−2k+1) Tn−2k+1(a1, a2n−2k+1−2s) Tn−2k+1(a2n−2k+1−2s+1, a2n−2k+1−s)
)
我們定義↓
= N(ϕ(n− 2k + 1, s)) + N(ABB(n − 2k + 1, s)).
定義 十 . 定義
• ABB(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• AAA(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• γ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
• ϕ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
定義 十九.
• ABB(m) =2
m−1
∑
s=1
N(ABB(m, s)).
• AAA(m) =2
m−1
∑
s=1
N(AAA(m, s)).
• γ(m) = 2
m−1
∑
s=1
N(γ(m, s)).
• ϕ(m) = 2
m−1
∑
s=1
N(ϕ(m, s)).
注意到 γ(m, s) 是矩陣,γ(m) 是數值,其他符號也是類似地。
引理 二十.
#
{Tn 的第III 類同色3-AP} =
2n
∑
−1−1 d=3N(MIII
(d))
=
∑
d∈S1
N(ABB(n, s))
∑ +
d∈S1
[N(γ(n
− 2k, s)) + N(AAA(n − 2k, s))]∑ +
d∈S1
[N(ϕ(n
− 2k + 1, s)) + N(ABB(n − 2k + 1, s))]=
ABB(n− 1)
+
⌊n
∑
−22 ⌋ k=1[γ(n
− 2k) + AAA(n − 2k)]+
⌊n−2
∑
2 ⌋ k=1[ϕ(n
− 2k + 1) + ABB(n − 2k + 1)].
3 計算 ABB(m), AAA(m), γ(m), ϕ(m) 的公式(巧解 個遞迴式)
定義 二十一. 定義
• ABB(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• AAB(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• AAA(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• ABA(m, s) =
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
.
• γ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
• θ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
• ϕ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
• δ(m, s) =
Tm(a2s+1, a2m) Tm(as+1, a2m−s)
Tm(a1, a2m−2s)
.
定義 二十二.
• ABB(m) =2
m−1
∑
s=1
N(ABB(m, s)).
• AAB(m) =2
m−1
∑
s=1
N(AAB(m, s)).
• AAA(m) =2
m−1∑
s=1
N(AAA(m, s)).
• ABA(m) =2
m−1
∑
s=1
N(ABA(m, s)).
• γ(m) = 2
m−1
∑
s=1
N(γ(m, s)).
• θ(m) =2
m−1
∑
s=1
N(θ(m, s)).
• ϕ(m) = 2
m−1∑
s=1
N(ϕ(m, s)).
• δ(m) =2
m−1
∑
s=1
N(δ(m, s)).
注意到 γ(m, s) 是矩陣,γ(m) 是數值,其他符號也是類似地。
引理 二十三.
ABB(n) AAB(n) AAA(n) ABA(n)
γ(n) θ(n) ϕ(n) δ(n)
=
0 0 0 2 1 0 0 0 0 0 2 0 0 1 0 0 1 1 0 0 0 0 1 0 1 1 0 0 0 0 0 1 0 0 2 2 2 0 0 0 0 2 2 0 0 2 0 0 2 1 0 1 0 0 2 0 2 1 0 1 0 0 0 2
·
ABB(n− 1) AAB(n− 1) AAA(n− 1) ABA(n− 1)
γ(n− 1) θ(n− 1) ϕ(n− 1) δ(n− 1)
+
2n−2
0 0 0 2n−2
0 0 0
Proof. 我們只證明第一個與第五個遞迴關係式,其餘的遞迴關係式證明方法皆類似。
ABB(m)
=
2∑m−1
s=1
N(ABB(m, s))
=
2∑m−1
s=1
N(
Tm(as+1, a2s) Tm(a1, as) Tm(a2m−s+1, a2m)
)
=
∑
s<2m−2
N(
Tm−1(as+1, a2s) Tm−1(a1, as) Tm−1(a2m−1−s+1, a2m−1)
) +
∑
s=2m−1
N(
Tm−1(a2m−2+1, a2m−1) Tm−1(a1, a2m−2) Tm−1(a2m−2+1, a2m−1)
) +
∑
s>2m−2
N(
[c|c]Tm−1(as+1, a2m+1) Tm−1(a1, a2s−2m−1) Tm−1(a1, a2m−1−s) Tm−1(a2m−1−s+1, as) Tm−1(a2m−1−s+1, a2m−2s) Tm−1(a2m−1−2s+1, a2m−1)
)
=
∑
s<2m−2
[N(ABA(m− 1, s))]
+ 2m−2
∑+
s>2m−2
[N(ABA(m− 1, 2m−2− s) + N(γ(m − 1, 2m−2− s))]
= 2· ABA(m − 1) + γ(m − 1) + 2m−2.