51nod1089(最长回文子串之manacher算法)

提示：扫一扫查出行
【扫一扫了解最新限行尾号】
复制提示

题目链接: https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1089

题意:中文题诶~

思路: 我前面做的那道回文子串的题目是枚举中间字符O(n^2)时间过的,不过这题字符串长度限制为1e5,O(n^2)肯定会超时啦;

有个叫 manacher 的算法是时间复杂度为 O(n), 本题就是 manacher 模板题啦;

我们先看一下 manacher 算法这个东东;

首先回文串匹配奇数长度和偶数长度操作是不一样的, 我们可以在每个字符的两边都加上一个特殊字符 '#', 那么无论原串是长度是奇数还是偶数都会编程奇数长度,并且它的所有回文子串也都变成了奇数长度. 因为对于长度为 len 的字符串 str 我们要使每个字符两边都有一个 '#' 字符,需加 len+1 个 '#' 字符, 那么改变后的字符串 s 长度为 2*lne+1, 其必为奇数;那么现在我们只需要处理奇数的情况了;

我们在用 vis 数组存储以i中心的最大回文串的回文半径(包括 s[i] 字符),那么i对应的原串 str 的回文串长度为 vis[i]-1. 因为str中该回文串长度为 2*vis[i]-1, 其中有 vis[i]个'#'字符嘛, 那么我们只要求出所有 vis[i], 那么答案也就知道了啦;

我们再维护两个变量 id 为当前最长回文子串的中心位置, mx 为其右边界位置;

那么对于当前 i (我们是从前往后求 vis 的, vis[id] 已知, vis[i]未知, 所以 i 一定是在 id 后面的) 我们可以分情况讨论:

对于 mx>i 的情况, 即 i 字符在当前最大回文串里面. 我们可以找到 i 关于 id 的对称位置 j=2*id-i , 若 vis[j] <= mx-i 即以 j 位置为中心的回文串最左端没有超出最大回文串的范围, 即以 j 为中心和以 i 为中心的最大回文子串都在以 id 为中心的那个回文串里面, 那么由回文串的对称性我们可以知道这两个字符串是一样的, 所以有 vis[i]=vis[j];

如果vis[j] > mx-i, 即以 j 为中心的最大回文子串超出了最大回文子串的范围, 那么 vis[i]>=mx-i, 对于更长的范围就需要我们一个个去匹配了啦;

对于 mx<i 的情况我们并没有什么信息可以利用, 所以需要重新匹配~

代码:

#include <bits/stdc++.h>
#define MAXN 200010
using namespace std; int main(void){
char s[MAXN], str[MAXN];
int vis[MAXN], j=;
memset(vis, , sizeof(vis));
scanf("%s", str);
for(int i=; str[i]!='\0'; i++){ //在每个字符的两边添加一个特殊字符'#'
s[j++]='#';
s[j++]=str[i];
}
s[j++]='#', vis[]=;
s[j++]='\0';
int len=j-, mx=, id=, ans=; //**id表示当前最大回文子串的中心位置, mx表示当前最大回文子串的最右端位置, ans表示当前已知的最大回文串长度
for(int i=; i<len; i++){
if(mx>i){ //i在当前最大回文串内的情况
vis[i]=min(mx-i, vis[*id-i]);
}else{
vis[i]=;
}
while(s[i-vis[i]]==s[i+vis[i]]&&i>=vis[i]){ //注意边界
vis[i]++; //**匹配没有可利用信息的字符
}
if(vis[i]+i>mx){ //**更新id和mx的值
id=i;
mx=i+vis[i];
}
ans=max(ans, vis[i]); //**更新ans
}
printf("%d\n", ans-);
return ;
}