HDU 3068 manacher算法(O(n)的最长回文串算法)

最新推荐文章于 2024-09-21 14:48:47 发布

humeay

最新推荐文章于 2024-09-21 14:48:47 发布

阅读量439

点赞数

分类专栏：字符串文章标签：算法 acm

本文链接：https://blog.csdn.net/qq_27925701/article/details/47991459

版权

字符串专栏收录该内容

5 篇文章 0 订阅

订阅专栏

宿舍的交换机坏了，因为之前还没开学，所以就没人修这里写图片描述 ,然后就拖到现在寝室才有网。然后这几天因为省赛选拔赛所以心好累
/**************************************************************

manacher是很久之前就听说的想学的，今天晚上因为选拔完了就闲了看的

数组:
s 字符原串
str 通过在每两个数组元素之间添加’#’形成的新串
p[i] 代表以str[i]为中心的回文串最大长度
元素：
mx 当前位置之前的回文串能达到最右边的位置
id 当前形成mx的回文串的中心点

/此处将末尾改为0纯粹是因为想跑FOR方便而已/
众所周知，回文串长度可以为奇，可以为偶，所以进行回文串判断过程中就会重复一些不必要的过程，manacher算法通过每两个元素之间添加新无关元素来实现长度为奇+偶=奇来统一长度减少时间复杂度

一般方法，跑到当前i节点然后向两边扩展，最坏时间复杂度为o(n^2)
manacher算法通过运用回文串的对称性来达到跑到i节点，可以从k长度开始扩展,k=与i点关于id点对称的j点的p[j],但是考虑到一种情况

这里写图片描述

所以,k=min(j-x,p[j])

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define rfor(i,a,b) for(i=a;i<=b;++i)
#define lfor(i,a,b) for(i=a;i>=b;--i)
#define sfor(i,a,h) for(i=h[a];i!=-1;i=e[i].next)
#define mem(a,b) memset(a,b,sizeof(a))
#define mec(a,b) memcpy(a,b,sizeof(b))
#define cheak(i) printf("%d ",i)
#define min(a,b) (a>b?b:a)
#define max(a,b) (a>b?a:b)
#define inf 0x7fffffff
typedef long long LL;
#define maxn 500005
char str[maxn],s[maxn];
int p[maxn];
int manacher(char* s,int len)
{
    int i,mx=0,l,id,ans=0;
    str[0]='@';str[1]='#';
    for(i=1;i<=len;++i)
    str[i<<1]=s[i],str[(i<<1)+1]='#';
    str[len*2+1]='#';str[len*2+2]=0;
    for(i=1;str[i];++i)
    {
        if(mx>i) p[i]=min(p[2*id-i],mx-i);
        else p[i]=1;
        while(str[i+p[i]]==str[i-p[i]])
        p[i]++;
        if(i+p[i]>mx)
        {
            mx=i+p[i];
            id=i;
        }
        ans=max(p[i],ans);
    }
    return ans-1;
}
int main()
{
    int len,i;
    while(~scanf("%s",s+1))
    {
        len=strlen(s+1);
        int ans=manacher(s,len);
        printf("%d\n",ans);
    }
    return 0;
}