题目链接:【hdu 3068】
输入有a~z组成的字符串,输出最长回文子串的长度
*Manacher算法模板 ==> 一个专用于求最长回文串的算法
1、转化字符串:abcba ==> $#a#b#c#b#a#(前面加$防止数组溢出)
2、s: $ # a # b # c # b # a #
p: 1 2 1 2 1 6 1 2 1 2 1
p[i]的含义:以s[i]为中心的回文串右扩张的最长长度,p[i] - 1就是转化前的字符串以i为中心的回文子串的长度
3、求p[],设id为回文子串的中心,mx为边界,则mx = id + p[id],i与j关于id对称
(1)id<i<mx
==>p[j] < mx-i
整张图是关于id对称的,所以p[i] = p[i]
==>p[j]>=mx-i
此时p[i] = mx-i
(2)i>=mx
这里的p[i]初始值只能是1,通过后续计算慢慢累加
void getp()
{
memset(p, 0, sizeof(p));
int id = 0, mx = 0;
for(int i=1; i<len; i++)
{
int j = 2*id-i;
p[i] = mx>i ? min(p[j], mx-i):1;
while(s[i+p[i]] == s[i-p[i]]) p[i]++;
if(i+p[i] > mx)
{
mx = i+p[i];
id = i;
}
}
}
所以本题的完整代码就是:
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;
#define inf 310010
char s[inf], str[inf];
int p[inf], len;
void gets()
{
int li = strlen(str);
len = 0;
s[len++] = '&';
s[len++] = '#';
for(int i=0; i<li; i++)
{
s[len++] = str[i];
s[len++] = '#';
}
}
void getp()
{
gets();
memset(p, 0, sizeof(p));
int id = 0, mx = 0;
for(int i=1; i<len; i++)
{
int j = 2*id-i;
p[i] = mx>i ? min(p[j], mx-i):1;
while(s[i+p[i]] == s[i-p[i]]) p[i]++;
if(i+p[i] > mx)
{
mx = i+p[i];
id = i;
}
}
}
int main()
{
while(~scanf("%s", str))
{
getp();
int ans=1;
for(int i=0; i<len; i++) ans = max(ans, p[i]-1);
printf("%d\n", ans);
}
return 0;
}