今天上午在实验室窝了一上午,水了一上午的题,
唯一有收获的就是新学习了一个最长回文子串的算法
Manacher's ALGORITHM: O(n)时间求字符串的最长回文子串
https://www.felix021.com/blog/read.php?2040
算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符,当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了(见下面的一个例子,回文串长度全为奇数了),然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P[id]记录的是以字符str[id]为中心的最长回文串,当以str[id]为第一个字符,这个最长回文串向右延伸了P[id]个字符。
原串: w a a b w s w f d
新串: # w # a # a # b # w # s # w # f # d #
辅助数组P:1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
这里有一个很好的性质,P[id]-1就是该回文子串在原串中的长度(包括‘#’)。如果这里不是特别清楚,可以自己拿出纸来画一画,自己体会体会。当然这里可能每个人写法不尽相同,不过我想大致思路应该是一样的吧。
好,我们继续。现在的关键问题就在于怎么在O(n)时间复杂度内求出P数组了。只要把这个P数组求出来,最长回文子串就可以直接扫一遍得出来了。
由于这个算法是线性从前往后扫的。那么当我们准备求P[i]的时候,i以前的P[j]我们是已经得到了的。我们用mx记在i之前的回文串中,延伸至最右端的位置+1处。同时用id这个变量记下取得这个最优mx时的id值。(注:为了防止字符比较的时候越界,我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’,故我的新串下标是从1开始的)
好,到这里,我们可以先贴一份代码了。
const int MAXN=100000+1000;
char S[MAXN*2];//改造后的串
char T[MAXN];//原始串
int p[MAXN*2];
void init()//改造原始串
{
int i,j;
S[0]='@',S[1]='#';
for(int i=0,j=2;T[i];i++,j+=2)
{
S[j]=T[i];
S[j+1]='#';
}
S[j]=0;
}
void manacher()
{
int id,i,mx=0;//mx指的是那个id能触及的最大位置 +1 的地方
for(int i=1;S[i];i++)
{
if(mx>i) p[i]=min(p[2*id-i] , mx-i);
else p[i]=1;
while(S[i-p[i]] == S[i+p[i]])
p[i]++;
if(i+p[i]>mx)
{
mx=i+p[i];
id=i;
}
}
}