KMP算法

最新推荐文章于 2023-06-28 23:43:43 发布
最新推荐文章于 2023-06-28 23:43:43 发布
阅读量387 收藏
点赞数
分类专栏： 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43845922/article/details/119295005
版权

KMP

预备知识

字符串的前缀，后缀和部分匹配值

前缀：指除最后一个字符以外，字符串的所有头部子串。
后缀：指除了第一个字符外，字符串的所有尾部子串。
部分匹配值：则为字符串的前缀和后缀的最长相等前后缀长度。

以 ababa 为例说明

  1. 'a’的前后缀均为空集，最长相等前后缀长度为0
    'ab’的前缀为{a}, 后缀为{b}，{a}∩{b}=∅，最长相等前后缀长度为0
  2. 'aba’的前缀为{a,ab}, 后缀为{ba,a}, {a,ab}∩{ba,a}={a},最长相等前后缀长度为1；
  3. 'abab’的前缀为{a,ab,aba}, 后缀为{bab,ab,b}, {a,ab,aba}∩{bab,ab,b}={ab},最长相等前后缀长度为2
  4. 'ababa’的前缀为{a,ab,aba,abab}, 后缀为{baba,aba,ba,a}, {a,ab,aba,abab}∩{baba,aba,ba,a}={a, aba},最长相等前后缀长度为3（即aba的长度）
  5. next[j] = 字符串S1…j-1的最长相等前后缀的长度 + 1

手算过程：

下面只是举例说明过程，按具体的某个next[ j ]来求解，其他位置的求解也类似。KMP的next值是自己和自己一位一位地移动后算出的每一位的next值，无需其他字符串。 规定 next[0]=0, next[1]=1, 其他位置可以按下面方法计算。
Eg1
在这里插入图片描述
求next[5]=1；第5个元素即e用红色标出
依次将串右移一位，因为假设第5位不匹配，则至少要右移一次。
下图中一共右移了4次，直到串的首字母 a 与 e同一列，说明此时在此过程中没有任何一个子串与 绿色的 bcd或bcd的后缀 有完全重合的
所以下一次则应和串的第一个字母比较，即得 next[5]=1
在这里插入图片描述
Eg2
在这里插入图片描述
求next[5]=2；第5个元素即 b 用红色标出
依次将串右移一位，因为假设第5位不匹配，则至少要右移一次。
下图中一共右移了3次，直到有一个子串与 绿色的 bca或bcd的后缀 有完全重合的, 也即第3次移动后，绿色的a和绿色的bca的后缀 a 完全重合，而第1,2次的粉色的均不与 绿色的 bcd或bcd的后缀 有完全重合的
所以下一次则应和串的第二个字母比较，即得 next[5]=2

这里因为第三次移动后就找找了，所以第四次的移动就不需要,下图中显示出来只是为了说明问题。
在这里插入图片描述
Eg3
在这里插入图片描述
求next[6]=4；第6个元素即 a 用红色标出
依次将串右移一位，因为假设第6位不匹配，则至少要右移一次。
下图中一共右移了2次，直到有一个子串与 绿色的 baba或baba的后缀 有完全重合的, 也即第2次移动后，绿色的aba和绿色的baba的后缀 aba 完全重合，而第1次的粉色的不与 绿色的 bcd 完全重合
所以下一次则应和串的第4个字母 b 比较，即得 next[5]=4

这里因为第2次移动后就找找了，所以余下几次的移动就不需要,下图中显示出来只是为了说明问题。

在这里插入图片描述
Eg4
在这里插入图片描述
求next[8]=7；第8个元素即 a 用红色标出
依次将串右移一位，因为假设第8位不匹配，则至少要右移一次。
下图中一共右移了1次，直到有一个子串与 绿色的 aaaaaa或aaaaaa的后缀 有完全重合的, 也即第1次移动后，绿色的 aaaaaa 和 绿色的 aaaaaa 完全重合
所以下一次则应和串的第二个字母比较，即得 next[8]=7

这里因为第1次就找找了，所以余下几次的移动就不需要，下图中显示出来只是为了说明问题。
在这里插入图片描述

KMP算法的实现

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

//串的顺序存储  从数组下标为1开始存储字符
#define MAXLEN 255//预定义最大串长为255
typedef struct{
    char ch[MAXLEN];//每个分量存储一个字符
    int length;//串的实际长度
}SString;

