KMP学习笔记

最新推荐文章于 2022-05-11 21:59:10 发布
最新推荐文章于 2022-05-11 21:59:10 发布
阅读量179 收藏 1
点赞数 1
分类专栏： 左神笔记
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45045689/article/details/107431278
版权

KMP算法是一种字符串匹配算法，可以在 O(n+m) 的时间复杂度内实现两个字符串的匹配。

我觉得KMP关键就两点：1.实现在匹配的过程中，主串不回退。   2.用最长前缀 等于 最长后缀的性质快速求出next[]数组。

 

我们先看关键点一：为什么要实现，在匹配的过程中，主串不回退？？

用暴力来解决   洛谷3375 【模板】KMP字符串匹配 这题试试。

 

题目描述

给出两个字符串 s1​ 和 s2​，若 s1​ 的区间 [l, r] 子串与 s2​ 完全相同，则称 s2​ 在 s1​ 中出现了，其出现位置为 L。
现在请你求出 s2​ 在 s1​ 中所有出现的位置。

定义一个字符串 s 的 border 为 s 的一个非 s 本身的子串 t，满足 t 既是 s 的前缀，又是 s 的后缀。
对于 s2​，你还需要求出对于其每个前缀 s′ 的最长 border t′ 的长度。

 

输入样例：

ABABABC
ABA

输出样例：

1
3
0 0 1

用暴力解法就没有next数组这一说了。。

#include<string>
#include<iostream>
using namespace std;

void bruteForce(string S, string P)
{
	int lenS = S.length(), lenP = P.length();
	int flag;

	for (int i = 0; i <= lenS - lenP; i++){
		flag = true;

		for (int j = 0; j < lenP; j++){
			if (S[i + j] != P[j]){
				flag = false;
				break;
			}
		}//inner for
		if (flag){
			cout << i + 1 << endl;
		}//if

	}//for

}


int main()
{
	string S, P;
	cin >> S >> P;
	
	bruteForce(S, P);

	return 0;
}

时间复杂度为O(N * M) ，太慢了。

 

我们用next数组跳过 那些一定不成功的比较。那什么是next数组呢？？(next数组的定义在各教材中有所区别。)

定义一：

　next数组是对于模式串而言的。P 的 next 数组定义为：next[i] 表示 P[0] ~ P[i] 这一个子串，使得 前k个字符恰等于后k个字符 的最大的k. 特别地，k不能取i+1（因为这个子串一共才 i+1 个字符，自己肯定与自己相等，就没有意义了）。

定义二：

   next 数组相当于 最大长度值 整体向右移动一位，然后初始值赋值为 -1。  (这个 最大长度值 指 定义一 中的 next数组)

 

因为next数组的定义不同，所以导致用next数组进行匹配 的代码 有点不同。很明显，该题要我们输出定义一中的next数组，不过，我更偏向用定义二中的next数组。

 

这张图很好的解释了用next数组快速匹配的核心：

图片来源：https://www.zhihu.com/question/21923021

 

#include<string>
#include<iostream>
using namespace std;
const int maxn = 1e6 + 5;

int* getNextArray(string P)
{
	static int next[maxn];
	next[0] = 0;  // 只有0号下标一个元素， 最长前缀 和 最长后缀 的长度 肯定为0

	int PLen = P.length();
	int i, j;

	for (i = 1; i < PLen; i++){ //枚举 模式串的下标
		for (j = i; j >= 0; j--){  //枚举 最长前缀 和 最长后缀 的长度

			if (P.substr(0, j).compare(P.substr(i - j + 1, j)) == 0){
				next[i] = j;
				break;
			}//if
		}//inner for

	}//extren for

	//整体向后移一位，同时第0为赋值为-1
	for (i = PLen; i >= 1; i--){
		next[i] = next[i - 1];
	}//for
	next[0] = -1;

	return next;
}


