中心拓展法o(n2)
遍历每一个可能的中心点,以该中心点为基础,向两边延伸,生成回文串,直到两末端字符不相等,此时得到以该点为中心的最长回文串。
Manacher算法
相关知识
时间复杂度为O(N),空间复杂度为O(N)
Manacher用于解决最长回文子串的线性方法,同时也可以顺便解决回文子串的计数问题,巧妙运用了字符串的对称性。
因为字符串可能是偶数长度,例如字符串baab,此时中心是在两个aa中间,此时枚举较为麻烦,所以我们可以加上一些#,此时字符串变为#b#a#a#b#,为奇数长度,同时,为了在后续的代码中避免处理数组越界的情况,我们可以把字符串变为^#b#a#a#b#$,即在两端位置添加与#不同的符号,此时循环体在遇到这两个字符时会自动退出
视频讲解
如有需要,可以观看此视频,加深理解http://【马拉车算法 | Coding Club】https://www.bilibili.com/video/BV1Sx4y1k7jG?vd_source=4b861fc9f463823ccf14e08ce49d7bf2
代码细节
用 A的下标 i 减去 A [ i ],再除以 2 ,就是最原始原字符串的开头下标了。
int i_index=2*c-i;
求解A[i]
我们用 C 表示回文串的中心,用 R 表示回文串的右边半径。所以
R = C + A[ i ]。
C 和 R 所对应的回文串是当前循环中 R 最靠右的回文串。
完整代码
#include <bits/stdc++.h>
using namespace std;
const int N=1e5+5;
int a[N];
int main(){
string s,t;
while(cin>>s){
int len=s.size();
int n=0;
//改变字符串,在字母空隙添加字符
t="^#";
for(int i=0;i<len;i++){
t+=s[i];
t+='#';
}
t+='$';
memset(a,0,sizeof(a));
int c=0,r=0; //c表示回文串的中心,r表示回文串的右边半径
for(int i=1;i<t.size();i++){
int i_index=2*c-i;//表示当前所求的第i个字符关于中心c对应的下标
if(i<r){
a[i]=min(r-i,a[i_index]);//防止超出回文串的右边半径即防止超出r
}
else{
a[i]=0;//i等于r时,a[i]赋值为0,通过中心拓展一步一步拓展
}
//中心拓展法
while(t[i-1-a[i]]==t[i+1+a[i]]){
a[i]++;
}
if(i+a[i]>r){
//当我们求出的a[i]的右边界大于r时,此时需要更新c和r为当前的回文串,因为必须保证i在r里面,所以要不断更新r
c=i;
r=i+a[i];
}
}
int maxx=0;
for(int i=1;i<t.size();i++){
maxx=max(maxx,a[i]);//求最长回文子串的长度
}
cout<<maxx<<endl;
}
return 0;
}
例题https://ac.nowcoder.com/acm/problem/14517