BZOJ 2565 最长回文串

最新推荐文章于 2019-08-04 16:04:43 发布

weixin_33738555

最新推荐文章于 2019-08-04 16:04:43 发布

阅读量115

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/hanyuweining/p/10321910.html

本文详细介绍了一种高效的回文自动机算法，该算法通过两次遍历字符串，分别从前向后和从后向前构建自动机，以计算字符串中所有回文子串的长度。文章附带了完整的C++代码实现，展示了算法的高效性和简洁性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门

回文自动机！

正着跑一遍记录以每个点作为回文子串的右端点的最大长度

倒过来跑一遍记录每个点作为左端点的最大长度

求个和就好啦

附代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 20021225
#define ll long long
#define mxn 100010
using namespace std;

struct node{int fa,ch[26],len;}t[mxn];
int poi,lt,f[mxn],s[mxn],n,len,ans;
char ch[mxn];
void init()
{
	memset(s,0,sizeof(s));
	for(int i=0;i<=poi;i++)	memset(t[i].ch,0,sizeof(t[i].ch)),t[i].fa=t[i].len=0;
	s[0] = -1; t[1].fa = t[0].fa =1;
	t[1].len = -1; t[0].len = 0;
	poi = lt = 1; n=0;
}
int id(char c){return c-'a';}
void extend(int c,int b,int id)
{
	//printf("%d\n",c);
	int p = lt; s[++n] = c;
	while(s[n-1-t[p].len]!=c)	p=t[p].fa;
	if(!t[p].ch[c])
	{
		int q = ++poi; t[q].len = t[p].len+2;
		int np = t[p].fa;
		while(s[n-1-t[np].len]!=c)	np = t[np].fa;
		t[q].fa = t[np].ch[c]; t[p].ch[c] = q;
	}
	lt = t[p].ch[c];
	if(b) f[id+1] = t[lt].len;
	else ans=max(f[id]+t[lt].len,ans);
}
int main()
{
	scanf("%s",ch); len = strlen(ch);
	init(); for(int i=0;i<len;i++)	extend(id(ch[i]),1,i);
	init(); for(int i=len-1;~i;i--) extend(id(ch[i]),0,i);
	printf("%d\n",ans);
	return 0;
}

这玩意跑的奇快无比。

转载于:https://www.cnblogs.com/hanyuweining/p/10321910.html