青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

KMP和擴展KMP

Posted on 2011-04-17 19:11 Mato_No1 閱讀(5837) 評論(2)  編輯 收藏 引用 所屬分類: 字符串匹配
KMP:給出兩個字符串A(稱為模板串)和B(稱為子串),長度分別為lenA和lenB,要求在線性時間內,對于每個A[i](0<=i<lenA),求出A[i]往前和B的前綴匹配的最大匹配長度,記為ex[i](或者說,ex[i]為滿足A[i-z+1..i]==B[0..z-1]的最大的z值)。KMP的主要目的是求B是不是A的子串,以及若是,B在A中所有出現的位置(當ex[i]=lenB時)。
【算法】
設next[i]為滿足B[i-z+1..i]==B[0..z-1]的最大的z值(也就是B的自身匹配)。設目前next[0..lenB-1]與ex[0..i-1]均已求出,要用它們來求ex[i]的值。
根據ex的定義,有A[i-1-ex[i-1]+1..i-1]==B[0..ex[i-1]-1],這時,若有A[i]==B[ex[i-1]],則可以直接得到ex[i]=ex[i-1]+1(因為i-1-ex[i-1]+1即i-ex[i-1],現在由于A[i]==B[ex[i-1]],可得A[i-ex[i-1]..i]==B[0..ex[i-1]],即A[i-ex[i-1]+1-1..i]==B[0..ex[i-1]+1-1],所以ex[i]=ex[i-1]+1)。若A[i]!=B[ex[i-1]]?
設j=next[ex[i-1]-1],則根據next定義得B[ex[i-1]-j..ex[i-1]-1]==B[0..j-1],又因為A[i-ex[i-1]..i-1]==B[0..ex[i-1]-1]得A[i-j..i-1]==B[ex[i-1]-j..ex[i-1]-1],這樣有A[i-j..i-1]==B[0..j-1]!也就是此時只需再比較A[i]與B[j]的值是否相等即可,若相等,可得ex[i]=j+1,若仍不相等,則更新j為next[j-1],繼續比較A[i]與B[j]是否相等……直到A[i]與B[j]相等或直到j==0時,A[i]仍不等于B[j],此時ex[i]=0。邊界:求ex[0]時,初始j(用來代替ex[i-1])為0。
現在還有一個問題,如何求next?顯然next就是以B自身為模板串,B為子串的“自身匹配”,用類似的辦法即可,唯一不同的是next[0]=lenB可以直接得到,求next[1]時,初始j(代替next[i-1])為0。
【核心代碼】
    lenA = strlen(A); lenB = strlen(B);
    next[
0= lenB;
    
int j = 0;
    re2(i, 
1, lenB) {
        
while (j && B[i] != B[j]) j = next[j - 1];
        
if (B[i] == B[j]) j++;
        next[i] 
= j;
    }
    j 
= 0;
    re(i, lenA) {
        
while (j && A[i] != B[j]) j = next[j - 1];
        
if (A[i] == B[j]) j++;
        ex[i] 
= j;
    }
擴展KMP:給出模板串A和子串B,長度分別為lenA和lenB,要求在線性時間內,對于每個A[i](0<=i<lenA),求出A[i..lenA-1]與B的最長公共前綴長度,記為ex[i](或者說,ex[i]為滿足A[i..i+z-1]==B[0..z-1]的最大的z值)。擴展KMP可以用來解決很多字符串問題,如求一個字符串的最長回文子串和最長重復子串。
【算法】
設next[i]為滿足B[i..i+z-1]==B[0..z-1]的最大的z值(也就是B的自身匹配)。設目前next[0..lenB-1]與ex[0..i-1]均已求出,要用它們來求ex[i]的值。
設p為目前A串中匹配到的最遠位置,k為讓其匹配到最遠位置的值(或者說,k是在0<=i0<i的所有i0值中,使i0+ex[i0]-1的值最大的一個,p為這個最大值,即k+ex[k]-1),顯然,p之后的所有位都是未知的,也就是目前還無法知道A[p+1..lenA-1]中的任何一位和B的任何一位是否相等。
根據ex的定義可得,A[k..p]==B[0..p-k],因為i>k,所以又有A[i..p]==B[i-k..p-k],設L=next[i-k],則根據next的定義有B[0..L-1]==B[i-k..i-k+L-1]。考慮i-k+L-1與p-k的關系:
(1)i-k+L-1<p-k,即i+L<=p。這時,由A[i..p]==B[i-k..p-k]可以得到A[i..i+L-1]==B[i-k..i-k+L-1],又因為B[0..L-1]==B[i-k..i-k+L-1]所以A[i..i+L-1]==B[0..L-1],這就說明ex[i]>=L。又由于next的定義可得,A[i+L]必然不等于B[L](否則A[i..i+L]==B[0..L],因為i+L<=p,所以A[i..i+L]==B[i-k..i-k+L],這樣B[0..L]==B[i-k..i-k+L],故next[i-k]的值應為L+1或更大),這樣,可以直接得到ex[i]=L!
(2)i+k-L+1>=p-k,即i+L>p。這時,首先可以知道A[i..p]和B[0..p-i]是相等的(因為A[i..p]==B[i-k..p-k],而i+k-L+1>=p-k,由B[0..L-1]==B[i-k..i-k+L-1]可得B[0..p-i]==B[i-k..p-k],即A[i..p]==B[0..p-i]),然后,對于A[p+1]和B[p-i+1]是否相等,目前是不知道的(因為前面已經說過,p是目前A串中匹配到的最遠位置,在p之后無法知道任何一位的匹配信息),因此,要從A[p+1]與B[p-i+1]開始往后繼續匹配(設j為目前B的匹配位置的下標,一開始j=p-i+1,每次比較A[i+j]與B[j]是否相等,直到不相等或者越界為止,此時的j值就是ex[i]的值)。在這種情況下,p的值必然會得到延伸,因此更新k和p的值。
邊界:ex[0]的值需要預先求出,然后將初始的k設為0,p設為ex[0]-1。
對于求next數組,也是“自身匹配”,類似KMP的方法處理即可。唯一的不同點也在邊界上:可以直接知道next[0]=lenB,next[1]的值預先求出,然后初始k=1,p=ex[1]。

需要嚴重注意的是,在上述的情況(2)中,本該從A[p+1]與B[p-i+1]開始匹配,但是,若p+1<i,也就是p-i+1<0(這種情況是有可能發生的,當ex[i-1]=0,且前面的ex值都沒有延伸到i及以后的時候)的話,需要將A、B的下標都加1(因為此時p必然等于i-2,如果A、B的下標用兩個變量x、y控制的話,x和y都要加1)!!

【核心代碼】
lenA = strlen(A); lenB = strlen(B);
    next[
0= lenB; next[1= lenB - 1;
    re(i, lenB
-1if (B[i] != B[i + 1]) {next[1= i; break;}
    
int j, k = 1, p, L;
    re2(i, 
2, lenB) {
        p 
= k + next[k] - 1; L = next[i - k];
        
if (i + L <= p) next[i] = L; else {
            j 
= p - i + 1;
            
if (j < 0) j = 0;
            
while (i + j < lenB && B[i + j] == B[j]) j++;
            next[i] 
= j; k = i;
        }
    }
    
int minlen = lenA <= lenB ? lenA : lenB; ex[0= minlen;
    re(i, minlen) 
if (A[i] != B[i]) {ex[0= i; break;}
    k 
= 0;
    re2(i, 
1, lenA) {
        p 
= k + ex[k] - 1; L = next[i - k];
        
if (i + L <= p) ex[i] = L; else {
            j 
= p - i + 1;
            
if (j < 0) j = 0;
            
while (i + j < lenA && j < lenB && A[i + j] == B[j]) j++;
            ex[i] 
= j; k = i;
        }
    }
【時間復雜度分析】
在KMP和擴展KMP中,不管是A串還是B串,其匹配位置都是單調遞增的,故總時間復雜度是線性的,都為O(lenA + lenB)(只是擴展KMP比KMP的常數更大一些)。
【應用】
KMP和擴展KMP在解決字符串問題中有大用。很多看上去很猥瑣的字符串問題,都可以歸結到這兩種算法之中。另外,這里的“字符串”可以延伸為一切類型的數組,而不僅僅是字符數組。

Feedback

# re: KMP和擴展KMP  回復  更多評論   

2011-12-14 22:34 by ~~
KMP算法里面next[0]應該是0才對吧?
next[i]為滿足B[i-z+1..i]==B[0..z-1]的最大的z值(也就是B的自身匹配)
對于B[0],i- 1是-1,不能向前延伸了~

# re: KMP和擴展KMP  回復  更多評論   

2015-08-02 14:21 by yajun
擴展KMP
if (i + L <= p) ex[i] = L; else {
相等的時候,需要向后匹配,因為可能會比L更長。
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            欧美一区2区视频在线观看 | 中日韩美女免费视频网址在线观看 | 亚洲精品一二三| 久久青草福利网站| 欧美在线看片| 老司机久久99久久精品播放免费| 久久婷婷色综合| 欧美黄色一区| 亚洲免费成人av| 亚洲午夜日本在线观看| 一区二区三区欧美日韩| 亚洲视频免费| 欧美中文字幕在线视频| 久久中文欧美| 欧美精品七区| 国产精品你懂得| 狠狠爱www人成狠狠爱综合网| 韩国av一区二区三区| 亚洲欧洲中文日韩久久av乱码| 亚洲精品国产系列| 亚洲一区影音先锋| 久久深夜福利| 亚洲精品老司机| 99国产精品99久久久久久| 亚洲视频观看| 久久婷婷国产综合国色天香| 欧美大香线蕉线伊人久久国产精品| 欧美日韩在线综合| 在线 亚洲欧美在线综合一区| av不卡在线看| 久久嫩草精品久久久精品一| 亚洲人成艺术| 久久夜精品va视频免费观看| 欧美小视频在线观看| 狠狠色综合一区二区| 亚洲天天影视| 欧美国产日韩在线| 亚洲欧美高清| 欧美精品亚洲一区二区在线播放| 久久久精品国产免费观看同学| 欧美高清在线一区二区| 国产一区二区三区精品欧美日韩一区二区三区 | 日韩天堂在线观看| 午夜精品一区二区三区在线视 | 激情亚洲网站| 一区二区三区欧美| 亚洲福利在线视频| 久久国产88| 欧美 日韩 国产一区二区在线视频 | 亚洲一区精品在线| 欧美夜福利tv在线| 91久久精品日日躁夜夜躁国产| 欧美久久久久久久| 久久久久一区二区三区四区| 亚洲区在线播放| 欧美专区在线| 亚洲人成精品久久久久| 欧美在线观看日本一区| 国产精品推荐精品| 亚洲一区二区在线视频| 亚洲毛片av| 欧美日韩成人综合| 亚洲精品国产精品乱码不99按摩 | 99热这里只有精品8| 美女性感视频久久久| 午夜精品久久久久影视| 国产精品影片在线观看| 午夜久久福利| 亚洲深夜福利视频| 国产精品久久久久久久免费软件 | 久久综合狠狠综合久久激情| 亚洲欧美清纯在线制服| 国产日韩欧美二区| 欧美日韩国产一区二区三区| 亚洲国产欧美国产综合一区| 亚洲日本在线视频观看| 欧美在线一区二区| 亚洲日韩成人| 中文国产一区| 国产一区二区三区四区五区美女 | 欧美精品乱码久久久久久按摩| 久久er精品视频| 欧美岛国激情| 欧美乱大交xxxxx| 欧美特黄一级| 亚洲午夜精品一区二区三区他趣| 欧美成人精品在线视频| 亚洲伦理久久| 亚洲欧美一区二区视频| 91久久精品国产91性色tv| 亚洲午夜精品在线| 国产精品免费一区豆花| 欧美在线观看网站| 久久久天天操| 99pao成人国产永久免费视频| 亚洲精品乱码久久久久久黑人| 欧美揉bbbbb揉bbbbb| 久久九九精品| 欧美成在线观看| 欧美一区二区三区日韩| 久久午夜精品| 亚洲一级黄色| 免费看的黄色欧美网站| 亚洲字幕在线观看| 久久综合国产精品台湾中文娱乐网| 亚洲精品一区二区网址| 亚洲综合成人在线| 91久久夜色精品国产九色| 亚洲免费观看视频| 怡红院精品视频| 亚洲性人人天天夜夜摸| 亚洲福利视频三区| 亚洲免费影视| 夜夜嗨av一区二区三区| 久久99伊人| 亚洲午夜伦理| 免费久久精品视频| 午夜精品在线视频| 欧美久久久久久久| 久久躁日日躁aaaaxxxx| 欧美三级不卡| 欧美大片免费观看在线观看网站推荐| 国产精品家教| 日韩特黄影片| 亚洲第一综合天堂另类专| 亚洲特级片在线| 一本色道久久综合亚洲精品不| 久久免费少妇高潮久久精品99| 午夜精品久久久久久久久| 欧美日韩国产成人| 亚洲国产精品ⅴa在线观看 | 欧美高清你懂得| 久久婷婷影院| 国语自产偷拍精品视频偷| 亚洲欧美成人一区二区在线电影| 最近中文字幕日韩精品 | 久久综合久久综合久久| 免费观看一区| 久久久久久穴| 欧美激情久久久久久| 亚洲午夜国产一区99re久久| 欧美亚洲午夜视频在线观看| 久久久在线视频| 欧美一区二区在线看| 亚洲七七久久综合桃花剧情介绍| 一区二区精品在线| 久久九九国产精品怡红院| 亚洲欧美美女| 亚洲片在线观看| 亚洲国产成人tv| 亚洲精品日产精品乱码不卡| 另类激情亚洲| 亚洲狼人综合| 欧美激情一区| 久久精品国产综合精品| 亚洲女同在线| 亚洲香蕉成视频在线观看| 欧美大片在线观看| 黄色日韩在线| 91久久国产综合久久蜜月精品| 亚洲一区二区三区中文字幕| 亚洲人成艺术| 国产三级欧美三级日产三级99| 国产女主播一区二区| 国产综合激情| 国产一区免费视频| 欧美日韩1080p| 亚洲国产一区二区a毛片| 老鸭窝91久久精品色噜噜导演| 亚洲韩国日本中文字幕| 加勒比av一区二区| 国内成人在线| 久久精品视频在线| 国产精品九九| 久久激情网站| 亚洲欧美日本国产专区一区| 欧美日韩国产一级| 欧美一级播放| 欧美a级一区| 亚洲一区观看| 国产午夜久久久久| 欧美精品久久久久久久久久| 亚洲香蕉成视频在线观看| 久久人体大胆视频| 在线一区二区三区四区| 国产精品素人视频| 久久久久综合网| 日韩视频免费看| 久久久久久9| 日韩视频免费观看高清完整版| 欧美日韩国产色站一区二区三区| 一区二区三区波多野结衣在线观看| 久久国产综合精品| 亚洲午夜一二三区视频| 国产一区在线看| 欧美久久婷婷综合色| 欧美一级淫片播放口| 宅男噜噜噜66国产日韩在线观看| 影音先锋一区| 国产裸体写真av一区二区|