首页 > 代码库 > 第四章:2.串 -- 串的模式匹配算法(KMP)

第四章:2.串 -- 串的模式匹配算法(KMP)

前言:  

 

目录:

  1.串类型的定义

  2.串的表示和实现

  3.串的模式匹配算法

  4.串操作应用举例

 

正文:

  串的模式匹配即,在给定主串S 中,搜索子串T 的位置,如果存在T 则返回其所在位置,否则返回 0

  串的模式匹配算法

  主串 S: a b c a b c d s v t

  子串 T: a b c d

 

 

  一、原始算法

    匹配一旦失败,子串即向右移动一个单位,直到完全匹配停止。

    第一次匹配:(注:红色代表不匹配(失配))

      S: a b c a b c a b c d s v t

        T: a b c d

    第二次匹配:

      S: a b c a b c d s v t

        T:    a

 

  代码实现:(在 上一节已经实现的堆分配存储结构的基础 上实现)    

//模式匹配(子串定位通常称为模式匹配)
//求主串S 第 pos 个字符之后,与模式串 T 相等的位置
int Index(HString S,HString T,int pos){
    int i=pos;
    int j=0;
    //没有到主串末尾
    while(i<S.length){
        if(S.ch[i]==T.ch[j]){

            //如果子串已经匹配结束,那么返回主串S 中开始匹配的位置
            if(j+1==T.length)        
                return i-j+1;

            //否则继续往后比较
            i++;
            j++;
        }else{
            //存在字符不相等,i指向主串S 当前开始匹配位置的下一位。 j 在子串中重头开始
            i=i-j+1;
            j=0;
        }
    }
    //匹配失败
    return 0;    
}

    运行结果:

      匹配位置: 7  、比较总次数:16次

  

  二、算法的第一次改进

    第一次匹配:

       下标: 0 1 2 3 

      S[i]: a b c a b c a b c d s v t

        T[j]: a b c d

    原始算法中,第一次匹配时,在i=3, j=3 时,a <> b。这时发现将 T 向右移动一个位置 和 自动两个位置都是没有必要的。这时我们在第二次匹配时,直接将子串 T 向右移动 3 个位置即可。

    第二次匹配:

      S: a b c a b c a b c d s v t

        T:         a b c d

    T之所以可以向右移动 3 个位置,是因为 在 已匹配的位置 a b c 中, b c 和 a 不相等,所有串不可能匹配起来。

    总结如下:

         当已匹配的串中(包含n个字符)

        1.除首个字符外,不存在字符和首字符相等,那么子串 T 应向右移动 n-1个位置

        2.如果存在 和首字符相等的元素 Tj,那么应该把子串 T 向右移动 j 个位置

    因此,一旦子串确定,我们就可以确定当子串T 中的 Tj 和主串S 中的Si 匹配失败时,子串应该移动的位置。此时,我们只需要把注意力放在子串T 上即可。

    以上面子串T 为例,设step[j]  为比较到T[j] 匹配失败时,子串T 应该向右移动的位置长度。

         下标j: 0 1 2 3

            T[j]: a b c d

        step[j]: 1 1 2 3

    综上所述,只要子串T 确定, step[j] 数组就可以被确定。那么推广到一般情况,我们通过函数构造 step数组的算法如下。

