Problem Description
给出一个只由小写英文字符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

题意:模板题,求最长回文子串的长度,字符串长度有11W,可以用manacher算法,注意输入的时候如果写成while(scanf("%s", a))会tle,应该写成while(scanf("%s", a) == 1)

manacher(马拉车)算法:复杂度O(n)

作用:给出一个字符串,求它的最长回文子串(的长度),求具体的串可以在模板上改一改

1.将字符串之间插入串中不可能有的字符,如将a a b a b a a b处理成 $ # a  # a # b # a # b # a # a # b #,

一个回文串是对称的,既然对称,就一定有那个“对称轴”,对于长度为奇数的回文,它的对称轴是一个字符,对于长度为偶数的回文串,它的对称轴是两个字符。为了方便处理,将所有回文串的对称轴变为一个字符,可以如上插入不相关字符,这样连上刚插入的字符,原串就变成了奇数长度。 最前面的  '$' 是为了防止溢出(我也不太懂)

2.设数组Mp[i]为以i为中心的回文串的半径长度,

i          0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17

ma[i]   $  #  a  #  a  #  b  #  a  #  b  #   a   #   a   #   b   #

mp[i]   0  0  1  2  1  0  3  1  7  ...

如此可以观察到mp[i]也代表 原串中 以i为中心的回文串的长度 ! 所以最长回文串的长度就是mp数组中的最大值了

3.求mp数组

用i将整个ma数组遍历一遍,在遍历的过程中维护一个已知的最长回文串,这个回文串要求它的右边界最靠右、右边界相同的话中心最靠左,

以下转载于https://www.felix021.com/blog/read.php?2040

该算法增加两个辅助变量(其实一个就够了,两个更清晰)id和mx,其中 id 为已知的 {右边界最大} 的回文子串的中心,mx则为id+mp[id],也就是这个子串的右边界。

然后可以得到一个非常神奇的结论,这个算法的关键点就在这里了:如果mx > i,那么P[i] >= MIN(P[2 * id - i], mx - i)。就是这个串卡了我非常久。实际上如果把它写得复杂一点,理解起来会简单很多:

//记j = 2 * id - i,也就是说 j 是 i 关于 id 的对称点(j = id - (i - id))
if (mx - i > mp[j]) 
    mp[i] = mp[j];
else /* mp[j] >= mx - i */
    mp[i] = mx - i; // mp[i] >= mx - i,取最小值,之后再匹配更新。

当然光看代码还是不够清晰,还是借助图来理解比较容易。

当 mx - i > P[j] 的时候,以S[j]为中心的回文子串包含在以S[id]为中心的回文子串中,由于 i 和 j 对称,以S[i]为中心的回文子串必然包含在以S[id]为中心的回文子串中,所以必有 P[i] = P[j],见下图。

当 P[j] >= mx - i 的时候,以S[j]为中心的回文子串不一定完全包含于以S[id]为中心的回文子串中,但是基于对称性可知,下图中两个绿框所包围的部分是相同的,也就是说以S[i]为中心的回文子串,其向右至少会扩张到mx的位置,也就是说 P[i] >= mx - i。至于mx之后的部分是否对称,就只能老老实实去匹配了。

对于 mx <= i 的情况,无法对 P[i]做更多的假设,只能P[i] = 1,然后再去匹配了。

于是代码如下:

//输入,并处理得到字符串s
int mp[1000], mx = 0, id = 0;
memset(mp, 0, sizeof(mp));
for (i = 1; s[i] != '\0'; i++) {
    p[i] = mx > i ? min(mp[2*id-i], mx-i) : 1;
    while (s[i + mp[i]] == s[i - mp[i]]) mp[i]++;
    if (i + mp[i] > mx) {
        mx = i + mp[i];
        id = i;
    }
}
//找出mp[i]中最大的

整个manacher算法就是这些,另外贴一下该题目的代码,也是manacher的模板,模板取自上海大学kuangbin

#include<cstdio>
#include<cstring>
#include<algorithm> #define clr(a, x) memset(a, x, sizeof(a))
typedef long long ll;
typedef long double ld; using namespace std; const int maxn = 110011;
char ma[maxn * 2];
int mp[maxn * 2]; void manacher(char s[], int len){
int l = 0;
ma[l++] = '$';
ma[l++] = '#';
for(int i = 0 ; i < len; i++){
ma[l++] = s[i];
ma[l++] = '#';
}
ma[l] = 0;
int mx = 0, id = 0;
for(int i = 0; i < l; i++){
mp[i] = mx > i ? min(mp[2*id - i],mx - i): 1;
while(ma[i + mp[i]] == ma[i - mp[i]])mp[i]++;
if(i + mp[i] > mx){
mx = i + mp[i];
id = i;
}
}
} int main(){
char s[maxn];
while(scanf("%s", s) == 1){
int len = strlen(s);
manacher(s, len);
int ans = 0;
for(int i = 0; i < 2 * len + 2; i++)
ans = max(ans, mp[i] -1);
printf("%d\n", ans);
}
return 0;
}