void KMP(string S, string P)
{
	int i1 = 0, i2 = 0;
	int SLen = S.length(), PLen = P.length();

	int *next = getNextArray(P);  //获得next数组

	while (i1 < SLen){
		if (S[i1] == P[i2]){
			i1++;
			i2++;
		}
		else if (i2){
			i2 = next[i2];
		}
		else if (next[i2] == -1){
			i1++;
		}

		if (i2 == PLen){
			cout << i1 - i2 + 1 << endl;
			i1--;
			i2 = next[i2 - 1];
		}//if
	}//while

	for (int i = 1; i <= PLen; i++){
		cout << next[i] << " ";
	}
	cout << endl;
}


int main()
{
	string S, P;
	cin >> S >> P;
	KMP(S, P);

	return 0;
}

 

这样求next数组，时间复杂度为O(m * m) 

我们来到了 KMP 的要点2 ：用最长前缀 等于 最长后缀的性质快速求出next[]数组。

毕竟是由数据状况 和 要解决的 问题出发 思考算法嘛。用两个例子来讲清楚怎么实现吧！

比如我想得到 下标为 10 的next数组，及在蓝色括号内的找最长前缀 和 最长后缀。可以知道下标为9 的next数组的值为4。如果在模式串中 i - 1 位置上的 字符 与 cn 位置上的字符相等。说明可以从i - 1 得到的最长前缀 和 最长后缀的基础上 追加一个 cn位置上的字符。在这个例子中，就是从原来abab 变为 ababc 。所以next 数组的值 为 next[i - 1] + 1。

 

如果在模式串中 i - 1 位置上的 字符 与 cn 位置上的字符不相等呢？？

我们就看next[cn] 上的值，并把这个值赋值给cn。   比如在下面这个例子中，在模式串中 i - 1 和 cn 的值不相等。而 next[cn] == 2 。 

看模式串下标为 cn 时 ，与下标为 i -1 上的值是否相等。如果相等，那么next[i] = cn。

如果还不相等，那么继续进行这个操作。如果cn 为 0 了，模式串下标为 cn  ，与下标为 i -1 上的值还不相等，那么next[i] = 0.

总结：

cn 为 i - 1位置上next数组最长前缀的后一个位置，i-1 为 i - 1 位置上next数组最长后缀的后一个位置。

 

#include<string>
#include<iostream>
using namespace std;
const int maxn = 1e6 + 5;

int* getNextArray(string P)
{
	static int next[maxn];
	next[0] = -1;
	int PLen = P.length();

	if (PLen == 1){
		return next;
	}//if

	next[1] = 0;
	int i = 2;
	int cn = 0;
	while (i <= PLen){  //得到的next[PLen] 是没有用到KMP中的
		if (P[i - 1] == P[cn]){
			next[i++] = ++cn;
		}//if
		else if (cn > 0){
			cn = next[cn];
		}//else if
		else{
			next[i++] = 0;
		}//else
	}//while

	return next;
}


void KMP(string S, string P)
{
	int i1 = 0, i2 = 0;
	int SLen = S.length(), PLen = P.length();

	int *next = getNextArray(P);  //获得next数组

	while (i1 < SLen){
		if (S[i1] == P[i2]){
			i1++;
			i2++;
		}
		else if (i2){
			i2 = next[i2];
		}
		else if (next[i2] == -1){
			i1++;
		}

		if (i2 == PLen){
			cout << i1 - i2 + 1 << endl;
			i1--;
			i2 = next[i2 - 1];
		}//if
	}//while

	for (int i = 1; i <= PLen; i++){
		cout << next[i] << " ";
	}
	cout << endl;
}


int main()
{
	//freopen("D:\\in.txt","r",stdin);
	string S, P;
	cin >> S >> P;
	KMP(S, P);

	return 0;
}

 

再来几道题加深印象：

剪花布条

 HDU - 2087 

题意：在一个主串中 ，找可以分为多少个模式串。模板题。

