【APIO2014】bzoj3676 回文串【解法二】

Description

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

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

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

解法一【manacher+后缀数组】见【这里】。
解法三【pam】见【这里】。
manacher找出来本质不同的回文串【O(n)个】然后到SAM上查询出现次数。

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;#define LL long longchar s[600010];int fa[600010][22],trans[600010][27],val[600010],cnt[600010],num[600010],ord[600010],pos[600010],f[600010],n,tot=1;int main(){    int p,q,np,nq,last=1,c,id,mx;    LL ans=1;    scanf("%s",s+1);    n=strlen(s+1);    for (int i=1;i<=n;i++)    {        p=last;        val[last=np=++tot]=val[p]+1;        c=s[i]-'a'+1;        while (p&&!trans[p][c])        {            trans[p][c]=np;            p=fa[p][0];        }        if (!p) fa[np][0]=1;        else        {            q=trans[p][c];            if (val[q]==val[p]+1) fa[np][0]=q;            else            {                val[nq=++tot]=val[p]+1;                fa[nq][0]=fa[q][0];                fa[q][0]=fa[np][0]=nq;                for (int j=1;j<=26;j++) trans[nq][j]=trans[q][j];                while (p&&trans[p][c]==q)                {                    trans[p][c]=nq;                    p=fa[p][0];                }            }        }    }    p=1;    for (int i=1;i<=n;i++)    {        p=trans[p][s[i]-'a'+1];        pos[i]=p;        cnt[p]++;    }    for (int i=1;i<=tot;i++) num[val[i]]++;    for (int i=1;i<=n;i++) num[i]+=num[i-1];    for (int i=1;i<=tot;i++) ord[num[val[i]]--]=i;    for (int i=tot;i;i--) cnt[fa[ord[i]][0]]+=cnt[ord[i]];    for (int k=1;(1<<k)<=tot;k++)        for (int i=1;i<=tot;i++)            fa[i][k]=fa[fa[i][k-1]][k-1];    for (int i=n;i;i--)    {        s[i*2-1]=s[i];        s[i*2-2]='$';    }    s[2*n]='$';    id=mx=0;    for (int i=1;i<=n*2-1;i++)    {        if (i<=mx) f[i]=min(mx-i,f[2*id-i]);        while (i-f[i]-1>=0&&s[i+f[i]+1]==s[i-f[i]-1]) f[i]++;        for (int j=mx-i+2;j<=f[i];j+=2)        {            p=pos[(i+1)/2+j/2];            for (int k=20;k>=0;k--)                if (val[fa[p][k]]>=j)                    p=fa[p][k];            ans=max(ans,(LL)cnt[p]*j);        }        if (i+f[i]>mx)        {            id=i;            mx=i+f[i];        }    }    printf("%lld\n",ans);}