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 = 4000000+10;
const int MAXM = 2000010;
const double eps = 1e-12;

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

int manacher()
{
    mem0(rad);
    for(int i=1,j=0,k,len=strlen(str)<<1;i<len;)
    {
        while (str[i-j-1]==str[i+j+1]) j++;  //扫描得出rad值
        rad[i]=j;
        for (k=1; 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,0);  //更新下一个rad值的初始值
    }
    int ans = 0;
    for(int i=1;i<strlen(str);i++)
    {
        int x ;
        if(str[i] != '#') x = rad[i]/2*2+1;
        else x = (rad[i]+1)/2*2;
        ans = max(ans, x);
    }
    return ans;
}

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