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

标签：

这篇博文对Manacher讲的还不错

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;
}

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

标签：

原文地址：http://blog.csdn.net/u013167299/article/details/45217957