1、为了规避奇偶串的不同情况,对原串进行一下处理:
原串 abcd —> #a#b#c#d#
原串 1234 —> #1#2#3#4#
即间隔插入分隔符 ‘#’ ;
这样以来,所有的回文子串都变成的奇数长度。
2、引入变量:p数组,mx, mid;
设str是对原串处理后的新串。
p[i]表示以str[i]为中心的最长回文子串的半径,不难推出原串的长度是p[i] - 1。
mx是已经计算过的所有回文串能到达的最右端的边界,mid是该回文串的中心。
接下来是求解p[i]的过程,(求p[i]时,p[1~i-1]已经计算完毕):
显然有mid < i,(因为p[mid]已经计算过)。
情况一:i < mx
此时我们可以利用回文串的对称性,找到i关于mid的对称点j(j = 2 * mid - i)。
p[j]有两种情况:
一是以j为中心的最大回文串被包含在以mid为中心的最大回文串中,由对称性得到p[i]至少也是p[j],也就是p[i] >= p[j],接下来就暴力求解不确定的部分。
二是以j为中心的最大回文串超出了以mid为中心的最大回文串的范围,由于超出的部分不在以mid为中心的最大回文串内,不具有对称性,但是在mid回文串以内的部分一定与i对应的位置相同,p[i]至少是mx - i,即p[i] >= mx - i,接下来不确定的外围部分暴力求解。
情况二:i>= mx
此时,没有任何其他信息可以辅助计算p[i],令p[i] = 1,然后暴力求解即可。
例题:
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<map>
#include<stack>
#include<queue>
#include<string>
#include<vector>
#include<algorithm>
#define ll long long
#define ull unsigned long long
#define mem(a) memset(a,0,sizeof(a))
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
#define inf 0x3f3f3f3f
#define sf(n) scanf("%d", &n)
#define endl '\n'
#define cyes cout<<"yes"<<endl
#define debug(x) cout<<"-----"<<x<<endl
using namespace std;
const int N = 1e6 + 5;
int n, m;
char s[N], str[N<<1];
int p[N<<1];
int manacher()
{
//由原串s得到新串str
str[0] = '$'; str[1] = '#';
for(int i = 1; i <= n; i++)
{
str[2 * i] = s[i];
str[2 * i + 1] = '#';
}
m = 2 * n + 1;
int mx = 0, mid, res = 0;
for(int i = 1; i <= m; i++)
{
p[i] = i < mx ? min(p[2 * mid - i], mx - i) : 1;
while(str[i + p[i]] == str[i - p[i]]) p[i]++; //向外扩展
if(i + p[i] > mx) //更新mx和mid
{
mx = i + p[i];
mid = i;
}
res = max(res, p[i] - 1); //对每一个回文子串取最大
}
return res;
}
int main()
{
int kase = 0;
while(scanf("%s", s + 1), s[1] != 'E')
{
n = strlen(s + 1);
printf("Case %d: %d\n", ++kase, manacher());
}
return 0;
}