1. 程式人生 > >【Manacher模板】HDU 3068——求最長迴文子串

【Manacher模板】HDU 3068——求最長迴文子串

直接做會超時,需要優化,網上通行的演算法是manacher演算法(具體原理還不是很明白),這裡可以當模板使。

// 原串最大長度N
// 返回最大回文字串 res
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<string>
#include<algorithm>
using namespace std;
const int N = 110009;

int rid[N<<2];
string manacher(char *s){
    string t, res;    
    int L = strlen(s);    
    // init
    t += '?';
    for(int i = 0; i < L; i++)
        t += '#', t += s[i];
    t += "#*";
    // getrid
    for(int i=1,j=0,k,End=(int)t.size();i<End;){
        while( t[i-j-1]==t[i+j+1] ) j++;
        rid[i] = j;
        for(k=1;k<=j&&(rid[i-k]!=rid[i]-k);k++) 
            rid[i+k] = min( rid[i-k], rid[i]-k );
        i += k;
        j = max( 0,j-k );
    }    
    // Max
    int m = 0, pos;
    for(int i = 1; i < (int)t.size(); i++){
        if( m < rid[i]*2+(t[i]!='#') )
            m = rid[i]*2+(t[i]!='#'), pos = i;
    }
    for(int i = pos-rid[pos]; i <= pos+rid[pos]; i++){
        if( t[i] != '#' ) res += t[i];    
    }
    return res;
}


int main()
{
    char s[N];
    while(scanf("%s", s)!=EOF)
    {
    	 string ans = manacher(s); 
	//    printf("MaxLength = %d\n", ans.size() );
	      printf("%d\n", ans.size() );
    		
    }
   return 0;
}