#pragma GCC optimize(3,"Ofast","inline")
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#include<queue>
#include<stack>
#include<list>
#include<map>
#include<set>
#include<cmath>
#include<sstream>
#include<cstdlib>
#include<bitset>
#include<climits>
#define F(i,s,t) for(int i=(s);i<=(t);i++)
#define D(i,s,t) for(int i=(s);i>=(t);i--)
#define dBug(i) printf("Value=%d\n",i)
#define ddBug(i,j) printf("Value=%d %d\n",i,j)
#define ed putchar('\n')
#define FO freopen("D:\\in.txt","r",stdin)
#define IOS cin.tie(0) ,cout.tie(0), cout.sync_with_stdio(0)
typedef long long ll;
//const int INF = 1 << 30;
//const double EPS = 1e-6;
//#define MX 102
//#define Mod 10000
using namespace std;
const int maxn = 1e6 + 5;

int* getNextArray(string P)
{
	static int next[maxn];
	next[0] = -1;
	int PLen = P.length();

	if (PLen == 1){
		return next;
	}

	next[1] = 0;
	int i = 2;
	int cn = 0;
	while (i <= PLen){
		if (P[i - 1] == P[cn]){
			next[i++] = ++cn;
		}//if
		else if (cn > 0){
			cn = next[cn];
		}//else if
		else{
			next[i++] = 0;
		}

	}//while

	return next;
}

int KMP(string S, string P)
{
	int i1 = 0, i2 = 0;
	int SLen = S.length(), PLen = P.length();
	int result = 0;

	int *next = getNextArray(P);

	while (i1 < SLen){
		if (S[i1] == P[i2]){
			i1++;
			i2++;
		}//if
		else if (i2){
			i2 = next[i2];
		}//else if
		else if (next[i2] == -1){
			i1++;
		}

		if (i2 == PLen){
			result++;
			i2 = 0;
		}//if
	}//while

	return result;
}


int main()
{
	string S, P;
	while (cin >> S){
		if (S.compare("#") == 0){
			break;
		}//if
		cin >> P;

		cout << KMP(S, P) << endl;
		S.clear();
		P.clear();
	}//while
	return 0;
}

 

Simpsons’ Hidden Talents

 HDU - 2594 

题意：

两个串，找到第一个串前缀 等于 第二个串的后缀 的子序列。

题解:

把这两个串连接起来，再求它的next数组，可是next数组的值不可以大于两个串的长度。

#pragma GCC optimize(3,"Ofast","inline")
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#include<queue>
#include<stack>
#include<list>
#include<map>
#include<set>
#include<cmath>
#include<sstream>
#include<cstdlib>
#include<bitset>
#include<climits>
#define F(i,s,t) for(int i=(s);i<=(t);i++)
#define D(i,s,t) for(int i=(s);i>=(t);i--)
#define dBug(i) printf("Value=%d\n",i)
#define ddBug(i,j) printf("Value=%d %d\n",i,j)
#define ed putchar('\n')
#define FO freopen("D:\\in.txt","r",stdin)
#define IOS cin.tie(0) ,cout.tie(0), cout.sync_with_stdio(0)
typedef long long ll;
//const int INF = 1 << 30;
//const double EPS = 1e-6;
//#define MX 102
//#define Mod 10000
using namespace std;
const int maxn = 5e4 + 5;

int* getNextArray(char* P)
{
	static int next[maxn * 2];
	next[0] = -1;
	int PLen = strlen(P);

	if (PLen == 1){
		return next;
	}

	next[1] = 0;
	int i = 2;
	int cn = 0;
	while (i <= PLen){
		if (P[i - 1] == P[cn]){
			next[i++] = ++cn;
		}//if
		else if (cn > 0){
			cn = next[cn];
		}//else if
		else{
			next[i++] = 0;
		}

	}//while

	return next;
}

int main()
{
	char s1[maxn], s2[maxn];
	int *Next;
	int s1L, s2L, ans;

	while (scanf("%s %s", s1, s2) != EOF){
		s1L = strlen(s1);
		s2L = strlen(s2);
		strcat(s1, s2);
		Next = getNextArray(s1);

		ans = Next[s1L + s2L];
		while (ans > s1L || ans > s2L){
			ans = Next[ans];
		}//inner while

		if (ans){
			for (int i = 0; i < ans; i++){
				printf("%c",s1[i]);
			}//for
			printf(" %d\n",ans);
		}//if
		else{
			printf("0\n");
		}
	}//while


	return 0;
}

 

 

