题目来源:hihicoder最长回文子串
题目中的题解没看懂,在网上看了别人的BLOG,看懂了这中转移方程 : p[i] >= min(mxlen-i,p[2*id-i] ) ;
在hihocoder中用string会超时,所以建议用char[]
下面代码用了5种方法
#include<iostream>
#include<cstring>
using namespace std;
#define MAX 2001007
string lps0(string s){
int pos=0,len=0;
//以i开始,以j结尾的子串 O(n^3) ,TLE
for(int i=0;i<s.length();++i)
for(int j=i;j<s.length();++j){
string sub=s.substr(i,j+1);
bool flag=true;
for(int k=0;flag&&k<sub.length()/2;++k)
if(sub[k] != sub[sub.length()-1-k])flag=false;
if(flag&&sub.length()>len)pos=i,len=sub.length();
}
return s.substr(pos,len);
}
string lps1(string s){
int pos=0,len=0;
//O(n^2) ,TLE
for(int i=0;i<s.length();++i){
int tpos=i;
//以i为中心
for(int j=i-1,k=i+1;j>=0&&k<s.length();--j,++k)
if(s[j] != s[k] ){
tpos=j+1;
break;
}else if(j == 0 || k == s.length()-1){
tpos=j;
break;
}
if((i-tpos)*2+1 > len){
len = (i-tpos)*2+1;
pos = tpos;
}
//不以i为中心
tpos=i+1;
for(int j=i,k=i+1;j>=0&&k<s.length();--j,++k)
if(s[j] != s[k] ){
tpos=j+1;
break;
}else if(j == 0 || k == s.length()-1){
tpos=j;
break;
}
if((i-tpos+1)*2 > len){
len = (i-tpos+1)*2;
pos = tpos;
}
}
return s.substr(pos,len);
}
char s[MAX],t[MAX];
//对上述方法的改进 AC
int lps11(){
int mlen=0,len=strlen(s);
for(int i=0;i<len;++i){
int b=i,e=i;
while(s[i] == s[e+1])++e;
i=e;
while(b>=0&&e<len&&s[b] == s[e])b--,e++;
if(e-b+1 > mlen)mlen=e-b+1;
}
return mlen-2;
}
int p[MAX];
int lps2(){
int len = strlen(s);
//初始化
t[0]='$';
for(int i=0;i<len;++i)
t[2*i+1]='#',t[2*i+2]=s[i];
len = 2*len+2;
t[len-1]='#',t[len]='\0';
// cout<<t<<endl;
int pos=0,mx=0,mlen=0,mpos;
for(int i=0;i<len;i++){
if(mx>i)p[i]=min(mx-i,p[2*pos-i]);
else p[i]=1;
while(t[i-p[i]] == t[i+p[i]])p[i]++; //可以不用担心越界访问的问题,因为t[0]='$',而t[length]='\0'
if(i+p[i] > mx){
mx = i+p[i];
pos = i;
}
if(mlen < p[i]){
mlen=p[i];
mpos=i;
}
}
// for(int j=mpos-mlen+1;j<mpos+mlen;++j)
// if(t[j] != '#')
// ans = ans + t[j];
return mlen-1;
}
//用string会超时 TLE
int lps3(string& s){
string t="$",ans="";
//初始化
for(int i=0;i<s.length();++i)
t = t + "#" + s[i];
t = t + "#";
//cout<<t<<endl;
int pos=0,mx=0,mlen=0,mpos;
for(int i=0;i<t.length();i++){
if(mx>i)p[i]=min(mx-i,p[2*pos-i]);
else p[i]=1;
while(t[i-p[i]] == t[i+p[i]])p[i]++; //可以不用担心越界访问的问题,因为t[0]='$',而t[length]='\0'
if(i+p[i] > mx){
mx = i+p[i];
pos = i;
}
if(mlen < p[i]){
mlen=p[i];
mpos=i;
}
}
// for(int j=mpos-mlen+1;j<mpos+mlen;++j)
// if(t[j] != '#')
// ans = ans + t[j];
return mlen-1;
}
int main(){
int N;
cin>>N;
while(N--){
cin>>s;
// sub = lps2(s);
// cout<<sub<<" ";
// cout<<sub.length()<<endl;
cout<<lps11()<<endl;
}
return 0;
}