HDU 3068(manacher)

给出一个只由小写英文字符a,b,c…y,z组成的字符串S,求S中最长回文串的长度.
回文就是正反读都是一样的字符串,如aba, abba等
Input
输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c…y,z组成的字符串S
两组case之间由空行隔开(该空行不用处理)
字符串长度len <= 110000
Output
每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度.
Sample Input

aaaa

abab

Sample Output

4
3

解题思路:
解决这类最长回文子串的问题一般使用manacher算法,它可以O(n)的解决这类问题。
先说一下思路叭,先使用一个数组来存以i为中心的最长的回文串的长度,一个变量mx存目前的回文串能够延伸到的最右边界,也就是说存在一个i使得以i为中心的回文串是能够包含到第mx个元素的,那么如果我们遍历到某一个j时,如果j<mx,那么说明j是被一个以i为中心的回文串包含的,那么如果j>i,j相对于i对称的另一个的回文串我们已经知道了,此时以j为中心的回文串的长度我们可以知道一部分,就是min(p[2*i-j],mx-j),然后我们再循环去找剩下的回文串。
因为我们会用到对称这个思路,如果一个字符串的长度为偶数,我们是不能处理的,因为我们找不到对称点,我们枚举的是以i为中心的回文串,如果是偶数,我们就找不到i,所以我们要预处理一下,我们在每个字符中间加入一个字符,比如说’#’,举个栗子:原字符串是:abcab,预处理之后就变成了:#a#b#c#a#b,我们再加一个防止越界的字符‘$’就变成了:$#a#b#c#a#b,然后我们就可以枚举每一位了。
核心代码:

int manacher()
{
    int len=init();
    int id=0,mx=0,maxl=0;
    for(int i=1;i<len;i++){
        if(i<mx){
            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(mx<i+p[i]){//为了减少while循环的次数,更多的执行第一个if语句
            mx=i+p[i];
            id=i;
        }
        maxl=max(maxl,p[i]-1);
    }
    return maxl;
}

那么会了manacher这道题就套板子就可以了。
代码:

#pragma GCC optimize(2)
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string>
#include <vector>
#include <set>
#include <map>
#include <stack>
#include <bitset>
#include <queue>
//#include <random>
#include <time.h>
using namespace std;
#define int long long
#define ull unsigned long long
#define ls root<<1
#define rs root<<1|1
const int maxn=5e5+10;
const int maxm=3e6+10;
//std::mt19937 rnd(time(NULL));
int p[maxn];
char str1[maxn],str[maxn];
int init()
{
    int k=2;
    str[0]='$',str[1]='#';
    for(int i=0;str1[i]!='\0';i++){
        str[k++]=str1[i];
        str[k++]='#';
    }
    str[k]='\0';
    return k;
}
int malache()
{
    int len=init();
    int id=0,mx=0,maxl=0;
    for(int i=1;i<len;i++){
        if(i<mx){
            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(mx<i+p[i]){
            mx=i+p[i];
            id=i;
        }
        maxl=max(maxl,p[i]-1);
    }
    return maxl;
}
signed main()
{
    while(~scanf("%s",str1)){
        printf("%lld\n",malache());
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值