//构造step数组
int *next;
void InitStep(HString T){
    int *step;
    step=(int *)malloc(T.length*sizeof(int));
    if(!step) exit(OVERFLOW);
    step[0]=1;
    int equalPos=0;        //记录和首字符相等的字符首次出现的位置,为0说明不存在
    for(int i=1;i<T.length;i++){
        if(T.ch[i]==T.ch[0]){
            //记录出现位置
            if(equalPos==0)
                equalPos=i;

            step[i]=i;
        }else{
            if(equalPos==0){
                step[i]=i;
            }else{
                step[i]=equalPos;
            }
        }
    }
    next=step;
}

 

  此时,只需将 原始的模式匹配算法中 步长增1 改为 增step[j] 如下

     //存在字符不相等,i指向主串S 当前开始匹配位置的下一位。 j 在子串中重头开始

    i=i-j+1; (改为:i=i-j+next[j]; )

      j=0;

  算法即得到改进。

  运行结果:

    匹配位置: 7  、比较总次数:12次

  当主串数据量增大时,此改进后的算法的优势会越来越明显。

 

  三、算法的第二次改进即(KMP算法)

    假设此时的主串和子串为:

      S[i]: a b c r a b c a b c d s v t

        T[j]: a b c r a b c d

 

    第一次匹配:

       下标: 0 1 2 3 4 5 6 7

      S[i]: a  b c r a b c  s a b c r a b c d s v t

        T[j]: a  b c r a b c d  

 

    如果根据第一次改进的算法将会得到step:      

         下标j: 0 1 2 3 4 5 6 7

            T[j]:  a b c  r a b c d

        step[j]: 1 1 2 3 4 4 4 4

    可得第一次匹配失败时,step[7]=4,那么子串T 向右移动4 个长度的位置开始第二次匹配, 从主串的第4个下标开始比较,即 S[4]=a 开始比较。

 

    第二次匹配:

      下标: 0 1 2 3 4 5 6 7

      S:   a  b c r a b c  s a b c r a b c d s v t

        T:              a b c  r a b c d

 

    但是我们发现:

      下标4,5,6 对应的a,b,c。这三次比较也是没有必要进行的。

      即第二次匹配,我们仅仅需要令 S[7]=s 和 T[3]=r 开始比较即可。

    分析:

      之所以可以从 子串T 的第4 个字符开始进行比较,而不需要比较前3 个字符。是因为T[0,1,2] 分别等于 S[4,5,6],又因为在第一次匹配过程中,S[7]失配的时候,已匹配字符序列为:T[0,1,2,3,4,5,6] = S[0,1,2,3,4,5,6],可知:            

            T[0,1,2,3,4,5,6] = S[0,1,2,3,4,5,6] 可知 T[4,5,6] = S[4,5,6]

            又T[0,1,2] = S[4,5,6]

            因此有:T[0,1,2] = S[4,5,6] =T[4,5,6]

      在子串T 中: a b c r a b c d 中失配字符 d前面的子串sub:a b c r a b c,其前缀等于后缀(即两端相等)。

      总结,在 S[i]  、T[j] 的匹配过程中,当 S[i] 不等于 T[j] 时,在已匹配的串 T[0,1,...,j-1]中 其前缀等于后缀,那么下一次匹配只需要将其前缀和 上次后缀对其即可,然后从匹配前缀的下一位置和主串继续比较,S[7]=d 和 T[3]=r进行比较。那么当前 主串下标 i=7 ,子串下标 j=3。我们设此时的子串下标为 k。

      那么此时有:

          T[0,1,...,k-1] = T[j-k,...,j-1] 长度为 k 且 0<k<j-1

               下标: 0 1 2 3 4 5 6 7

                S[i]: a  b c r a b c  s a b c r a b c d s v t

           第一次匹配: T[j]: a  b c r a b c d 

        第二次匹配: T[j]:           a  b c r a b c d 

    这时,我们发现,只要子串确定,我们就可以确定当 T 在 j 出失配的时候,S[i] 应该继续和下标为k 的T 继续比较。

    推广到一般情况:

    设当 主串S 和子串T 失配的时候,下标分别为 i 和 j。

      1、当已匹配的串存在前缀 = 后缀,那么我们保持 i 不变,让S[i]  和 T[k]继续往下比较即可。

      2、当已匹配的串不存在 前缀 = 后缀的情况

             如果T[j] 不是T 的第一个字符,那么 S[i] 应该和 T 的第一个元素开始比较,此时 k=0

             如果T[j] 是T 的第一个字符, 那么 i 后移一位,和 T的第一元素开始比较,即S[++j] 和 T[0] 开始比较,我们规定此时 k = -1

     如上所述,我们设主串S[i]  和子串T[j]失配时,pos[j]=k。

     此时模式匹配算法实现如下:

int KMP(HString S,HString T,int position){
    int i=position;
    int j=0;
    o=0;
    //没有到主串末尾
    while(i<S.length){
        o++;
        if(S.ch[i]==T.ch[j]){

            //如果子串已经匹配结束,那么返回主串S 中开始匹配的位置
            if(j+1==T.length)        
                return i-j+1;

            //否则继续往后比较
            i++;
            j++;
        }else{
            //i,j失配
            if(pos[j]==-1){
                i++;
                j=0;
            }else{
                j=pos[j];
            }            
        }
    }
    //匹配失败
    return 0;    
}

 

    如上模式匹配算法已经实现,现在需要做的就是来实现 pos[j] 数组。

    pos:      

         下标j:   0 1 2 3 4 5 6 7

            T[j]:   a b c  r a b c d

         pos[j]: -1 0 0 0 0 1 2 3

    假设 pos[j]=k 成立,则对于已匹配部分:T[0,...,j-1] 有    

      如果T[j] = T[k],即T[0,...,k-1,k] = T[j-k,...,j-1,j]  那么 pos[j+1] = k+1 =pos[j] + 1

      如果T[j] <> T[k]

          此时可把求 pos 函数的问题看成是一个模式匹配问题,整个串及时主串又是子串,当前的匹配过程如下:

          T[j-k,...,j-1], T[j]

          T[0,......,k-1],T[k]

          则当 T[j] <> T[k]时,即让模式串右移至 第 pos[k] 个字符和主串中的第 j 个字符相比较。

          若 T[ pos[k] ] = T[j],那么 pos[j+1] = k+1 =pos[j] + 1pos[j+1]= pos[k]+1

          若 T[ pos[k] ] <> T[j] ,那么如上所述以此类推,直到两者相等,或者不存在子串 前后缀相匹配 即(pos[k]=-1或0)

    代码实现:

 

//构造pos数组
int *pos;
void InitPos(HString T){
    int *p;
    p=(int *)malloc(T.length*sizeof(int));
    if(!p) exit(OVERFLOW);
    p[0]=-1;
    int i=1;        
    int k=-1;        
    while(i<T.length){
        if(k==-1||T.ch[i-1]==T.ch[k]){            
            p[i]=k+1;
            k=p[i];
            ++i;
        }else{
            k=p[k];
        }
    }
    pos=p;
}

 

  运行结果: position=9,比较次数为 18次

  技术分享

  用原始Index 匹配算法测试当前的主串、子串,运行结果:position=9,比较次数为 26次

  用第一次的改进匹配算法测试当前的主串、子串,运行结果:position=9,比较次数为 21次

 

第四章:2.串 -- 串的模式匹配算法(KMP)