- 字符串算法 — Manacher算法
首先介绍基础入门知识,以下这部分来着一贴吧,由于是很久之前看的,最近才整理一下,发现没有保存链接,请原创楼主见谅。
//首先:大家都知道什么叫回文串吧,这个算法要解决的就是一个字符串中最长的回文子串有多长。 //这个算法可以在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
1 #include<iostream> 2 #include<cstring> 3 using namespace std; 4 5 const int N=1000010;; 6 7 char s1[N]; 8 char s[N*2]; 9 int p[N*2]; 10 int n; 11 12 int min(int a,int b) 13 { 14 if(a>b) 15 return b; 16 else 17 return a; 18 } 19 20 int max(int a,int b) 21 { 22 if(a>b) 23 return a; 24 else 25 return b; 26 } 27 28 void Manacher() 29 { 30 int i; 31 int mx=0; 32 int id; 33 for(i=1;i<n;i++) 34 { 35 if(mx>i) 36 p[i]=min(p[2*id-i],mx-i); 37 else 38 p[i]=1; 39 for( ;s[i+p[i]]==s[i-p[i]];p[i]++); 40 if(p[i]+i>mx) 41 { 42 mx=p[i]+i; 43 id=i; 44 } 45 } 46 } 47 48 int main() 49 { 50 int count=1; 51 while(scanf("%s",s1)!=EOF,strcmp(s1,"END")) 52 { 53 int ans=1; 54 int len=strlen(s1); 55 n=len*2; 56 s[0]='$'; 57 // s[1]='#'; 58 for(int i=0;i<=len;i++) 59 { 60 s[2*i+1]='#'; 61 s[2*i+2]=s1[i]; 62 } 63 memset(p,0,sizeof(p)); 64 Manacher(); 65 for(i=0;i<n;i++) 66 ans=max(ans,p[i]); 67 printf("Case %d: %d\n",count++,ans-1); 68 } 69 return 0; 70 }