LCS 演算法:Javascript 最長公共子序列

發表於2018-01-23

最長公共子序列(Longest Common Subsequence LCS)是從給定的兩個序列X和Y中取出儘可能多的一部分字元,按照它們在原序列排列的先後次序排列得到。LCS問題的演算法用途廣泛,如在軟體不同版本的管理中,用LCS演算法找到新舊版本的異同處;在軟體測試中,用LCS演算法對錄製和回放的序列進行比較,在基因工程領域,用LCS演算法檢查患者DNA連與鍵康DNA鏈的異同;在防抄襲系統中,用LCS演算法檢查論文的抄襲率。LCS演算法也可以用於程式程式碼相似度度量,人體執行的序列檢索,視訊段匹配等方面,所以對LCS演算法進行研究具有很高的應用價值。

基本概念

  1. 子序列(subsequence): 一個特定序列的子序列就是將給定序列中零個或多個元素去掉後得到的結果(不改變元素間相對次序)。例如序列的子序列有:等。
  2. 公共子序列(common subsequence): 給定序列X和Y,序列Z是X的子序列,也是Y的子序列,則Z是X和Y的公共子序列。例如X=[A,B,C,B,D,A,B],Y=[B,D,C,A,B,A[,那麼序列Z=[B,C,A]為X和Y的公共子序列,其長度為3。但Z不是X和Y的最長公共子序列,而序列[B,C,B,A]和[B,D,A,B]也均為X和Y的最長公共子序列,長度為4,而X和Y不存在長度大於等於5的公共子序列。對於序列[A,B,C]和序列[E,F,G]的公共子序列只有空序列[]。
  3. 最長公共子序列:給定序列X和Y,從它們的所有公共子序列中選出長度最長的那一個或幾個。
  4. 子串: 將一個序列從最前或最後或同時刪掉零個或幾個字元構成的新系列。區別與子序列,子序列是可以從中間摳掉字元的。cnblogs這個字串中子序列有多少個呢?很顯然有27個,比如其中的cb,cgs等等都是其子序列

給一個圖再解釋一下:

image_1c3ofq79gg4f1u7sf3o40gls29

我們可以看出子序列不見得一定是連續的,連續的是子串。

問題分析

我們還是從一個矩陣開始分析,自己推匯出狀態遷移方程。

首先,我們把問題轉換成前端夠為熟悉的概念,不要序列序列地叫了,可以認為是陣列或字串。一切從簡,我們就估且認定是兩個字串做比較吧。

我們重點留意”子序列“的概念,它可以刪掉多個或零個,也可以全部幹掉。這時我們的第一個子序列為空字串(如果我們的序列不是字串,我們還可以 )!這個真是千萬要注意到!許多人就是看不懂《演算法導論》的那個圖表,還有許多部落格的作者不懂裝懂。我們總是從左到右比較,當然了第一個字串,由於作為矩陣的高,就垂直放置了。

X “” B D C A B A
“”
A
B
C
D
A
B

假令X = “ABCDAB”, Y=”BDCABA”,各自取出最短的序列,也就是空字串與空字串比較。LCS的方程解為一個數字,那麼這個表格也只能填數字。兩個空字串的公同區域的長度為0.

X “” B D C A B A
“” 0
A
B
C
D
A
B

然後我們X不動,繼續讓空字串出陣,Y讓“B”出陣,很顯然,它們的公共區域的長度為0. Y換成其他字元, D啊,C啊, 或者, 它們的連續組合DC、 DDC, 情況沒有變, 依然為0. 因此第一行都為0. 然後我們Y不動,Y只出空字任串,那麼與上面的分析一樣,都為0,第一列都是0.

X “” B D C A B A
“” 0 0 0 0 0 0 0
A 0
B 0
C 0
D 0
A 0
B 0

LCS問題與揹包問題有點不一樣,揹包問題還可以設定-1行,而最長公共子序列因為有空子序列的出現,一開始就把左邊與上邊固定死了。

然後我們再將問題放大些,這次雙方都出一個字元,顯然只有兩都相同時,才有存在不為空字串的公共子序列,長度也理解數然為1。

A為”X”, Y為”BDCA”的子序列的任意一個

X “” B D C A B A
“” 0 0 0 0 0 0 0
A 0 0 0 0 1
B 0
C 0
D 0
A 0
B 0

繼續往右填空,該怎麼填?顯然,LCS不能大於X的長度,Y的從A字串開始的子序列與B的A序列相比,怎麼也能等於1。

X “” B D C A B A
“” 0 0 0 0 0 0 0
A 0 0 0 0 1 1 1
B 0
C 0
D 0
A 0
B 0

如果X只從派出前面個字元A,B吧,亦即是“”,“A”, “B”, “AB”這四種組合,前兩個已經解說過了。那我們先看B,${X_1} == ${Y_0}, 我們得到一個新的公共子串了,應該加1。為什麼呢?因為我們這個矩陣是一個狀態表,從左到右,從上到下描述狀態的遷移過程,並且這些狀態是基於已有狀態累加出來的。現在我們需要確認的是,現在我們要填的這個格子的值與它周圍已經填好的格子的值是存在何種關係。目前,資訊太少,就是一個孤點,直接填1。

X “” B D C A B A
“” 0 0 0 0 0 0 0
A 0 0 0 0 1 1 1
B 0 1
C 0
D 0
A 0
B 0

然後我們讓Y多出一個D做幫手,{“”,A,B,AB} vs {“”,B,D,BD},顯然,繼續填1. 一直填到Y的第二個B之前,都是1。 因為到BDCAB時,它們有另一個公共子序列,AB。

X “” B D C A B A
“” 0 0 0 0 0 0 0
A 0 0 0 0 1 1 1
B 0 1 1 1 1 2
C 0
D 0
A 0
B 0

到這一步,我們可以總結一些規則了,之後就是通過計算驗證我們的想法,加入新的規則或限定條件來完善。

LCS 演算法:Javascript 最長公共子序列Y將所有字元派上去,X依然是2個字元,經仔細觀察,還是填2.

看五行,X再多派一個C,ABC的子序列集合比AB的子序列集合大一些,那麼它與Y的B子序列集合大一些,就算不大,就不能比原來的小。顯然新增的C不能成為戰力,不是兩者的公共字元,因此值應該等於AB的子序列集合。

× “” B D C A B A
“” 0 0 0 0 0 0 0
A 0 0 0 0 1 1 1
B 0 1 1 1 1 2 2
C 0 1
D 0
A 0
B 0

LCS 演算法:Javascript 最長公共子序列

並且我們可以確定,如果兩個字串要比較的字元不一樣,那麼要填的格子是與其左邊或上邊有關,那邊大就取那個。

如果比較的字元一樣呢,稍安毋躁,剛好X的C要與Y的C進行比較,即ABC的子序列集合{“”,A,B,C,AB,BC,ABC}與BDC的子序列集合{””,B,D,C,BD,DC,BDC}比較,得到公共子串有“”,B,D 。這時還是與之前的結論一樣,當字元相等時,它對應的格子值等於左邊與右邊與左上角的值,並且左邊,上邊,左上邊總是相等的。這些奧祕需要更嚴格的數學知識來論證。

我們現在用下面的方程來繼續填表了。

image_1c3vdljfs1g5p1s5lpppcln1s88m

程式實現

LCS可以進一步簡化,只要通過挪位置,省去新陣列的生成

列印一個LCS

我們再給出列印函式,先看如何列印一個。我們從右下角開始尋找,一直找到最上一行終止。因此目標字串的構建是倒序。為了避免使用stringBuffer這樣麻煩的中間量,我們可以通過遞迴實現,每次執行程式時,只返回一個字串,沒有則返回一個空字串, 以printLCS(x,y,...) + str[i]相加,就可以得到我們要求的字串。

我們再寫出一個方法,來驗證我們得到的字串是否真正的LCS字串。作為一個已經工作的人,不能寫的程式碼像在校生那樣,不做單元測試就放到線上讓別人踩坑。

使用:

image_1c3v0n1qnhq1sqa36min2goi16

列印全部LCS

思路與上面差不多,我們注意一下,在LCS方法有一個Math.max取值,這其實是整合了三種情況,因此可以分叉出三個字串。我們的方法將返回一個es6集合物件,方便自動去掉。然後每次都用新的集合合併舊的集合的字任串。

使用:

image_1c3v1al221m5l1cdlmui48vd831j

空間優化

使用滾動陣列:

危險的遞迴解法

str1的一個子序列相應於下標序列{1, 2, …, m}的一個子序列,因此,str1共有${2^m}$個不同子序列(str2亦如此,如為${2^n}$),因此複雜度達到驚人的指數時間(${2^m * 2^n}$)。

參考連結

  1. http://blog.csdn.net/hrn1216/…
  2. https://segmentfault.com/a/11…
  3. https://www.cnblogs.com/ider/…
  4. http://www.cppblog.com/mysile…

相關文章