每找到一个回文串,就在所有的串中查找出现了多少次
因为暴力跳非常的慢,所以用倍增优化
f[i][j]表示从第i个节点向上跳2^j步到哪里
每次查询都是从末尾节点开始,倍增找到最后一个长度大于等于p的节点manacher算法证明了本质不同的回文串只有O(n)个,复杂度O(nlogn)
原来manacher可以直接求偶数长度的回文串呀
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 300010
using namespace std;
int ch[2*maxn][26],size[2*maxn],fa[20][2*maxn],mx[2*maxn];
char s[maxn];
int p[maxn],pos[maxn],q[2*maxn],v[maxn];
long long ans;
int n,m,tot=1,last=1;
void insert(int x,int id)
{
int p=last,np=last=++tot;
mx[np]=mx[p]+1;pos[id]=np;size[np]=1;
while (p && !ch[p][x]) ch[p][x]=np,p=fa[0][p];
if (!p) fa[0][np]=1;
else
{
int q=ch[p][x];
if (mx[q]==mx[p]+1) fa[0][np]=q;
else
{
int nq=++tot;mx[nq]=mx[p]+1;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
fa[0][nq]=fa[0][q];
fa[0][np]=fa[0][q]=nq;
while (p && ch[p][x]==q) ch[p][x]=nq,p=fa[0][p];
}
}
}
void bfs()
{
for (int i=1;i<=tot;i++) v[mx[i]]++;
for (int i=1;i<=n;i++) v[i]+=v[i-1];
for (int i=tot;i>=1;i--) q[v[mx[i]]--]=i;
for (int i=tot;i>=1;i--) size[fa[0][q[i]]]+=size[q[i]];
for (int j=1;j<=19;j++)
for (int i=1;i<=tot;i++)
fa[j][i]=fa[j-1][fa[j-1][i]];
}
void query(int x,int d)
{
for (int i=19;i>=0;i--)
if (mx[fa[i][x]]>=d) x=fa[i][x];
ans=max(ans,1ll*size[x]*d);
}
int main()
{
scanf("%s",s+1);n=strlen(s+1);
for (int i=1;i<=n;i++) insert(s[i]-'a',i);
s[0]='-';s[n+1]='+';
bfs();
int mx=0,id=0;
for (int i=1;i<=n;i++)
{
if (mx>i) p[i]=min(mx-i,p[2*id-i]); else p[i]=1,query(pos[i],1);
while (s[i+p[i]]==s[i-p[i]]) p[i]++,query(pos[i+p[i]-1],2*p[i]-1);
if (i+p[i]>mx) mx=i+p[i],id=i;
}
mx=0,id=0;
for (int i=1;i<=n;i++)
{
if (mx>i) p[i]=min(mx-i,p[2*id-i]); else p[i]=0;
while (s[i+p[i]+1]==s[i-p[i]]) p[i]++,query(pos[i+p[i]],2*p[i]);
if (i+p[i]>mx) mx=i+p[i],id=i;
}
printf("%lld\n",ans);
return 0;
}