Bil!
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
KMP 算法学习笔记.pdf
05-16
kmp算法
kmp-学习理解记录
01-20
本片是在回忆kmp的时候写下的，最近老师再讲这个，又发现之前写的都忘的差不多了在这里从新写一下笔记，纯小白文字： 基础暴力： 再看kmp的时候我们先来讲讲常规的模式串的暴力匹配，对于常规的一个文本串S以及一个...
2022.2.9
qq_63765733的博客
02-09 947
上午 学kmp，做题没做出来 下午 学长讲课 做kmp题 晚上 哈希 题目描述 给出两个字符串s_1s1​和s_2s2​，若s_1s1​的区间[l, r][l,r]子串与s_2s2​完全相同，则称s_2s2​在s_1s1​中出现了，其出现位置为ll。 现在请你求出s_2s2​在s_1s1​中所有出现的位置。 定义一个字符串ss的 border 为ss的一个非ss本身的子串tt，满足tt既是ss的前缀，又是ss的后缀。 对于s_2s2...
KMP算法详解
weixin_30273931的博客
10-08 98
背景： KMP算法(Knuth-Morris-Pratt Algorithm)是一种改进的字符串匹配算法，由D.E.Knuth，J.H.Morris和V.R.Pratt同时发现，因此人们称它为克努特——莫里斯——普拉特操作（简称KMP算法）。KMP算法的关键是利用匹配失败后的信息，尽量减少模式串与主串的匹配次数以达到快速匹配的目的。时间复杂度O(m+n)。 思路： ...
KMP——杨子曰算法
杨子曰：代码是神奇的
03-11 483
KMP——杨子曰算法 半年过后我又回来了 先给一道题： 如题，给出两个字符串s1和s2，其中s2为s1的子串，求出s2在s1中所有出现的位置。 （字符串长度&amp;amp;amp;amp;amp;amp;lt;=1000000） 大佬说：“哟，这道题的暴力好优啊！！” 就像下面这样指针i，j从头开始 哦，一样，i++，j++ 嗯，又一样，i++，j++ 继续 什么，不一样了，没事j从头来，i从第二个出发 一开始就不一样了，那再向前一格 …...
【C语言】编一个程序,将两个字符串s1和s2比较，若s1＞s2，输出一个正数；若s1=s2.输出0；若s1＜s2,输出一个负数。
魚肉的博客
05-11 8706
题目：编一个程序,将两个字符串s1和s2比较，若s1>s2，输出一个正数;若s1=s2，输出0;若s1<s2，输出一个负数。不要用strcpy函数。两个字符串用gets函数读入。输出的正数或负数的绝对值应是相比较的两个字符串相应字符的ASCII 码的差值。例如"A"与"C"相比，由于"A"<"C",应输出负数，同时由于'A'与'C'的ASCII码差值为2.因此应输出“-2”。同理: "And 和"Aid"比较，根据第2个字符比较结果,"n"比"i"的ASCII码大5.因此应输出“5” 解
编写一个程序，将两个字符串s1和s2比较，如果s1 > s2,输出一个正数；s1 = s2，输出0,；s1 < s2输出一个负数。不要使用strcmp函数。
热门推荐
桔子
10-17 2万+
两个字符串用gets函数读入。输出的正数或者负数的绝对值应是相比较的两个字符串相应字符的ASCII码的差值。例如“A”与“C”相比，由于“A” < “C”,应该输出负数，由于“A”与”C”的ASCII嘛差值为2，因此应该输出“-2”。同理：“And”和“Aid”比较，根据第二个字符比较的结果，”n”比”i”大5，因此输出5#include <stdio.h> #include <windows.h>
note:learning note 学习笔记
05-07
个人学习笔记 计算机学科基础 - 数据结构和算法 数据结构 哈希表 优先队列 后缀数组 红黑树 线段树 链表 栈 二叉树 图 算法 DFS、BFS 回溯法 排序 二分搜索 DP 贪心 KMP 目录 建设中... 掌握科班必备知识 c 语言 ...
algorithm：算法和数据结构学习笔记
02-12
主要处理边缘情况时间/空间复杂度分析专题排序链表题LRU 二分位运算栈和含量堆二叉树贪心并查集图动态规划单调栈滑动窗口和双指针完美洗牌问题斐波那契问题和扩展蓄水池算法计算表达式首要树Manacher算法KMP算法...
高级java笔试题-Lookoop:学习笔记
06-03
c++笔记归纳 Data Structures and Algorithm Analysis 数据结构与一些算法，来自算法导论，数据结构与算法分析-C语言描述，C Primer Plus, 数据结构-python描述，博客 ADT链表(c)-抽象链表实现 geometry（c++）-...
KMP算法学习笔记
jerryzhou的博客
09-13 423
参考1简书：https://www.jianshu.com/p/18410598c061（DFA_KMP) 参考2：https://segmentfault.com/a/1190000008575379（BruteForce、KMP） 参考2原文：https://subetter.com/algorithms-and-mathematics/kmp-algorithm.html ...
[算法]字符串的boader
codeingdog的专栏
07-24 390
定义: 我们定义, 一个字符串s, 它的border是既是它的前缀又是它的后缀, 且不等于s本身的串. 然后我们定义真border是所有border中的最长者. 这样我们发现每个border 其实是KMP算法中, 字符串S的next数组. ...
kmp学习笔记
蒟蒻QMQ
11-13 169
没错，博主又开始学算法了。。。话说我以前不是叫自己猴子吗？ 还是那句话，因为是个人的感受与体验，所以写上原创。 模版题 kmp简介 平常我们字符串匹配，暴力O(nm)O(nm)O(nm)是不是感觉自己已经可以称霸全世界了？ 不是！还有STL STL的大佬就不要看我这种蒟蒻的博客了。 当然，时间慢自然会有大佬去开发一种算法来加速。 就是他们，研究出了kmp看毛片算法。 sto,sto,sto km...
[KMP]【学习笔记】
weixin_33810302的博客
02-04 61
POJ3461 Oulipo   KMP裸题 出现几次   关于KMP 字符串从0开始，所以p[i]就是地i+1个字符 f[i]是失配函数，表示已经匹配了i个字符，i+1(就是p[i])失配转移到哪里 令j=f[i]，就是说以位置i-1结尾的后缀包括了0...j-1这个前缀，再检查p[j]＝＝s[i]（即j+1  i+1） // // main.cpp // p...
（填坑计划）KMP算法——利用pattern字符串信息优化匹配过程
一只菜鸡的博客
11-26 446
其实在此之前我已经看过好几回KMP教程了，各种讲解都挺详细的。但是由于平时实在用不上，导致一段时间不看就忘了，大概的思想还是记得的。昨天下午手撸了一面草稿纸，从匹配过程到next数组的建立，然后动手写code，然后找了一道模板题提交上去发现自己写的KMP是没有问题的。遂想记录一篇博客，把KMP算法的思路给记录下来，以后再遇到KMP，应该也能在一段时间内推导出来。 问题描述 首先KMP算法解决的是字符串匹配问题，给定模式串pattern与字符串s，需要找到s中是否有子串与pattern完全相同。那么显然这
[学习笔记]kmp学习笔记
Edwin VanCleef的博客
09-27 188
kmp学习笔记 问题： 给定两个字符串，将两个字符串进行匹配 原理： n^2暴力？那还要这个算法有啥用…… 先看一下简便的匹配方法： 这时候我们发现在D位置是失配，如果按照暴力方式的话我们要从位置A开始继续匹配 但是由于前面一定是匹配上的，我们就可以从C开始。 对于这种失配后跳到的节点我们把它叫做当前节点的失配。（个人的理解是
怎样能快速理解学习KMP算法
最新发布
03-31
2. 学习KMP算法的实现步骤：KMP算法的实现可以分为两步：首先，预处理模式串，生成next数组；然后，在文本串中匹配模式串，利用next数组来进行匹配。 3. 掌握next数组的原理和计算方法：next数组是KMP算法的核心之...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值