题目意思:
题目链接:http://hihocoder.com/contest/hiho1/problem/1
给出一个字符串(n<=10^6),求出这个字符串的最长回文字串。输入有多组数据,只需要输出对应的结果即可。
样例输入
3
abababa
aaaabaa
acacdas
样例输出
7
5
3
解析:当我们看到n<=10^6时,证明一般的模拟算法肯定是过不去的,对于时间复杂度O(n^3),O(n^2)的算法肯定是要TL的,因此需要提到Manacher算法,在O(n)时间求出最长回文字串。看一下小hi的提示:
提示三:
小Ho答道:“我想想,如果以第5个字符为中心的最长回文子串的长度是5的话,这就告诉了我[3, 7]这一段是一个回文子串,所以呢?”
小Hi继续提示道:“假设这时候你想要计算以第6个字符为中心的最长回文子串的长度,你有没有什么已知的信息了?”
小Ho边想边说道:“唔,首先第6个字符和第4个字符是一样的,第7个字符和第3个字符是一样的,而第5个字符本身就肯定和第5个字符一样,那么如果[3, 5]这一段是回文子串的话,那么[5, 7]这一段肯定也是回文子串。
Manacher算法:引用文章
AC代码:
/**
*Author: xiaoran
*Solution:
*
*/
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#define LL long long
using namespace std;
const int MAXN=1000005;
char s[MAXN];
bool huiWen1(char *s,int low,int high){//枚举起始位置
while(low<high){
if(s[low]!=s[high]) return false;//不是回文
else low++,high--;
}
return true;//是回文
}
int huiWen2(char *s,int mid){//枚举中心位置mid
int len=strlen(s);
//奇数
int max1=1;
int k=max(mid,len-mid-1);
for(int i=1;i<=k;i++){
if(s[mid-i]==s[mid+i]) max1+=2;
else break;
}
//偶数
int max2=0;
k=max(mid,len-mid-1);
for(int i=0;i<k;i++){
if(s[mid-i]==s[mid+i+1]) max2+=2;
else break;
}
return max(max1,max2);
}
char Ma[MAXN*2];
int Mp[MAXN*2];
//Mp[i]表示第i个字符的回文半径,即rad[i]尽可能大
//且满足Ma[i-Mp[i],i-1]=Ma[i+1,i+Mp[i]]
//Manacher算法。时间O(n),空间O(n)最长回文串
int Manacher(char *s,int len){
int k=0;
//初始化字符串,将所有字符串转化为奇数长度
Ma[k++]='$';
Ma[k++]='#';
for(int i=0;i<len;i++){
Ma[k++]=s[i];
Ma[k++]='#';
}
Ma[k]='\0';//字符串结束
//其中id表示最大回文子串中心的位置,mx则为id+P[id],也就是最大回文子串的边界
int mx=0,id=0;
for(int i=0;i<k;i++){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]]) Mp[i]++;
if(i+Mp[i]>mx){
mx=i+Mp[i];
id=i;
}
}
}
int main()
{
//freopen("E:/input.txt","r",stdin);
//freopen("E:/output.txt","w",stdout);
int t;
scanf("%d",&t);
while(t--){
scanf("%s",s);
int len=strlen(s);
int lenm=1,maxlen=len;
int ok=0;
/**超时O(n^3)
for(int i=0;i<len;i++){
for(int j=i+1;j<len;j++){
if(huiWen1(s,i,j)&&lenm<j-i+1){
lenm=j-i+1;
}
}
}
cout<<lenm<<endl;
**/
/**
//从大到小枚举长度
while(!ok){//还是超时,O(n^3)系数可能会小一些
for(int i=0;i+maxlen<=len;i++){
if(huiWen1(s,i,i+maxlen-1)){
ok=1;//找发最长回文
break;
}
}
--maxlen;//遍历下一个长度
}
cout<<maxlen+1<<endl;
**/
/**
//枚举中心位置(考虑奇偶)O(n^2)超时
for(int i=1;i<len-1;i++){
lenm=max(lenm,huiWen2(s,i));
}
cout<<lenm<<endl;
**/
Manacher(s,len);
int ans=0;
for(int i=0;i<2*len+2;i++){
ans=max(ans,Mp[i]-1);
}
printf("%d\n",ans);
}
return 0;
}
注:参考文章和kuangbin的代码。(小猫的重学,请大神见谅!)