简介
在字符串的题目中,有时会遇上 回文串 这样一个名词
顾名思义,回文串 就是 正读和反读都一样的字符串
而 最长回文子串 ,就是在一个字符串的所有子串中,是回文串且长度最长的那个
求最长回文子串最普通的方法是 O(N2) ,即枚举一个点往两边扩展
但是在有些题目中,N 却十分的大
那么我们就要用到 时间空间复杂度都是 O(N) 的 Manacher 算法
用法
在处理回文串时,我们常常会被 中间字符是一个还是两个 这样的问题困扰
但是在机智的 Manacher 算法 中,这个问题得到了完美的解决
在每两个字符中间插入一个不会出现的分隔符(如:#)
之后在头尾插入一个还是没有出现的分隔符(如:*)来防止 While 出界
这样处理起来就方便很多了!
设读入的字符串为 s[i] ,
记录 p[i] 表示 以 s[i] 为中心往两边扩展的最大长度
观察可知,实际的回文串长度即为当前的 p[i]−1
再记录一个数 id , p[id]+id 表示在 i 位置前所有回文串中能延伸到的最右端的位置
如下图:
算法核心就是:
if(p[id]+id>i) p[i]=min(p[id∗2−i],p[id]+id−i); else p[i]=1;当之前所有回文串中能延伸到的最右端覆盖过 i 时,则取最小值,否则
p[i]=1 ,及自己本身这样不断维护 p[i] 和 id ,就能在 O(N) 内求出 最长回文子串 了!
至于为什么时间是线性的,由于最有端 p[id]+id 最多只能移动 N 次,
有效移动的操作就严格线性啦!!
下面附上模板:
void Manacher()
{
scanf("%s",s);
int len=strlen(s);
for(int i=len;i>=0;i--)
{
int k=i*2+1;
s[k+1]=s[i],s[k]='#';
}//插入分隔符
len*=2;
s[ans=id=0]='*';
for(int i=2;i<=len;i++)
{
if(p[id]+id>i) p[i]=min(p[id*2-i],p[id]+id-i); else p[i]=1;
while(s[i-p[i]]==s[i+p[i]]) p[i]++;
if(p[i]+i>p[id]+id) id=i;
if(p[i]>ans) ans=p[i];
}//处理、计算
}
- 注释: s[i],p[i],id 如题意义, ans 表示 最长回文子串 的长度,而 len 是原串长度
总结
这个 Manacher 算法效率极高,时间空间都是 O(N) 线性的
再者代码极短,所以使用起来十分方便,应多多使用!!!