[APIO2014] [Uoj103] [Bzoj3676] Palindromes回文串 [Manacher，后缀数组]

用Manacher算法枚举回文子串，每次在后缀数组排序后的后缀数组中二分，因为用某一后缀和其他子串分别求匹配的长度，匹配长度在排序后该后缀的两侧具有单调性（匹配长度为min{H[x]|i<=x<=j}，所以对于查询min(H[x])用ST表O(n)预处理，O(1)查询即可。Manacher时间复杂度O(n),后缀数组复杂度O(nlogn)，总复杂度O(nlogn)。注意二分时的边界条件！

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <algorithm>

using namespace std;

long long    Ans;
int    A[],B[],U[];
int    SA[],*Rank,H[],Tmp[];
int    lg2[],ST[][],p[];
char    str[];

void    Get_H(const int n)
{
    int    i,j,k=;
    for(i=;i<n;H[Rank[i++]]=k)
        for(k?k--:k,j=SA[Rank[i]-];str[i+k]==str[j+k];++k);
    for(i=;i<=n;++i)lg2[i]=lg2[i>>]+;
    for(i=;i<=n;++i)ST[i][]=H[i];
    for(j=;(<<j)<=n;++j)
    {
        for(i=;i+(<<j)-<=n;++i)
        {
            ST[i][j]=min(ST[i][j-],ST[i+(<<(j-))][j-]);
        }
    }
    return ;
}

int    Query(const int l,const int r)
{
    int    temp=lg2[r-l+];
    return min(ST[l][temp],ST[r-(<<temp)+][temp]);
}

bool    cmp(const int * s,const int a,const int b,const int l)
{
    return s[a]==s[b] && s[a+l]==s[b+l];
}

int*    Get_SA(const int n,int    m)
{
    int    i,j,_p,*x=A,*y=B;
    for(i=;i<m;++i)U[i]=;
    for(i=;i<n;++i)U[x[i]=str[i]]++;
    for(i=;i<m;++i)U[i]+=U[i-];
    for(i=n-;i>=;--i)SA[--U[x[i]]]=i;

    for(j=,_p=;_p<n;m=_p,j<<=)
    {
        for(_p=,i=n-j;i<n;++i)y[_p++]=i;
        for(i=;i<n;++i)if(SA[i]>=j)y[_p++]=SA[i]-j;
        for(i=;i<n;++i)Tmp[i]=x[y[i]];
        for(i=;i<m;++i)U[i]=;
        for(i=;i<n;++i)U[Tmp[i]]++;
        for(i=;i<m;++i)U[i]+=U[i-];
        for(i=n-;i>=;--i)SA[--U[Tmp[i]]]=y[i];
        for(swap(x,y),_p=,x[SA[]]=,i=;i<n;++i)
            x[SA[i]]=cmp(y,SA[i-],SA[i],j)?_p-:_p++;
    }
    return x;
}

long long    Calc(int l,int r,const int n)
{
    l=(l-)>>,r=(r-)>>;
    int    pos=Rank[l],L,R,temp=;

    L=,R=pos;

    while(L<R-)
    {
        int    mid=L+((R-L)>>);
        if(Query(mid+,pos)>=r-l+)R=mid;
        else    L=mid;
    }

    temp+=pos-R;
    L=pos,R=n+;

    while(L<R-)
    {
        int    mid=L+((R-L)>>);
        if(Query(pos+,mid)>=r-l+)L=mid;
        else    R=mid;
    }
    temp+=L-pos;
    return (long long)temp*(r-l+);
}

void    Manacher(const int n)
{
    int    i,pos=;
    for(i=n-;i>=;--i)
    {
        str[i+i+]=str[i];
        str[i+i+]='#';
    }
    str[]='^';str[n<<|]='#';str[(n+)<<]='$';
    for(i=;i<=(n<<|);++i)
    {
        if(p[pos]+pos>i)
            p[i]=min(p[(pos<<)-i],p[pos]+pos-i);
        else    p[i]=;
        while(str[i-p[i]]==str[i+p[i]])
        { if(i+p[i]>p[pos]+pos) Ans=max(Ans,Calc(i-p[i],i+p[i],n)); p[i]++; }
        if(pos+p[pos]<i+p[i])pos=i;
    }

    return ;
}

int main()
{
    int    n;

    scanf("%s",str);

    n=strlen(str);
    Rank=Get_SA(n+,);
    Get_H(n);
    Manacher(n);

    printf("%lld",Ans);

    return ;
}