最長迴文子串LPS（Hihocoder1032）

題目來源：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;
}