最长回文子串--轻松理解Manacher算法
最长回文子串这个问题的Manacher算法,看了很多博客,好不容易理解了,做一下记录。
这个算法的核心就是:将已经查询过的子字符串的最右端下标保存下来,在计算下标为i的回文字符串时,不需要从左右相邻的地方开始比较遍历,而是从某个初始值开始。
那么求这个初值就是该算法的关键。
1.字符串的初始化
先将字符串的每两个字符之间插入标识符,如“#”,然后在头尾也插入,插入什么符号这个其实影响不大。我是在头部和尾部也插入的“#”。这一步是为了让对称轴都在字符串数组中。
实例字符串 abada --> #a#b#a#d#a#
2.计算回文长度的初始
下标i从1 ~ length-1 开始遍历。同时借助2个辅助量 mid 和 maxRight ,maxRight 是用来存储我们扫描过的字符串的最右端的下标。mid为扫描最右端时的对称轴下标。注意,最右端的回文串不一定是最长的回文串。我们只是用maxRight来标识扫描长度。
我们用p[i]来存储回文字符串的单边长度。即以下标i的字符为对称轴的回文串的最右端与i的差值+1。回文串长度/2+1。
示例:
注意:maxRight对应下标时需要-1去对应。因为每次从center开始扩张时,结束条件是maxRight时center不能成为回文中心。
核心代码为这一句:
p[i] = maxRight > i ? Math.min(p[2*center-i],maxRight-i) : 1 ;
可能这样比较难以理解。我们可以结合以上的例子,来慢慢分析。
str[center] - str[maxRight-1]为我们的最右端回文字符串。我们维护这个字符串的位置。
当我们处理str[i]这个字符的回文长度时。
if ( maxRight > i){
//如果当前的位置已经被扫描过了(maxRight为我们扫描过的最右端)
//center肯定在左边,因为是随着i++扫描来更新center的
//假设j是i关于center的对称下标,即i+j=2*center
//string下标示意i图
//0----j---center---i----maxRight------length-1
int j = 2*center-i;
//算法核心:
//由于i j 关于center对称 而 center到maxRight在center~【center-maxRight】也是相同的
//另外由于j的最长回文串长度我们在计算i的时候已经算好了,那么可以认为i的周围也有这么p[j]长度的对称回文串
//center左右对称 j左右长度p[j]字符串对称 那么i左右也有对称字符串
//那么此时分为2种情况
//第一种 ---(j-p[j])--j--(j+p[i])-----center---(i-p[j])----i---(i+p[j])----maxRight
//第二种 ------(i-p[j])----center------i---maxRight----(i+p[j])
if ( p[j] + i<maxRight){
//第一种情况,已扫描范围内就已经找到最大值了
p[i] = p[j];
}else {
//第二张情况,maxRight太小,后面的可能相等,需要左右扩张来判断
p[i] = maxRight - i ;
}
}else {
//0-------maxRight--i
//p[i]还没有扫描呢,因此初始长度就是1
p[i] = 1;
}
上面的核心分析过程简化后就是上面的那一行代码。
核心代码理解了,后续的就简单了,从p[i]的初始值开始,向左右扩张,判断左右边界是不是相等,并同时更新center和maxRight的值。最后扫描p[]数组,其中最大的就是回文子串的长度,下标就是扩充了#的字符串下标。
最后,贴上源代码,大家看一看就可以理解了。
public static String longestPalindrome(String s) {
if ( s.length() <= 1 )return s ;
StringBuilder builder = new StringBuilder();
builder.append("#");
for (int i = 0; i <s.length() ; i++) {
builder.append(s.charAt(i));
builder.append("#");
}
int center = 0 ,maxRight = 0 ,len = builder.length();
int[] p = new int[len];
for (int i = 1; i < len-1; i++) {
//确定已扫描过的字符串中 最长回文串的初始值
//减少的时间复杂度 就在这里
//maxRight都被扫描,初值就不用从0开始了。
p[i] = maxRight > i ? Math.min(p[2*center-i],maxRight-i) : 1 ;
//确定p[i]的初始值后继续扩张判断回文
while (i+p[i] <len && i-p[i]>=0 &&builder.charAt(i+p[i])==builder.charAt(i-p[i])){
p[i]++;
}
//更新maxRight center
if ( i + p[i] > maxRight){
maxRight = i + p[i] ;
center = i ;
}
}
//遍历p数组求极大值
int mid= 0 , maxLength = 0 ;
for (int i = 1; i < builder.length(); i++) {
if ( p[i] > maxLength){
maxLength =p[i];
mid = i;
}
}
//分割字符串,消去"#"
return builder.toString().substring(mid-maxLength+1,mid+maxLength-1).replace("#","");
}
最长回文子串--轻松理解Manacher算法的更多相关文章
- 1089 最长回文子串 V2(Manacher算法)
1089 最长回文子串 V2(Manacher算法) 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注 回文串是指aba.abba.cccbccc.aaaa ...
- 51Nod 1089:最长回文子串 V2(Manacher算法)
1089 最长回文子串 V2(Manacher算法) 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注 回文串是指aba.abba.cccbccc.aaa ...
- 51NOD 1088 最长回文子串&1089 最长回文子串 V2(Manacher算法)
回文串是指aba.abba.cccbccc.aaaa这种左右对称的字符串. 输入一个字符串Str,输出Str里最长回文子串的长度. Input 输入Str(Str的长度 <= 1000(第二题要 ...
- [51Nod1089] 最长回文子串 V2(Manacher算法)
1089 最长回文子串 V2(Manacher算法) 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 回文串是指aba.abba.cccbccc.aaaa这种左右对称 ...
- hdu5371 最长回文子串变形(Manacher算法)
pid=5371">http://acm.hdu.edu.cn/showproblem.php? pid=5371 Problem Description Hotaru Ichijou ...
- 【回文字符串】 最长回文子串O(N) Manacher算法
原理讲的清晰:Manacher's ALGORITHM: O(n)时间求字符串的最长回文子串 注意: ①动态生命P[]和newStr数组后,不要忘记delete[] //其实这是基本的编码习惯 ②最终 ...
- 最长回文子串问题-Manacher算法
转:http://blog.csdn.net/dyx404514/article/details/42061017 Manacher算法 算法总结第三弹 manacher算法,前面讲了两个字符串相算法 ...
- 51nod 1089 最长回文子串 V2(Manacher算法)
回文串是指aba.abba.cccbccc.aaaa这种左右对称的字符串. 输入一个字符串Str,输出Str里最长回文子串的长度. 收起 输入 输入Str(Str的长度 <= 100000) ...
- 【51NOD-0】1089 最长回文子串 V2(Manacher算法)
[算法]回文树 #include<cstdio> #include<algorithm> #include<cstring> using namespace std ...
随机推荐
- javascript(6)
javascript的冒泡排序 var arr=[5,0,-56,900,12]; //大的排序次数(arr.length-1) var flag=false; for(var i=0;i<=a ...
- Linux下查看端口占用情况
用启动服务的账号登录,然后运行命令: lsof -i:<端口号> 也可使用命令: netstat -apn|grep <端口号> 找到进程号以后,再使用以下命令查看详细信息: ...
- mac快速正确的安装 Ruby, Rails 运行环境
Mac OS X 任意 Linux 发行版本(Ubuntu,CentOS, Redhat, ArchLinux ...) 强烈新手使用 Ubuntu 省掉不必要的麻烦! 以下代码区域,带有 $ 打头的 ...
- install docker
摘要: 我的环境是:CentOS-7-x86_64-Minimal-1511.iso , 也可参考docker官网文档,来安装, url : https://docs.docker.com/engin ...
- codeforces 558B B. Amr and The Large Array(水题)
题目链接: B. Amr and The Large Array time limit per test 1 second memory limit per test 256 megabytes in ...
- 第十四章-MySQL
1 安装 MySQL常见的版本 GA: 广泛使用的版本 RC: 最接近正式版本 Alpha和Bean: 内测版本和公测版本 有两种安装方式: 安装包和压缩包 1) 安装msi文件 2) 解压zip文件 ...
- (Nginx + Gunicorn) 时 nginx 400 request line is too large (4360 4094)
查看nginx下面两个参数 值是否满足 client_header_buffer_size 512k;large_client_header_buffers 4 512k; 满足依然出现 如果ngin ...
- bzoj 4753 [Jsoi2016]最佳团体——0/1分数规划
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4753 0/1分数规划裸题. #include<iostream> #includ ...
- VijosP1180:选课
描述 学校实行学分制.每门的必修课都有固定的学分,同时还必须获得相应的选修课程学分.学校开设了N(N<300)门的选修课程,每个学生可选课程的数量M是给定的.学生选修了这M门课并考核通过就能获得 ...
- POJ1986(LCA应用:求两结点之间距离)
Distance Queries Time Limit: 2000MS Memory Limit: 30000K Total Submissions: 11304 Accepted: 3985 ...