关于manacher
  • 板块学术版
  • 楼主冷笑叹秋萧
  • 当前回复2
  • 已保存回复2
  • 发布时间2021/11/19 08:37
  • 上次更新2023/11/4 00:11:07
查看原帖
关于manacher
37885
冷笑叹秋萧楼主2021/11/19 08:37

以下是manachermanacher算法的板子

#include<bits/stdc++.h>
#define N 30000005
#define R register int
using namespace std;
char s[N];
int p[N],ans;
int main()
{
	char ch=getchar();s[0]='@';int len=1;s[1]='#';
	while (ch>='a' && ch<='z') s[++len]=ch,s[++len]='#',ch=getchar();s[len+1]='*';
	int mx=0,mid=0;
	for (R i=1;i<=len;++i)
	{
		if (i<mx) p[i]=min(p[(mid<<1)-i],mx-i);else p[i]=1;
		while (s[i+p[i]]==s[i-p[i]]) ++p[i];
		if (i+p[i]>mx) mx=i+p[i],mid=i;ans=max(ans,p[i]-1);
	}
	printf("%d\n",ans);
	return 0;
}

以字符串aabbccbbdaabbccbbd为例

跑出来的回文半径如下:

aabbccbbd*a*a*b*b*c*c*b*b*d*

1232123212721232121

也就是说当前这个回文中心的位置也会产生一个贡献,那代码第12行的mximx-i不应该是mxi+1mx-i+1吗。以及第16行的i+p[i]i+p[i]不应该是i+p[i]1i+p[i]-1吗?(虽然跑出来的结果是一样的,因为往后也会更新到,但是这样的打法本质上是有问题的?)

2021/11/19 08:37
加载中...