exkmp略解
阿新 • • 發佈:2017-11-27
最長 因此 更新 code () 並且 mark 我們 -a
推導
ext[i]
表示母串s[i..lens]
和子串t[1..lent]
的最長公共前綴。
nxt[i]
表示t[i..lent]
和t[1..lent]
的最長公共前綴。
假設ext[1..k]
已經算好,現在想求ext[k]
。
我們先找出1..k
中,i+ext[i]-1
的最大值p
,並且記此時的i
為a
。
可知s[a..p]=t[1..p-a+1]
,自然有s[k+1..p]=t[k-a+2..p-a+1]
,記l=nxt[k-a+2]
。
出現兩種情況:
k+l<p
可想而知,ext[k+1]=l
k+l>=p
這樣會發現有未知情況,因此就再度匹配,直至適配,然後更新a
。
int lena, lenb, nxt[100005 ], ext[100005];
char a[100005], b[100005];
void mknxt(){
int j=0, po;
nxt[0] = lenb;
while(j+1<lenb && b[j]==b[j+1]) j++;
nxt[1] = j;
po = 1;
for(int i=2; i<len; i++){
if(i+nxt[i-po]<nxt[po]+po)
nxt[i] = nxt[i-po];
else{
j = max(nxt[po] + po - i, 0 );
while(i+j<lenb && b[i+j]==b[j]) j++;
nxt[i] = j;
po = i;
}
}
}
void exkmp(){
int j=0, po;
while(a[j]==b[j] && j<lena && j<lenb)
j++;
ext[0] = j;
po = 0;
for(int i=1; i<lena; i++){
if (nxt[i-po]+i<ext[po]+po)
ex[i] = nxt[i-po];
else{
j = max(ex[po] + po - i, 0);
while(i+j<lena && j<lenb && a[i+j]==b[j]) j++;
ext[i] = j;
po = i;
}
}
}
exkmp略解