【算法学习】KMP 算法

合集 - 【算法学习】(35)

1.【算法学习】排序2024-08-272.【算法学习】Manacher 马拉车2024-08-29

3.【算法学习】KMP 算法2024-09-01

4.LCA 最近公共祖先（树链和倍增）这次真有树链了！！！2024-08-305.线段覆盖问题2024-08-296.【算法学习】学换根dp有感2024-09-207.二分图最大匹配2024-09-178.【算法学习】01BFS2024-09-179.洛谷 P1892 [BOI2003] 团伙 种类并查集 扩展域并查集2024-09-1110.【算法学习】高斯消元法2024-09-1011.贝叶斯公式2024-09-0712.背包2024-09-0513.【算法学习】模拟退火2024-09-0514.【算法学习】基环树2024-10-2915.【算法学习】树链部分2024-10-2916.【算法学习】莫队2024-10-2917.【算法学习】分块九讲2024-10-2818.平衡树2024-10-2319.圆方树2024-10-2220.【算法学习】点分治2024-10-1121.公式2024-10-0922.【算法学习】笛卡尔树2024-10-0423.【算法学习】悬线法2024-10-0424.欧几里得算法与 EX2024-09-2325.【算法学习】逆元与求解2024-09-2326.【算法学习】费马定理2024-09-2327.三分2024-09-2328.裴蜀定理2024-09-2329.【算法学习】欧拉函数φ2024-09-2330.【算法学习】二维转一维问题2024-11-0331.【算法学习】扫描线2024-11-0132.【算法学习】矩阵乘法2024-10-3133.【算法学习】同余最短路2024-11-0434.【算法学习】组合数学2024-11-0535.【算法学习】反悔贪心2024-11-06

收起

学习笔记

我认为我这个算法可能无法讲明白，而且工作量巨大，所以为了让你快速学会我推荐学习下列笔记。

学习笔记1

学习笔记2

学习笔记3

放个模板

#include <bits/stdc++.h>
using namespace std;
const int N=1e6+10;
char a[N],b[N];
int lena,lenb;
int nxt[N];
 
int main(){
    ios::sync_with_stdio(false);
    cin>>a+1;
    cin>>b+1;
    lena=strlen(a+1);
    lenb=strlen(b+1);
    
    int j=0;
    
    for(int i=2;i<=lenb;i++){
    	while(j&&b[i]!=b[j+1]){
    		j=nxt[j];
		}
		if(b[j+1]==b[i]){
			j++;
		}
		nxt[i]=j;
	}
    j=0;
    for(int i=1;i<=lena;i++){
    	while(j&&a[i]!=b[j+1]){
    		j=nxt[j];
		}
		if(b[j+1]==a[i]){
			j++;
		}
		if(j==lenb){
			cout<<i-lenb+1<<"\n";
		}
	}
    
   	for(int i=1;i<=lenb;i++){
   		cout<<nxt[i]<<" ";
	} 
    
    return 0;
}

例题

感觉经过上述的学习，你一定有所收获吧（如果没有的话还是菜就多练吧），所以接下来我会举出一些题目，应该会对你的学习有些帮助。

1.【模板】KMP（洛谷P3375）

如果你学 KMP 请先会做对这道题，这将会检验你的代码能力，请务必做对。

2.字符串中的最大值（51nod P1277）

这道题是对 KMP nxt数组 很好的应用，题意是求其所有前缀中，字符长度与出现次数的乘积的最大值。

我们可知 nxt 数组指示了当前模式串在该位置失配时，应该将模式串的哪一位与此位对齐，所以我们可以从后向前统计每个前缀的出现次数 \(res[]\)，并将 \(res[nxt[i]]+=res[i]\) ，最后再从前向后计算最大值。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N=1e6+10;
ll res[N];
char b[N];
int nxt[N];
int lenb;
 
int main(){
    ios::sync_with_stdio(false);
	
	cin>>b+1;
	
	lenb=strlen(b+1);
	
	int j=0;
	
	for(int i=2;i<=lenb;i++){
		while(j&&b[i]!=b[j+1]){
			j=nxt[j];
		}
		if(b[i]==b[j+1]){
			j++;
		}	
		nxt[i]=j;
	}
	ll ans=0;
	
	for(int i=lenb;i>=1;i--){
		res[i]++;
		if(nxt[i]!=0){
			res[nxt[i]]+=res[i];
		}
	}
	
	for(int i=1;i<=lenb;i++){
//		cout<<res[i]<<" ";
		ans=max(ans,res[i]*i);
	} 
	cout<<ans;
    return 0;
}

2.逆向工程（51nod P3118）

nxt 第 \(i\) 位不等于 0 时可以求出，若为 0 时，需不断向前找 nxt 将遇到的数字打标记，然后从小到大找第一个没打过标记的字母即可。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N=1e5+10;
int n;
int nxt[N];
char b[N];
int a[30];
 
int main(){
    ios::sync_with_stdio(false);
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>nxt[i];
	}
	b[1]='a';
	for(int i=2;i<=n;i++){
		if(nxt[i]!=0){
			b[i]=b[nxt[i]];
		}
		else{
			int j=nxt[i-1]+1;
			while(1){
				a[b[j]-'a']=1;
				if(j==1){
					break;
				}
				j=nxt[j-1]+1;
			}
			int flag=0;
			for(int k=0;k<26;k++){
				if(a[k]==1){
					a[k]=0;
					continue;
				}
				if(a[k]==0&&flag==0){
					b[i]=k+'a';
					flag=1;
				}
			}
		}
	}
	for(int i=1;i<=n;i++){
		cout<<b[i];
	}
    return 0;
}