最近学习了马拉车求回文串长度的算法,感觉这个比较难理解,放的时间稍微一长,今天一看,发现就快忘记了,所以写一篇博客总结一下。
算法复杂度为o(n),实际上只是初始化串一遍和走过串一遍,效率比字符串哈希和两层for循环好的多。
具体详情请参看注释,适合之前已经了解马拉车算法,现在复习时看看。
洛谷p3805
题目描述
给出一个只由小写英文字符a,b,c…y,z组成的字符串S,求S中最长回文串的长度.
字符串长度为n
输入格式
一行小写英文字符a,b,c…y,z组成的字符串S
输出格式
一个整数表示答案
#include <bits/stdc++.h>
using namespace std;
char a[11000000+100],b[11000000*2+200];
int mp[11000000*2+200];//里面最终存的是在这个串的这个位置的回文串长度+1
int main(){
scanf("%s",a);
int len = strlen(a),cnt = 1;
b[0] = '$';
b[cnt++]='#';
for(int i = 0;i < len;i++){
b[cnt++] = a[i];
b[cnt++] = '#';
}
b[cnt] = '0';
// 初始化加#为了统一奇数串和偶数串的做法,加$和0为了不用判断边界条件,
// 切结,这三个字符串不能在原字符串中出现
int mx = 0,id = 0;
//mx的意义在于该串现阶段最远到达的位置,id是用于计算现在走到的位置的那个点通过现阶段的中点位置的对称位置,
//可以找到下标在哪里,通过那一个位置串的回文串大小既可以得到现在位置可以到的最远位置。
for(int i = 0;i < cnt;i++){
if(mx > i){//判断现在遍历到的位置i 是否已经超过这个串判断过的最远距离,如果没有
mp[i] = min(mp[2*id-i],mx-i);//通过id计算出对称点回文串长度,和现在已经走到最远位置做比较,
//取小的,如果对称点的回文串长度大于现在已经遍历到的最远距离,则不能确定在没遍历到那个位置是否匹配,
//需要在下面循环中判断,如果是mx-i大的话,说明这个正好可以直接用对称点那里匹配的长度
}
else{//当mx<=i时,就不在范围之内了,因为mx当时设定的就是未匹配上的第一个字符
mp[i] = 1; //如果超过了,只能在下面老老实实匹配
}
while(b[i+mp[i]] == b[i-mp[i]])//老老实实匹配的过程
mp[i]++;
if(i+mp[i]>mx){//判断在这次更新之后,这个串到达的最远位置是否发生改变,发生了就更改mx和id的值
mx = i+ mp[i];
id = i;
}
}
printf("%d\n",*max_element(mp+1,mp+cnt)-1);
return 0;
}