//KMP算法求next值
void get_next(SString T, int next[]){
    int i=1,j=0;//i依次指向 T 的每一位,
    //j表示假设T.ch[i]不匹配，则下一次应和串的第j个字母比较
    /*这里不容易理解：
    只有满足T.ch[i]==T.ch[j]，(为了说明问题，记m=++i; n=++j)
    而后经过++i,++j,  才会执行到next[m]=n
    说明前面的n-1个字符是完全重合的，
    则当下一次T和其他所给的主串S进行匹配时，
    若到m时出现不匹配情况
    则下一次应和串T的第n个字符比较，即next[m]=n;
    */
    next[1]=0;//规定next[1]为0
    while(i<T.length){//@@@
        if(j==0||T.ch[i]==T.ch[j]){//
            ++i;
            ++j;
            next[i]=j;
        }
        else{
            j=next[j];
        }
    }
}

void PrintString(SString S){
    printf("\nString is:");
    for(int i=1;i<=S.length;i++){
        printf("%c",S.ch[i]);
    }
    printf("\n");
}

int Index_KMP(SString S,SString T){
    int i=1,j=1;
    int next[255];
    get_next(T,next);
    printf("\nnext array is:\n");
    for(int k=1;k <= T.length;k++){//next数组从1开始存数
        printf("%d ",next[k]);
    }
    printf("\n");
    while(i<=S.length && j<=T.length){
        if(j==0 || S.ch[i]==T.ch[j]){
        //j=0说明此时S.ch[i]和模式串T的任何一个字符都不匹配，所以i,j都要向后移
            ++i;++j;//继续比较后继字符
        }
        else{
            j=next[j];//S.ch[i]!=T.ch[j]说明下一次要和T.ch[ next[j] ]进行比较
        }
    }
    if(j>T.length)
        return i-T.length;//匹配成功
    else
        return 0;//匹配失败
}


int main(){
    SString S1,S2;
    //char str1[40]="hello?How are you world";
    //char str2[40]="How are you";
    char str1[40]="aaabaaaab";
    char str2[40]="aaaab";

    for(int i=0;i<strlen(str1);i++){
        S1.ch[i+1]=str1[i];
    }
    S1.length=strlen(str1);


    for(int i=0;i<strlen(str2);i++){
        S2.ch[i+1]=str2[i];
    }
    S2.length=strlen(str2);
    PrintString(S1);
    PrintString(S2);
    int x = Index_KMP(S1,S2);
    printf("\nindex is %d\n",x);
    return 0;
}

在这里插入图片描述
在这里插入图片描述

KMP算法的进一步优化


上图中S[4]和T[4]失配，如果用next数组，因为next[4]=3,则下一个需要用T[3]和S[4]比较，而T[3]为a，也失配，同理于next[3],next[2],next[1], 可以发现T的前三个字符和当前失配的a字符即T[4]是相同的，则再和S[4]相比较，必然会导致继续失配，这样的比较毫无意义，所以就产生了nextval数组。

如果出现上述情况即T[j]=T[ next[j] ],则需要再次递归，将next[j]修正为next[ next[j] ],
直到两者不相等为止，更新后的数组命为nextval.

Eg1
在这里插入图片描述

  1. j=2, T[2]=b, next[2]=1; b!= T[next[2]] 即(T[1] : a), 所以next[2]不变
  2. j=3, T[3]=a, next[3]=1; a = T[next[3]] 即(T[1] : a), 所以next[3]应修改为 next[next[3]]=next[1]=0
  3. j=4, T[4]=b, next[4]=2; b = T[next[4]] 即(T[2] : b), 所以next[4]应修改为 next[next[4]]=next[2]=1
  4. j=5, T[5]=a, next[5]=3; b = T[next[5]] 即(T[3] : b),所以next[5]应修改为next[next[5]]=next[3], 而next[3]在第二步中已经修改为0，故next[5]也修改为最新的next[3]即0
  5. j=6, T[6]=a, next[6]=4; a != T[next[6]] 即(T[4] : b), 所以next[6]不变
  6. 同理next[7]也不变

Eg2
在这里插入图片描述
下面代码可以运行出上面两个示例的 next 和 nextval 结果

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

//串的顺序存储  从数组下标为1开始存储字符
#define MAXLEN 255//预定义最大串长为255
typedef struct{
    char ch[MAXLEN];//每个分量存储一个字符
    int length;//串的实际长度
}SString;

//KMP算法求next值
void get_next(SString T, int next[]){
    int i=1,j=0;//i依次指向 T 的每一位,
    //j表示假设T.ch[i]不匹配，则下一次应和串的第j个字母比较
    /*这里不容易理解：
    只有满足T.ch[i]==T.ch[j]，(为了说明问题，记m=++i; n=++j)
    而后经过++i,++j,  才会执行到next[m]=n
    说明前面的n-1个字符是完全重合的，
    则当下一次T和其他所给的主串S进行匹配时，
    若到m时出现不匹配情况
    则下一次应和串T的第n个字符比较，即next[m]=n;
    */
    next[1]=0;//规定next[1]为0
    while(i<T.length){//@@@
        if(j==0||T.ch[i]==T.ch[j]){//
            ++i;
            ++j;
            next[i]=j;
        }
        else{
            j=next[j];
        }
    }
}



