【APIO2014】【BZOJ3676】回文串

Description

考虑一个只包含小写拉丁字母的字符串s。我们定义s的一个子串t的“出
现值”为t在s中的出现次数乘以t的长度。请你求出s的所有回文子串中的最
大出现值。

Input

输入只有一行，为一个只包含小写字母(a -z)的非空字符串s。

Output

输出一个整数，为逝查回文子串的最大出现值。

Sample Input

【样例输入l】

abacaba

【样例输入2]

www

Sample Output

【样例输出l】

7

【样例输出2]

4

HINT

一个串是回文的，当且仅当它从左到右读和从右到左读完全一样。

在第一个样例中，回文子串有7个：a，b，c，aba，aca，bacab，abacaba，其中：

● a出现4次，其出现值为4：1：1=4

● b出现2次，其出现值为2：1：1=2

● c出现1次，其出现值为l：1：l=l

● aba出现2次，其出现值为2：1：3=6

● aca出现1次，其出现值为1=1：3=3

●bacab出现1次，其出现值为1：1：5=5

● abacaba出现1次，其出现值为1：1：7=7

故最大回文子串出现值为7。

【数据规模与评分】

数据满足1≤字符串长度≤300000。

Source

为了做这题特地学了Manacher…
用的是黄学长的做法所以我就不再赘述了..
黄学长题解传送门

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 300010
#define LL long long
using namespace std;
int n;
LL ans;
int p[MAXN],c[MAXN],tmp[MAXN<<1],Pow[20]={1};
char ch[MAXN];
struct sam
{
    int last,cnt,p,q,np,nq;
    int a[MAXN<<1][26],fa[MAXN<<1],len[MAXN<<1],right[MAXN<<1],pos[MAXN<<1],deep[MAXN<<1];
    int Fa[MAXN<<1][20];
    sam()   {last=++cnt;}
    void insert(int c,int id)
    {
        p=last;np=last=++cnt;len[np]=len[p]+1;right[np]=1;pos[id]=last;
        while (!a[p][c]&&p) a[p][c]=np,p=fa[p];
        if (!p) fa[np]=1;
        else
        {
            q=a[p][c];
            if (len[q]==len[p]+1)   fa[np]=q;
            else
            {
                nq=++cnt;len[nq]=len[p]+1;memcpy(a[nq],a[q],sizeof(a[q]));
                fa[nq]=fa[q];fa[q]=fa[np]=nq;
                while (a[p][c]==q)  a[p][c]=nq,p=fa[p];
            }
        }
    }
    void init()
    {
        for (int i=1;i<=cnt;i++)    c[len[i]]++;
        for (int i=1;i<=n;i++)  c[i]+=c[i-1];
        for (int i=cnt;i;i--)   tmp[c[len[i]]--]=i;
        for (int i=cnt;i;i--)   right[fa[tmp[i]]]+=right[tmp[i]];
        for (int i=1;i<=cnt;i++)
        {
            int x=tmp[i];
            deep[x]=deep[fa[x]]+1;Fa[x][0]=fa[x];
            for (int j=1;Pow[j]<=deep[x];j++)   Fa[x][j]=Fa[Fa[x][j-1]][j-1];
        }
    }
    void update(int l,int r)
    {
        int mid=pos[r];
        for (int i=19;i>=0;i--)
        {
            int x=Fa[mid][i];
            if (len[x]>=r-l+1)  mid=x;
        }
        ans=max(ans,(LL)(right[mid])*(r-l+1));
    }
}sam;
void manacher()
{
    int maxn=0,id=0;
    for (int i=1;i<=n;i++)
    {
        if (maxn>i) p[i]=min(p[2*id-i-1],maxn-i);
        else    p[i]=0;
        while (ch[i+p[i]+1]==ch[i-p[i]])    ++p[i],sam.update(i-p[i]+1,i+p[i]);
        if (p[i]+i>maxn)    maxn=p[i]+i,id=i;
    }
    maxn=0;
    for (int i=1;i<=n;i++)
    {
        if (maxn>i) p[i]=min(p[2*id-i],maxn-i-1);
        else    p[i]=1,sam.update(i-p[i]+1,i+p[i]-1);
        while (ch[i+p[i]]==ch[i-p[i]])  ++p[i],sam.update(i-p[i]+1,i+p[i]-1);
        if (p[i]+i>maxn)    maxn=p[i]+i,id=i;
    }
}
int main()
{
    for (int i=1;i<=19;i++) Pow[i]=Pow[i-1]<<1;
    scanf("%s",ch+1);n=strlen(ch+1);
    for (int i=1;i<=n;i++)  sam.insert(ch[i]-'a',i);
    sam.init();ch[0]='+';ch[n+1]='*';manacher();
    cout<<ans<<endl;
}