#1032 : 最長回文子串 （Manacher算法）

這篇博文對Manacher講的還不錯

#1032 : 最長回文子串

時間限制:1000ms

單點時限:1000ms

記憶體限制:64MB

描述

小Hi和小Ho是一對好朋友，出生在資訊化社會的他們對程式設計産生了莫大的興趣，他們約定好互相幫助，在程式設計的學習道路上一同前進。

這一天，他們遇到了一連串的字元串，于是小Hi就向小Ho提出了那個經典的問題：“小Ho，你能不能分别在這些字元串中找到它們每一個的最長回文子串呢？”

小Ho奇怪的問道：“什麼叫做最長回文子串呢？”

小Hi回答道：“一個字元串中連續的一段就是這個字元串的子串，而回文串指的是12421這種從前往後讀和從後往前讀一模一樣的字元串，是以最長回文子串的意思就是這個字元串中最長的身為回文串的子串啦~”

小Ho道：“原來如此！那麼我該怎麼得到這些字元串呢？我又應該怎麼告訴你我所計算出的最長回文子串呢？

小Hi笑着說道：“這個很容易啦，你隻需要寫一個程式，先從标準輸入讀取一個整數N（N<=30)，代表我給你的字元串的個數，然後接下來的就是我要給你的那N個字元串（字元串長度<=10^6)啦。而你要告訴我你的答案的話，隻要将你計算出的最長回文子串的長度按照我給你的順序依次輸出到标準輸出就可以了！你看這就是一個例子。”

提示一 提示二 提示三 提示四

樣例輸入

3

abababa

aaaabaa

acacdas

樣例輸出

7

5

3

#include<bits/stdc++.h>
using namespace std;
template<class T>inline T read(T&x)
{
    char c;
    while((c=getchar())<=32)if(c==EOF)return 0;
    bool ok=false;
    if(c=='-')ok=true,c=getchar();
    for(x=0; c>32; c=getchar())
        x=x*10+c-'0';
    if(ok)x=-x;
    return 1;
}
template<class T> inline T read_(T&x,T&y)
{
    return read(x)&&read(y);
}
template<class T> inline T read__(T&x,T&y,T&z)
{
    return read(x)&&read(y)&&read(z);
}
template<class T> inline void write(T x)
{
    if(x<0)putchar('-'),x=-x;
    if(x<10)putchar(x+'0');
    else write(x/10),putchar(x%10+'0');
}
template<class T>inline void writeln(T x)
{
    write(x);
    putchar('\n');
}
//-------ZCC IO template------
const int maxn=1e6+1000;
const double inf=999999999;
#define lson (rt<<1),L,M
#define rson (rt<<1|1),M+1,R
#define M ((L+R)>>1)
#define For(i,t,n) for(int i=(t);i<(n);i++)
typedef long long  LL;
typedef double DB;
typedef pair<int,int> P;
#define bug printf("---\n");
#define mod  100000000

char s[maxn];
char a[maxn*2];
int d[maxn*2];
int init()
{
    int n,len=strlen(s);
    a[0]='$';
    a[1]='#';
    for(int i=0;i<len;i++)
    {
        a[i*2+2]=s[i];
        a[i*2+3]='#';
    }
    n=len*2+2;
    a[n]=0;
    return n;
}
void Manacher(int len)
{
    int pos=0,mx=0;
    for(int i=1;i<len;i++)
    {
        if(mx>i)
            d[i]=min(d[pos*2-i],mx-i);
        else
            d[i]=1;
        while(a[i+d[i]]==a[i-d[i]])d[i]++;
        if(d[i]+i>mx)
        {
            mx=d[i]+i;
            pos=i;
        }
    }
}

int main()
{
    int T,n;
    read(T);
    while(T--)
    {
scanf("%s",s);
        n=init();Manacher(n);
        writeln(*max_element(d,d+n)-1);
    }
    return 0;
}