void PrintString(SString S){
    printf("\nString is:");
    for(int i=1;i<=S.length;i++){
        printf("%c",S.ch[i]);
    }
    printf("\n");
}

//KMP算法求nextval值
void get_nextval(SString T, int nextval[]){
    int i=1,j=0;//i依次指向 T 的每一位,
    nextval[1]=0;//规定next[1]为0
    while(i<T.length){//@@@
        if(j==0||T.ch[i]==T.ch[j]){//
            ++i;
            ++j;
            if(T.ch[i]!=T.ch[j]){
                nextval[i]=j;
            }
            else{//T.ch[i]==T.ch[j]
//则若遇到T.ch[i]!=S.ch[i']时，也必有T.ch[j]!=S.ch[i'],
/*求next[i]就是为了T和S在匹配过程中发生了i之前的都匹配
而T.ch[i]!=S.ch[i']，为了避免 i’ 回溯，所以要看next[i]的值，
而当T.ch[i]==T.ch[j]，也必有T.ch[j]!=S.ch[i'],
也即有nextval[i]=nextval[j]
*/
                nextval[i]=nextval[j];
            }
        }
        else{
            j=nextval[j];
        }
    }
}

int Index_KMP(SString S,SString T){
    int i=1,j=1;
    int nextval[255];


    int next[255];
    get_next(T,next);
    printf("\nnext array is:\n");
    for(int k=1;k <= T.length;k++){//next数组从1开始存数
        printf("%d ",next[k]);
    }
    printf("\n");


    get_nextval(T,nextval);
    printf("\nnextval array is:\n");
    for(int k=1;k <= T.length;k++){//next数组从1开始存数
        printf("%d ",nextval[k]);
    }
    printf("\n");

    while(i<=S.length && j<=T.length){
        if(j==0 || S.ch[i]==T.ch[j]){
        //j=0说明此时S.ch[i]和模式串T的任何一个字符都不匹配，所以i,j都要向后移
            ++i;++j;//继续比较后继字符
        }
        else{
            j=nextval[j];//S.ch[i]!=T.ch[j]说明下一次要和T.ch[ next[j] ]进行比较
        }
    }
    if(j>T.length)
        return i-T.length;//匹配成功
    else
        return 0;//匹配失败
}


int main(){
    SString S1,S2;
    //char str1[40]="hello?How are you world";
    //char str2[40]="How are you";
    char str1[40]="aaabaaaab";
    //char str2[40]="aaaab";
    //char str2[40]="ababaaaba";
    char str2[40]="aaaaaaaab";
    for(int i=0;i<strlen(str1);i++){
        S1.ch[i+1]=str1[i];
    }
    S1.length=strlen(str1);


    for(int i=0;i<strlen(str2);i++){
        S2.ch[i+1]=str2[i];
    }
    S2.length=strlen(str2);
   // PrintString(S1);
    PrintString(S2);
    int x = Index_KMP(S1,S2);
   // printf("\nindex is %d\n",x);
    return 0;
}

对应改进KMP算法 示例中的Eg1
对应改进KMP算法 示例中的Eg2

