ACM -- 算法小结(八)字符串算法之Manacher算法
2013-09-26 14:23
288 查看
字符串算法 -- Manacher算法
首先介绍基础入门知识,以下这部分来着一贴吧,由于是很久之前看的,最近才整理一下,发现没有保存链接,请原创楼主见谅。
核心代码如下:
接下来到了小试牛刀的时刻啦~~ poj 3974 ,hdoj3068
首先介绍基础入门知识,以下这部分来着一贴吧,由于是很久之前看的,最近才整理一下,发现没有保存链接,请原创楼主见谅。
//首先:大家都知道什么叫回文串吧,这个算法要解决的就是一个字符串中最长的回文子串有多长。 //这个算法可以在O(n)的时间复杂度内既线性时间复杂度的情况下,求出以每个字符为中心的最长回文有多长, //这个算法有一个很巧妙的地方,它把奇数的回文串和偶数的回文串统一起来考虑了。 //这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是 //充分利用了字符匹配的特殊性,避免了大量不必要的重复匹配。 //算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符,当然这个分隔符要在原串中没有出现过。 //一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了 //(见下面的一个例子,回文串长度全为奇数了),然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。 //P[id]记录的是以字符str[id]为中心的最长回文串,当以str[id]为第一个字符,这个最长回文串向右延伸了P[id]个字符。 // 原串: w aa bwsw f d // 新串: # w # a # a # b # w # s # w # f # d # //辅助数组P 1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1 // 这里有一个很好的性质,P[id]-1就是该回文子串在原串中的长度(包括‘#’)。如果这里不是特别清楚, //可以自己拿出纸来画一画,自己体会体会。当然这里可能每个人写法不尽相同,不过我想大致思路应该是一样的吧。 // 好,我们继续。现在的关键问题就在于怎么在O(n)时间复杂度内求出P数组了。只要把这个P数组求出来,最长回文子串就可以直接扫一遍得出来了。 // 由于这个算法是线性从前往后扫的。那么当我们准备求P[i]的时候,i以前的P[j]我们是已经得到了的。 //我们用mx记在i之前的回文串中,延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。 //(注:为了防止字符比较的时候越界,我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’,故我的新串下标是从1开始的) //好,到这里,我们可以先贴一份代码了。
核心代码如下:
void pk() { int i; int mx = 0; int id; for(i=1; i<n; i++) { if( mx > i ) p[i] = MIN( p[2*id-i], mx-i ); else p[i] = 1; for(; str[i+p[i]] == str[i-p[i]]; p[i]++) ; if( p[i] + i > mx ) { mx = p[i] + i; id = i; } } }
接下来到了小试牛刀的时刻啦~~ poj 3974 ,hdoj3068
#include<iostream> #include<cstring> using namespace std; const int N=1000010;; char s1 ; char s[N*2]; int p[N*2]; int n; int min(int a,int b) { if(a>b) return b; else return a; } int max(int a,int b) { if(a>b) return a; else return b; } void Manacher() { int i; int mx=0; int id; for(i=1;i<n;i++) { if(mx>i) p[i]=min(p[2*id-i],mx-i); else p[i]=1; for( ;s[i+p[i]]==s[i-p[i]];p[i]++); if(p[i]+i>mx) { mx=p[i]+i; id=i; } } } int main() { int count=1; while(scanf("%s",s1)!=EOF,strcmp(s1,"END")) { int ans=1; int len=strlen(s1); n=len*2; s[0]='$'; // s[1]='#'; for(int i=0;i<=len;i++) { s[2*i+1]='#'; s[2*i+2]=s1[i]; } memset(p,0,sizeof(p)); Manacher(); for(i=0;i<n;i++) ans=max(ans,p[i]); printf("Case %d: %d\n",count++,ans-1); } return 0; }
相关文章推荐
- ACM -- 算法小结(五)字符串算法之Sunday算法
- [ACM] [算法基础] POJ 2503 在线翻译Babelfish(二分查找 字符串)
- ACM -- 算法小结(一)利用数组存放实现排序
- ACM中常用算法—-字符串
- ACM -- 算法小结(三)反转句子顺序与反转单词组成顺序
- ACM -- 算法小结(二)错排公式的应用
- ACM经典算法之字符串处理:字符串查找
- ACM -- 算法小结(四)KMP(POJ3461)
- 字符串模式匹配算法小结
- ACM经典算法之字符串处理:字符串截取
- 一些C语言中字符串的算法问题解决实例小结
- ACM经典算法之字符串处理:字符串替换
- ACM 算法竞赛 八大基础算法 模拟法 字符串 递归分治 排序 枚举 贪心 递推
- ACM -- 算法小结(六)逆波兰表达式
- ACM -- 算法小结(七)Phone list解题报告
- ACM -- 算法小结(九)DP之Humble numbers
- ACM经典算法之字符串处理:字符串替换
- ACM经典算法之字符串处理
- 算法练习5. Longest Palindromic Substring 最长回文子字符串(Manacher算法)
- 字符串专题小结:Manacher算法求最长回文串