HDU 3068 最长回文 (Manacher最长回文串)的更多相关文章

  1. hdu 3068 最长回文 manacher算法(视频)

    感悟: 首先我要Orz一下qsc,我在网上很难找到关于acm的教学视频,但偶然发现了这个,感觉做的很好,链接:戳戳戳 感觉这种花费自己时间去教别人的人真的很伟大. manacher算法把所有的回文都变 ...

  2. HDU - 3068 最长回文(manacher)

    HDU - 3068 最长回文 Time Limit: 2000MS   Memory Limit: 32768KB   64bit IO Format: %I64d & %I64u Subm ...

  3. hdu 3068 最长回文 (Manacher算法求最长回文串)

    参考博客:Manacher算法--O(n)回文子串算法 - xuanflyer - 博客频道 - CSDN.NET 从队友那里听来的一个算法,O(N)求得每个中心延伸的回文长度.这个算法好像比较偏门, ...

  4. HDU 3068 最长回文( Manacher模板题 )

    链接:传送门 思路:Manacher模板题,寻找串中的最长回文子串 /***************************************************************** ...

  5. HDU 3068 (Manacher) 最长回文

    求一个字符串的最长子串,Manacher算法是一种O(n)的算法,很给力! s2[0] = '$',是避免在循环中对数组越界的检查. 老大的代码: http://www.cnblogs.com/Big ...

  6. 最长回文子串(百度笔试题和hdu 3068)

    版权所有.所有权利保留. 欢迎转载,转载时请注明出处: http://blog.csdn.net/xiaofei_it/article/details/17123559 求一个字符串的最长回文子串.注 ...

  7. hdu 3068 最长回文(manachar求最长回文子串)

    题目连接:hdu 3068 最长回文 解题思路:通过manachar算法求最长回文子串,如果用遍历的话绝对超时. #include <stdio.h> #include <strin ...

  8. hdu_3068 最长回文(Manacher算法)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3068 最长回文 Time Limit: 4000/2000 MS (Java/Others)    M ...

  9. Manacher 算法(hdu 3068 && hdu 3294)

    今天打算补前晚 BC 的第二题,发现要用到能在 O(n) 时间求最大回文子串长度的 Manacher 算法,第一次听,于是便去百度了下,看了大半天,总算能看懂了其思想,至于他给出的代码模板我没能完全看 ...

  10. 【BZOJ-2342】双倍回文 Manacher + 并查集

    2342: [Shoi2011]双倍回文 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1799  Solved: 671[Submit][Statu ...

随机推荐

  1. Qt中设置窗口图标

    转:https://blog.csdn.net/weiren2006/article/details/7438028 1.通过qtcreator新建一个文件filename.qrc,将图片添加到fil ...

  2. 拥有 GitHub 开源项目的小伙伴,免费申请 JetBrains 全家桶的全流程详解

    工欲善其事,必先利其器.如果您想要学习 Java.PHP.Ruby.Python.JavaScript.Objective-C..NET 中的任何一种开发技术,国际知名且屡获殊荣的 JetBrains ...

  3. python中方法调用和函数调用的区别

    函数调用: 传几个参数,就会有几个实参方法调用: 默认传递一个参数self,至少要定义一个形参

  4. win10内存泄漏怎么办

    我的电脑莫名的内存占用越来越高,但是实际内存很低 查了很多资料 下边是我找到的最有效的方法: 原因:是系统自带的网络数据监控和和Killer网卡的监控程序冲突,导致 非页面缓存无法释放. 解决方法:解 ...

  5. 关于爬虫的日常复习(8)—— 实战:request+正则爬取猫眼榜单top100

  6. 2019CSP初赛游记

    Day 0 作为一个初三的小蒟蒻…… 对于J+S两场比赛超级紧张的…… 教练发的神奇的模拟卷…… 我基本不会…… 就这样吧…… Day 1 Morning 不知道怎么就进了考场…… 周围坐的全是同学( ...

  7. Java入门 - 语言基础 - 11.switch_case

    原文地址:http://www.work100.net/training/java-switch-case.html 更多教程:光束云 - 免费课程 switch_case 序号 文内章节 视频 1 ...

  8. 利用Springmvc的AbstractXlsxView下载Excel文件

    设计一个模型,在针对多中数据类型进行拓展 public abstract class ExcelView extends AbstractXlsxView { public CellStyle cel ...

  9. springcloud初次zuul超时报错com.netflix.zuul.exception.ZuulException:Forwarding error

    报错如下 com.netflix.zuul.exception.ZuulException:Forwarding error Caused by: com.netflix.hystrix.except ...

  10. 将jar包安装到本地仓库

    通过cmd切换到apache maven 的bin目录 mvn install:install-file -DgroupId=com.antgroup.zmxy -DartifactId=zmxy-s ...