Manacher算法O(n)求最长回文

最新推荐文章于 2021-08-23 23:27:31 发布

FightingForFinal

最新推荐文章于 2021-08-23 23:27:31 发布

阅读量281

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/daniel_csdn/article/details/51490942

版权

数论专栏收录该内容

22 篇文章

订阅专栏

本文介绍了一种使用Manacher算法高效求解字符串中最长回文子串的方法。通过预处理字符串并在每个位置上扩展以找到回文中心，最终确定最长回文子串的长度。

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

证明见此处
 51nod1089最长回文子串 V2（Manacher算法）

#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<string>
#include<sstream>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<deque>
#include<stack>
using namespace std;
const int maxn=100000;

int p[maxn*2+10];

int manacher(string &s){
    memset(p,0,sizeof(p));
    int id=0,mx=0;
    for(int i=0;i<s.size();i++){
        p[i]=mx>i?min(p[2*id-i],mx-i):1;

        while(i-p[i]>=0&&i+p[i]<=s.size()&&s[i-p[i]]==s[i+p[i]])p[i]++;

        if(mx<i+p[i])mx=i+p[i],id=i;
    }

    //寻找长回文的中心
    int Max=0,Maxi=0;
    for(int i=0;i<s.size();i++){
        if(Max<p[i])Max=p[i],Maxi=i;
    }

    //根据最长回文开头是否是#来算出长度
    int b=Maxi-Max+1;
    int e=Maxi+Max-1;
    return s[b]=='#'?(e-b)/2:(e-b)/2+1;

//    for(int i=b;i<=e;i++){//遍历一遍输出
//        if(s[i]!='#')cout<<s[i];
//    }
//    cout<<endl;
}

int main(){
    string s,ts;
    cin>>ts;
    s+=ts[0];
    for(int i=1;i<ts.size();i++){
        s+='#';        //不要写成s=s+'#'+ts[i],会超时
        s+=ts[i];
    }
    cout<<manacher(s);
}