题目参考:http://hihocoder.com/problemset/problem/1032
参考链接:http://blog.sina.com.cn/s/blog_70811e1a01014esn.html
大体思路其实就是找出一个中心点,判断中心点两端的回文串半径是多少;
但由于找中心点的方法只适用于奇数长的回文串,故可以在每两个字符间插入一个间隔符来帮助结算;
用rd[i]表示以经过填充后的字符串里的第i个字符为中心,它的回文串长度;
可以得知,在【i-rd[i],i-1】这一段和【i+1,i+rd[i]】的字符串呈回文状;
然后对于k在【1,rd[i]】的范围内,找出rd[i+k]的最小值(或是一个准确的值);
根据回文串的对称关系,可以将rd[i+k]分成两种情况
1、rd[i]-k!=rd[i-k]时,rd[i+k]=min(rd[i]-k,rd[i-k]),该式子证明可以在上文链接中看到,大致是根据第i个字符的左边k个字符和右边k个字符的对称性来得出答案的;
2、rd[i]-k==rd[i-k]时,只能判断出rd[i+k]的最小值为rd[i]-k,但并不知道rd[i+k]的准确值,因为此时应打断k的循环,将i移至i+k再作一次回文串长度的判断;
算法的时间复杂度分析:
虽然实现的代码里看上去是两层循环,但是实际上,当遍历到 i 的时候,能保证 i 后面的字符被遍历不超过3次,因此时间复杂度应该是O(n)。
AC代码:
#include <bits/stdc++.h>
using namespace std;
char s[1000006];
char road[2000006];
int rd[2000006];
int main()
{
int T;
cin>>T;
while(T--)
{
cin>>s;
int len=strlen(s);
memset(road,0,sizeof(road));
for(int i=0 ; i<len ; i++)
{
road[i*2]='#';
road[i*2+1]=s[i];
}
road[2*len]='#';
road[2*len+1]='\0';
len= 2*len+1;
int j=0;
int k;
int res=0;
for(int i=1 ; i<len ; )
{
while(i+j+1<len&&road[i-j-1]==road[i+j+1])
j++;
rd[i]=j;
res=max(rd[i],res);
for(k=1; k<=j&&rd[i-k]!=rd[i]-k ; k++)
{
rd[i+k]=min(rd[i-k],rd[i]-k);
}
i+=k;
j=max(j-k,0);
}
cout<<res<<endl;
}
return 0;
}