最长回文子串之Manacher算法【模板】

最新推荐文章于 2022-10-18 11:04:44 发布

原创最新推荐文章于 2022-10-18 11:04:44 发布 · 228 阅读

0 ·

CC 4.0 BY-SA版权

算法专栏收录该内容

356 篇文章

订阅专栏

本文详细介绍马拉车(Manacher)算法的实现原理与步骤，通过C和C++代码示例展示了如何寻找字符串中最长的回文子串。适用于需要解决回文相关问题的开发者。

模板

C语言

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>

using namespace std;

const int N=1000010;
int p[2*N];
char s[2*N];

int Manacher()
{
	int mx=0,id=0,l=strlen(s+1);
	memset(p,0,sizeof(p));
	s[0]='@';
	for(int i=l;i>=1;i--) s[2*i]=s[i];
	l=2*l+1;//处理后字符串的长度 
	for(int i=1;i<=l;i+=2) s[i]='#';
	s[l+1]='%';
	for(int i=1;i<=l;i++)
	{
		p[i]=mx>i?min(p[2*id-i],mx-i):1;
		while(s[i+p[i]]==s[i-p[i]]) p[i]++;
		if(i+p[i]>mx)
		{
			id=i;
			mx=i+p[i];
		}
	}
	int ans=0;
	for(int i=1;i<=l;i++)
	  ans=max(ans,p[i]);
	return ans-1;
	
}

int main()
{
	int cnt=1;
	while(~scanf("%s",s+1))
	{
		printf("Case %d: %d\n",cnt++,Manacher());
	 } 
	 return 0;
}

C++

#include<iostream>
#include<string>
#include<cstring>

using namespace std;

const int N=100000;
int p[2*N+10];

int Manacher(string s)
{
	if(s=="")
	  return 0;
	int mx=0,id=0,len=s.length();
	memset(p,0,sizeof(p));
	for(int k=0;k<=len*2+1;k+=2)
	  s.insert(k,"#");
	for(int i=0;s[i];i++){
		p[i]=mx>i?min(p[2*id-i],mx-i):1;
		while((i+p[i]<s.length())&&(i-p[i]>=0)&&s[i+p[i]]==s[i-p[i]])
		  p[i]++;
		if(i+p[i]>mx){
			mx=i+p[i];
			id=i;
		}
	}
	int ans=1;
	for(int i=0;i<s.length();i++)
	  if(ans<p[i])
	    ans=p[i];
	return ans-1;
}

int main()
{
	string s;
	cin>>s;
	cout<<Manacher(s)<<endl;
	return 0;
}