学习manacher(最长公共回文串算法)
给定一个字符串求出其中最长个公共回文串.
举列子: abab -->回文串长度为2
以前的算法诸如: 扩展kmp求法过于麻烦,看到有一篇博文(http://leetcode.com/2011/11/longest-palindromic-substring-part-ii.html),写了一个关于这样的算法,按耐不住自己内心的小激动,就去学了下,于是将自己学习的一点点的理解记录下来:
它的处理方法是: 对于这样一个字符串: abab;
a b a b 我们对于每一个字符串的进行逐个匹配
1 2 2 1
对于奇数长度的字符串 ababa -->回文串的长度为:
a b a b a ( 这是对于每一个字符串进行模式匹配) 之间复杂度为 O(n*m)
1 2 3 2 1
而现在我们学习的一种算法是时间复杂度为O(n)
算法的思想是: 设置一个匹配半径ra , 然后分为三种情况进行讨论: (当然讨论主要是缩减不必要的重复匹配)
由于已经有好些博主写的较为详细,自己就不在胡乱填鸭了.....
以下引用的是 这位楼主的博客 (http://blog.sina.com.cn/s/blog_70811e1a01014esn.html)
首先,在字符串s中,用rad[i]表示第i个字符的回文半径,即rad[i]尽可能大,且满足: s[i-rad[i],i-1]=s[i+1,i+rad[i]] 很明显,求出了所有的rad,就求出了所有的长度为奇数的回文子串. 至于偶数的怎么求,最后再讲. 假设现在求出了rad[1..i-1],现在要求后面的rad值,并且通过前面的操作,得知了当前字符i的rad值至少为j.现在通过试图扩大j来扫描, 求出了rad[i].再假设现在有个指针k,从1循环到rad[i],试图通过某些手段来求出[i+1,i+rad[i]]的rad值. 根据定义,黑色的部分是一个回文子串,两段红色的区间全等. 因为之前已经求出了rad[i-k],所以直接用它.有3种情况:
①rad[i]-k<rad[i-k] 如图,rad[i-k]的范围为青色.因为黑色的部分是回文的,且青色的部分超过了黑色的部分,所以rad[i+k]肯定至少为rad[i]-k,即橙色 的部分.那橙色以外的部分就不是了吗?这是肯定的.因为如果橙色以外的部分也是回文的,那么根据青色和红色部分的关系,可以证明黑色部分再往外延伸一点也 是一个回文子串,这肯定不可能,因此rad[i+k]=rad[i]-k.为了方便下文,这里的rad[i+k]=rad[i]- k=min(rad[i]-k,rad[i-k]).
②rad[i]-k>rad[i-k] 如图,rad[i-k]的范围为青色.因为黑色的部分是回文的,且青色的部分在黑色的部分里面,根据定义,很容易得出:rad[i+k]=rad[i-k].为了方便下文,这里的rad[i+k]=rad[i-k]=min(rad[i]-k,rad[i-k]). 根据上面两种情况,可以得出结论:当rad[i]-k!=rad[i-k]的时候,rad[i+k]=min(rad[i]-k,rad[i-k]). 注意:当rad[i]-k==rad[i-k]的时候,就不同了,这是第三种情况:
如图,通过和第一种情况对比之后会发现,因为青色的部分没有超出黑色的部分,所以即使橙色的部分全等,也无法像第一种情况一样引出矛盾,因此橙色的部分是 有可能全等的,但是,根据已知的信息,我们不知道橙色的部分是多长,因此就把i指针移到i+k的位置,j=rad[i-k](因为它的rad值至少为 rad[i-k]),等下次循环的时候再做了. 整个算法就这样. 至于时间复杂度为什么是O(n),我已经证明了,但很难说清楚.所以自己体会吧. 上文还留有一个问题,就是这样只能算出奇数长度的回文子串,偶数的就不行.怎么办呢?有一种直接但比较笨的方法,就是做两遍(因为两个程序是差不多的,只 是rad值的意义和一些下标变了而已).但是写两个差不多的程序是很痛苦的,而且容易错.所以一种比较好的方法就是在原来的串中每两个字符之间加入一个特 殊字符,再做.如:aabbaca,把它变成a#a#b#b#a#c#a,这样的话,无论原来的回文子串长度是偶数还是奇数,现在都变成奇数了
于是有这种思想我们不难写出这样的代码:
代码:
1 #include<stdio.h>
2 #include<string.h>
3 #include<stdlib.h>
4 #define maxn 1000000
5 char str[maxn];
6 int ra[maxn];
7 int min(int a,int b)
8 {
9 return a<b?a:b;
10 }
11 void init(char s[],int len){
12 int en=(len<<1)+1;
13 memset(ra,0,sizeof(int)*(en+1));
14 while(en>0){
15 if(en&1) str[en]='#';
16 else str[en]=str[--len];
17 en--;
18 }
19 s[en]='$';
20 }
21 int manacher (int len ){
22 len=len*2-1;
23 int ce=0,mx=0,ans=0;
24 for(int i=1;i<=len;i++){
25 if(mx>i) ra[i]=min(ra[(ce<<1)-i],mx-i); //位置在回文串的范围之内
26 else ra[i]=1;
27 while(str[i+ra[i]]==str[i-ra[i]]) //以i点为中心进行模式匹配
28 ra[i]++;
29 if(i+ra[i]>mx){
30 mx=i+ra[i];
31 ce=i;
32 }
33 if((i&1)&&ans<ra[i])
34 ans=ra[i];
35 }
36 memset(str,0,sizeof(str));
37 return ans;
38 }
39 int main(){
40 while(scanf("%s",str)!=EOF){
41 int len=strlen(str);
42 init(str,len);
43 printf("%dn",manacher(len)-1);
44 }
45 return 0;
46 }
- 一分钟教你在博客园中制作自己的动态云球形标签页
- Python标准库06 子进程 (subprocess包)
- 摩拜、美团也伸手共享汽车拉!众多品牌里面,这三家最牛!
- Android Studio快捷键每日一练(2)
- Linux进程间通信
- Android Studio快捷键每日一练(1)
- Android开发中的全屏背景显示方案
- Android中的FragmentManager的问题
- Premiere Pro & After Effects插件开发调试方法
- RegQueryValueEx正确使用方法
- 区块链+医疗五类应用前景广阔,英美已有企业试水
- Linux进程关系
- 基于OBS的插件开发总结
- Linux从程序到进程
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- python实现线性回归之岭回归
- 操作系统实验之存储管理
- MySQL EXPLAIN 的使用
- mybatis文件映射之关联查询初探(一)
- python实现线性回归之弹性网回归
- 【原创】python倒排索引之查找包含某主题或单词的文件
- python实现逻辑回归
- Linux文件管理参考
- CloudBase Framework丨第一个 Deno 部署工具是如何打造的?
- 关于null通过+" ",String.ValueOf转换为字符串的问题!!!
- Java实现尺取法
- 【自然语言处理】利用朴素贝叶斯进行新闻分类(自己处理数据)
- mybatis文件映射之select操作返回Map
- shm进程间通信失败了!!!
- Vue3.0 beta源码学习笔记(二)