POJ3974 Palindrome (manacher算法)

题目大意就是说在给定的字符串里找出一个长度最大的回文子串。

才开始接触到manacher，不过这个算法的确很强大，这里转载了一篇有关manacher算法的讲解，可以去看看：地址

神器：

 #include <map>
 #include <set>
 #include <stack>
 #include <queue>
 #include <cmath>
 #include <ctime>
 #include <vector>
 #include <cstdio>
 #include <cctype>
 #include <cstring>
 #include <cstdlib>
 #include <iostream>
 #include <algorithm>
 using namespace std;
 #define INF 0x3f3f3f3f
 #define inf (-((LL)1<<40))
 #define lson k<<1, L, mid
 #define rson k<<1|1, mid+1, R
 #define mem0(a) memset(a,0,sizeof(a))
 #define mem1(a) memset(a,-1,sizeof(a))
 #define mem(a, b) memset(a, b, sizeof(a))
 #define FOPENIN(IN) freopen(IN, "r", stdin)
 #define FOPENOUT(OUT) freopen(OUT, "w", stdout)

 template<class T> T CMP_MIN(T a, T b) { return a < b; }
 template<class T> T CMP_MAX(T a, T b) { return a > b; }
 template<class T> T MAX(T a, T b) { return a > b ? a : b; }
 template<class T> T MIN(T a, T b) { return a < b ? a : b; }
 template<class T> T GCD(T a, T b) { return b ? GCD(b, a%b) : a; }
 template<class T> T LCM(T a, T b) { return a / GCD(a,b) * b;    }

 //typedef __int64 LL;
 typedef long long LL;
 const int MAXN = +;
 const int MAXM = ;
 const double eps = 1e-;

 char str[MAXN];
 int rad[MAXN];

 int manacher()
 {
     mem0(rad);
     for(int i=,j=,k,len=strlen(str)<<;i<len;)
     {
         while (str[i-j-]==str[i+j+]) j++;  //扫描得出rad值
         rad[i]=j;
         for (k=; k<=j && rad[i-k]!=rad[i]-k; k++) rad[i+k]=min(rad[i-k],rad[i]-k);  //k指针扫描
         i+=k;  //i跳到下一个需要计算rad值的位置
         j=max(j-k,);  //更新下一个rad值的初始值
     }
     int ans = ;
     for(int i=;i<strlen(str);i++)
     {
         int x ;
         if(str[i] != '#') x = rad[i]/*+;
         else x = (rad[i]+)/*;
         ans = max(ans, x);
     }
     return ans;
 }

 int main()
 {
     int T = ;
     while(gets(str)!=NULL && str[] != 'E')
     {
         int len = strlen(str);
         for(int i=len-;i>=;i--)
         {
             str[(i<<)+] = '#';
             str[(i<<)+] = str[i];
         }
         str[] = '?'; str[len<<] = '*';
         //printf("%s\n", str);
         int ans = manacher();
         printf("Case %d: %d\n", ++T, ans);
         mem0(str);
     }
     return ;
 }