P3805 【模板】manacher算法

木落淮南，雨晴雲夢

于 2021-03-28 10:18:13 发布

阅读量145

点赞数

CC 4.0 BY-SA版权

分类专栏：洛谷 Manacher

本文链接：https://blog.csdn.net/lcl1234567890lcl/article/details/115276988

洛谷同时被 2 个专栏收录

29 篇文章

订阅专栏

Manacher

1 篇文章

订阅专栏

本文介绍了一种字符串处理算法——Manacher's Algorithm，用于高效地找到重复字符子串。通过C++实现，它优化了字符串匹配过程，适用于长度庞大的输入。核心在于利用已知信息加速动态规划，提高了查找最长回文子串的效率。

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

链接

//#pragma GCC optimize(2)
//#pragma GCC optimize(3,"Ofast","inline")
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <queue>
#include <deque>
#include <algorithm>
#include <cmath>
#include <vector>
#include <stack>
#include <bitset>
#include <set>
#include <map>
#define inf 0x7fffffff
#define ll long long
//#define int long long
//#define double long double
#define eps 1e-5
//#define mod 1e9+7
using namespace std;
const int mod=10007;
const int M=10007;
const int N=5*1e7+5;//空间最大限制 4e8
char s[N];
char s1[N];
int len;
int p[N];
void init()
{
	len=2;
//	memset(s,0,strlen(s));
	s[0]='@';
	s[1]='#';
	int sz=strlen(s1);
	for(int i=0;i<sz;i++)
	{
		s[len++]=s1[i];
		s[len++]='#';
	}
	s[len]='\0';
}
int manacher()
{
	int ans=-1,id=0,mx=0;
//	cout<<len<<endl;
	memset(p,0,sizeof(p));
	for(int i=1;i<len;i++)
	{
		if(i<mx)  p[i]=min(p[id*2-i],mx-i);
		else p[i]=1;
		while(s[i-p[i]]==s[i+p[i]])  p[i]++;
		if(p[i]+i>mx)
		{
			mx=p[i]+i;
			id=i;
		}
		ans=max(ans,p[i]-1);
	}
	return ans;
}
signed main()
{
//	ios::sync_with_stdio(false);
    while(scanf("%s",s1)!=EOF)
    {
    	init();
    	printf("%d\n",manacher());
    }
    return 0;
}