为啥不能修改昵称啊
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
KMP算法KMP算法KMP算法KMP算法
07-19
KMP算法
KMP 算法实例详解
01-01
KMP算法，是由Knuth，Morris，Pratt共同提出的模式匹配算法，其对于任何模式和目标序列，都可以在线性时间内完成匹配查找，而不会发生退化，是一个非常优秀的模式匹配算法。 分析：KMP模板题、KMP的关键是求出next的...
KMP 算法详解
Hello World!
04-22 1万+
KMP算法是一种高效的字符串匹配算法，算法名称取自于三位共同发明人名字的首字母组合。该算法的主要使用场景就是在字符串（也叫主串）中的模式串（也叫字串）定位问题，常见的有“求子串出现的起始位置”、“求子串的出现次数”等。
详解kmp算法
qq_51950769的博客
02-05 6566
1.朴素模式匹配 1.1基本思想： 从主串test 和子串pattern 的第一个字符开始，将两字符串的字符一一比对，如果出现某个字符不匹配，主串回溯到第二个字符，子串回溯到第一个字符再进行一一比对。如果出现某个字符不匹配，主串回溯到第三个字符，子串回溯到第一个字符再进行一一比对…循环一直到子串字符全部匹配成功。 1.2例如： 第一次：不匹配 第二次：主串从第2个字符，模式串从头开始，直到b与c不匹配 第三次：主串从第3个位置开始，模式串从头开始，a与c...
图解KMP算法，带你彻底吃透KMP
qq_43869106的博客
01-23 1万+
KMP算法一直是一个比较难以理解的算法，本篇文章主要根据《大话数据结构》中关于KMP算法的讲解，结合自己的思考，对于KMP算法进行一个比较详细的解释。
（详细图解）KMP算法（C语言）------可算领悟了这个由三位前辈研究的算法
在夜里行走的博客
10-23 4403
字数可能较多，希望能看下去
kmp算法演示
11-20
kmp算法求Next值演示，可以自己多研究一下，上面详细的讲解了KMP算法求Next值的具体方法
kmp算法分析
06-12
kmp算法分析，针对kmp算法进行了简要分析，对于初学者能有一定帮助，如果内容有误请联系我更改。
史上最详细的KMP算法教程，看这一篇就够了
⚡
02-04 8627
前缀函数、next数组、nextval数组，BF算法，KMP算法，你想要知道的都在这了！
数据结构KMP算法详解
最新发布
强哥的微博
06-28 2539
KMP算法的改进可以简述为： 如果a位字符与它next值指向的b位字符相等，则该a位的nextval就指向b位的nextval值，如果不等，则该a位的nextval值就是它自己a位的next值。我们一开始想到的是暴力求解，我们是将子串和主串逐一匹配，如果第一个字符相等就继续匹配第二个字符，直到子串与主串全都匹配成功，就返回子串的位置，一旦其中某两个字符匹配不成功，主串就回到开始匹配字符的下一字符，而子串回到第一字符。这个例子中当‘b’与‘c’不匹配时应该‘b’与’c’前一位的‘a’比,这显然是不匹配的。
KMP算法详解及实现
Koikoi12的博客
11-04 969
KMP算法的核心，是一个被称为部分匹配表(Partial Match Table)的数组。对于字符串“abababca”，它的PMT如下表所示： 字符串的前缀和后缀定义： 如果字符串A和B，存在A=BS，其中S是任意的非空字符串，那就称B为A的前缀。例如，”Harry”的前缀包括{”H”, ”Ha”, ”Har”, ”Harr”}，我们把所有前缀组成的集合，称为字符串的前缀集合。 同样可以定义后缀A=SB， 其中S是任意的非空字符串，那就称B为A的后缀，例如，”Potter”的后缀包括{”otter”, ”
KMP 算法
七秒钟的记忆
03-25 275
KMP 算法感觉还是用的动态规划啊
串文档之KMP
二木成林
05-15 480
字符串的 KMP 算法
leetcode 028.实现strStr()，即查找重复字符串（KMP算法）
学无止境,勤则可达；志存高远,恒亦能成！
04-20 1421
前言 本题是经典的字符串单模匹配的模型，因此可以使用字符串匹配算法解决，常见的字符串匹配算法包括暴力匹配、Knuth-Morris-Pratt 算法、Boyer-Moore 算法、Sunday 算法等，本文 ## 前言 本题是经典的字符串单模匹配的模型，因此可以使用字符串匹配算法解决，常见的字符串匹配算法包括暴力匹配、Knuth-Morris-Pratt 算法、Boyer-Moore 算法、Sunday 算法等，本文将讲解朴素解法（暴力匹配）和KMP算法 因为哈希方法可能出现哈希值相等但是字符串不相等的情况
KMP算法最浅显理解——一看就明白
热门推荐
好记性不如烂笔头的专栏
02-07 31万+
说明 KMP算法看懂了觉得特别简单，思路很简单，看不懂之前，查各种资料，看的稀里糊涂，即使网上最简单的解释，依然看的稀里糊涂。 我花了半天时间，争取用最短的篇幅大致搞明白这玩意到底是啥。 这里不扯概念，只讲算法过程和代码理解： KMP算法求解什么类型问题 字符串匹配。给你两个字符串，寻找其中一个字符串是否包含另一个字符串，如果包含，返回包含的起始位置。 如下面两个字符串： c...
[数据结构与算法-21]KMP算法
学哥斌的博客
04-21 658
文章目录1、应用场景2、暴力匹配算法2.1 案例展示2.2 时间复杂度分析2.3 代码实现2.4 KMP算法对暴力匹配的改进3、 KMP算法3.1 算法原理 1、应用场景   Knuth-Morris-Pratt 字符串查找算法，简称为 KMP算法，常用于在一个文本串 S 内查找一个模式串 P 的出现位置。是用来代替暴力匹配算法的有效算法。   如果找到，即返回第一个字符的索引。 2、暴力匹配算...
使用kmp算法
05-29
好的，下面我来为您讲解使用KMP算法解决字符串匹配问题的方法。 KMP算法是一种比较高效的字符串匹配算法，它的核心思想是利用已经匹配过的信息，尽可能减少不必要的字符比较次数，从而提高匹配效率。具体实现时，...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值