KMP算法初步

KMP算法太难了！！！为了理解花了好长时间。

 

想找个通俗点的比喻一时半会也想不出来，先强记吧...

学习视频：1.帮你把KMP算法学个通透！（理论篇）_哔哩哔哩_bilibili

                  2.帮你把KMP算法学个通透！（求next数组代码篇）_哔哩哔哩_bilibili

模板题：P3375 【模板】KMP字符串匹配 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

模板代码：

#include<bits/stdc++.h>
using namespace std;
char s1[2000005], s2[2000005];//s1目标串,s2模式串
int Next[2000005], length_s1, length_s2;
void getNext()
{
	int len = 0, last;//len表示前缀串最后一位的序号同时也表示前(包括)last位的最长公共前后缀长度，last表示后缀串最后一位的序号
	for (last = 1; last < length_s2; last++)
	{
		while (len > 0 && s2[len] != s2[last])
		{
			len = Next[len - 1];
		}
		if (s2[len] == s2[last])
		{
			len++;
			Next[last] = len;
		}
	}
}
void KMP()
{
	int pos_s2 = 0;//指向模式串与目标串比较的点
	for (int i = 0; i < length_s1; i++)
	{
		if (s1[i] == s2[pos_s2] && pos_s2 == length_s2 - 1)
		{
			cout << i - length_s2 + 2 << endl;//输出模式串在目标串中第一次出现的位置
			pos_s2 = Next[pos_s2];
			continue;
		}
		else if (s1[i] == s2[pos_s2])
		{
			pos_s2++;
			continue;
		}
		while (s1[i] != s2[pos_s2] && pos_s2)
		{
			pos_s2 = Next[pos_s2 - 1];
			i--;//保持比较位置不变
		}
	}
}
int main(void)
{
	cin >> s1;
	cin >> s2;
	length_s1 = strlen(s1);
	length_s2 = strlen(s2);
	getNext();
	KMP();
	for (int i = 0; i < length_s2; i++)//输出Next数组
	{
		cout << Next[i] << ' ';
	}
	return 0;
}

上述代码是第一次学习(2021.8.15)KMP时写的，学得不是很透彻，只得了60分。

---------------------------------------------------------------------------------------------------------------------------------

再次学习：图解字符串匹配的KMP算法 (qq.com)

2021.10.25 再次学习KMP算法，按照自己的想法写的代码，经过调试AC

#include<bits/stdc++.h>
using namespace std;
char s1[2000005], s2[2000005];//s1目标串,s2主串
int Next[2000005], length_s1, length_s2;
void getNext()
{
	for (int i = 0; i < length_s1; i++)//求目标串第i位的部分匹配值
	{
		for (int j = i; j >= 1; j--)//假设目标串第i位的最大部分匹配值为j
		{
			int a = 0, flag = 0;
			for (int k = i - j + 1; k <= i ; k++)
			{
				if (s1[a] != s1[k])//s1[a]是前缀,s[k]是后缀
				{
					flag = 1;
					break;
				}
				a++;
			}
			if (flag == 0)
			{
				Next[i] = j;
				break;
			}
		}
	}
}
void KMP()
{
	int p_s1 = 0, p_s2 = 0;
	while (1)
	{
		if (s1[p_s1] == s2[p_s2])
		{
			p_s1++;
			p_s2++;
		}
		if (p_s1 == length_s1)
		{
			cout << p_s2 - p_s1 + 1 << endl;
			if (p_s1 >= 1)
				p_s1 = Next[p_s1 - 1];
		}
		if(p_s2 == length_s2)
			break;
		if(s1[p_s1] != s2[p_s2])
		{
			if (p_s1 == 0)//s1的第一个字符与主串不一样，主串指针+1
				p_s2++;
			if(p_s1 >= 1)
				p_s1 = Next[p_s1 - 1];
		}
	}
}
int main(void)
{
	cin >> s2;
	cin >> s1;
	length_s1 = strlen(s1);
	length_s2 = strlen(s2);
	getNext();
	KMP();
	for (int i = 0; i < length_s1; i++)//输出Next数组
		cout << Next[i] << ' ';
	return 0;
}

例题：HDU 2087 剪花布条

Problem Description

一块花布条，里面有些图案，另有一块直接可用的小饰条，里面也有一些图案。对于给定的花布条和小饰条，计算一下能从花布条中尽可能剪出几块小饰条来呢？

Input

输入中含有一些数据，分别是成对出现的花布条和小饰条，其布条都是用可见ASCII字符表示的，可见的ASCII字符有多少个，布条的花纹也有多少种花样。花纹条和小饰条不会超过1000个字符长。如果遇见#字符，则不再进行工作。

Output

输出能从花纹布中剪出的最多小饰条个数，如果一块都没有，那就老老实实输出0，每个结果之间应换行。

Sample Input

abcde a3

aaaaaa aa

#

Sample Outpu

0

3

#include<bits/stdc++.h>
using namespace std;
char s1[1005], s2[1005];//s1目标串,s2模式串
int Next[1005], length_s1, length_s2, ans;
void getNext()
{
	memset(Next, 0, sizeof(Next));
	int len = 0, last;//len表示前缀串最后一位的序号同时也表示前(包括)last位的最长公共前后缀长度，last表示后缀串最后一位的序号
	for (last = 1; last < length_s2; last++)
	{
		while (len > 0 && s2[len] != s2[last])
		{
			len = Next[len - 1];
		}
		if (s2[len] == s2[last])
		{
			len++;
			Next[last] = len;
		}
	}
}
void KMP()
{
	int pos_s2 = 0;
	for (int i = 0; i < length_s1; i++)
	{
		if (s1[i] == s2[pos_s2] && pos_s2 == length_s2 - 1)
		{
			ans++;
			pos_s2 = 0;//注意模式串不能重叠比较
			continue;
		}
		else if (s1[i] == s2[pos_s2])
		{
			pos_s2++;
			continue;
		}
		while (s1[i] != s2[pos_s2] && pos_s2)
		{
			pos_s2 = Next[pos_s2 - 1];
			i--;
		}
	}
}
int main(void)
{
	while (1)
	{
		cin >> s1;
		if (s1[0] == '#')
			break;
		cin >> s2;
		length_s1 = strlen(s1);
		length_s2 = strlen(s2);
		getNext();
		KMP();
		cout<<ans<<endl;
		ans = 0;
	}
	return 0;
}

一道深入理解next数组的题：Problem - 3336 (hdu.edu.cn)

我理解不了呀...