BZOJ3676 APIO2014回文串（manacher+后缀自动机）

　　由于本质不同的回文子串数量是O(n)的，考虑在对于每个回文子串在第一次找到它时对其暴力统计。可以发现manacher时若右端点移动则找到了一个新回文串。注意这样会漏掉串长为1的情况，特判一下。

　　现在问题变为统计一个子串的出现次数。可以用SA，二分乱搞一下即可。这里使用SAM。以parent树上表示该子串的节点为起点，用倍增往上跳，找到深度最小的满足len限制的点就好了，出现次数就是其right集合的大小。

　　uojAC，luoguRE一个点，bzojMLE……

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 600010
int n,last,cnt,len[N],fail[N],son[N][],size[N],pos[N>>],p[N];
long long ans=;
char s[N];
namespace tree
{
    int p[N],t=,fa[N][];
    struct data{int to,nxt;
    }edge[N];
    void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
    void dfs(int k)
    {
        for (int i=p[k];i;i=edge[i].nxt)
        {
            fa[edge[i].to][]=k;
            dfs(edge[i].to);
            size[k]+=size[edge[i].to];
        }
    }
    void build()
    {
        for (int i=;i<=cnt;i++) addedge(fail[i],i);
        fa[][]=;dfs();
        for (int j=;j<;j++)
            for (int i=;i<=cnt;i++)
            fa[i][j]=fa[fa[i][j-]][j-];
    }
    int calc(int l,int r)
    {
        int x=pos[r];
        for (int j=;~j;j--) if (len[fa[x][j]]>=r-l+) x=fa[x][j];
        return size[x];
    }
}
using tree::calc;
void ins(int c,int n)
{
    int x=++cnt,p=last;last=x;len[x]=n;size[x]=;pos[n]=x;
    while (!son[p][c]&&p) son[p][c]=x,p=fail[p];
    if (!p) fail[x]=;
    else
    {
        int q=son[p][c];
        if (len[p]+==len[q]) fail[x]=q;
        else
        {
            int y=++cnt;len[y]=len[p]+;
            memcpy(son[y],son[q],sizeof(son[q]));
            fail[y]=fail[q];fail[q]=fail[x]=y;
            while (son[p][c]==q) son[p][c]=y,p=fail[p];
        }
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3676.in","r",stdin);
    freopen("bzoj3676.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    last=cnt=;
    scanf("%s",s+);n=strlen(s+);
    for (int i=;i<=n;i++) ins(s[i]-,i);
    tree::build();
    for (int i=n;i>=;i--) s[i*-]=s[i];
    for (int i=;i<n;i++) s[i<<]='$';
    int x=;
    for (int i=;i<=n;i++) ans=max(ans,1ll*calc(i,i));
    for (int i=;i<n*;i++)
    {
        if (x+p[x]>i) p[i]=min(x+p[x]-i,p[x-(i-x)]);
        while (i-p[i]->=&&i+p[i]+<n*&&s[i+p[i]+]==s[i-p[i]-])
        {
            p[i]++;
            if (s[i+p[i]]!='$') ans=max(ans,1ll*((i+p[i]>>)-(i-p[i]>>)+)*calc((i-p[i]>>)+,(i+p[i]>>)+));
        }
        if (i+p[i]>x+p[x]) x=i;
    }
    cout<<ans;
    return ;
}