首次学习manacher算法,感谢分享优秀博文的博主:
http://blog.csdn.net/xingyeyongheng/article/details/9310555
读完后,我对其有了大体的认识。
经过manacher算法处理后的字符串,其内容完全改变。比如字符串212:
设p[i]是i位置处向外拓展的半径。
after the change, it’s
type | contents |
---|---|
changed string | *#2#1#2#0 |
index (i) | 012345678 |
p[i] | 0021412 |
manacher算法的核心代码:
(下面代码主要用于说明)
int manacher(char *s, int p[]){
memset(p,0,sizeof(p));
int len=strlen(s);
for(int i=len;i>=0;i--){
s[i+i+2]=s[i];
s[i+i+1]='#';
}
s[0]='*';
len<<=1;
len++;
int l=0,maxn=0;
for(int i=2;i<len;i++){
if(p[l]+l>i) p[i]=min(p[l]+l-i,p[2*l-i]);
else p[i]=1;
while(s[i-p[i]]==s[i+p[i]]) p[i]++;
if(l+p[l]<i+p[i]) l=i;
if(maxn<p[i]) maxn=p[i];
}
return maxn-1;
}
接下来献上三道基础题:
POj 3974 Palindrome
http://poj.org/problem?id=3974
求出串中回文子串的最大长度。
分析:manacher最简单的应用。
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N=1e6+10;
char s[2*N+10];
int p[2*N+10];
int main(int argc, char *argv[]) {
int ca=1;
while(~scanf("%s",s)){
if(strcmp(s,"END")==0) break;
memset(p,0,sizeof(p));
int len=strlen(s);
for(int i=len;i>=0;i--){
s[i+i+2]=s[i];
s[i+i+1]='#';
}
s[0]='*';
len<<=1;
len++;
int l=0,maxn=0;
for(int i=2;i<len;i++){
if(p[l]+l>i) p[i]=min(p[l]+l-i,p[2*l-i]);
else p[i]=1;
while(s[i-p[i]]==s[i+p[i]]) p[i]++;
if(l+p[l]<i+p[i]) l=i;
if(maxn<p[i]) maxn=p[i];
}
printf("Case %d: %d\n",ca++,maxn-1);
}
return 0;
}
hdu 3294 Girls’ research
http://acm.hdu.edu.cn/showproblem.php?pid=3294
大意:先将整个字符串转化成原本的串(由给出的密码字符确定),然后由manacher算法求出最长的回文字符串,要求输出起点,终点,内容。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=2e5+10,M=2*N+10;
char s[M],key;
int p[M];
int main()
{
//freopen("cin.txt","r",stdin);
int ca=1;
while(~scanf("%c%s",&key,s)){
getchar();
int d=key-'a';
for(int i=0;s[i];i++) {
s[i]-='a';
s[i]=(s[i]-d+26)%26;
s[i]+='a'; //cout<<s[i];
}
memset(p,0,sizeof(p));
int len=strlen(s);
for(int i=len;i>=0;i--){
s[i+i+2]=s[i];
s[i+i+1]='#';
}
s[0]='*';
len<<=1; len++;
int l=0,maxn=0,ld=0,rd=0;
for(int i=2;i<len;i++){
int tl=0,tr=0;
if(p[l]+l>i) p[i]=min(p[l]+l-i,p[2*l-i]);
else p[i]=1;
while(s[i-p[i]]==s[i+p[i]]) {
tl=i-p[i]; tr=i+p[i];
p[i]++;
}
if(l+p[l]<p[i]+i) l=i;
if(maxn<p[i]) {
maxn=p[i];
ld=tl; rd=tr;
}
}
if(maxn-1>=2){
while(s[ld]<97||s[ld]>122) ld++;
while(s[rd]<97||s[rd]>122) rd--;
printf("%d %d\n",ld/2-1,rd/2-1);
for(int i=ld;i<=rd;i++) if(s[i]<=122&&s[i]>=97) printf("%c",s[i]);
puts("");
}
else puts("No solution!");
}
return 0;
}
hdu 4513 吉哥系列故事――完美队形II
http://acm.hdu.edu.cn/showproblem.php?pid=4513
大意:求出完美串——回文子串,且满足H[1] <= H[2] <= H[3] …. <= H[mid] 山峰。
分析:中间过程处理一下,根据manacher算法得到的字符串的特点来更新山峰回文序列。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=1e5+10;
int s[2*N+10];
int p[2*N+10];
int main() {
//freopen("cin.txt","r",stdin);
int t,n;
cin>>t;
while(t--){
scanf("%d",&n);
for(int i=0;i<n;i++) scanf("%d",&s[i]);
memset(p,0,sizeof(p));
int len=n;
for(int i=len;i>=0;i--){
s[i+i+2]=s[i];
s[i+i+1]=1;
}
s[0]=1000;
len<<=1;
len++;
s[len+1]=1001;
int l=0,maxn=0;
for(int i=2;i<len;i++){
if(p[l]+l>i) p[i]=min(p[l]+l-i,p[2*l-i]);
else p[i]=1;
int tt=s[i];
if(tt<50){
if(s[i-p[i]]==s[i+p[i]]) {
tt=s[i+p[i]];
p[i]++;
}
else continue;
}
while(s[i+p[i]]<50) p[i]++; // very important
while(s[i-p[i]]==s[i+p[i]] && s[i+p[i]]<=tt)
{
tt=s[i+p[i]];
p[i]+=2;
}
if(p[l]+l<p[i]+i) l=i;
if(maxn<p[i]) maxn=p[i];
}
//for(int i=0;i<=len+2;i++) cout<<s[i]<<" "; cout<<endl;
//for(int i=0;i<=len+2;i++) cout<<p[i]<<" "; cout<<endl;
printf("%d\n",maxn-1);
}
return 0;
}