【bzoj3676】[Apio2014]回文串 后缀自动机+倍增+manacher

每找到一个回文串，就在所有的串中查找出现了多少次

因为暴力跳非常的慢，所以用倍增优化

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;
}